T cells made with messenger RNA (mRNA) encoding chimeric antigen receptor (CAR) offer a safe alternative to those transduced with viral CARs by mitigating the side effects of constitutively active T cells. Previous studies have shown that mRNA CAR T cells are transiently effective but lack persistence and potency across tumor types. It was hypothesized that the efficacy of mRNA CARs could be improved by utilizing recent advancements in RNA technology, such as incorporating a modified nucleoside, 1-methylpseudouridine, into the mRNA and applying a novel purification method using RNase III to eliminate dsRNA contaminants. T cells electroporated with nucleoside-modified and purified mRNA encoding CD19 CAR showed an initial twofold increase in CAR surface expression, as well as a twofold improvement in cytotoxic killing of leukemia cells that persisted up to 5 days. T cells generated with nucleoside-modified and purified CAR mRNA also showed reduced expression of checkpoint regulators and a differential pattern of genetic activation compared to those made with conventional mRNA. In vivo studies using a leukemia mouse model revealed that the most robust 100-fold suppression of leukemic burden was achieved using T cells electroporated with purified mRNAs, regardless of their nucleoside modification. The results provide a novel approach to generate mRNA for clinical trials, and poise mRNA CAR T cells for increased efficacy during testing as new CAR targets emerge.

We present an extensive assessment of mutation burden through sequencing analysis of >81,000 tumors from pediatric and adult patients, including tumors with hypermutation caused by chemotherapy, carcinogens, or germline alterations. Hypermutation was detected in tumor types not previously associated with high mutation burden. Replication repair deficiency was a major contributing factor. We uncovered new driver mutations in the replication-repair-associated DNA polymerases and a distinct impact of microsatellite instability and replication repair deficiency on the scale of mutation load. Unbiased clustering, based on mutational context, revealed clinically relevant subgroups regardless of the tumors' tissue of origin, highlighting similarities in evolutionary dynamics leading to hypermutation. Mutagens, such as UV light, were implicated in unexpected cancers, including sarcomas and lung tumors. The order of mutational signatures identified previous treatment and germline replication repair deficiency, which improved management of patients and families. These data will inform tumor classification, genetic testing, and clinical trial design.

Prolonged exit from quiescence by hematopoietic stem cells (HSCs) progressively impairs their homeostasis in the bone marrow through an unidentified mechanism. We show that Rb proteins, which are major enforcers of quiescence, maintain HSC homeostasis by positively regulating thrombopoietin (Tpo)-mediated Jak2 signaling. Rb family protein inactivation triggers the progressive E2f-mediated transactivation of Socs3, a potent inhibitor of Jak2 signaling, in cycling HSCs. Aberrant activation of Socs3 impairs Tpo signaling and leads to impaired HSC homeostasis. Therefore, Rb proteins act as a central hub of quiescence and homeostasis by coordinating the regulation of both cell cycle and Jak2 signaling in HSCs.

Metastatic breast cancer is a leading cause of cancer-related deaths in women worldwide. Patients with triple negative breast cancer (TNBCs), a highly aggressive tumor subtype, have a particularly poor prognosis. Multiple reports demonstrate that altered content of the multicopy mitochondrial genome (mtDNA) in primary breast tumors correlates with poor prognosis. We earlier reported that mtDNA copy number reduction in breast cancer cell lines induces an epithelial-mesenchymal transition associated with metastasis. However, it is unknown whether the breast tumor subtypes (TNBC, Luminal and HER2+) differ in the nature and amount of mitochondrial defects and if mitochondrial defects can be used as a marker to identify tumors at risk for metastasis. By analyzing human primary tumors, cell lines and the TCGA dataset, we demonstrate a high degree of variability in mitochondrial defects among the tumor subtypes and TNBCs, in particular, exhibit higher frequency of mitochondrial defects, including reduced mtDNA content, mtDNA sequence imbalance (mtRNR1:ND4), impaired mitochondrial respiration and metabolic switch to glycolysis which is associated with tumorigenicity. We identified that genes involved in maintenance of mitochondrial structural and functional integrity are differentially expressed in TNBCs compared to non-TNBC tumors. Furthermore, we identified a subset of TNBC tumors that contain lower expression of epithelial splicing regulatory protein (ESRP)-1, typical of metastasizing cells. The overall impact of our findings reported here is that mitochondrial heterogeneity among TNBCs can be used to identify TNBC patients at risk of metastasis and the altered metabolism and metabolic genes can be targeted to improve chemotherapeutic response.

