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.

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.

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.

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 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.

Molecular diagnosis of retinoblastoma by circulating tumor DNA analysis.

PURPOSE: The analysis of circulating tumor DNA (ctDNA), a fraction of total cell-free DNA (cfDNA), might be of special interest in retinoblastoma patients. Because the accessibility to tumor tissue is very limited in these patients, either for histopathological diagnosis of suspicious intraocular masses (biopsies are proscribed) or for somatic RB1 studies and genetic counseling (due to current successful conservative approaches), we aim to validate the detection of ctDNA in plasma of non-hereditary retinoblastoma patients by molecular analysis of RB1 gene. EXPERIMENTAL DESIGN: In a cohort of 19 intraocular unilateral non-hereditary retinoblastoma patients for whom a plasma sample was available at diagnosis, we performed high-deep next-generation sequencing (NGS) of RB1 in cfDNA. Two different bioinformatics/statistics approaches were applied depending on whether the somatic RB1 status was available or not. RESULTS: Median plasma sample volume was 600 µL [100-1000]; median cfDNA plasma concentration was 119 [38-1980] and 27 [11-653] ng/mL at diagnosis and after complete remission, respectively. In the subgroup of patients with known somatic RB1 alterations (n = 11), seven of nine somatic mutations were detected (median allele fraction: 6.7%). In patients without identified somatic RB1 alterations (n = 8), six candidate variants were identified for seven patients. CONCLUSIONS: Despite small tumor size, blood-ocular barrier, poor ctDNA blood release and limited plasma sample volumes, we confirm that it is possible to detect ctDNA with high-deep NGS in plasma from patients with intraocular non-hereditary retinoblastoma. This may aid in diagnosis of suspicious cases, family genetic counseling or follow-up of residual intraocular disease.

Frequency and Prognostic Impact of ALK Amplifications and Mutations in the European Neuroblastoma Study Group (SIOPEN) High-Risk Neuroblastoma Trial (HR-NBL1).

PURPOSE: In neuroblastoma (NB), the ALK receptor tyrosine kinase can be constitutively activated through activating point mutations or genomic amplification. We studied ALK genetic alterations in high-risk (HR) patients on the HR-NBL1/SIOPEN trial to determine their frequency, correlation with clinical parameters, and prognostic impact. MATERIALS AND METHODS: Diagnostic tumor samples were available from 1,092 HR-NBL1/SIOPEN patients to determine ALK amplification status (n = 330), ALK mutational profile (n = 191), or both (n = 571). RESULTS: Genomic ALK amplification (ALKa) was detected in 4.5% of cases (41 out of 901), all except one with MYCN amplification (MNA). ALKa was associated with a significantly poorer overall survival (OS) (5-year OS: ALKa [n = 41] 28% [95% CI, 15 to 42]; no-ALKa [n = 860] 51% [95% CI, 47 to 54], [P < .001]), particularly in cases with metastatic disease. ALK mutations (ALKm) were detected at a clonal level (> 20% mutated allele fraction) in 10% of cases (76 out of 762) and at a subclonal level (mutated allele fraction 0.1%-20%) in 3.9% of patients (30 out of 762), with a strong correlation between the presence of ALKm and MNA (P < .001). Among 571 cases with known ALKa and ALKm status, a statistically significant difference in OS was observed between cases with ALKa or clonal ALKm versus subclonal ALKm or no ALK alterations (5-year OS: ALKa [n = 19], 26% [95% CI, 10 to 47], clonal ALKm [n = 65] 33% [95% CI, 21 to 44], subclonal ALKm (n = 22) 48% [95% CI, 26 to 67], and no alteration [n = 465], 51% [95% CI, 46 to 55], respectively; P = .001). Importantly, in a multivariate model, involvement of more than one metastatic compartment (hazard ratio [HR], 2.87; P < .001), ALKa (HR, 2.38; P = .004), and clonal ALKm (HR, 1.77; P = .001) were independent predictors of poor outcome. CONCLUSION: Genetic alterations of ALK (clonal mutations and amplifications) in HR-NB are independent predictors of poorer survival. These data provide a rationale for integration of ALK inhibitors in upfront treatment of HR-NB with ALK alterations.

Multimodal analysis of cell-free DNA whole-genome sequencing for pediatric cancers with low mutational burden.

