CAR-T cell expansion platforms yield distinct T cell differentiation states.

With investigators looking to expand engineered T cell therapies such as CAR-T to new tumor targets and patient populations, a variety of cell manufacturing platforms have been developed to scale manufacturing capacity using closed and/or automated systems. Such platforms are particularly useful for solid tumor targets, which typically require higher CAR-T cell doses. Although T cell phenotype and function are key attributes that often correlate with therapeutic efficacy, how manufacturing platforms influence the final CAR-T cell product is currently unknown. We compared 4 commonly used T cell manufacturing platforms (CliniMACS Prodigy, Xuri W25 rocking platform, G-Rex gas-permeable bioreactor, static bag culture) using identical media, stimulation, culture length, and donor starting material. Selected CD4+CD8+ cells were transduced with lentiviral vector incorporating a CAR targeting FGFR4, a promising target for pediatric sarcoma. We observed significant differences in overall expansion over the 14-day culture; bag cultures had the highest capacity for expansion while the Prodigy had the lowest (481-fold versus 84-fold, respectively). Strikingly, we also observed considerable differences in the phenotype of the final product, with the Prodigy significantly enriched for CCR7+CD45RA+ naïve/stem central memory (Tn/scm)-like cells at 46% compared to bag and G-Rex with 16% and 13%, respectively. Gene expression analysis also showed that Prodigy CAR-Ts are more naïve, less cytotoxic and less exhausted than bag, G-Rex, and Xuri CAR-Ts, and pointed to differences in cell metabolism that were confirmed via metabolic assays. We hypothesized that dissolved oxygen level, which decreased substantially during the final 3 days of the Prodigy culture, may contribute to the observed differences in T cell phenotype. By culturing bag and G-Rex cultures in 1% O2 from day 5 onward, we could generate >60% Tn/scm-like cells, with longer time in hypoxia correlating with a higher percentage of Tn/scm-like cells. Intriguingly, our results suggest that oxygenation is responsible, at least in part, for observed differences in T cell phenotype among bioreactors and suggest hypoxic culture as a potential strategy prevent T cell differentiation during expansion. Ultimately, our study demonstrates that selection of bioreactor system may have profound effects not only on the capacity for expansion, but also on the differentiation state of the resulting CAR-T cells.

Dentithecamides A-H, Diacylated Zoanthoxanthin Derivatives with PAX3-FOXO1 Inhibitory Activity from the Hydroid Dentitheca habereri.

Chemical investigation of the marine hydroid Dentitheca habereri led to the identification of eight new diacylated zoanthoxanthin alkaloids, named dentithecamides A-H (1-8), along with three previously reported analogues, zoamides B-D (9-11). The structures of compounds 1-11 were elucidated by spectroscopic and spectrometric analyses, including IR, HRESIMS, and NMR experiments, and by comparison with literature data. Compounds 1-11 are the first zoanthoxanthin alkaloids to be reported from a hydroid. Dentithecamides A (1) and B (2) along with zoamides B-D (9-11), which all share a conformationally mobile cycloheptadiene core, inhibited PAX3-FOXO1 regulated transcriptional activity and thus provided a structural framework for the potential development of more potent PAX3-FOXO1 inhibitors.

SRAssembler: Selective Recursive local Assembly of homologous genomic regions.

BACKGROUND: The falling cost of next-generation sequencing technology has allowed deep sequencing across related species and of individuals within species. Whole genome assemblies from these data remain high time- and resource-consuming computational tasks, particularly if best solutions are sought using different assembly strategies and parameter sets. However, in many cases, the underlying research questions are not genome-wide but rather target specific genes or sets of genes. We describe a novel assembly tool, SRAssembler, that efficiently assembles only contigs containing potential homologs of a gene or protein query, thus enabling gene-specific genome studies over large numbers of short read samples. RESULTS: We demonstrate the functionality of SRAssembler with examples largely drawn from plant genomics. The workflow implements a recursive strategy by which relevant reads are successively pulled from the input sets based on overlapping significant matches, resulting in virtual chromosome walking. The typical workflow behavior is illustrated with assembly of simulated reads. Applications to real data show that SRAssembler produces homologous contigs of equivalent quality to whole genome assemblies. Settings can be chosen to not only assemble presumed orthologs but also paralogous gene loci in distinct contigs. A key application is assembly of the same locus in many individuals from population genome data, which provides assessment of structural variation beyond what can be inferred from read mapping to a reference genome alone. SRAssembler can be used on modest computing resources or used in parallel on high performance computing clusters (most easily by invoking a dedicated Singularity image). CONCLUSIONS: SRAssembler offers an efficient tool to complement whole genome assembly software. It can be used to solve gene-specific research questions based on large genomic read samples from multiple sources and would be an expedient choice when whole genome assembly from the reads is either not feasible, too costly, or unnecessary. The program can also aid decision making on the depth of sequencing in an ongoing novel genome sequencing project or with respect to ultimate whole genome assembly strategies.

