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.

Multimodal analysis for human ex vivo studies shows extensive molecular changes from delays in blood processing.

Multi-omic profiling of human peripheral blood is increasingly utilized to identify biomarkers and pathophysiologic mechanisms of disease. The importance of these platforms in clinical and translational studies led us to investigate the impact of delayed blood processing on the numbers and state of peripheral blood mononuclear cells (PBMC) and on the plasma proteome. Similar to previous studies, we show minimal effects of delayed processing on the numbers and general phenotype of PBMC up to 18 hours. In contrast, profound changes in the single-cell transcriptome and composition of the plasma proteome become evident as early as 6 hours after blood draw. These reflect patterns of cellular activation across diverse cell types that lead to progressive distancing of the gene expression state and plasma proteome from native in vivo biology. Differences accumulating during an overnight rest (18 hours) could confound relevant biologic variance related to many underlying disease states.

The COVID-19 immune landscape is dynamically and reversibly correlated with disease severity.

BACKGROUNDDespite a rapidly growing body of literature on coronavirus disease 2019 (COVID-19), our understanding of the immune correlates of disease severity, course, and outcome remains poor.METHODSUsing mass cytometry, we assessed the immune landscape in longitudinal whole-blood specimens from 59 patients presenting with acute COVID-19 and classified based on maximal disease severity. Hospitalized patients negative for SARS-CoV-2 were used as controls.RESULTSWe found that the immune landscape in COVID-19 formed 3 dominant clusters, which correlated with disease severity. Longitudinal analysis identified a pattern of productive innate and adaptive immune responses in individuals who had a moderate disease course, whereas those with severe disease had features suggestive of a protracted and dysregulated immune response. Further, we identified coordinate immune alterations accompanying clinical improvement and decline that were also seen in patients who received IL-6 pathway blockade.CONCLUSIONThe hospitalized COVID-19 negative cohort allowed us to identify immune alterations that were shared between severe COVID-19 and other critically ill patients. Collectively, our findings indicate that selection of immune interventions should be based in part on disease presentation and early disease trajectory due to the profound differences in the immune response in those with mild to moderate disease and those with the most severe disease.FUNDINGBenaroya Family Foundation, the Leonard and Norma Klorfine Foundation, Glenn and Mary Lynn Mounger, and the National Institutes of Health.

SBML Level 3: an extensible format for the exchange and reuse of biological models.

Systems biology has experienced dramatic growth in the number, size, and complexity of computational models. To reproduce simulation results and reuse models, researchers must exchange unambiguous model descriptions. We review the latest edition of the Systems Biology Markup Language (SBML), a format designed for this purpose. A community of modelers and software authors developed SBML Level 3 over the past decade. Its modular form consists of a core suited to representing reaction-based models and packages that extend the core with features suited to other model types including constraint-based models, reaction-diffusion models, logical network models, and rule-based models. The format leverages two decades of SBML and a rich software ecosystem that transformed how systems biologists build and interact with models. More recently, the rise of multiscale models of whole cells and organs, and new data sources such as single-cell measurements and live imaging, has precipitated new ways of integrating data with models. We provide our perspectives on the challenges presented by these developments and how SBML Level 3 provides the foundation needed to support this evolution.

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.

Integrative network modeling reveals mechanisms underlying T cell exhaustion.

Failure to clear antigens causes CD8+ T cells to become increasingly hypo-functional, a state known as exhaustion. We combined manually extracted information from published literature with gene expression data from diverse model systems to infer a set of molecular regulatory interactions that underpin exhaustion. Topological analysis and simulation modeling of the network suggests CD8+ T cells undergo 2 major transitions in state following stimulation. The time cells spend in the earlier pro-memory/proliferative (PP) state is a fixed and inherent property of the network structure. Transition to the second state is necessary for exhaustion. Combining insights from network topology analysis and simulation modeling, we predict the extent to which each node in our network drives cells towards an exhausted state. We demonstrate the utility of our approach by experimentally testing the prediction that drug-induced interference with EZH2 function increases the proportion of pro-memory/proliferative cells in the early days post-activation.

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.