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.

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

Introduction: Typical medullary breast carcinoma (MBC) has recently been recognized to be part of the basal-like carcinoma spectrum, a feature in agreement with the high rate of TP53 mutations previously reported in MBCs. The present study was therefore designed to identify phenotypic and genetic alterations that distinguish MBCs from basal-like carcinomas (BLC). Methods: Expression levels of estrogen receptor (ER), progesterone receptor (PR), ERBB2, TP53, cytokeratins (KRTs) 5/6, 14, 8/18, epidermal growth factor receptor and KIT, as well as TP53 gene sequence and high-density array comparative genomic hybridization (CGH) profiles, were assessed and compared in a series of 33 MBCs and 26 BLCs. Results: All tumors were negative for ER, PR and ERBB2. KRTs 5/6 were more frequently expressed in MBCs (94%) than in BLCs (56%) (p = 0.0004). TP53 mutations were disclosed in 20/26 MBCs (77%) and 20/24 BLCs (83%). Array CGH analysis showed that a higher number of gains (95 regions) and losses (34 regions) was observed in MBCs than in BLCs (36 regions of gain; 13 regions of losses). In addition, gains of 1q and 8q, and losses of X were found to be common to the two groups, whereas gains of 10p (53% of the cases), 9p (30.8% of the cases) and 16q (25.8% of the cases), and losses of 4p (34.8% of the cases), and amplicons of 1q, 8p, 10p and 12p were the genetic alterations found to characterize MBC. Conclusion: Our study has revealed that MBCs are part of the basal-like group and share common genomic alterations with BLCs, the most frequent being 1q and 8q gains and X losses; however, MBCs are a distinct entity within the basal-like spectrum, characterized by a higher rate of KRT 5/6 expression, a higher rate of gains and losses than BLCs, recurrent 10p, 9p and 16q gains, 4p losses, and 1q, 8p, 10p and 12p amplicons. Our results thus contribute to a better understanding of the heterogeneity in basal-like breast tumors and provide potential diagnostic tools.

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.