Accelerating cures for children with cancer remains an immediate challenge as a result of extensive oncogenic heterogeneity between and within histologies, distinct molecular mechanisms evolving between diagnosis and relapsed disease, and limited therapeutic options. To systematically prioritize and rationally test novel agents in preclinical murine models, researchers within the Pediatric Preclinical Testing Consortium are continuously developing patient-derived xenografts (PDXs)-many of which are refractory to current standard-of-care treatments-from high-risk childhood cancers. Here, we genomically characterize 261 PDX models from 37 unique pediatric cancers; demonstrate faithful recapitulation of histologies and subtypes; and refine our understanding of relapsed disease. In addition, we use expression signatures to classify tumors for TP53 and NF1 pathway inactivation. We anticipate that these data will serve as a resource for pediatric oncology drug development and will guide rational clinical trial design for children with cancer.

With the increasing rates of exome and whole genome sequencing, the ability to classify large sets of germline sequencing variants using up-to-date American College of Medical Genetics - Association for Molecular Pathology (ACMG-AMP) criteria is crucial. Here, we present Automated Germline Variant Pathogenicity (AutoGVP), a tool that integrates germline variant pathogenicity annotations from ClinVar and sequence variant classifications from a modified version of InterVar (PVS1 strength adjustments, removal of PP5/BP6). This tool facilitates large-scale, clinically-focused classification of germline sequence variants in a research setting.

Craniopharyngiomas are epithelial tumors that typically arise in the suprasellar region of the brain. Patients experience substantial clinical sequelae from both extension of the tumors and therapeutic interventions that damage the optic chiasm, the pituitary stalk and the hypothalamic area. Using whole-exome sequencing, we identified mutations in CTNNB1 (β-catenin) in nearly all adamantinomatous craniopharyngiomas examined (11/12, 92%) and recurrent mutations in BRAF (resulting in p.Val600Glu) in all papillary craniopharyngiomas (3/3, 100%). Targeted genotyping revealed BRAF p.Val600Glu in 95% of papillary craniopharyngiomas (36 of 39 tumors) and mutation of CTNNB1 in 96% of adamantinomatous craniopharyngiomas (51 of 53 tumors). The CTNNB1 and BRAF mutations were clonal in each tumor subtype, and we detected no other recurrent mutations or genomic aberrations in either subtype. Adamantinomatous and papillary craniopharyngiomas harbor mutations that are mutually exclusive and clonal. These findings have important implications for the diagnosis and treatment of these neoplasms.

The genetic aetiology of sporadic neuroblastoma is still largely unknown. We have identified diverse neuroblastoma susceptibility loci by genomewide association studies (GWASs); however, additional SNPs that likely contribute to neuroblastoma susceptibility prompted this investigation for identification of additional variants that are likely hidden among signals discarded by the multiple testing corrections used in the analysis of genomewide data. There is evidence suggesting the CDKN1B, coding for the cycle inhibitor p27Kip1, is involved in neuroblastoma. We thus assess whether genetic variants of CDKN1B are associated with neuroblastoma. We imputed all possible genotypes across CDKN1B locus on a discovery case series of 2101 neuroblastoma patients and 4202 genetically matched controls of European ancestry. The most significantly associated rs34330 was analysed in an independent Italian cohort of 311 cases and 709 controls. In vitro functional analysis was carried out in HEK293T and in neuroblastoma cell line SHEP-2, both transfected with pGL3-CDKN1B-CC or pGL3-CDKN1B-TT constructs. We identified an association of the rs34330 T allele (-79C/T) with the neuroblastoma risk (Pcombined = 0.002; OR = 1.17). The risk allele (T) of this single nucleotide polymorphism led to a lower transcription rate in cells transfected with a luciferase reporter driven by the polymorphic p27Kip1 promoter (P < 0.05). Three independent sets of neuroblastoma tumours carrying -79TT genotype showed a tendency towards lower CDKN1B mRNA levels. Our study shows that a functional variant, associated with a reduced CDKN1B gene transcription, influences neuroblastoma susceptibility.

