Imaging characterization of paediatric tumours with the neurotrophic tyrosine receptor kinase fusion transcript

OBJECTIVES

The neurotrophic tyrosine receptor kinase (NTRK) fusion transcript (FT) is a major genetic landmark of infantile fibrosarcoma (IFS) and cellular congenital mesoblastic nephroma (cCMN) but is also described in other tumours. The recent availability of NTRK-targeted drugs enhances the need for better identification. We aimed to describe the anatomic locations and imaging features of tumours with NTRK-FT in children.

CASE SERIES

Imaging characteristics of NTRK-FT tumours of 41 children (median age: 4 months; 63% <1 year old; range: 0-188) managed between 2001 and 2019 were retrospectively analysed. The tumours were located in the soft tissues (n = 24, including 19 IFS), kidneys (n = 9, including 8 cCMN), central nervous system (CNS) (n = 5), lung (n = 2), and bone (n = 1). The tumours were frequently deep-located (93%) and heterogeneous (71%) with necrotic (53%) or haemorrhagic components (29%). Although inconstant, enlarged intratumoural vessels were a recurrent finding (70%) with an irregular distribution (63%) in the most frequent anatomical locations.

CONCLUSION

Paediatric NTRK-FT tumours mainly occur in infants with very variable histotypes and locations. Rich and irregular intra-tumoural vascularization are recurrent findings.

ADVANCES IN KNOWLEDGE

Apart from IFS of soft tissues and cCMN of the kidneys, others NTRK-FT tumours locations have to be known, as CNS tumours. Better knowledge of the imaging characteristics may help guide the pathological and biological identification.

Oncogenic ETS fusions promote DNA damage and proinflammatory responses via pericentromeric RNAs in extracellular vesicles

Aberrant expression of the E26 transformation-specific (ETS) transcription factors characterizes numerous human malignancies. Many of these proteins, including EWS:FLI1 and EWS:ERG fusions in Ewing sarcoma (EwS) and TMPRSS2:ERG in prostate cancer (PCa), drive oncogenic programs via binding to GGAA repeats. We report here that both EWS:FLI1 and ERG bind and transcriptionally activate GGAA-rich pericentromeric heterochromatin. The respective pathogen-like HSAT2 and HSAT3 RNAs, together with LINE, SINE, ERV, and other repeat transcripts, are expressed in EwS and PCa tumors, secreted in extracellular vesicles (EVs), and are highly elevated in plasma of patients with EwS with metastatic disease. High human satellite 2 and 3 (HSAT2,3) levels in EWS:FLI1- or ERG-expressing cells and tumors were associated with induction of G2/M checkpoint, mitotic spindle, and DNA damage programs. These programs were also activated in EwS EV-treated fibroblasts, coincident with accumulation of HSAT2,3 RNAs, proinflammatory responses, mitotic defects, and senescence. Mechanistically, HSAT2,3-enriched cancer EVs induced cGAS-TBK1 innate immune signaling and formation of cytosolic granules positive for double-strand RNAs, RNA-DNA, and cGAS. Hence, aberrantly expressed ETS proteins derepress pericentromeric heterochromatin, yielding pathogenic RNAs that transmit genotoxic stress and inflammation to local and distant sites. Monitoring HSAT2,3 plasma levels and preventing their dissemination may thus improve therapeutic strategies and blood-based diagnostics.

The pelvic organs receive no parasympathetic innervation

The pelvic organs (bladder, rectum, and sex organs) have been represented for a century as receiving autonomic innervation from two pathways - lumbar sympathetic and sacral parasympathetic - by way of a shared relay, the pelvic ganglion, conceived as an assemblage of sympathetic and parasympathetic neurons. Using single-cell RNA sequencing, we find that the mouse pelvic ganglion is made of four classes of neurons, distinct from both sympathetic and parasympathetic ones, albeit with a kinship to the former, but not the latter, through a complex genetic signature. We also show that spinal lumbar preganglionic neurons synapse in the pelvic ganglion onto equal numbers of noradrenergic and cholinergic cells, both of which therefore serve as sympathetic relays. Thus, the pelvic viscera receive no innervation from parasympathetic or typical sympathetic neurons, but instead from a divergent tail end of the sympathetic chains, in charge of its idiosyncratic functions.

Evaluation of the role of AXL in fusion-positive pediatric rhabdomyosarcoma identifies the small-molecule inhibitor bemcentinib (BGB324) as potent chemosensitizer

Rhabdomyosarcoma (RMS) is a highly aggressive pediatric cancer with features of skeletal muscle differentiation. Over 80% of the high-risk patients ultimately fail to respond to chemotherapy treatment, leading to limited therapeutic options and dismal prognostic rates. The lack of response - and subsequent tumor recurrence - is driven in part by stem cell-like cells, the tumor subpopulation that is enriched after treatment, and characterized by expression of the AXL receptor tyrosine kinase (AXL). AXL mediates survival, migration and therapy resistance in several cancer types; however, its function in RMS remains unclear. In this study, we investigated the role of AXL in RMS tumorigenesis, migration and chemotherapy response, and whether targeting of AXL with small molecule inhibitors could potentiate the efficacy of chemotherapy. We show that AXL is expressed in a heterogeneous manner in patient-derived xenografts (PDXs), primary cultures and cell line models of RMS, consistent with its stem-cell-state selectivity. By generating a CRISPR/Cas9 AXL knock-out and overexpressing models, we show that AXL contributes to the migratory phenotype of RMS, but not to chemotherapy resistance. Instead, pharmacological blockade with the AXL inhibitors bemcentinib (BGB324), cabozantinib and NPS-1034 rapidly killed RMS cells in an AXL-independent manner, and augmented the efficacy of the chemotherapeutics vincristine and cyclophosphamide. In vivo administration of the combination of bemcentinib and vincristine exerted strong anti-tumoral activity in a rapidly progressing PDX mouse model, significantly reducing tumor bruden compared to single-agent treatment. Collectively, our data identify bemcentinib as a promising drug to improve chemotherapy efficacy in RMS patients.

Fusion-negative rhabdomyosarcoma 3D organoids to predict effective drug combinations: A proof-of-concept on cell death inducers

Rhabdomyosarcoma (RMS) is the main form of pediatric soft-tissue sarcoma. Its cure rate has not notably improved in the last 20 years following relapse, and the lack of reliable preclinical models has hampered the design of new therapies. This is particularly true for highly heterogeneous fusion-negative RMS (FNRMS). Although methods have been proposed to establish FNRMS organoids, their efficiency remains limited to date, both in terms of derivation rate and ability to accurately mimic the original tumor. Here, we present the development of a next-generation 3D organoid model derived from relapsed adult and pediatric FNRMS. This model preserves the molecular features of the patients' tumors and is expandable for several months in 3D, reinforcing its interest to drug combination screening with longitudinal efficacy monitoring. As a proof-of-concept, we demonstrate its preclinical relevance by reevaluating the therapeutic opportunities of targeting apoptosis in FNRMS from a streamlined approach based on transcriptomic data exploitation.

PAX3-FOXO1 uses its activation domain to recruit CBP/P300 and shape RNA Pol2 cluster distribution

Activation of oncogenic gene expression from long-range enhancers is initiated by the assembly of DNA-binding transcription factors (TF), leading to recruitment of co-activators such as CBP/p300 to modify the local genomic context and facilitate RNA-Polymerase 2 (Pol2) binding. Yet, most TF-to-coactivator recruitment relationships remain unmapped. Here, studying the oncogenic fusion TF PAX3-FOXO1 (P3F) from alveolar rhabdomyosarcoma (aRMS), we show that a single cysteine in the activation domain (AD) of P3F is important for a small alpha helical coil that recruits CBP/p300 to chromatin. P3F driven transcription requires both this single cysteine and CBP/p300. Mutants of the cysteine reduce aRMS cell proliferation and induce cellular differentiation. Furthermore, we discover a profound dependence on CBP/p300 for clustering of Pol2 loops that connect P3F to its target genes. In the absence of CBP/p300, Pol2 long range enhancer loops collapse, Pol2 accumulates in CpG islands and fails to exit the gene body. These results reveal a potential novel axis for therapeutic interference with P3F in aRMS and clarify the molecular relationship of P3F and CBP/p300 in sustaining active Pol2 clusters essential for oncogenic transcription.

Simple droplet microfluidics platform for drug screening on cancer spheroids