PdumBase: a transcriptome database and research tool for Platynereis dumerilii and early development of other metazoans.

BACKGROUND: The marine polychaete annelid Platynereis dumerilii has recently emerged as a prominent organism for the study of development, evolution, stem cells, regeneration, marine ecology, chronobiology and neurobiology within metazoans. Its phylogenetic position within the spiralian/ lophotrochozoan clade, the comparatively high conservation of ancestral features in the Platynereis genome, and experimental access to any stage within its life cycle, make Platynereis an important model for elucidating the complex regulatory and functional molecular mechanisms governing early development, later organogenesis, and various features of its larval and adult life. High resolution RNA-seq gene expression data obtained from specific developmental stages can be used to dissect early developmental mechanisms. However, the potential for discovery of these mechanisms relies on tools to search, retrieve, and compare genome-wide information within Platynereis, and across other metazoan taxa. RESULTS: To facilitate exploration and discovery by the broader scientific community, we have developed a web-based, searchable online research tool, PdumBase, featuring the first comprehensive transcriptome database for Platynereis dumerilii during early stages of development (2 h ~ 14 h). Our database also includes additional stages over the P. dumerilii life cycle and provides access to the expression data of 17,213 genes (31,806 transcripts) along with annotation information sourced from Swiss-Prot, Gene Ontology, KEGG pathways, Pfam domains, TmHMM, SingleP, and EggNOG orthology. Expression data for each gene includes the stage, the normalized FPKM, the raw read counts, and information that can be leveraged for statistical analyses of differential gene expression and the construction of genome-wide co-expression networks. In addition, PdumBase offers early stage transcriptome expression data from five further species as a valuable resource for investigators interested in comparing early development in different organisms. To understand conservation of Platynereis gene models and to validate gene annotation, most Platynereis gene models include a comprehensive phylogenetic analysis across 18 species representing diverse metazoan taxa. CONCLUSIONS: PdumBase represents the first online resource for the early developmental transcriptome of Platynereis dumerilii. It serves as a research platform for discovery and exploration of gene expression during early stages, throughout the Platynereis life cycle, and enables comparison to other model organisms. PdumBase is freely available at http://pdumbase.gdcb.iastate.edu .

The asymmetric cell division machinery in the spiral-cleaving egg and embryo of the marine annelid Platynereis dumerilii.