Sequencing of cell-free DNA in the blood of cancer patients (liquid biopsy) provides attractive opportunities for early diagnosis, assessment of treatment response, and minimally invasive disease monitoring. To unlock liquid biopsy analysis for pediatric tumors with few genetic aberrations, we introduce an integrated genetic/epigenetic analysis method and demonstrate its utility on 241 deep whole-genome sequencing profiles of 95 patients with Ewing sarcoma and 31 patients with other pediatric sarcomas. Our method achieves sensitive detection and classification of circulating tumor DNA in peripheral blood independent of any genetic alterations. Moreover, we benchmark different metrics for cell-free DNA fragmentation analysis, and we introduce the LIQUORICE algorithm for detecting circulating tumor DNA based on cancer-specific chromatin signatures. Finally, we combine several fragmentation-based metrics into an integrated machine learning classifier for liquid biopsy analysis that exploits widespread epigenetic deregulation and is tailored to cancers with low mutation rates. Clinical associations highlight the potential value of cfDNA fragmentation patterns as prognostic biomarkers in Ewing sarcoma. In summary, our study provides a comprehensive analysis of circulating tumor DNA beyond recurrent genetic aberrations, and it renders the benefits of liquid biopsy more readily accessible for childhood cancers.

Low-frequency variation near common germline susceptibility loci are associated with risk of Ewing sarcoma.

BACKGROUND: Ewing sarcoma (EwS) is a rare, aggressive solid tumor of childhood, adolescence and young adulthood associated with pathognomonic EWSR1-ETS fusion oncoproteins altering transcriptional regulation. Genome-wide association studies (GWAS) have identified 6 common germline susceptibility loci but have not investigated low-frequency inherited variants with minor allele frequencies below 5% due to limited genotyped cases of this rare tumor. METHODS: We investigated the contribution of rare and low-frequency variation to EwS susceptibility in the largest EwS genome-wide association study to date (733 EwS cases and 1,346 unaffected controls of European ancestry). RESULTS: We identified two low-frequency variants, rs112837127 and rs2296730, on chromosome 20 that were associated with EwS risk (OR = 0.186 and 2.038, respectively; P-value < 5×10-8) and located near previously reported common susceptibility loci. After adjusting for the most associated common variant at the locus, only rs112837127 remained a statistically significant independent signal (OR = 0.200, P-value = 5.84×10-8). CONCLUSIONS: These findings suggest rare variation residing on common haplotypes are important contributors to EwS risk. IMPACT: Motivate future targeted sequencing studies for a comprehensive evaluation of low-frequency and rare variation around common EwS susceptibility loci.

Study of chromatin remodeling genes implicates SMARCA4 as a putative player in oncogenesis in neuroblastoma.

In neuroblastoma (NB), genetic alterations in chromatin remodeling (CRGs) and epigenetic modifier genes (EMGs) have been described. We sought to determine their frequency and clinical impact. Whole exome (WES)/whole genome sequencing (WGS) data and targeted sequencing (TSCA®) of exonic regions of 33 CRGs/EMGs were analyzed in tumor samples from 283 NB patients, with constitutional material available for 55 patients. The frequency of CRG/EMG variations in NB cases was then compared to the Genome Aggregation Database (gnomAD). The sequencing revealed SNVs/small InDels or focal CNAs of CRGs/EMGs in 20% (56/283) of all cases, occurring at a somatic level in 4 (7.2%), at a germline level in 12 (22%) cases, whereas for the remaining cases, only tumor material could be analyzed. The most frequently altered genes were ATRX (5%), SMARCA4 (2.5%), MLL3 (2.5%) and ARID1B (2.5%). Double events (SNVs/small InDels/CNAs associated with LOH) were observed in SMARCA4 (n = 3), ATRX (n = 1) and PBRM1 (n = 1). Among the 60 variations, 24 (8.4%) targeted domains of functional importance for chromatin remodeling or highly conserved domains but of unknown function. Variations in SMARCA4 and ATRX occurred more frequently in the NB as compared to the gnomAD control cohort (OR = 4.49, 95%CI: 1.63-9.97, p = 0.038; OR 3.44, 95%CI: 1.46-6.91, p = 0.043, respectively). Cases with CRG/EMG variations showed a poorer overall survival compared to cases without variations. Genetic variations of CRGs/EMGs with likely functional impact were observed in 8.4% (24/283) of NB. Our case-control approach suggests a role of SMARCA4 as a player of NB oncogenesis.

Circulating tumor DNA analysis enables molecular characterization of pediatric renal tumors at diagnosis.