3D in vitro biological systems are progressively replacing 2D systems to increase the physiological relevance of cellular studies. Microfluidics-based approaches can be powerful tools towards such biomimetic systems, but often require high-end complicated and expensive processes and equipment for microfabrication. Herein, a drug screening platform is proposed, minimizing technicality and manufacturing steps. It provides an alternate way of spheroid generation in droplets in tubes. Droplet microfluidics then elicit multiple droplets merging events at programmable times, to submit sequentially the spheroids to chemotherapy and to reagents for cytotoxicity screening. After a comprehensive study of tumorogenesis within the droplets, the system is validated for drug screening (IC50) with chemotherapies in cancer cell lines as well as cells from a patient-derived-xenografts (PDX). As compared to microtiter plates methods, our system reduces the initial number of cells up to 10 times and opens new avenues towards primary tumors drug screening approaches.

Adaptive nanopore sequencing to determine pathogenicity of BRCA1 exonic duplication

BRCA1 and BRCA2 are tumour suppressor genes that have been characterised as predisposition genes for the development of hereditary breast and ovarian cancers among other malignancies. The molecular diagnosis of this predisposition syndrome is based on the detection of inactivating variants of any type in those genes. But in the case of structural variants, functional consequences can be difficult to assess using standard molecular methods, as the precise resolution of their sequence is often impossible with short-read next generation sequencing techniques. It has been recently demonstrated that Oxford Nanopore long-read sequencing technology can accurately and rapidly provide genetic diagnoses of Mendelian diseases, including those linked to pathogenic structural variants. Here, we report the accurate resolution of a germline duplication event of exons 18-20 of BRCA1 using Nanopore sequencing with adaptive sampling target enrichment. This allowed us to classify this variant as pathogenic within a short timeframe of 10 days. This study provides a proof-of-concept that nanopore adaptive sampling is a highly efficient technique for the investigation of structural variants of tumour suppressor genes in a clinical context.

Imaging and multi-omics datasets converge to define different neural progenitor origins for ATRT-SHH subgroups

Atypical teratoid rhabdoid tumors (ATRT) are divided into MYC, TYR and SHH subgroups, suggesting diverse lineages of origin. Here, we investigate the imaging of human ATRT at diagnosis and the precise anatomic origin of brain tumors in the Rosa26-CreERT2::Smarcb1flox/flox model. This cross-species analysis points to an extra-cerebral origin for MYC tumors. Additionally, we clearly distinguish SHH ATRT emerging from the cerebellar anterior lobe (CAL) from those emerging from the basal ganglia (BG) and intra-ventricular (IV) regions. Molecular characteristics point to the midbrain-hindbrain boundary as the origin of CAL SHH ATRT, and to the ganglionic eminence as the origin of BG/IV SHH ATRT. Single-cell RNA sequencing on SHH ATRT supports these hypotheses. Trajectory analyses suggest that SMARCB1 loss induces a de-differentiation process mediated by repressors of the neuronal program such as REST, ID and the NOTCH pathway.

Sequential genomic analysis using a multisample/multiplatform approach to better define rhabdomyosarcoma progression and relapse

The genomic spectrum of rhabdomyosarcoma (RMS) progression from primary to relapse is not fully understood. In this pilot study, we explore the sensitivity of various targeted and whole-genome NGS platforms in order to assess the best genomic approach of using liquid biopsy in future prospective clinical trials. Moreover, we investigate 35 paired primary/relapsed RMS from two contributing institutions, 18 fusion-positive (FP-RMS) and 17 fusion-negative RMS (FN-RMS) by either targeted DNA or whole exome sequencing (WES). In 10 cases, circulating tumor DNA (ctDNA) from multiple timepoints through clinical care and progression was analyzed for feasibility of liquid biopsy in monitoring treatment response/relapse. ctDNA alterations were evaluated using a targeted 36-gene custom RMS panel at high coverage for single-nucleotide variation and fusion detection, and a shallow whole-genome sequencing for copy number variation. FP-RMS have a stable genome with relapse, with common secondary alterations CDKN2A/B, MYCN, and CDK4 present at diagnosis and impacting survival. FP-RMS lacking major secondary events at baseline acquire recurrent MYCN and AKT1 alterations. FN-RMS acquire a higher number of new alterations, most commonly SMARCA2 missense mutations. ctDNA analyses detect pathognomonic variants in all RMS patients within our collection at diagnosis, regardless of type of alterations, and confirmed at relapse in 86% of FP-RMS and 100% FN-RMS. Moreover, a higher number of fusion reads is detected with increased disease burden and at relapse in patients following a fatal outcome. These results underscore patterns of tumor progression and provide rationale for using liquid biopsy to monitor treatment response.

A biobank of pediatric patient-derived-xenograft models in cancer precision medicine trial MAPPYACTS for relapsed and refractory tumors

Pediatric patients with recurrent and refractory cancers are in most need for new treatments. This study developed patient-derived-xenograft (PDX) models within the European MAPPYACTS cancer precision medicine trial (NCT02613962). To date, 131 PDX models were established following heterotopical and/or orthotopical implantation in immunocompromised mice: 76 sarcomas, 25 other solid tumors, 12 central nervous system tumors, 15 acute leukemias, and 3 lymphomas. PDX establishment rate was 43%. Histology, whole exome and RNA sequencing revealed a high concordance with the primary patient's tumor profile, human leukocyte-antigen characteristics and specific metabolic pathway signatures. A detailed patient molecular characterization, including specific mutations prioritized in the clinical molecular tumor boards are provided. Ninety models were shared with the IMI2 ITCC Pediatric Preclinical Proof-of-concept Platform (IMI2 ITCC-P4) for further exploitation. This PDX biobank of unique recurrent childhood cancers provides an essential support for basic and translational research and treatments development in advanced pediatric malignancies.

Transmorph: a unifying computational framework for modular single-cell RNA-seq data integration

Data integration of single-cell RNA-seq (scRNA-seq) data describes the task of embedding datasets gathered from different sources or experiments into a common representation so that cells with similar types or states are embedded close to one another independently from their dataset of origin. Data integration is a crucial step in most scRNA-seq data analysis pipelines involving multiple batches. It improves data visualization, batch effect reduction, clustering, label transfer, and cell type inference. Many data integration tools have been proposed during the last decade, but a surge in the number of these methods has made it difficult to pick one for a given use case. Furthermore, these tools are provided as rigid pieces of software, making it hard to adapt them to various specific scenarios. In order to address both of these issues at once, we introduce the transmorph framework. It allows the user to engineer powerful data integration pipelines and is supported by a rich software ecosystem. We demonstrate transmorph usefulness by solving a variety of practical challenges on scRNA-seq datasets including joint datasets embedding, gene space integration, and transfer of cycle phase annotations. transmorph is provided as an open source python package.

Cell-Free DNA Extracted from CSF for the Molecular Diagnosis of Pediatric Embryonal Brain Tumors

BACKGROUND

Liquid biopsies are revolutionary tools used to detect tumor-specific genetic alterations in body fluids, including the use of cell-free DNA (cfDNA) for molecular diagnosis in cancer patients. In brain tumors, cerebrospinal fluid (CSF) cfDNA might be more informative than plasma cfDNA. Here, we assess the use of CSF cfDNA in pediatric embryonal brain tumors (EBT) for molecular diagnosis.

METHODS

The CSF cfDNA of pediatric patients with medulloblastoma (n = 18), ATRT (n = 3), ETMR (n = 1), CNS NB FOXR2 (n = 2) and pediatric EBT NOS (n = 1) (mean cfDNA concentration 48 ng/mL; range 4-442 ng/mL) and matched tumor genomic DNA were sequenced by WES and/or a targeted sequencing approach to determine single-nucleotide variations (SNVs) and copy number alterations (CNA). A specific capture covering transcription start sites (TSS) of genes of interest was also used for nucleosome footprinting in CSF cfDNA.

RESULTS

15/25 CSF cfDNA samples yielded informative results, with informative CNA and SNVs in 11 and 15 cases, respectively. For cases with paired tumor and CSF cfDNA WES (n = 15), a mean of 83 (range 1-160) shared SNVs were observed, including SNVs in classical medulloblastoma genes such as SMO and KMT2D. Interestingly, tumor-specific SNVs (mean 18; range 1-62) or CSF-specific SNVs (mean 5; range 0-25) were also observed, suggesting clonal heterogeneity. The TSS panel resulted in differential coverage profiles across all 112 studied genes in 7 cases, indicating distinct promoter accessibility.