BACKGROUND: Over one third of all animal phyla utilize a mode of early embryogenesis called 'spiral cleavage' to divide the fertilized egg into embryonic cells with different cell fates. This mode is characterized by a series of invariant, stereotypic, asymmetric cell divisions (ACDs) that generates cells of different size and defined position within the early embryo. Astonishingly, very little is known about the underlying molecular machinery to orchestrate these ACDs in spiral-cleaving embryos. Here we identify, for the first time, cohorts of factors that may contribute to early embryonic ACDs in a spiralian embryo. RESULTS: To do so we analyzed stage-specific transcriptome data in eggs and early embryos of the spiralian annelid Platynereis dumerilii for the expression of over 50 candidate genes that are involved in (1) establishing cortical domains such as the partitioning defective (par) genes, (2) directing spindle orientation, (3) conveying polarity cues including crumbs and scribble, and (4) maintaining cell-cell adhesion between embryonic cells. In general, each of these cohorts of genes are co-expressed exhibiting high levels of transcripts in the oocyte and fertilized single-celled embryo, with progressively lower levels at later stages. Interestingly, a small number of key factors within each ACD module show different expression profiles with increased early zygotic expression suggesting distinct regulatory functions. In addition, our analysis discovered several highly co-expressed genes that have been associated with specialized neural cell-cell recognition functions in the nervous system. The high maternal contribution of these 'neural' adhesion complexes indicates novel general adhesion functions during early embryogenesis. CONCLUSIONS: Spiralian embryos are champions of ACD generating embryonic cells of different size with astonishing accuracy. Our results suggest that the molecular machinery for ACD is already stored as maternal transcripts in the oocyte. Thus, the spiralian egg can be viewed as a totipotent yet highly specialized cell that evolved to execute fast and precise ACDs during spiral cleaving stages. Our survey identifies cohorts of factors in P. dumerilii that are candidates for these molecular mechanisms and their regulation, and sets the stage for a functional dissection of ACD in a spiral-cleaving embryo.

Defining the transcriptomic landscape of the developing enteric nervous system and its cellular environment.

BACKGROUND: Motility and the coordination of moving food through the gastrointestinal tract rely on a complex network of neurons known as the enteric nervous system (ENS). Despite its critical function, many of the molecular mechanisms that direct the development of the ENS and the elaboration of neural network connections remain unknown. The goal of this study was to transcriptionally identify molecular pathways and candidate genes that drive specification, differentiation and the neural circuitry of specific neural progenitors, the phox2b expressing ENS cell lineage, during normal enteric nervous system development. Because ENS development is tightly linked to its environment, the transcriptional landscape of the cellular environment of the intestine was also analyzed. RESULTS: Thousands of zebrafish intestines were manually dissected from a transgenic line expressing green fluorescent protein under the phox2b regulatory elements [Tg(phox2b:EGFP) w37 ]. Fluorescence-activated cell sorting was used to separate GFP-positive phox2b expressing ENS progenitor and derivatives from GFP-negative intestinal cells. RNA-seq was performed to obtain accurate, reproducible transcriptional profiles and the unbiased detection of low level transcripts. Analysis revealed genes and pathways that may function in ENS cell determination, genes that may be identifiers of different ENS subtypes, and genes that define the non-neural cellular microenvironment of the ENS. Differential expression analysis between the two cell populations revealed the expected neuronal nature of the phox2b expressing lineage including the enrichment for genes required for neurogenesis and synaptogenesis, and identified many novel genes not previously associated with ENS development. Pathway analysis pointed to a high level of G-protein coupled pathway activation, and identified novel roles for candidate pathways such as the Nogo/Reticulon axon guidance pathway in ENS development. CONCLUSION: We report the comprehensive gene expression profiles of a lineage-specific population of enteric progenitors, their derivatives, and their microenvironment during normal enteric nervous system development. Our results confirm previously implicated genes and pathways required for ENS development, and also identify scores of novel candidate genes and pathways. Thus, our dataset suggests various potential mechanisms that drive ENS development facilitating characterization and discovery of novel therapeutic strategies to improve gastrointestinal disorders.

A transcriptional blueprint for a spiral-cleaving embryo.

BACKGROUND: The spiral cleavage mode of early development is utilized in over one-third of all animal phyla and generates embryonic cells of different size, position, and fate through a conserved set of stereotypic and invariant asymmetric cell divisions. Despite the widespread use of spiral cleavage, regulatory and molecular features for any spiral-cleaving embryo are largely uncharted. To address this gap we use RNA-sequencing on the spiralian model Platynereis dumerilii to capture and quantify the first complete genome-wide transcriptional landscape of early spiral cleavage. RESULTS: RNA-sequencing datasets from seven stages in early Platynereis development, from the zygote to the protrochophore, are described here including the de novo assembly and annotation of ~17,200 Platynereis genes. Depth and quality of the RNA-sequencing datasets allow the identification of the temporal onset and level of transcription for each annotated gene, even if the expression is restricted to a single cell. Over 4000 transcripts are maternally contributed and cleared by the end of the early spiral cleavage phase. Small early waves of zygotic expression are followed by major waves of thousands of genes, demarcating the maternal to zygotic transition shortly after the completion of spiral cleavages in this annelid species. CONCLUSIONS: Our comprehensive stage-specific transcriptional analysis of early embryonic stages in Platynereis elucidates the regulatory genome during early spiral embryogenesis and defines the maternal to zygotic transition in Platynereis embryos. This transcriptome assembly provides the first systems-level view of the transcriptional and regulatory landscape for a spiral-cleaving embryo.