Circulating tumor DNA (ctDNA) is a powerful tool for the molecular characterization of cancer. The most frequent pediatric kidney tumors (KT) are Wilms' tumors (WT), but other diagnoses may occur. According to the SIOP strategy, in most countries pediatric KT have a presumptive diagnosis of WT if they are clinically and radiologically compatible. The histologic confirmation is established after post-chemotherapy nephrectomy. Thus, there is a risk for a small fraction of patients to receive neoadjuvant chemotherapy that is not adapted to the disease. The aim of this work is to perform molecular diagnosis of pediatric KT by tumor genetic characterization based on the analysis of ctDNA. We analyzed ctDNA extracted from plasma samples of 18 pediatric patients with KT by whole-exome sequencing and compared the results to their matched tumor and germline DNA. Copy number alterations (CNAs) and single nucleotide variations (SNVs) were analyzed. We were able to detect tumor cell specific genetic alterations-CNAs, SNVs or both-in ctDNA in all patients except in one (for whom the plasma sample was obtained long after nephrectomy). These results open the door to new applications for the study of ctDNA with regards to the molecular diagnosis of KT, with a possibility of its usefulness for adapting the treatment early after diagnosis, but also for disease monitoring and follow up.

Genome-wide association study identifies multiple new loci associated with Ewing sarcoma susceptibility.

Ewing sarcoma (EWS) is a pediatric cancer characterized by the EWSR1-FLI1 fusion. We performed a genome-wide association study of 733 EWS cases and 1346 unaffected individuals of European ancestry. Our study replicates previously reported susceptibility loci at 1p36.22, 10q21.3 and 15q15.1, and identifies new loci at 6p25.1, 20p11.22 and 20p11.23. Effect estimates exhibit odds ratios in excess of 1.7, which is high for cancer GWAS, and striking in light of the rarity of EWS cases in familial cancer syndromes. Expression quantitative trait locus (eQTL) analyses identify candidate genes at 6p25.1 (RREB1) and 20p11.23 (KIZ). The 20p11.22 locus is near NKX2-2, a highly overexpressed gene in EWS. Interestingly, most loci reside near GGAA repeat sequences and may disrupt binding of the EWSR1-FLI1 fusion protein. The high locus to case discovery ratio from 733 EWS cases suggests a genetic architecture in which moderate risk SNPs constitute a significant fraction of risk.

Genomic Amplifications and Distal 6q Loss: Novel Markers for Poor Survival in High-risk Neuroblastoma Patients.

Background: Neuroblastoma is characterized by substantial clinical heterogeneity. Despite intensive treatment, the survival rates of high-risk neuroblastoma patients are still disappointingly low. Somatic chromosomal copy number aberrations have been shown to be associated with patient outcome, particularly in low- and intermediate-risk neuroblastoma patients. To improve outcome prediction in high-risk neuroblastoma, we aimed to design a prognostic classification method based on copy number aberrations. Methods: In an international collaboration, normalized high-resolution DNA copy number data (arrayCGH and SNP arrays) from 556 high-risk neuroblastomas obtained at diagnosis were collected from nine collaborative groups and segmented using the same method. We applied logistic and Cox proportional hazard regression to identify genomic aberrations associated with poor outcome. Results: In this study, we identified two types of copy number aberrations that are associated with extremely poor outcome. Distal 6q losses were detected in 5.9% of patients and were associated with a 10-year survival probability of only 3.4% (95% confidence interval [CI] = 0.5% to 23.3%, two-sided P = .002). Amplifications of regions not encompassing the MYCN locus were detected in 18.1% of patients and were associated with a 10-year survival probability of only 5.8% (95% CI = 1.5% to 22.2%, two-sided P < .001). Conclusions: Using a unique large copy number data set of high-risk neuroblastoma cases, we identified a small subset of high-risk neuroblastoma patients with extremely low survival probability that might be eligible for inclusion in clinical trials of new therapeutics. The amplicons may also nominate alternative treatments that target the amplified genes.

QuantumClone: clonal assessment of functional mutations in cancer based on a genotype-aware method for clonal reconstruction.