The HERBY trial was a phase II open-label, randomized, multicenter trial evaluating bevacizumab (BEV) in addition to temozolomide/radiotherapy in patients with newly diagnosed non-brainstem high-grade glioma (HGG) between the ages of 3 and 18 years. We carried out comprehensive molecular analysis integrated with pathology, radiology, and immune profiling. In post-hoc subgroup analysis, hypermutator tumors (mismatch repair deficiency and somatic POLE/POLD1 mutations) and those biologically resembling pleomorphic xanthoastrocytoma ([PXA]-like, driven by BRAF_V600E or NF1 mutation) had significantly more CD8+ tumor-infiltrating lymphocytes, and longer survival with the addition of BEV. Histone H3 subgroups (hemispheric G34R/V and midline K27M) had a worse outcome and were immune cold. Future clinical trials will need to take into account the diversity represented by the term "HGG" in the pediatric population.

Pediatric low-grade gliomas (PLGGs) are frequently associated with activating BRAF gene fusions, such as KIAA1549-BRAF, that aberrantly drive the mitogen activated protein kinase (MAPK) pathway. Although RAF inhibitors (RAFi) have been proven effective in BRAF-V600E mutant tumors, we have previously shown how the KIAA1549-BRAF fusion can be paradoxically activated by RAFi. While newer classes of RAFi, such as PLX8394, have now been shown to inhibit MAPK activation by KIAA1549-BRAF, we sought to identify alternative MAPK pathway targeting strategies using clinically relevant MEK inhibitors (MEKi), along with potential escape mechanisms of acquired resistance to single-agent MAPK pathway therapies. We demonstrate effectiveness of multiple MEKi against diverse BRAF-fusions with novel N-terminal partners, with trametinib being the most potent. However, resistance to MEKi or PLX8394 develops via increased RTK expression causing activation of PI3K/mTOR pathway in BRAF-fusion expressing resistant clones. To circumvent acquired resistance, we show potency of combinatorial targeting with trametinib and everolimus, an mTOR inhibitor (mTORi) against multiple BRAF-fusions. While single-agent mTORi and MEKi PLGG clinical trials are underway, our study provides preclinical rationales for using MEKi and mTORi combinatorial therapy to stave off or prevent emergent drug-resistance in BRAF-fusion driven PLGGs.

We report a comprehensive proteogenomics analysis, including whole-genome sequencing, RNA sequencing, and proteomics and phosphoproteomics profiling, of 218 tumors across 7 histological types of childhood brain cancer: low-grade glioma (n = 93), ependymoma (32), high-grade glioma (25), medulloblastoma (22), ganglioglioma (18), craniopharyngioma (16), and atypical teratoid rhabdoid tumor (12). Proteomics data identify common biological themes that span histological boundaries, suggesting that treatments used for one histological type may be applied effectively to other tumors sharing similar proteomics features. Immune landscape characterization reveals diverse tumor microenvironments across and within diagnoses. Proteomics data further reveal functional effects of somatic mutations and copy number variations (CNVs) not evident in transcriptomics data. Kinase-substrate association and co-expression network analysis identify important biological mechanisms of tumorigenesis. This is the first large-scale proteogenomics analysis across traditional histological boundaries to uncover foundational pediatric brain tumor biology and inform rational treatment selection.

Despite advancements in molecular and histopathologic characterization of pediatric low-grade gliomas (pLGGs), there remains significant phenotypic heterogeneity among tumors with similar categorizations. We hypothesized that an unsupervised machine learning approach based on radiomic features may reveal distinct pLGG imaging subtypes.

A more complete understanding of aberrant oncogenic signaling in neuroblastoma, a malignancy of the developing sympathetic nervous system, is paramount to improving patient outcomes. Recently, we identified LIN28B as an oncogenic driver in high-risk neuroblastoma. Here, we identify the oncogene RAN as a LIN28B target and show regional gain of chromosome 12q24 as an additional somatic alteration resulting in increased RAN expression. We show that LIN28B influences RAN expression by promoting RAN Binding Protein 2 expression and by directly binding RAN mRNA. Further, we demonstrate a convergence of LIN28B and RAN signaling on Aurora kinase A activity. Collectively, these findings demonstrate that LIN28B-RAN-AURKA signaling drives neuroblastoma oncogenesis, suggesting that this pathway may be amenable to therapeutic targeting.