CONCLUSION

CSF cfDNA sequencing yielded informative results in 60% (15/25) of all cases, with informative results in 83% (15/18) of all cases analyzed by WES. These results pave the way for the implementation of these novel approaches for molecular diagnosis and minimal residual disease monitoring.

Neurotrophic tropomyosin receptor kinase (NTRK) fusion positive tumors: a historical cohort analysis

BACKGROUND

NTRK gene fusions have been identified in various tumors; some requiring aggressive therapy and sometimes new TRK inhibitors (TRKi). We aimed to describe a national, unselected, retrospective, multicenter cohort.

RESEARCH DESIGN AND METHODS

Patients were identified through the French sarcoma diagnostic laboratory at Institut Curie through samples analyzed by RT-qPCR or whole-transcriptome sequencing.

RESULTS

From 2001 to 2019, 65 NTRK fusion tumors were identified within 2120 analyses (3.1%): 58 by RNA sequencing (including 20 after RT-qPCR analysis) and 7 exclusively by RT-qPCR. Of the 61 patients identified, 37 patients had infantile soft tissue or kidney fibrosarcomas (IFS), 15 other mesenchymal (Other-MT) and nine central nervous system (CNS) tumors. They encompassed 14 different tumor types with variable behaviors. Overall, 53 patients had surgery (3 mutilating), 38 chemotherapy (20 alkylating agents/anthracycline), 11 radiotherapy, two 'observation strategy' and 13 received TRKi. After a median follow-up of 61.0 months [range, 2.5-226.0], 10 patients died. Five-year overall survival is, respectively, 91.9% [95%CI, 83.5-100.0], 61.1% [95%CI, 34.2-100.0] and 64.8% [95%CI, 39.3-100.0] for IFS, Other-MT, and CNS groups.

CONCLUSIONS

NTRK-fusion positive tumors are rare but detection is improved through RNA sequencing. TRKi could be considered at diagnosis for CNS NTRK-fusion positive tumors, some IFS, and Other-MT.

TRIAL REGISTRATION

Not adapted.

Impact of age on survival according to molecular tumor findings in children and adolescents with soft-tissue and bone sarcoma: The BIOSCA project

BACKGROUND

Adolescents (15-19 years) with sarcoma are known to have significantly worse survival than children (0-14 years). One possible reason may be that the adolescent sarcomas exhibit specific biological characteristics resulting in differences in clinical presentation and treatment resistance behaviors. The BIOSCA project aims to further explore these age-related differences in survival accounting for molecular tumor characteristic in children and adolescents with sarcoma.

METHODS

A retrospective national population-based observational study with documented somatic genetic analyses was conducted between 2011 and 2016 of all patients aged from 0 to 17 years with a diagnosis of sarcoma using the National Registry of Childhood Cancers Database.

RESULTS

A total of 1637 children (0-9years: 40%), preadolescents (10-14years: 35%) and adolescents (15-17 years: 25%) with a diagnosis of bone (N = 845) or soft-tissue (N = 792) sarcoma were included. Adolescents had significantly worse outcome for undifferentiated small round cell sarcoma (USRCS), alveolar rhabdomyosarcoma (ARMS), and epithelioid sarcoma. Five-year overall survivals were worse among CIC-rearranged USRCS cases (47% [95%CI:21-69]) as compared to other USRCS, and PAX3::FOXO1 ARMS patients (44% [95%CI:32-55]) as compared to other ARMS. Adjusting for stage and genomic-profiling status, adolescents with USRCS were 1.6-fold more likely to die than children (P = 0.05), while the difference in survival between age of ARMS patients was weaken. Indeed, the prevalence of PAX3::FOXO1 increased significantly with age.

CONCLUSION

Age was an independent prognostic factor of outcome only in patients with USRCS, while the association between age and survival of patients with ARMS could be partly explained by differences in prevalence of PAX3::FOXO1.

Refined high-content imaging-based phenotypic drug screening in zebrafish xenografts

Zebrafish xenotransplantation models are increasingly applied for phenotypic drug screening to identify small compounds for precision oncology. Larval zebrafish xenografts offer the opportunity to perform drug screens at high-throughput in a complex in vivo environment. However, the full potential of the larval zebrafish xenograft model has not yet been realized and several steps of the drug screening workflow still await automation to increase throughput. Here, we present a robust workflow for drug screening in zebrafish xenografts using high-content imaging. We established embedding methods for high-content imaging of xenografts in 96-well format over consecutive days. In addition, we provide strategies for automated imaging and analysis of zebrafish xenografts including automated tumor cell detection and tumor size analysis over time. We also compared commonly used injection sites and cell labeling dyes and show specific site requirements for tumor cells from different entities. We demonstrate that our setup allows us to investigate proliferation and response to small compounds in several zebrafish xenografts ranging from pediatric sarcomas and neuroblastoma to glioblastoma and leukemia. This fast and cost-efficient assay enables the quantification of anti-tumor efficacy of small compounds in large cohorts of a vertebrate model system in vivo. Our assay may aid in prioritizing compounds or compound combinations for further preclinical and clinical investigations.

Reversible transitions between noradrenergic and mesenchymal tumor identities define cell plasticity in neuroblastoma

Noradrenergic and mesenchymal identities have been characterized in neuroblastoma cell lines according to their epigenetic landscapes and core regulatory circuitries. However, their relationship and relative contribution in patient tumors remain poorly defined. We now document spontaneous and reversible plasticity between the two identities, associated with epigenetic reprogramming, in several neuroblastoma models. Interestingly, xenografts with cells from each identity eventually harbor a noradrenergic phenotype suggesting that the microenvironment provides a powerful pressure towards this phenotype. Accordingly, such a noradrenergic cell identity is systematically observed in single-cell RNA-seq of 18 tumor biopsies and 15 PDX models. Yet, a subpopulation of these noradrenergic tumor cells presents with mesenchymal features that are shared with plasticity models, indicating that the plasticity described in these models has relevance in neuroblastoma patients. This work therefore emphasizes that intrinsic plasticity properties of neuroblastoma cells are dependent upon external cues of the environment to drive cell identity.

Abstract 3525: VGLL2-NCOA2 leverages developmental programs for pediatric sarcomagenesis

Clinical sequencing efforts are rapidly identifying sarcoma gene fusions that have not been functionally validated for their transformation capacity and biological activity. An example is the new fusion of transcriptional coactivators, VGLL2-NCOA2, found in infantile rhabdomyosarcoma. To delineate VGLL2-NCOA2 tumorigenic mechanisms and identify therapeutic vulnerabilities, we implemented a cross-species comparative oncology approach with zebrafish, mouse allograft, and patient samples. We found that in our transgenic zebrafish and mouse allograft models, VGLL2-NCOA2 is sufficient to generate mesenchymal tumors that display features of immature skeletal muscle. Zebrafish and mouse VGLL2-NCOA2 tumors are consistent with the human disease histologically and transcriptionally, and express diagnostic rhabdomyosarcoma markers. We found a shared molecular feature of these tumors was a suppression of myogenic differentiation pathways and reactivation of developmental programs. To identify impacted developmental programs, we transcriptionally clustered zebrafish VGLL2-NCOA2 tumors with the trajectory of early zebrafish development. A subset of VGLL2-NCOA2 zebrafish tumors cluster with embryonic somitogenesis and pinpoint VGLL2-NCOA2 developmental targets, including a RAS family GTPase, ARF6. In VGLL2-NCOA2 zebrafish, mouse allograft, and patient tumors, ARF6 is highly expressed and is absent from mature skeletal muscle. ARF6 knockout in our C2C12-VGLL2-NCOA2 cell culture model suppresses VGLL2-NCOA2 soft agar colony formation capacity, suggesting a genetic cooperating event in the disease. Our data indicate that VGLL2-NCOA2 is an oncogene which leverages developmental programs for tumorigenesis, and that the reactivation or persistence of ARF6 could represent a therapeutic opportunity.