Motivation: In cancer, clonal evolution is assessed based on information coming from single nucleotide variants and copy number alterations. Nonetheless, existing methods often fail to accurately combine information from both sources to truthfully reconstruct clonal populations in a given tumor sample or in a set of tumor samples coming from the same patient. Moreover, previously published methods detect clones from a single set of variants. As a result, compromises have to be done between stringent variant filtering [reducing dispersion in variant allele frequency estimates (VAFs)] and using all biologically relevant variants. Results: We present a framework for defining cancer clones using most reliable variants of high depth of coverage and assigning functional mutations to the detected clones. The key element of our framework is QuantumClone, a method for variant clustering into clones based on VAFs, genotypes of corresponding regions and information about tumor purity. We validated QuantumClone and our framework on simulated data. We then applied our framework to whole genome sequencing data for 19 neuroblastoma trios each including constitutional, diagnosis and relapse samples. We confirmed an enrichment of damaging variants within such pathways as MAPK (mitogen-activated protein kinases), neuritogenesis, epithelial-mesenchymal transition, cell survival and DNA repair. Most pathways had more damaging variants in the expanding clones compared to shrinking ones, which can be explained by the increased total number of variants between these two populations. Functional mutational rate varied for ancestral clones and clones shrinking or expanding upon treatment, suggesting changes in clone selection mechanisms at different time points of tumor evolution. Availability and implementation: Source code and binaries of the QuantumClone R package are freely available for download at https://CRAN.R-project.org/package=QuantumClone. Contact: gudrun.schleiermacher@curie.fr or valentina.boeva@inserm.fr. Supplementary information: Supplementary data are available at Bioinformatics online.

Whole-Exome Sequencing of Cell-Free DNA Reveals Temporo-spatial Heterogeneity and Identifies Treatment-Resistant Clones in Neuroblastoma.

Purpose: Neuroblastoma displays important clinical and genetic heterogeneity, with emergence of new mutations at tumor progression.Experimental Design: To study clonal evolution during treatment and follow-up, an innovative method based on circulating cell-free DNA (cfDNA) analysis by whole-exome sequencing (WES) paired with target sequencing was realized in sequential liquid biopsy samples of 19 neuroblastoma patients.Results: WES of the primary tumor and cfDNA at diagnosis showed overlap of single-nucleotide variants (SNV) and copy number alterations, with 41% and 93% of all detected alterations common to the primary neuroblastoma and cfDNA. CfDNA WES at a second time point indicated a mean of 22 new SNVs for patients with progressive disease. Relapse-specific alterations included genes of the MAPK pathway and targeted the protein kinase A signaling pathway. Deep coverage target sequencing of intermediate time points during treatment and follow-up identified distinct subclones. For 17 seemingly relapse-specific SNVs detected by cfDNA WES at relapse but not tumor or cfDNA WES at diagnosis, deep coverage target sequencing detected these alterations in minor subclones, with relapse-emerging SNVs targeting genes of neuritogenesis and cell cycle. Furthermore a persisting, resistant clone with concomitant disappearance of other clones was identified by a mutation in the ubiquitin protein ligase HERC2Conclusions: Modelization of mutated allele fractions in cfDNA indicated distinct patterns of clonal evolution, with either a minor, treatment-resistant clone expanding to a major clone at relapse, or minor clones collaborating toward tumor progression. Identification of treatment-resistant clones will enable development of more efficient treatment strategies. Clin Cancer Res; 24(4); 939-49. ©2017 AACR.

Genomic Profiles of Neuroblastoma Associated With Opsoclonus Myoclonus Syndrome.

Opsoclonus myoclonus syndrome (OMS), often called "dancing eyed syndrome," is a rare neurological condition associated with neuroblastoma in the majority of all childhood cases. Genomic copy number profiles have shown to be of prognostic significance for neuroblastoma patients. The aim of this retrospective multicenter study was to analyze the genomic copy number profiles of tumors from children with neuroblastoma presenting with OMS at diagnosis. In 44 cases of neuroblastoma associated with OMS, overall genomic profiling by either array-comparative genomic hybridization or single nucleotide polymorphism array proved successful in 91% of the cases, distinguishing tumors harboring segmental chromosome alterations from those with numerical chromosome alterations only. A total of 23/44 (52%) tumors showed an segmental chromosome alterations genomic profile, 16/44 (36%) an numerical chromosome alterations genomic profile, and 1 case displayed an atypical profile (12q amplicon). No recurrently small interstitial copy number alterations were identified. With no tumor relapse nor disease-related deaths, the overall genomic profile was not of prognostic impact with regard to the oncological outcome in this series of patients. Thus, the observation of an excellent oncological outcome, even for those with an unfavorable genomic profile of neuroblastoma, supports the hypothesis that an immune response might be involved in tumor control in these patients with OMS.

Radiogenomics of neuroblastomas: Relationships between imaging phenotypes, tumor genomic profile and survival.