Lineage specific transcription factor waves reprogram neuroblastoma from self-renewal to differentiation.

Temporal regulation of super-enhancer (SE) driven transcription factors (TFs) underlies normal developmental programs. Neuroblastoma (NB) arises from an inability of sympathoadrenal progenitors to exit a self-renewal program and terminally differentiate. To identify SEs driving TF regulators, we use all-trans retinoic acid (ATRA) to induce NB growth arrest and differentiation. Time-course H3K27ac ChIP-seq and RNA-seq reveal ATRA coordinated SE waves. SEs that decrease with ATRA link to stem cell development (MYCN, GATA3, SOX11). CRISPR-Cas9 and siRNA verify SOX11 dependency, in vitro and in vivo. Silencing the SOX11 SE using dCAS9-KRAB decreases SOX11 mRNA and inhibits cell growth. Other TFs activate in sequential waves at 2, 4 and 8 days of ATRA treatment that regulate neural development (GATA2 and SOX4). Silencing the gained SOX4 SE using dCAS9-KRAB decreases SOX4 expression and attenuates ATRA-induced differentiation genes. Our study identifies oncogenic lineage drivers of NB self-renewal and TFs critical for implementing a differentiation program.

KDM3B inhibitors disrupt the oncogenic activity of PAX3-FOXO1 in fusion-positive rhabdomyosarcoma.

Fusion-positive rhabdomyosarcoma (FP-RMS) is an aggressive pediatric sarcoma driven primarily by the PAX3-FOXO1 fusion oncogene, for which therapies targeting PAX3-FOXO1 are lacking. Here, we screen 62,643 compounds using an engineered cell line that monitors PAX3-FOXO1 transcriptional activity identifying a hitherto uncharacterized compound, P3FI-63. RNA-seq, ATAC-seq, and docking analyses implicate histone lysine demethylases (KDMs) as its targets. Enzymatic assays confirm the inhibition of multiple KDMs with the highest selectivity for KDM3B. Structural similarity search of P3FI-63 identifies P3FI-90 with improved solubility and potency. Biophysical binding of P3FI-90 to KDM3B is demonstrated using NMR and SPR. P3FI-90 suppresses the growth of FP-RMS in vitro and in vivo through downregulating PAX3-FOXO1 activity, and combined knockdown of KDM3B and KDM1A phenocopies P3FI-90 effects. Thus, we report KDM inhibitors P3FI-63 and P3FI-90 with the highest specificity for KDM3B. Their potent suppression of PAX3-FOXO1 activity indicates a possible therapeutic approach for FP-RMS and other transcriptionally addicted cancers.

Piperacetazine Directly Binds to the PAX3::FOXO1 Fusion Protein and Inhibits Its Transcriptional Activity.

The tumor-specific chromosomal translocation product, PAX3::FOXO1, is an aberrant fusion protein that plays a key role for oncogenesis in the alveolar subtype of rhabdomyosarcoma (RMS). PAX3::FOXO1 represents a validated molecular target for alveolar RMS and successful inhibition of its oncogenic activity is likely to have significant clinical applications. Even though several PAX3::FOXO1 function-based screening studies have been successfully completed, a directly binding small-molecule inhibitor of PAX3::FOXO1 has not been reported. Therefore, we screened small-molecule libraries to identify compounds that were capable of directly binding to PAX3::FOXO1 protein using surface plasmon resonance technology. Compounds that directly bound to PAX3::FOXO1 were further evaluated in secondary transcriptional activation assays. We discovered that piperacetazine can directly bind to PAX3::FOXO1 protein and inhibit fusion protein-derived transcription in multiple alveolar RMS cell lines. Piperacetazine inhibited anchorage-independent growth of fusion-positive alveolar RMS cells but not embryonal RMS cells. On the basis of our findings, piperacetazine is a molecular scaffold upon which derivatives could be developed as specific inhibitors of PAX3::FOXO1. These novel inhibitors could potentially be evaluated in future clinical trials for recurrent or metastatic alveolar RMS as novel targeted therapy options. SIGNIFICANCE: RMS is a malignant soft-tissue tumor mainly affecting the pediatric population. A subgroup of RMS with worse prognosis harbors a unique chromosomal translocation creating an oncogenic fusion protein, PAX3::FOXO1. We identified piperacetazine as a direct inhibitor of PAX3::FOXO1, which may provide a scaffold for designing RMS-specific targeted therapy.