Angiocentric gliomas are pediatric low-grade gliomas (PLGGs) without known recurrent genetic drivers. We performed genomic analysis of new and published data from 249 PLGGs, including 19 angiocentric gliomas. We identified MYB-QKI fusions as a specific and single candidate driver event in angiocentric gliomas. In vitro and in vivo functional studies show that MYB-QKI rearrangements promote tumorigenesis through three mechanisms: MYB activation by truncation, enhancer translocation driving aberrant MYB-QKI expression and hemizygous loss of the tumor suppressor QKI. To our knowledge, this represents the first example of a single driver rearrangement simultaneously transforming cells via three genetic and epigenetic mechanisms in a tumor.

Tumor-associated macrophages (TAMs) play an important role in tumor immunity and comprise of subsets that have distinct phenotype, function, and ontology. Transcriptomic analyses of human medulloblastoma, the most common malignant pediatric brain cancer, showed that medulloblastomas (MBs) with activated sonic hedgehog signaling (SHH-MB) have significantly more TAMs than other MB subtypes. Therefore, we examined MB-associated TAMs by single-cell RNA sequencing of autochthonous murine SHH-MB at steady state and under two distinct treatment modalities: molecular-targeted inhibitor and radiation. Our analyses reveal significant TAM heterogeneity, identify markers of ontologically distinct TAM subsets, and show the impact of brain microenvironment on the differentiation of tumor-infiltrating monocytes. TAM composition undergoes dramatic changes with treatment and differs significantly between molecular-targeted and radiation therapy. We identify an immunosuppressive monocyte-derived TAM subset that emerges with radiation therapy and demonstrate its role in regulating T cell and neutrophil infiltration in MB.

Pancreatic ductal adenocarcinoma (PDAC) is the third leading cause of cancer death in the US. Despite multiple large-scale genetic sequencing studies, identification of predictors of patient survival remains challenging. We performed a comprehensive assessment and integrative analysis of large-scale gene expression datasets, across multiple platforms, to enable discovery of a prognostic gene signature for patient survival in pancreatic cancer. PDAC RNA-Sequencing data from The Cancer Genome Atlas was stratified into Survival+ (>2-year survival) and Survival-(<1-year survival) cohorts (n = 47). Comparisons of RNA expression profiles between survival groups and normal pancreatic tissue expression data from the Gene Expression Omnibus generated an initial PDAC specific prognostic differential expression gene list. The candidate prognostic gene list was then trained on the Australian pancreatic cancer dataset from the ICGC database (n = 103), using iterative sampling based algorithms, to derive a gene signature predictive of patient survival. The gene signature was validated in 2 independent patient cohorts and against existing PDAC subtype classifications. We identified 707 candidate prognostic genes exhibiting differential expression in tumor versus normal tissue. A substantial fraction of these genes was also found to be differentially methylated between survival groups. From the candidate gene list, a 5-gene signature (ADM, ASPM, DCBLD2, E2F7, and KRT6A) was identified. Our signature demonstrated significant power to predict patient survival in two distinct patient cohorts and was independent of AJCC TNM staging. Cross-validation of our gene signature reported a better ROC AUC (≥ 0.8) when compared to existing PDAC survival signatures. Furthermore, validation of our signature through immunohistochemical analysis of patient tumor tissue and existing gene expression subtyping data in PDAC, demonstrated a correlation to the presence of vascular invasion and the aggressive squamous tumor subtype. Assessment of these genes in patient biopsies could help further inform risk-stratification and treatment decisions in pancreatic cancer.

