PURPOSE: Liquid biopsies can be used to investigate tumor-derived DNA, circulating in the cell-free DNA (cfDNA) pool in blood. We aimed to develop a droplet digital polymerase chain reaction (ddPCR) assay detecting hypermethylation of tumor suppressor gene RASSF1A as a simple standard test to detect various pediatric tumor types in small volume blood samples and to evaluate this test for monitoring treatment response of patients with high-risk neuroblastoma.METHODS: We developed a ddPCR assay to sensitively detect tumor-derived hypermethylated RASSF1A DNA in liquid biopsies. We tested this assay in plasma of 96 patients with neuroblastoma, renal tumors, rhabdomyosarcoma, or Hodgkin lymphoma at diagnosis and in cerebrospinal fluid of four patients with brain tumors. We evaluated the presence of hypermethylated RASSF1A in plasma samples during treatment and follow-up in 47 patients with neuroblastoma treated according to high-risk protocol and correlated results with blood mRNA–based and bone marrow mRNA–based minimal residual disease detection and clinical outcomes.RESULTS: The total cfDNA level was significantly higher in patients with metastatic neuroblastoma and nephroblastoma compared with healthy adult and pediatric controls. Hypermethylated RASSF1A was present in 41 of 42 patients with metastatic neuroblastoma and in all patients with nephroblastoma, with the median percentage of 69% and 21% of total RASSF1A, respectively. Hypermethylated RASSF1A levels decreased during therapy and recurred at relapse.CONCLUSION: Our findings demonstrate the value of ddPCR-based detection of hypermethylated RASSF1A as a circulating molecular tumor marker in neuroblastoma. Our preliminary investigation of RASSF1A hypermethylation detection in circulating cfDNA of other pediatric tumor entities demonstrates potential as a pan-tumor marker, but requires investigation in larger cohorts to evaluate its use and limitations.

Paediatric tumours are often characterised by the presence of recurrent DNA copy number alterations (CNAs). These DNA copy number profiles, obtained from a tissue biopsy, can aid in the correct prognostic classification and therapeutic stratification of several paediatric cancer entities (e.g. MYCN amplification in neuroblastoma) and are part of the routine diagnostic practice. Liquid biopsies (LQBs) offer a potentially safer alternative for such invasive tumour tissue biopsies and can provide deeper insight into tumour heterogeneity. The robustness and reliability of LQB CNA analyses was evaluated. We performed retrospective CNA profiling using shallow whole-genome sequencing (sWGS) on paired plasma circulating cell-free DNA (cfDNA) and tissue DNA samples from routinely collected samples from paediatric patients (n = 128) representing different tumour entities, including osteosarcoma, Ewing sarcoma, rhabdomyosarcoma, Wilms tumour, brain tumours and neuroblastoma. Overall, we observed a good concordance between CNAs in tissue DNA and cfDNA. The main cause of CNA discordance was found to be low cfDNA sample quality (i.e. the ratio of cfDNA (700 bp)). Furthermore, CNAs were observed that were present in cfDNA and not in tissue DNA, or vice-versa. In neuroblastoma samples, no false-positives or false-negatives were identified for the detection of the prognostic marker MYCN amplification. In future prospective studies, CNA analysis on LQBs that are of sufficient quality can serve as a complementary assay for CNA analysis on tissue biopsies, as either cfDNA or tissue DNA can contain CNAs that cannot be identified in the other biomaterial.

MRNA RT-qPCR is shown to be a very sensitive technique to detect minimal residual disease (MRD) in patients with neuroblastoma. Multiple mRNA markers are known to detect heterogeneous neuroblastoma cells in bone marrow (BM) or blood from patients. However, the limited volumes of BM and blood available can hamper the detection of multiple markers. To make optimal use of these samples, we developed a multiplex RT-qPCR for the detection of MRD in neuroblastoma. GUSB and PHOX2B were tested as single markers. The adrenergic markers TH, GAP43, CHRNA3 and DBH and mesenchymal markers POSTN, PRRX1 and FMO3 were tested in multiplex. Using control blood and BM, we established new thresholds for positivity. Comparison of multiplex and singleplex RT-qPCR results from 21 blood and 24 BM samples from neuroblastoma patients demonstrated a comparable sensitivity. With this multiplex RT-qPCR, we are able to test seven different neuroblastoma mRNA markers, which overcomes tumor heterogeneity and improves sensitivity of MRD detection, even in those samples of low RNA quantity. With resources and time being saved, reduction in sample volume and consumables can assist in the introduction of MRD by RT-qPCR into clinical practice.