A SNAI2/CTCF Interaction is Required for NOTCH1 Expression in Rhabdomyosarcoma.

Rhabdomyosarcoma (RMS) is a pediatric malignancy of the muscle with characteristics of cells blocked in differentiation. NOTCH1 is an oncogene that promotes self-renewal and blocks differentiation in the fusion negative-RMS sub-type. However, how NOTCH1 expression is transcriptionally maintained in tumors is unknown. Analyses of SNAI2 and CTCF chromatin binding and HiC analyses revealed a conserved SNAI2/CTCF overlapping peak downstream of the NOTCH1 locus marking a sub-topologically associating domain (TAD) boundary. Deletion of the SNAI2-CTCF peak showed that it is essential for NOTCH1 expression and viability of FN-RMS cells. Reintroducing constitutively activated NOTCH1-ΔE in cells with the SNAI2-CTCF peak deleted restored cell-viability. Ablation of SNAI2 using CRISPR/Cas9 reagents resulted in the loss of majority of RD and SMS-CTR FN-RMS cells. However, the few surviving clones that repopulate cultures have recovered NOTCH1. Cells that re-establish NOTCH1 expression after SNAI2 ablation are unable to differentiate robustly as SNAI2 shRNA knockdown cells; yet, SNAI2-ablated cells continued to be exquisitely sensitive to ionizing radiation. Thus, we have uncovered a novel mechanism by which SNAI2 and CTCF maintenance of a sub-TAD boundary promotes rather than represses NOTCH1 expression. Further, we demonstrate that SNAI2 suppression of apoptosis post-radiation is independent of SNAI2/NOTCH1 effects on self-renewal and differentiation.

Predicting Molecular Subtype and Survival of Rhabdomyosarcoma Patients Using Deep Learning of H&E Images: A Report from the Children's Oncology Group.

PURPOSE: Rhabdomyosarcoma (RMS) is an aggressive soft-tissue sarcoma, which primarily occurs in children and young adults. We previously reported specific genomic alterations in RMS, which strongly correlated with survival; however, predicting these mutations or high-risk disease at diagnosis remains a significant challenge. In this study, we utilized convolutional neural networks (CNN) to learn histologic features associated with driver mutations and outcome using hematoxylin and eosin (H&E) images of RMS. EXPERIMENTAL DESIGN: Digital whole slide H&E images were collected from clinically annotated diagnostic tumor samples from 321 patients with RMS enrolled in Children's Oncology Group (COG) trials (1998-2017). Patches were extracted and fed into deep learning CNNs to learn features associated with mutations and relative event-free survival risk. The performance of the trained models was evaluated against independent test sample data (n = 136) or holdout test data. RESULTS: The trained CNN could accurately classify alveolar RMS, a high-risk subtype associated with PAX3/7-FOXO1 fusion genes, with an ROC of 0.85 on an independent test dataset. CNN models trained on mutationally-annotated samples identified tumors with RAS pathway with a ROC of 0.67, and high-risk mutations in MYOD1 or TP53 with a ROC of 0.97 and 0.63, respectively. Remarkably, CNN models were superior in predicting event-free and overall survival compared with current molecular-clinical risk stratification. CONCLUSIONS: This study demonstrates that high-risk features, including those associated with certain mutations, can be readily identified at diagnosis using deep learning. CNNs are a powerful tool for diagnostic and prognostic prediction of rhabdomyosarcoma, which will be tested in prospective COG clinical trials.