PURPOSE: This study investigated relationships between neuroblastomas (NBs) imaging phenotypes, tumor genomic profile and patient outcome. PATIENTS AND METHODS: This IRB-approved retrospective observational study included 133 NB patients (73 M, 60 F; median age 15 months, range 0-151) treated in a single institution between 1998 and 2012. A consensus review of imaging (CT-scan, MRI) categorized tumors according to both the primarily involved compartment (i.e., neck, chest, abdomen or pelvis) and the sympathetic anatomical structure the tumors rose from (i.e., cervical, paravertebral or periarterial chains, or adrenal gland). Tumor shape, volume and image-defined surgical risk factors (IDRFs) at diagnosis were recorded. Genomic profiles were assessed using array-based comparative genomic hybridization and divided into three groups: "numerical-only chromosome alterations" (NCA), "segmental chromosome alterations" (SCA) and "MYCN amplification" (MNA). Statistical analyses included Kruskal-Wallis, Chi2 and Fisher's exact tests and the Kaplan-Meier method with log-rank tests and Cox model for univariate and multivariate survival analyses. RESULTS: A significant association between the sympathetic structure origin of tumors and genomic profiles was demonstrated. NBs arising from cervical sympathetic chains were all NCA. Paravertebral NBs were NCA or SCA in 75% and 25%, respectively and none were MNA. Periarterial NBs were NCA, SCA or MNA in 33%, 56% and 11%, respectively. Adrenal NBs were NCA, SCA or MNA in 16%, 36% and 48%, respectively. Among MNA NBs, 92% originated from the adrenal gland. The sympathetic anatomical classification was significantly better correlated to overall survival than the compartmental classification (P < .0003). The tumor volume of MNA NBs was significantly higher than NCA or SCA NBs (P < .0001). Patients with initial volume less than 160 mL had significantly better overall survival (P < .009). A "single mass" pattern was significantly more frequent in NCA NBs (P = .0003). The number of IDRFs was significantly higher in MNA NBs (P < .0001). CONCLUSION: Imaging phenotypes of neuroblastomas, including tumor origin along the sympathetic system, correlate with tumor genomic profile and patient outcome.

Heterogeneity of neuroblastoma cell identity defined by transcriptional circuitries.

Neuroblastoma is a tumor of the peripheral sympathetic nervous system, derived from multipotent neural crest cells (NCCs). To define core regulatory circuitries (CRCs) controlling the gene expression program of neuroblastoma, we established and analyzed the neuroblastoma super-enhancer landscape. We discovered three types of identity in neuroblastoma cell lines: a sympathetic noradrenergic identity, defined by a CRC module including the PHOX2B, HAND2 and GATA3 transcription factors (TFs); an NCC-like identity, driven by a CRC module containing AP-1 TFs; and a mixed type, further deconvoluted at the single-cell level. Treatment of the mixed type with chemotherapeutic agents resulted in enrichment of NCC-like cells. The noradrenergic module was validated by ChIP-seq. Functional studies demonstrated dependency of neuroblastoma with noradrenergic identity on PHOX2B, evocative of lineage addiction. Most neuroblastoma primary tumors express TFs from the noradrenergic and NCC-like modules. Our data demonstrate a previously unknown aspect of tumor heterogeneity relevant for neuroblastoma treatment strategies.

DNA methylation heterogeneity defines a disease spectrum in Ewing sarcoma.

Developmental tumors in children and young adults carry few genetic alterations, yet they have diverse clinical presentation. Focusing on Ewing sarcoma, we sought to establish the prevalence and characteristics of epigenetic heterogeneity in genetically homogeneous cancers. We performed genome-scale DNA methylation sequencing for a large cohort of Ewing sarcoma tumors and analyzed epigenetic heterogeneity on three levels: between cancers, between tumors, and within tumors. We observed consistent DNA hypomethylation at enhancers regulated by the disease-defining EWS-FLI1 fusion protein, thus establishing epigenomic enhancer reprogramming as a ubiquitous and characteristic feature of Ewing sarcoma. DNA methylation differences between tumors identified a continuous disease spectrum underlying Ewing sarcoma, which reflected the strength of an EWS-FLI1 regulatory signature and a continuum between mesenchymal and stem cell signatures. There was substantial epigenetic heterogeneity within tumors, particularly in patients with metastatic disease. In summary, our study provides a comprehensive assessment of epigenetic heterogeneity in Ewing sarcoma and thereby highlights the importance of considering nongenetic aspects of tumor heterogeneity in the context of cancer biology and personalized medicine.