Citation Format: Sarah Watson, Collette A. LaVigne, Lin Xu, Didier Surdez, Joanna Cyrta, Delia Calderon, Matthew V. Cannon, Matthew R. Kent, Katherine M. Silvius, Jack P. Kucinski, Emma N. Harrison, Whitney Murchison, Dinesh Rakheja, Franck Tirode, Olivier Delattre, James F. Amatruda, Genevieve C. Kendall. VGLL2-NCOA2 leverages developmental programs for pediatric sarcomagenesis. [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 3525.

Abstract 96: Integrated molecular analysis of human dedifferentiated liposarcoma identifies a population of liposarcoma progenitors vulnerable to TGF beta inhibition

Introduction: Dedifferentiated liposarcomas (DDLPS) are the most frequent high-grade soft tissue sarcoma in adults. From a pathological point of view, these tumors are composed of high-grade undifferentiated tumor cells (DD), often showing an abrupt transition from a compartment of well-differentiated adipocytic tumor cells (WD). Tumor cells from both WD and DD compartments are characterized by recurrent MDM2 amplification, but their cellular origin and the molecular mechanisms associated with dedifferentiation are poorly understood.

Methods: We performed an integrated molecular analysis of tumors collected from 11 patients undergoing surgery for primary untreated DDLPS. DDLPS tumors were analyzed by single-cell RNA sequencing (scRNAseq) and bulk RNA sequencing on paired WD and DD samples from the same tumors. Results were validated in vitro and in vivo in an additional cohort of human tumors, patient-derived xenografts and DDLPS cell lines.

Results: Through RNA-sequencing of 102,753 individual cells from 11 primary DDLPS lesions, major cell clusters were identified based on unsupervised clustering of gene expression profiles and canonical markers. They include 31 tumor microenvironment clusters and 11 tumor cell clusters. A cluster of tumor cells from the WD compartment is characterized by signatures of early adipocytic progenitors, previously identified as TGFβ-dependent, DPP4-positive stromal progenitors. We show that these cells harbor specifically the truncal genomic alterations of the cancer, with further subclonal mutations identifiable in both WD and DD compartments of DDLPS. Furthermore, these cells have multipotent properties and their differentiation towards the adipocytic lineage is inhibited by TGFβ. Treatment of DD tumor cells with TGFβ inhibitors restores their adipocytic phenotype in vitro and in vivo.

Conclusion: We provide the first single-cell atlas of human DDLPS tumor and microenvironment and identify a population of adipocytic tumor progenitors at the origin of both WD and DD compartments. This study provides rationale for the development of therapeutic strategies based on TGFβ inhibition in advanced DDLPS.

Citation Format: Nadège Gruel, Chloé Quignot, Julien Vibert, Sylvie Bonvalot, Sophie El Zein, Dimitri Tzanis, Sylvain Baulande, Odette Mariani, Joshua Waterfall, Olivier Delattre, Sarah Watson. Integrated molecular analysis of human dedifferentiated liposarcoma identifies a population of liposarcoma progenitors vulnerable to TGF beta inhibition [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 96.

PAX-FOXO1 fusion status in children and adolescents with alveolar rhabdomyosarcoma: Impact on clinical, pathological, and survival features

BACKGROUND

Alveolar rhabdomyosarcoma (ARMS) is an aggressive pediatric cancer and cases with fusion PAX3-FOXO1 and PAX7-FOXO1 seem to have a poor prognosis. The aim is to evaluate whether PAX-FOXO1 alterations influence clinical outcome in childhood and adolescence population with ARMS.

PROCEDURE

A population-based study was conducted between 2011 and 2016 in patients less than 17 years with a diagnosis of ARMS. Overall survival (OS) depending on fusion status with clinical factors was analyzed.

RESULTS

Out of 111 ARMS patients recorded in the French National Childhood Cancer Registry during the 2011-2016 period, 61% expressed PAX3-FOXO1, 15% expressed PAX7-FOXO1, 13% were FOXO1 fusion-positive without PAX specification, and 7% were PAX-FOXO1 negative (n = 4 missing data). Compared to patients with PAX7-FOXO1 positive ARMS, those with PAX3-FOXO1 positive tumor were significantly older (10-17 years: 57.4% vs. 29.4%), and had more often a metastatic disease (54.4% vs. 23.5%). Poorer 5-year OS for patients with PAX3-FOXO1 and PAX not specified FOXO1-positive tumor were observed (44.0% [32.0-55.4] and 35.7% [13.1-59.4], respectively). After adjustment for stage at diagnosis, patients with positive tumor for PAX3-FOXO1 were 3.6-fold more likely to die than those with positive tumor for PAX7-FOXO1.

CONCLUSION

At the population level, PAX3-FOXO1 was associated with a significant higher risk of death compared to PAX7-FOXO1-positive and PAX-FOXO1-negative tumors, and could explain poorer 5-year OS observed in adolescence population diagnosed with ARMS. A continuous risk score derived from the combination of clinical parameters with PAX3-FOXO1 fusion status represents a robust approach to improving current risk-adapted therapy for ARMS.

Targeted long-read sequencing of the Ewing sarcoma 6p25.1 susceptibility locus identifies germline-somatic interactions with EWSR1-FLI1 binding

Ewing sarcoma (EwS) is a rare bone and soft tissue malignancy driven by chromosomal translocations encoding chimeric transcription factors, such as EWSR1-FLI1, that bind GGAA motifs forming novel enhancers that alter nearby expression. We propose that germline microsatellite variation at the 6p25.1 EwS susceptibility locus could impact downstream gene expression and EwS biology. We performed targeted long-read sequencing of EwS blood DNA to characterize variation and genomic features important for EWSR1-FLI1 binding. We identified 50 microsatellite alleles at 6p25.1 and observed that EwS-affected individuals had longer alleles (>135 bp) with more GGAA repeats. The 6p25.1 GGAA microsatellite showed chromatin features of an EWSR1-FLI1 enhancer and regulated expression of RREB1, a transcription factor associated with RAS/MAPK signaling. RREB1 knockdown reduced proliferation and clonogenic potential and reduced expression of cell cycle and DNA replication genes. Our integrative analysis at 6p25.1 details increased binding of longer GGAA microsatellite alleles with acquired EWSR-FLI1 to promote Ewing sarcomagenesis by RREB1-mediated proliferation.

Single-cell profiling of alveolar rhabdomyosarcoma reveals RAS pathway inhibitors as cell-fate hijackers with therapeutic relevance

Rhabdomyosarcoma (RMS) is a group of pediatric cancers with features of developing skeletal muscle. The cellular hierarchy and mechanisms leading to developmental arrest remain elusive. Here, we combined single-cell RNA sequencing, mass cytometry, and high-content imaging to resolve intratumoral heterogeneity of patient-derived primary RMS cultures. We show that the aggressive alveolar RMS (aRMS) subtype contains plastic muscle stem-like cells and cycling progenitors that drive tumor growth, and a subpopulation of differentiated cells that lost its proliferative potential and correlates with better outcomes. While chemotherapy eliminates cycling progenitors, it enriches aRMS for muscle stem-like cells. We screened for drugs hijacking aRMS toward clinically favorable subpopulations and identified a combination of RAF and MEK inhibitors that potently induces myogenic differentiation and inhibits tumor growth. Overall, our work provides insights into the developmental states underlying aRMS aggressiveness, chemoresistance, and progression and identifies the RAS pathway as a promising therapeutic target.

High-content drug screening in zebrafish xenografts reveals high efficacy of dual MCL-1/BCL-XL inhibition against Ewing sarcoma

Ewing sarcoma is a pediatric bone and soft tissue cancer with an urgent need for new therapies to improve disease outcome. To identify effective drugs, phenotypic drug screening has proven to be a powerful method, but achievable throughput in mouse xenografts, the preclinical Ewing sarcoma standard model, is limited. Here, we explored the use of xenografts in zebrafish for high-throughput drug screening to discover new combination therapies for Ewing sarcoma. We subjected xenografts in zebrafish larvae to high-content imaging and subsequent automated tumor size analysis to screen single agents and compound combinations. We identified three drug combinations effective against Ewing sarcoma cells: Irinotecan combined with either an MCL-1 or an BCL-XL inhibitor and in particular dual inhibition of the anti-apoptotic proteins MCL-1 and BCL-XL, which efficiently eradicated tumor cells in zebrafish xenografts. We confirmed enhanced efficacy of dual MCL-1/BCL-XL inhibition compared to single agents in a mouse PDX model. In conclusion, high-content screening of small compounds on Ewing sarcoma zebrafish xenografts identified dual MCL-1/BCL-XL targeting as a specific vulnerability and promising therapeutic strategy for Ewing sarcoma, which warrants further investigation towards clinical application.

Alarm communication predates eusociality in termites

Termites (Blattodea: Isoptera) have evolved specialized defensive strategies for colony protection. Alarm communication enables workers to escape threats while soldiers are recruited to the source of disturbance. Here, we study the vibroacoustic and chemical alarm communication in the wood roach Cryptocercus and in 20 termite species including seven of the nine termite families, all life-types, and all feeding and nesting habits. Our multidisciplinary approach shows that vibratory alarm signals represent an ethological synapomorphy of termites and Cryptocercus. In contrast, chemical alarms have evolved independently in several cockroach groups and at least twice in termites. Vibroacoustic alarm signaling patterns are the most complex in Neoisoptera, in which they are often combined with chemical signals. The alarm characters correlate to phylogenetic position, food type and hardness, foraging area size, and nesting habits. Overall, species of Neoisoptera have developed the most sophisticated communication system amongst termites, potentially contributing to their ecological success.

Combination Therapies Targeting Alk-Aberrant Neuroblastoma in Preclinical Models

BACKGROUND

ALK activating mutations are identified in approximately 10% of newly diagnosed neuroblastomas and ALK amplifications in a further 1-2% of cases. Lorlatinib, a third generation ALK inhibitor, will soon be given alongside induction chemotherapy for children with ALK-aberrant neuroblastoma. However, resistance to single agent treatment has been reported and therapies that improve the response duration are urgently required. We studied the preclinical combination of lorlatinib with chemotherapy, or with the MDM2 inhibitor, idasanutlin, as recent data has suggested that ALK inhibitor resistance can be overcome through activation of the p53-MDM2 pathway.

AIMS

To study the preclinical activity of ALK inhibitors alone and combined with chemotherapy or idasanutlin.

METHODS

We compared different ALK inhibitors in preclinical models prior to evaluating lorlatinib in combination with chemotherapy or idasanutlin. We developed a triple chemotherapy (CAV: cyclophosphamide, doxorubicin and vincristine) in vivo dosing schedule and applied this to both neuroblastoma genetically engineered mouse models (GEMM) and patient derived xenografts (PDX).

RESULTS

Lorlatinib in combination with chemotherapy was synergistic in immunocompetent neuroblastoma GEMM. Significant growth inhibition in response to lorlatinib was only observed in the ALK-amplified PDX model with high ALK expression. In this PDX lorlatinib combined with idasanutlin resulted in complete tumor regression and significantly delayed tumor regrowth.

CONCLUSION

In our preclinical neuroblastoma models, high ALK expression was associated with lorlatinib response alone or in combination with either chemotherapy or idasanutlin. The synergy between MDM2 and ALK inhibition warrants further evaluation of this combination as a potential clinical approach for children with neuroblastoma.

Staging of newly diagnosed Ewing sarcoma: Results of bone marrow aspiration and biopsy versus (18)FDG-PET/CT imaging for bone marrow involvement

BACKGROUND

Ewing sarcoma (ES) is an aggressive bone or extraosseous tumour with an unfavourable prognosis when bone marrow metastases are present at diagnosis. The gold standard diagnosis for bone marrow (BM) involvement is cytological and pathological analysis through bone marrow aspiration and biopsy (BMAB). Several recent studies suggest that these invasive and painful procedures could be replaced by 18F-fluorodeoxyglucose-positron emission tomography/computed tomography ((18)FDG-PET/CT), as this nuclear imaging technique is highly sensitive at detecting bone and extraosseous metastases of ES.

METHODS

In order to study the precision of (18)FDG-PET/CT in the evaluation of bone marrow metastases at diagnosis, we compared the imaging results with cytological/histological analyses performed on BM samples. We retrospectively studied 180 patients with ES recorded at the Léon Bérard Centre over the past 10 years, who were evaluated by (18)FDG-PET/CT and BMAB at diagnosis.

RESULTS

Of the 180 patients, 13 displayed marrow metastases by cytological/histological examination, and only one of these did not have (18)FDG-PET/CT signs of bone marrow involvement, whereas the 167 remaining patients without marrow metastasis all had a negative (18)FDG-PET/CT, except for one. Hence, the sensitivity and specificity of (18)FDG-PET/CT in these patients was 92.3% and 99.4%, respectively. The overall survival at five years of all patients was 67.4% but decrease to 38.5% in the group with bone marrow metastases.

CONCLUSION

Given the results presented herein the bone sarcoma group of the French Sarcoma Group suggests that invasive BMAB no longer be systematically performed for the staging at the diagnosis of ES.

Systematic multi-omics cell line profiling uncovers principles of Ewing sarcoma fusion oncogene-mediated gene regulation

Ewing sarcoma (EwS) is characterized by EWSR1-ETS fusion transcription factors converting polymorphic GGAA microsatellites (mSats) into potent neo-enhancers. Although the paucity of additional mutations makes EwS a genuine model to study principles of cooperation between dominant fusion oncogenes and neo-enhancers, this is impeded by the limited number of well-characterized models. Here we present the Ewing Sarcoma Cell Line Atlas (ESCLA), comprising whole-genome, DNA methylation, transcriptome, proteome, and chromatin immunoprecipitation sequencing (ChIP-seq) data of 18 cell lines with inducible EWSR1-ETS knockdown. The ESCLA shows hundreds of EWSR1-ETS-targets, the nature of EWSR1-ETS-preferred GGAA mSats, and putative indirect modes of EWSR1-ETS-mediated gene regulation, converging in the duality of a specific but plastic EwS signature. We identify heterogeneously regulated EWSR1-ETS-targets as potential prognostic EwS biomarkers. Our freely available ESCLA (http://r2platform.com/escla/) is a rich resource for EwS research and highlights the power of comprehensive datasets to unravel principles of heterogeneous gene regulation by chimeric transcription factors.

Molecular Characterization of Circulating Tumor DNA in Pediatric Rhabdomyosarcoma: A Feasibility Study

PURPOSE

Rhabdomyosarcomas (RMS) are rare neoplasms affecting children and young adults. Efforts to improve patient survival have been undermined by a lack of suitable disease markers. Plasma circulating tumor DNA (ctDNA) has shown promise as a potential minimally invasive biomarker and monitoring tool in other cancers; however, it remains underexplored in RMS. We aimed to determine the feasibility of identifying and quantifying ctDNA in plasma as a marker of disease burden and/or treatment response using blood samples from RMS mouse models and patients.

METHODS

We established mouse models of RMS and applied quantitative polymerase chain reaction (PCR) and droplet digital PCR (ddPCR) to detect ctDNA within the mouse plasma. Potential driver mutations, copy-number alterations, and DNA breakpoints associated with PAX3/7-FOXO1 gene fusions were identified in the RMS samples collected at diagnosis. Patient-matched plasma samples collected from 28 patients with RMS before, during, and after treatment were analyzed for the presence of ctDNA via ddPCR, panel sequencing, and/or whole-exome sequencing.

RESULTS

Human tumor-derived DNA was detectable in plasma samples from mouse models of RMS and correlated with tumor burden. In patients, ctDNA was detected in 14/18 pretreatment plasma samples with ddPCR and 7/7 cases assessed by sequencing. Levels of ctDNA at diagnosis were significantly higher in patients with unfavorable tumor sites, positive nodal status, and metastasis. In patients with serial plasma samples (n = 18), fluctuations in ctDNA levels corresponded to treatment response.

CONCLUSION

Comprehensive ctDNA analysis combining high sensitivity and throughput can identify key molecular drivers in RMS models and patients, suggesting potential as a minimally invasive biomarker. Preclinical assessment of treatments using mouse models and further patient testing through prospective clinical trials are now warranted.

Abstract IA009: Oncogenic transcription factors drive highly tumor-specific LincRNA

Many cancers are characterized by gene fusions encoding oncogenic chimeric transcription factors (TFs) such as EWS::FLI1 in Ewing sarcoma (EwS). EWS::FLI1 induces a complete reprogramming of mesenchymal stem cells, the putative cell of origin of Ewing sarcoma, through binding and subsequent generation of thousands of neoenhancers at GGAA microsatellite repeats. We show that EWS::FLI1 induces the robust expression of a specific set of novel spliced and polyadenylated transcripts within otherwise transcriptionally silent regions of the genome. These neo-genes (NGs) are virtually undetectable in large collections of normal tissues or non-EwS tumors and can be silenced by CRISPR interference at regulatory EWS::FLI1-bound microsatellites. Ribosome profiling and proteomics further show that some NGs are translated into highly EwS-specific peptides. More generally, hundreds of NGs can be detected in diverse cancers characterized by chimeric TFs. Altogether, we identify the transcription, processing, and translation of novel, specific, highly expressed multi-exonic transcripts from otherwise silent regions of the genome as a new activity of aberrant TFs in cancer. These neotranscripts, which can be classified as LincRNAs, or neopeptides may have specific functions in the tumor cells. The latter may also constitute ideal neoepitopes as they are tumor-specific, clonally expressed in sarcoma cells, induced by the driver genetic alterations and public for a given sarcoma.

Citation Format: Olivier Delattre, Julien Vibert, Olivier Saulnier, Joshua Waterfall. Oncogenic transcription factors drive highly tumor-specific LincRNA [abstract]. In: Proceedings of the AACR Special Conference: Sarcomas; 2022 May 9-12; Montreal, QC, Canada. Philadelphia (PA): AACR; Clin Cancer Res 2022;28(18_Suppl):Abstract nr IA009.

Failure of human rhombic lip differentiation underlies medulloblastoma formation

Medulloblastoma (MB) comprises a group of heterogeneous paediatric embryonal neoplasms of the hindbrain with strong links to early development of the hindbrain^1-4. Mutations that activate Sonic hedgehog signalling lead to Sonic hedgehog MB in the upper rhombic lip (RL) granule cell lineage^5-8. By contrast, mutations that activate WNT signalling lead to WNT MB in the lower RL^9,10. However, little is known about the more commonly occurring group 4 (G4) MB, which is thought to arise in the unipolar brush cell lineage^3,4. Here we demonstrate that somatic mutations that cause G4 MB converge on the core binding factor alpha (CBFA) complex and mutually exclusive alterations that affect CBFA2T2, CBFA2T3, PRDM6, UTX and OTX2. CBFA2T2 is expressed early in the progenitor cells of the cerebellar RL subventricular zone in Homo sapiens, and G4 MB transcriptionally resembles these progenitors but are stalled in developmental time. Knockdown of OTX2 in model systems relieves this differentiation blockade, which allows MB cells to spontaneously proceed along normal developmental differentiation trajectories. The specific nature of the split human RL, which is destined to generate most of the neurons in the human brain, and its high level of susceptible EOMES⁺KI67⁺ unipolar brush cell progenitor cells probably predisposes our species to the development of G4 MB.

Single-cell transcriptomics reveals shared immunosuppressive landscapes of mouse and human neuroblastoma

BACKGROUND

High-risk neuroblastoma is a pediatric cancer with still a dismal prognosis, despite multimodal and intensive therapies. Tumor microenvironment represents a key component of the tumor ecosystem the complexity of which has to be accurately understood to define selective targeting opportunities, including immune-based therapies.

METHODS

We combined various approaches including single-cell transcriptomics to dissect the tumor microenvironment of both a transgenic mouse neuroblastoma model and a cohort of 10 biopsies from neuroblastoma patients, either at diagnosis or at relapse. Features of related cells were validated by multicolor flow cytometry and functional assays.

RESULTS

We show that the immune microenvironment of MYCN-driven mouse neuroblastoma is characterized by a low content of T cells, several phenotypes of macrophages and a population of cells expressing signatures of myeloid-derived suppressor cells (MDSCs) that are molecularly distinct from the various macrophage subsets. We document two cancer-associated fibroblasts (CAFs) subsets, one of which corresponding to CAF-S1, known to have immunosuppressive functions. Our data unravel a complex content in myeloid cells in patient tumors and further document a striking correspondence of the microenvironment populations between both mouse and human tumors. We show that mouse intratumor T cells exhibit increased expression of inhibitory receptors at the protein level. Consistently, T cells from patients are characterized by features of exhaustion, expressing inhibitory receptors and showing low expression of effector cytokines. We further functionally demonstrate that MDSCs isolated from mouse neuroblastoma have immunosuppressive properties, impairing the proliferation of T lymphocytes.

CONCLUSIONS

Our study demonstrates that neuroblastoma tumors have an immunocompromised microenvironment characterized by dysfunctional T cells and accumulation of immunosuppressive cells. Our work provides a new and precious data resource to better understand the neuroblastoma ecosystem and suggest novel therapeutic strategies, targeting both tumor cells and components of the microenvironment.

Abstract 6245: 3D-models of pediatric bone sarcomas for personalized therapeutic screening

Osteosarcoma (OS) and Ewing sarcoma (ES) are the most common bone cancers in children. They are rare cancers and thus difficult to study due to scarcity of patient material, large genomic instability and a wide histological heterogeneity (in OS) or a lack of satisfactory transgenic animal model and availability of preclinical tests (in ES). There is a dire need for new models and novel therapeutic approaches. Although patient-derived xenografts (PDXs) may recapitulate human tumor biology and predict drug response, propagating PDXs in mice limits its use as a drug-testing platform.

We have established and standardized ES and OS spheroid culture and developed a semi-automated drug-screening platform in tumor spheroids. We established several robust techniques for spheroid formation, with clear pathophysiological gradients, but without central necroses at the onset of drug treatment. We performed RNA-seq comparing spheroid transcription profiles to 2D culture and observed dramatic changes in overall expression patterns. We observed upregulation of genes shown to correlate with poor prognosis in OS patients. We saw upregulation of processes associated with regulation of cell migration, negative regulation of proliferation and modulation of the extracellular matrix (ECM). In addition to ES spheroid models, we created bioprinted 3D-models of ES cell lines and of cells obtained from ES PDXs, using extrusion bioprinting techniques (where cells are encapsulated within the cross-linked polymers, thus allowing homogeneous distribution and high cell density). PDX-derived cells were kept in liquid culture and as 3D-bioprinted constructs, while their transcription profiles were compared with the initial PDX. The mevalonate pathway was the most overrepresented in all ES 3D-models, consistent with predominant upregulation of this metabolic pathway integral to tumor growth and progression. After 15 days in 3D-bioprinted culture, we observed pronounced upregulation of genes involved in ECM signaling, suggesting that the construct promoted in vivo-like tumor-ECM interactions, without further promoting main proliferation and cell survival pathways, which was observed in liquid culture. Furthermore, we showed potential for combinatorial treatment with statins and confirmed feasibility of drug testing in patient-derived 3D models.

Finally, as our spheroid models showed upregulation of many processes involved in metastasis (genes associated with invasion, migration, angiogenesis and hypoxia), we focused on lung as the most common site of metastasis in ES and OS patients. We are thus establishing mixed airway organoid/tumoroid cultures, to investigate further the lung metastatic niche, with a goal to provide proof of concept for patient-specific 3D-models of lung metastatic tumors to guide personalized drug selection for patients with advanced disease.

Citation Format: Branka Radic Sarikas, Mathias Ilg, Marica Markovic, Caterina Sturtzel, Eva Scheuringer, Justine Zulini, Martin Metzelder, Florian Halbritter, Martin Distel, Didier Surdez, Olivier Delattre, Aleksandr Ovsianikov, Heinrich Kovar. 3D-models of pediatric bone sarcomas for personalized therapeutic screening [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 6245.

MEDB-84. The French experience of ELP1-related medulloblastomas

Medulloblastoma (MB), the most frequent embryonic tumor of the cerebellum is classified into four molecular subgroups (WNT group, SHH group, group 3 and group 4). Although the vast majority of MB are sporadic, predisposing genetic diseases have been described in rare WNT MB and more frequently in the SHH group. In a recent pediatric series of SHH-MB, germline alterations of the ELP1 gene have been described in 14% of cases, making this gene the most frequent genetic predisposition in MB. We have investigated the potential interest of ELP1 immunostaining on a large cohort of 132 MB. A complete loss of ELP1 staining was observed in 12 SHH MB (among 57 total SHH MB: 21%). The loss of ELP1 immunostaining was well correlated with the presence of a bi-allelic alteration of the gene except for one case for which the MB had a loss of ELP1 protein expression demonstrated by immunohistochemistry (IHC) and confirmed by whole proteome analysis, although no obvious genetic alteration in the coding sequence of ELP1 could be found. Molecular analysis of a large “molecular” cohort of 266 MB from French centers for which somatic ELP1 was sequenced allows to identify 12 additional MB with bi-allelic ELP1 genetic alterations. Our results demonstrate the benefit of the ELP1 IHC as an accurate and reliable tool to screen ELP1-deficient MB. This new immunohistochemical tool will now be advantageously used to screen SHH MB upfront for genetic alteration in ELP1, and will subsequently help orientating these patients towards genetic counseling.

PATH-12. A new powerful tool to simultaneously deliver DNA methylation profile and driver fusions detection based on Nanopore Sequencing

The fast advances in molecular biology have led to the discovery of new genetic and epigenetic driver events that play a key role in the pathogenesis of central nervous system (CNS) tumors. These data showed the importance of molecular biology, and in particular methylation profiling and fusion detection, to complete the histological examination. Exploring all these complex genomic and epigenomic alterations with the technologies implemented in the routine diagnosis requires massive laboratory equipment. In addition, it could take weeks before getting data that could be important to adapt the diagnosis, prognosis, and to ultimately guide the treatment’s choice. We describe the innovative potential of a nanopore sequencing device to obtain fast and multiple information to characterize CNS tumors. Adaptive-sampling is a software-controlled enrichment unique method that allows simultaneous methylation profiling, copy number landscape assessment, and fusion gene detection. Thereby, this technology brings comprehensive molecular information that could be very important for CNS tumor classification in a very short delay (below one week) and with a relative low cost compared to the combined other molecular analysis needed to get the same information. We report ten cases of diverse CNS tumors with already known oncogenic structural variants such as YAP1-MAMLD1 fused ependymoma, CNS neuroblastoma with FOXR2 activation, CNS high-grade neuroepithelial tumor with BCOR and MEN1 alteration, and atypical teratoid rhabdoid tumors (ATRT). For these tumors, using adaptive sampling with the nanopore sequencing, we were able to detect fusion genes as driver events, while obtaining at the same time the methylation profile to corroborate the tumor classification. Our results with this combined approach, never described before, demonstrate its feasibility as a new extremely powerful tool and show remarkable perspectives for this technology to bring a fast and comprehensive tumor molecular characterization.

ATRT-11. Anatomico-biological correlations define a new layer for ATRT molecular subgroups pointing to potential lineages of origin

Atypical teratoid rhabdoid tumors (ATRT) are divided in three molecular subgroups, the so-called MYC, TYR and SHH subgroups. This heterogeneity suggests some diversity in the cells of origin, which remain hypothetical thus far. A careful radiological review of 55 MRI at diagnosis was performed in parallel with a careful analysis of mouse tumor origin in the Rosa26-CreERT2::Smarcbflox/flox model. Methylation, bulk RNAseq and scRNAseq analyses were integrated to these anatomic information to highlight potential origin for each molecularly and anatomically defined subgroups. We demonstrated that mouse Myc-ATRTs derive from extra-parenchymal meningeal areas, a finding consistant with many human MYC ATRT being clearly of intra-cranial extra-axial origin. Although this finding could point to a neural crest origin, transcriptomic features fail to unravel any lineage-specific signature. We also defined a distinct supra-tentorial SHH ATRT subgroup, characterized both in mouse and Humans by neural features pointing to the ganglionic eminence progenitors as the candidate origin. Finally we identified a distinct infra-tentorial SHH ATRT subgroup, not observed in mice, with hindbrain/midbrain boundary progenitor signature. scRNAseq from human SHH infra-tentorial tumors consistently suggest a dedifferentiation process involving the Notch pathway in the oncogenic transformation of hindbrain/midbrain neural progenitors.

INSP-15. ITCC-P4: A sustainable platform of molecularly well-characterized PDX models of pediatric cancers for high throughputin vivo testing

Thanks to state-of-the-art molecular profiling techniques we by now have a much better understanding of pediatric cancers and what is driving them. On the other hand, we have also realized that pediatric cancers are much more heterogeneous than previously thought. Many new types and subtypes of pediatric cancers have been identified with distinct molecular and clinical characteristics. However, for many if not most of these new types and subtypes there is no specific treatment available, yet. In order to develop specific treatment protocols and to increase survival rates for pediatric cancer patients further, both at diagnosis and relapse/metastasis, we need a large collection of well-characterized preclinical models representing all the different types and subtypes. These models can be used for preclinical drug testing to prioritize the pediatric development of anticancer drugs that would be best targeting pediatric tumor biology. The ITCC-P4 consortium, which is a collaboration between many academic centers across Europe, several companies involved in in vivo preclinical testing, and ten pharmaceutical companies, started in 2017 with the overall aim to establish a sustainable platform of &gt;400 molecularly well-characterized PDX models of high-risk pediatric cancers and to use them for in vivo testing of novel mechanism-of-action based treatments. Currently, 340 models have been fully established, including 87 brain tumor models and 253 non-brain tumor models, together representing many different tumor types both from primary and relapsed/metastatic disease. Out of these 340 models, 252 have been fully molecularly characterized, most of them together with their matching original tumors, and almost of all these models are currently being subjected to in vivo testing using three standard of care drugs and six novel mechanism-of-action based drugs. In this presentation, an update on the current status of the ITCC-P4 platform and the data we collectively have generated thus far will be presented.

ATRT-21. Contribution of germline mosaic alterations of SMARCB1 in rhabdoid tumor predisposition syndrome

Rhabdoid tumors are rare and aggressive tumors that usually arise in very young children. They are characterized by a bi-allelic inactivation of the SMARCB1 gene. Although the majority of alterations of SMARCB1 are acquired in tumors, a heterozygous germline alteration is seen in one third of patients and defines the rhabdoid tumors predisposition syndrome. Penetrance is almost complete and the vast majority of germline alterations of SMARCB1 are acquired de novo but rare familial cases with a healthy carrier have also been described. Since the advent of more sensitive molecular analysis technologies such as next-generation sequencing (NGS), the number of mosaicisms of genes involved in genetic diseases discovered from blood samples has increased considerably. The aim of our study was to explore the mosaicisms of SMARCB1 in the blood 1/ of children with rhabdoid tumors with at least one alteration of SMARCB1 previously identified in the tumor but not found in the blood with old-fashioned low-sensitivity technologies and 2/ in parents of children with heterozygous germline alteration of SMARCB1. We analyzed a custom NGS panel which covers the SMARCB1 gene with an average depth of 1.500X on blood samples of 111 children with rhabdoid tumors and 32 parents collected at the Institut Curie since 1999. The mosaicism rate found in index cases was 11.7% (13/111) and 3.1% (1/32) in parents. The variant allele frequency vary from 0.8% to 12.9%. Our results also indicate to be cautious about the possible confounding effect of circulating tumor DNA. This hitherto underestimated SMARCB1 mosaicism rate should motivate an optimization of the genetic counseling as well as the oncological monitoring of these children and thus have a significant medical impact given the catastrophic prognosis of rhabdoid tumors.

TBIO-06. Study of cell free DNA extracted from CSF enables molecular classification of embryonal brain tumors

Liquid biopsies are revolutionary tools to detect tumour-specific genetic alterations in body fluids. Here, we assess whether the circulating tumor DNA (ctDNA) in Cerebral Spinal Fluid (CSF) could be used for the tumor genetic profiling of pediatric embryonal brain tumors (EBT). Cell free DNA extracted from CSF from 4-5 lumbar puncture droplets from patients with Medulloblastoma (n=18), ATRT (n=3), ETMR(n=1), FOXR2 EBT (n=2) or other pediatric EBT (n=1) (cfDNA mean quantity 48 ng/ml ; range 5.6 - 442 ng/ml) and matched genomic DNA from primary tumors were sequenced by WES (Illumina 100PE) using Nimblegen Medexome Capture. SNVs/mutations were called using GATK-UnifiedGenotyper, GATK-HaplotypeCaller and Samtools. Copy Number profiles were generated with CNVkits. 10/13 cfDNA WES yielded satisfactory depth (&gt;10x). A mean of 466 (range 93-945) SNVs were detected in the primary tumor and 474 (range 18-922) in the CSF. A mean of 416 (range 18-872) commons SNVs were observed between the cfDNA and the primary tumor, comprising classical medulloblastoma genes such as SMO or MLL2. Interestingly, several SNVs were observed either in the tumor only (mean 50 ; range 3-115) or in CSF only (mean 58 ; range 0-148) suggesting a clonal heterogeneity. For 5 cases, Copy Number profiles were also available, allowing the detection of MYCN amplification or 19q13 miRNA cluster amplification. Altogether, we demonstrate the feasibility of WES on CSF with a low input of ctDNA. These results may pave the way for new tumor monitoring tools.

ETMR-03. Intra- and extra-cranial BCOR-ITD tumours are separate entities within the BCOR-rearranged family

BCOR-ITD tumours form an emerging family of aggressive entities with an internal tandem duplication (ITD) in the last exon of the BCOR gene. The family includes cerebral tumours, termed central nervous system BCOR-ITD (CNS BCOR-ITD), and sarcomatous types described in the kidney as clear cell sarcoma of the kidney (CCSK), in the endometrium as high-grade endometrial stromal sarcoma (HG-ESS), in bone, and in soft tissue as undifferentiated round cell sarcoma (URCS) or primitive myxoid mesenchymal tumour of infancy (PMMTI). Based on a series of 33 retrospective cases, including 10 CNS BCOR-ITD and 23 BCOR-ITD sarcomas, we interrogated the homogeneity of the entity regarding clinical, radiological and histopathological findings, and molecular signatures. Whole transcriptomic sequencing and DNA methylation profiling were used for unsupervised clustering. Histopathological review revealed marked differences between CNS BCOR-ITD and BCOR-ITD sarcomas. These two groups were consistently segregated by unsupervised clustering of expression (n=22) and DNA methylation (n=21) data. Proximity between the two groups may result from common somatic changes within key pathways directly related to the novel activity of the ITD itself. Conversely, comparison of gene signatures with single-cell RNAseq atlases suggests that the distinction between BCOR-ITD sarcomas and CNS BCOR-ITD may result from differences in cells of origin.

Republication: Targeting PI3KC2β Impairs Proliferation and Survival in Acute Leukemia, Brain Tumours and Neuroendocrine Tumours

BACKGROUND

Eight human catalytic phosphoinositide 3-kinase (PI3K) isoforms exist which are subdivided into three classes. While class I isoforms have been well-studied in cancer, little is known about the functions of class II PI3Ks.

MATERIALS AND METHODS

The expression pattern and functions of the class II PI3KC2β isoform were investigated in a panel of tumour samples and cell lines.

RESULTS

Overexpression of PI3KC2β was found in subsets of tumours and cell lines from acute myeloid leukemia (AML), glioblastoma multiforme (GBM), medulloblastoma (MB), neuroblastoma (NB), and small cell lung cancer (SCLC). Specific pharmacological inhibitors of PI3KC2β or RNA interference impaired proliferation of a panel of human cancer cell lines and primary cultures. Inhibition of PI3KC2β also induced apoptosis and sensitised the cancer cells to chemotherapeutic agents.

CONCLUSION

Together, these data show that PI3KC2β contributes to proliferation and survival in AML, brain tumours and neuroendocrine tumours, and may represent a novel target in these malignancies.

Upregulation of the Mevalonate Pathway through EWSR1-FLI1/EGR2 Regulatory Axis Confers Ewing Cells Exquisite Sensitivity to Statins

Simple Summary

The objective of this project was to search for new dependencies in Ewing sarcoma, a deadly disease for which new therapeutic approaches are urgently needed. A pharmacological screening of off-patent approved drugs (FDA agency) and the investigation of downstream targets of EGR2 were performed. The two approaches showed the MVA pathway as a major dependency in Ewing sarcoma and statin, an inhibitor of this pathway, as a potential new therapeutic agent for the treatment of Ewing sarcoma.

Abstract

Ewing sarcoma (EwS) is an aggressive primary bone cancer in children and young adults characterized by oncogenic fusions between genes encoding FET-RNA-binding proteins and ETS transcription factors, the most frequent fusion being EWSR1-FLI1. We show that EGR2, an Ewing-susceptibility gene and an essential direct target of EWSR1-FLI1, directly regulates the transcription of genes encoding key enzymes of the mevalonate (MVA) pathway. Consequently, Ewing sarcoma is one of the tumors that expresses the highest levels of mevalonate pathway genes. Moreover, genome-wide screens indicate that MVA pathway genes constitute major dependencies of Ewing cells. Accordingly, the statin inhibitors of HMG-CoA-reductase, a rate-limiting enzyme of the MVA pathway, demonstrate cytotoxicity in EwS. Statins induce increased ROS and lipid peroxidation levels, as well as decreased membrane localization of prenylated proteins, such as small GTP proteins. These metabolic effects lead to an alteration in the dynamics of S-phase progression and to apoptosis. Statin-induced effects can be rescued by downstream products of the MVA pathway. Finally, we further show that statins impair tumor growth in different Ewing PDX models. Altogether, the data show that statins, which are off-patent, well-tolerated, and inexpensive compounds, should be strongly considered in the therapeutic arsenal against this deadly childhood disease.

Oncogenic chimeric transcription factors drive tumor-specific transcription, processing, and translation of silent genomic regions

Many cancers are characterized by gene fusions encoding oncogenic chimeric transcription factors (TFs) such as EWS::FLI1 in Ewing sarcoma (EwS). Here, we find that EWS::FLI1 induces the robust expression of a specific set of novel spliced and polyadenylated transcripts within otherwise transcriptionally silent regions of the genome. These neogenes (NGs) are virtually undetectable in large collections of normal tissues or non-EwS tumors and can be silenced by CRISPR interference at regulatory EWS::FLI1-bound microsatellites. Ribosome profiling and proteomics further show that some NGs are translated into highly EwS-specific peptides. More generally, we show that hundreds of NGs can be detected in diverse cancers characterized by chimeric TFs. Altogether, this study identifies the transcription, processing, and translation of novel, specific, highly expressed multi-exonic transcripts from otherwise silent regions of the genome as a new activity of aberrant TFs in cancer.

The European MAPPYACTS Trial: Precision Medicine Program in Pediatric and Adolescent Patients with Recurrent Malignancies

ABSTRACT

MAPPYACTS (NCT02613962) is an international prospective precision medicine trial aiming to define tumor molecular profiles in pediatric patients with recurrent/refractory malignancies in order to suggest the most adapted salvage treatment. From February 2016 to July 2020, 787 patients were included in France, Italy, Ireland, and Spain. At least one genetic alteration leading to a targeted treatment suggestion was identified in 436 patients (69%) with successful sequencing; 10% of these alterations were considered "ready for routine use." Of 356 patients with follow-up beyond 12 months, 107 (30%) received one or more matched targeted therapies-56% of them within early clinical trials-mainly in the AcSé-ESMART platform trial (NCT02813135). Overall, matched treatment resulted in a 17% objective response rate, and of those patients with ready for routine use alterations, it was 38%. In patients with extracerebral tumors, 76% of actionable alterations detected in tumor tissue were also identified in circulating cell-free DNA (cfDNA).

SIGNIFICANCE

MAPPYACTS underlines the feasibility of molecular profiling at cancer recurrence in children on a multicenter, international level and demonstrates benefit for patients with selected key drivers. The use of cfDNA deserves validation in prospective studies. Our study highlights the need for innovative therapeutic proof-of-concept trials that address the underlying cancer complexity. This article is highlighted in the In This Issue feature, p. 1171.

Intra- and extra-cranial BCOR-ITD tumours are separate entities within the BCOR-rearranged family

BCOR-ITD tumours form an emerging family of aggressive entities with an internal tandem duplication (ITD) in the last exon of the BCOR gene. The family includes cerebral tumours, termed central nervous system BCOR-ITD (CNS BCOR-ITD), and sarcomatous types described in the kidney as clear cell sarcoma of the kidney (CCSK), in the endometrium as high-grade endometrial stromal sarcoma, and in the bone and soft tissue as undifferentiated round cell sarcoma or primitive myxoid mesenchymal tumour of infancy. Based on a series of 33 retrospective cases, including 10 CNS BCOR-ITD and 23 BCOR-ITD sarcomas, we interrogated the homogeneity of the entity regarding clinical, radiological, and histopathological findings, and molecular signatures. Whole-transcriptomic sequencing and DNA methylation profiling were used for unsupervised clustering. BCOR-ITD tumours mostly affected young children with a median age at diagnosis of 2.1 years (range 0-62.4). Median overall survival was 3.9 years and progression-free survival was 1.4 years. This dismal prognosis is shared among tumours in all locations except CCSK. Histopathological review revealed marked differences between CNS BCOR-ITD and BCOR-ITD sarcomas. These two groups were consistently segregated by unsupervised clustering of expression (n = 22) and DNA methylation (n = 21) data. Proximity between the two groups may result from common somatic changes within key pathways directly related to the novel activity of the ITD itself. Conversely, comparison of gene signatures with single-cell RNA-Seq atlases suggests that the distinction between BCOR-ITD sarcomas and CNS BCOR-ITD may result from differences in cells of origin.