A stem cell epigenome is associated with primary nonresponse to CD19 CAR T cells in pediatric acute lymphoblastic leukemia.

CD19 chimeric antigen receptor T-cell therapy (CD19-CAR) has changed the treatment landscape and outcomes for patients with pre-B-cell acute lymphoblastic leukemia (B-ALL). Unfortunately, primary nonresponse (PNR), sustained CD19+ disease, and concurrent expansion of CD19-CAR occur in 20% of the patients and is associated with adverse outcomes. Although some failures may be attributable to CD19 loss, mechanisms of CD19-independent, leukemia-intrinsic resistance to CD19-CAR remain poorly understood. We hypothesize that PNR leukemias are distinct compared with primary sensitive (PS) leukemias and that these differences are present before treatment. We used a multiomic approach to investigate this in 14 patients (7 with PNR and 7 with PS) enrolled in the PLAT-02 trial at Seattle Children's Hospital. Long-read PacBio sequencing helped identify 1 PNR in which 47% of CD19 transcripts had exon 2 skipping, but other samples lacked CD19 transcript abnormalities. Epigenetic profiling discovered DNA hypermethylation at genes targeted by polycomb repressive complex 2 (PRC2) in embryonic stem cells. Similarly, assays of transposase-accessible chromatin-sequencing revealed reduced accessibility at these PRC2 target genes, with a gain in accessibility of regions characteristic of hematopoietic stem cells and multilineage progenitors in PNR. Single-cell RNA sequencing and cytometry by time of flight analyses identified leukemic subpopulations expressing multilineage markers and decreased antigen presentation in PNR. We thus describe the association of a stem cell epigenome with primary resistance to CD19-CAR therapy. Future trials incorporating these biomarkers, with the addition of multispecific CAR T cells targeting against leukemic stem cell or myeloid antigens, and/or combined epigenetic therapy to disrupt this distinct stem cell epigenome may improve outcomes of patients with B-ALL.

Preclinical development of a chimeric antigen receptor T cell therapy targeting FGFR4 in rhabdomyosarcoma.

Pediatric patients with relapsed or refractory rhabdomyosarcoma (RMS) have dismal cure rates, and effective therapy is urgently needed. The oncogenic receptor tyrosine kinase fibroblast growth factor receptor 4 (FGFR4) is highly expressed in RMS and lowly expressed in healthy tissues. Here, we describe a second-generation FGFR4-targeting chimeric antigen receptor (CAR), based on an anti-human FGFR4-specific murine monoclonal antibody 3A11, as an adoptive T cell treatment for RMS. The 3A11 CAR T cells induced robust cytokine production and cytotoxicity against RMS cell lines in vitro. In contrast, a panel of healthy human primary cells failed to activate 3A11 CAR T cells, confirming the selectivity of 3A11 CAR T cells against tumors with high FGFR4 expression. Finally, we demonstrate that 3A11 CAR T cells are persistent in vivo and can effectively eliminate RMS tumors in two metastatic and two orthotopic models. Therefore, our study credentials CAR T cell therapy targeting FGFR4 to treat patients with RMS.

An optimized bicistronic chimeric antigen receptor against GPC2 or CD276 overcomes heterogeneous expression in neuroblastoma.

Chimeric antigen receptor (CAR) T cell therapies targeting single antigens have performed poorly in clinical trials for solid tumors due to heterogenous expression of tumor-associated antigens (TAAs), limited T cell persistence, and T cell exhaustion. Here, we aimed to identify optimal CARs against glypican 2 (GPC2) or CD276 (B7-H3), which were highly but heterogeneously expressed in neuroblastoma (NB), a lethal extracranial solid tumor of childhood. First, we examined CAR T cell expansion in the presence of targets by digital droplet PCR. Next, using pooled competitive optimization of CAR by cellular indexing of transcriptomes and epitopes by sequencing (CITE-Seq), termed P-COCC, we simultaneously analyzed protein and transcriptome expression of CAR T cells to identify high-activity CARs. Finally, we performed cytotoxicity assays to identify the most effective CAR against each target and combined the CARs into a bicistronic "OR" CAR (BiCisCAR). BiCisCAR T cells effectively eliminated tumor cells expressing GPC2 or CD276. Furthermore, the BiCisCAR T cells demonstrated prolonged persistence and resistance to exhaustion when compared with CARs targeting a single antigen. This study illustrated that targeting multiple TAAs with BiCisCAR may overcome heterogenous expression of target antigens in solid tumors and identified a potent, clinically relevant CAR against NB. Moreover, our multimodal approach integrating competitive expansion, P-COCC, and cytotoxicity assays is an effective strategy to identify potent CARs among a pool of candidates.