PURPOSE: Total cell-free DNA (cfDNA) and tumor-derived cfDNA (ctDNA) can be used to study tumor-derived genetic aberrations. We analyzed the diagnostic and prognostic potential of cfDNA and ctDNA, obtained from pediatric patients with rhabdomyosarcoma. METHODS: cfDNA was isolated from diagnostic plasma samples from 57 patients enrolled in the EpSSG RMS2005 study. To study the diagnostic potential, shallow whole genome sequencing (shWGS) and cell-free reduced representation bisulphite sequencing (cfRRBS) were performed in a subset of samples and all samples were tested using droplet digital polymerase chain reaction to detect methylated RASSF1A ( RASSF1A-M). Correlation with outcome was studied by combining cfDNA RASSF1A-M detection with analysis of our rhabdomyosarcoma-specific RNA panel in paired cellular blood and bone marrow fractions and survival analysis in 56 patients. RESULTS: At diagnosis, ctDNA was detected in 16 of 30 and 24 of 26 patients using shallow whole genome sequencing and cfRRBS, respectively. Furthermore, 21 of 25 samples were correctly classified as embryonal by cfRRBS. RASSF1A-M was detected in 21 of 57 patients. The presence of RASSF1A-M was significantly correlated with poor outcome (the 5-year event-free survival [EFS] rate was 46.2% for 21 RASSF1A-M ‒positive patients, compared with 84.9% for 36 RASSF1A-M ‒negative patients [ P < .001]). RASSF1A-M positivity had the highest prognostic effect among patients with metastatic disease. Patients both negative for RASSF1A-M and the rhabdomyosarcoma-specific RNA panel (28 of 56 patients) had excellent outcome (5-year EFS 92.9%), while double-positive patients (11/56) had poor outcome (5-year EFS 13.6%, P < .001). CONCLUSION: Analyzing ctDNA at diagnosis using various techniques is feasible in pediatric rhabdomyosarcoma and has potential for clinical use. Measuring RASSF1A-M in plasma at initial diagnosis correlated significantly with outcome, particularly when combined with paired analysis of blood and bone marrow using a rhabdomyosarcoma-specific RNA panel.

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.

Purpose: Survival of children with rhabdomyosarcoma that suffer from recurrent or progressive disease is poor. Identifying these patients upfront remains challenging, indicating a need for improvement of risk stratification. Detection of tumor-derived mRNA in bone marrow (BM) and peripheral blood (PB) using reverse-transcriptase qPCR (RT-qPCR) is a more sensitive method to detect disseminated disease. We identified a panel of genes to optimize risk stratification by RT-qPCR. Experimental Design: Candidate genes were selected using gene expression data from rhabdomyosarcoma and healthy hematologic tissues, and a multiplexed RT-qPCR was developed. Significance of molecular disease was determined in a cohort of 99 Dutch patients with rhabdomyosarcoma (72 localized and 27 metastasized) treated according to the European pediatric Soft tissue sarcoma Study Group (EpSSG) RMS2005 protocol. Results: We identified the following 11 rhabdomyosarcoma markers: ZIC1, ACTC1, MEGF10, PDLIM3, SNAI2, CDH11, TMEM47, MYOD1, MYOG, and PAX3/7-FOXO1. RT-qPCR was performed for this 11-marker panel on BM and PB samples from the patient cohort. Five-year event-free survival (EFS) was 35.5% [95% confidence interval (CI), 17.5%–53.5%] for the 33/99 RNA-positive patients, versus 88.0% (95% CI, 78.9%–97.2%) for the 66/99 RNA-negative patients (P < 0.0001). Five-year overall survival (OS) was 54.8% (95% CI, 36.2%–73.4%) and 93.7% (95% CI, 86.6%–100.0%), respectively (P < 0.0001). RNA panel positivity was negatively associated with EFS (Hazard Ratio = 9.52; 95% CI, 3.23–28.02), whereas the RMS2005 risk group stratification was not, in the multivariate Cox regression model. Conclusions: This study shows a strong association between PCR-based detection of disseminated disease at diagnosis with clinical outcome in pediatric patients with rhabdomyosarcoma, also compared with conventional risk stratification. This warrants further validation in prospective trials as additional technique for risk stratification.