We developed a computational pipeline designed to use RNA sequencing (n = 136) and gene expression profiling (n = 250) data from neuroblastoma tumors to identify cell surface proteins predicted to be highly expressed in MYCN amplified neuroblastomas and with little or no expression in normal human tissues. We then performed ChIP-seq in the MYCN amplified cell lines KELLY, NB-1643, and NGP to identify gene promoters that are occupied by MYCN protein to define the intersection with the differentially-expressed gene list. We initially identified 116 putative immunotherapy targets with predicted transmembrane domains, with the most significant differentially-expressed of these being the calmodulin kinase-like vesicle-associated gene (CAMKV, p = 2 × 10-6). CAMKV encodes a protein that binds calmodulin in the presence of calcium, but lacks the kinase activity of other calmodulin kinase family members. We confirmed that CAMKV is selectively expressed in 7/7 MYCN amplified neuroblastoma cell lines and showed that the transcription of CAMKV is directly controlled by MYCN. From membrane fractionation and immunohistochemistry, we verified that CAMKV is membranous in MYCN amplified neuroblastoma cell lines and patient-derived xenografts. Finally, immunohistochemistry showed that CAMKV is not expressed on normal tissues outside of the central nervous system. Together, these data demonstrate that CAMKV is a differentially-expressed cell surface protein that is transcriptionally regulated by MYCN, making it a candidate for targeting with antibodies or antibody-drug conjugates that do not cross the blood brain barrier.

Pediatric solid and central nervous system tumors are the leading cause of cancer-related death among children. Identifying new targeted therapies necessitates the use of pediatric cancer models that faithfully recapitulate the patient's disease. However, the generation and characterization of pediatric cancer models has significantly lagged behind adult cancers, underscoring the urgent need to develop pediatric-focused cell line resources. Herein, we establish a single-site collection of 261 cell lines, including 224 pediatric cell lines representing 18 distinct extracranial and brain childhood tumor types. We subjected 182 cell lines to multi-omics analyses (DNA sequencing, RNA sequencing, DNA methylation), and in parallel performed pharmacological and genetic CRISPR-Cas9 loss-of-function screens to identify pediatric-specific treatment opportunities and biomarkers. Our work provides insight into specific pathway vulnerabilities in molecularly defined pediatric tumor classes and uncovers biomarker-linked therapeutic opportunities of clinical relevance. Cell line data and resources are provided in an open access portal.

To elucidate the pathogenesis of vein of Galen malformations (VOGMs), the most common and most severe of congenital brain arteriovenous malformations, we performed an integrated analysis of 310 VOGM proband-family exomes and 336,326 human cerebrovasculature single-cell transcriptomes. We found the Ras suppressor p120 RasGAP (RASA1) harbored a genome-wide significant burden of loss-of-function de novo variants (2042.5-fold, p = 4.79 x 10-7). Rare, damaging transmitted variants were enriched in Ephrin receptor-B4 (EPHB4) (17.5-fold, p = 1.22 x 10-5), which cooperates with p120 RasGAP to regulate vascular development. Additional probands had damaging variants in ACVRL1, NOTCH1, ITGB1, and PTPN11. ACVRL1 variants were also identified in a multi-generational VOGM pedigree. Integrative genomic analysis defined developing endothelial cells as a likely spatio-temporal locus of VOGM pathophysiology. Mice expressing a VOGM-specific EPHB4 kinase-domain missense variant (Phe867Leu) exhibited disrupted developmental angiogenesis and impaired hierarchical development of arterial-capillary-venous networks, but only in the presence of a "second-hit" allele. These results illuminate human arterio-venous development and VOGM pathobiology and have implications for patients and their families.

Aberrant activation of the mammalian target of rapamycin complex 1 (mTORC1) is a common molecular event in a variety of pathological settings, including genetic tumor syndromes, cancer, and obesity. However, the cell-intrinsic consequences of mTORC1 activation remain poorly defined. Through a combination of unbiased genomic, metabolomic, and bioinformatic approaches, we demonstrate that mTORC1 activation is sufficient to stimulate specific metabolic pathways, including glycolysis, the oxidative arm of the pentose phosphate pathway, and de novo lipid biosynthesis. This is achieved through the activation of a transcriptional program affecting metabolic gene targets of hypoxia-inducible factor (HIF1alpha) and sterol regulatory element-binding protein (SREBP1 and SREBP2). We find that SREBP1 and 2 promote proliferation downstream of mTORC1, and the activation of these transcription factors is mediated by S6K1. Therefore, in addition to promoting protein synthesis, mTORC1 activates specific bioenergetic and anabolic cellular processes that are likely to contribute to human physiology and disease.