Interaction between SNAI2 and MYOD enhances oncogenesis and suppresses differentiation in Fusion Negative Rhabdomyosarcoma.

Rhabdomyosarcoma (RMS) is an aggressive pediatric malignancy of the muscle, that includes Fusion Positive (FP)-RMS harboring PAX3/7-FOXO1 and Fusion Negative (FN)-RMS commonly with RAS pathway mutations. RMS express myogenic master transcription factors MYOD and MYOG yet are unable to terminally differentiate. Here, we report that SNAI2 is highly expressed in FN-RMS, is oncogenic, blocks myogenic differentiation, and promotes growth. MYOD activates SNAI2 transcription via super enhancers with striped 3D contact architecture. Genome wide chromatin binding analysis demonstrates that SNAI2 preferentially binds enhancer elements and competes with MYOD at a subset of myogenic enhancers required for terminal differentiation. SNAI2 also suppresses expression of a muscle differentiation program modulated by MYOG, MEF2, and CDKN1A. Further, RAS/MEK-signaling modulates SNAI2 levels and binding to chromatin, suggesting that the differentiation blockade by oncogenic RAS is mediated in part by SNAI2. Thus, an interplay between SNAI2, MYOD, and RAS prevents myogenic differentiation and promotes tumorigenesis.

BAF complexes drive proliferation and block myogenic differentiation in fusion-positive rhabdomyosarcoma.

Rhabdomyosarcoma (RMS) is a pediatric malignancy of skeletal muscle lineage. The aggressive alveolar subtype is characterized by t(2;13) or t(1;13) translocations encoding for PAX3- or PAX7-FOXO1 chimeric transcription factors, respectively, and are referred to as fusion positive RMS (FP-RMS). The fusion gene alters the myogenic program and maintains the proliferative state while blocking terminal differentiation. Here, we investigated the contributions of chromatin regulatory complexes to FP-RMS tumor maintenance. We define the mSWI/SNF functional repertoire in FP-RMS. We find that SMARCA4 (encoding BRG1) is overexpressed in this malignancy compared to skeletal muscle and is essential for cell proliferation. Proteomic studies suggest proximity between PAX3-FOXO1 and BAF complexes, which is further supported by genome-wide binding profiles revealing enhancer colocalization of BAF with core regulatory transcription factors. Further, mSWI/SNF complexes localize to sites of de novo histone acetylation. Phenotypically, interference with mSWI/SNF complex function induces transcriptional activation of the skeletal muscle differentiation program associated with MYCN enhancer invasion at myogenic target genes, which is recapitulated by BRG1 targeting compounds. We conclude that inhibition of BRG1 overcomes the differentiation blockade of FP-RMS cells and may provide a therapeutic strategy for this lethal childhood tumor.

Tumor Mutation Burden, Expressed Neoantigens and the Immune Microenvironment in Diffuse Gliomas.

BACKGROUND: A consistent correlation between tumor mutation burden (TMB) and tumor immune microenvironment has not been observed in gliomas as in other cancers. METHODS: Driver germline and somatic mutations, TMB, neoantigen, and immune cell signatures were analyzed using whole exome sequencing (WES) and transcriptome sequencing of tumor and WES of matched germline DNA in a cohort of 66 glioma samples (44 IDH-mutant and 22 IDH-wildtype). RESULTS: Fourteen samples revealed a hypermutator phenotype (HMP). Eight pathogenic (P) or likely pathogenic (LP) germline variants were detected in 9 (19%) patients. Six of these 8 genes were DNA damage repair genes. P/LP germline variants were found in 22% of IDH-mutant gliomas and 12.5% of IDH-wildtype gliomas (p = 0.7). TMB was correlated with expressed neoantigen but showed an inverse correlation with immune score (R = -0.46, p = 0.03) in IDH-wildtype tumors and no correlation in IDH-mutant tumors. The Antigen Processing and Presentation (APP) score correlated with immune score and was surprisingly higher in NHMP versus HMP samples in IDH-wildtype gliomas, but higher in HMP versus NHMP in IDH-mutant gliomas. CONCLUSION: TMB was inversely correlated with immune score in IDH-wildtype gliomas and showed no correlation in IDH-mutant tumors. APP was correlated with immune score and may be further investigated as a biomarker for response to immunotherapy in gliomas. Studies of germline variants in a larger glioma cohort are warranted.

Genomic Classification and Clinical Outcome in Rhabdomyosarcoma: A Report From an International Consortium.

PURPOSE: Rhabdomyosarcoma is the most common soft tissue sarcoma of childhood. Despite aggressive therapy, the 5-year survival rate for patients with metastatic or recurrent disease remains poor, and beyond PAX-FOXO1 fusion status, no genomic markers are available for risk stratification. We present an international consortium study designed to determine the incidence of driver mutations and their association with clinical outcome. PATIENTS AND METHODS: Tumor samples collected from patients enrolled on Children's Oncology Group trials (1998-2017) and UK patients enrolled on malignant mesenchymal tumor and RMS2005 (1995-2016) trials were subjected to custom-capture sequencing. Mutations, indels, gene deletions, and amplifications were identified, and survival analysis was performed. RESULTS: DNA from 641 patients was suitable for analyses. A median of one mutation was found per tumor. In FOXO1 fusion-negative cases, mutation of any RAS pathway member was found in > 50% of cases, and 21% had no putative driver mutation identified. BCOR (15%), NF1 (15%), and TP53 (13%) mutations were found at a higher incidence than previously reported and TP53 mutations were associated with worse outcomes in both fusion-negative and FOXO1 fusion-positive cases. Interestingly, mutations in RAS isoforms predominated in infants < 1 year (64% of cases). Mutation of MYOD1 was associated with histologic patterns beyond those previously described, older age, head and neck primary site, and a dismal survival. Finally, we provide a searchable companion database (ClinOmics), containing all genomic variants, and clinical annotation including survival data. CONCLUSION: This is the largest genomic characterization of clinically annotated rhabdomyosarcoma tumors to date and provides prognostic genetic features that refine risk stratification and will be incorporated into prospective trials.

Genomic and Transcriptomic Analysis of Relapsed and Refractory Childhood Solid Tumors Reveals a Diverse Molecular Landscape and Mechanisms of Immune Evasion.

Children with treatment-refractory or relapsed (R/R) tumors face poor prognoses. As the genomic underpinnings driving R/R disease are not well defined, we describe here the genomic and transcriptomic landscapes of R/R solid tumors from 202 patients enrolled in Beat Childhood Cancer Consortium clinical trials. Tumor mutational burden (TMB) was elevated relative to untreated tumors at diagnosis, with one-third of tumors classified as having a pediatric high TMB. Prior chemotherapy exposure influenced the mutational landscape of these R/R tumors, with more than 40% of tumors demonstrating mutational signatures associated with platinum or temozolomide chemotherapy and two tumors showing treatment-associated hypermutation. Immunogenomic profiling found a heterogenous pattern of neoantigen and MHC class I expression and a general absence of immune infiltration. Transcriptional analysis and functional gene set enrichment analysis identified cross-pathology clusters associated with development, immune signaling, and cellular signaling pathways. While the landscapes of these R/R tumors reflected those of their corresponding untreated tumors at diagnosis, important exceptions were observed, suggestive of tumor evolution, treatment resistance mechanisms, and mutagenic etiologies of treatment. SIGNIFICANCE: Tumor heterogeneity, chemotherapy exposure, and tumor evolution contribute to the molecular profiles and increased mutational burden that occur in treatment-refractory and relapsed childhood solid tumors.