Top left: schematic of patient enrollment, sample collection in the operating room, radiology, hospital and oncology clinic, plasma and CSF processing, DNA extraction, sequencing and analysis. Top right: a brain diagram with ring charts of diagnosis and age by tumor location: 4 supratentorial, 23 infratentorial and 1 spinal tumor. Bottom: pie charts of surgical resection extent, patient status, chemotherapy or immunotherapy, and radiation for the 28 patients.

Children with brain tumors are monitored mostly by MRI, and a scan cannot always show early enough that a tumor is growing back. Circulating tumor DNA (ctDNA) is the small fraction of cell-free DNA (cfDNA) in blood and cerebrospinal fluid (CSF) that comes from the tumor. It could add a molecular readout, but pediatric brain tumors shed little of it, and laboratories define a positive sample in different ways. In our new preprint with colleagues at Phoenix Children’s, led by Shea Gallus and Sharvari Mankame, we combine three ways of detecting ctDNA into one score.

We followed 28 children with medulloblastoma, diffuse midline glioma, high-grade glioma, pilocytic astrocytoma or ependymoma. The clinical team collected 140 plasma and 27 CSF samples during clinic visits, MRI appointments and surgery, alongside each child’s routine care.

We sequenced each tumor, each matched germline sample and each cfDNA sample. Each cfDNA sample then received a score from three readouts: tumor-informed mutations from targeted sequencing, copy number tumor fraction from shallow whole genome sequencing, and the ratio of short to long cfDNA fragments. We set each threshold using the matched germline samples, and we call a sample ctDNA-positive when the composite score reaches 2.

Top: a scoring table. Targeted sequencing scores 0 with no mutations, 1 for a non-coding mutation below 1% VAF, and 2 for a pathogenic mutation at any VAF or a non-coding mutation at 1% VAF or more. Copy number tumor fraction scores 0, 1 or 2 at cutoffs of 3% and 5%, relative to the matched germline sample. Fragmentomics scores 1 when the short to long fragment ratio exceeds 0.381. Bottom: stacked bars of total score for each cfDNA sample, grouped by patient, for medulloblastoma, DMG and high-grade glioma, and lower-grade tumors, with a dashed line at the ctDNA-positive threshold of 2.
The scoring rules (top) and the composite score of each cfDNA sample by patient and tumor type (bottom). Green, blue and orange show the contributions of mutations, copy number and fragment size.

The three readouts together identified 38 ctDNA-positive samples out of 159. Mutations alone identified 26, and copy number alone identified 20. Two children with medulloblastoma had no mutations we could follow, and in both, copy number profiling detected ctDNA in later samples. Fragment size did not add positive samples on its own, but it often confirmed a call made by the other two readouts.

CSF was the more sensitive fluid. We detected ctDNA in 70% of CSF samples and in 17% of plasma samples. Plasma is far easier to collect at routine visits, however, and it still carried prognostic information. Children with ctDNA-positive plasma after diagnosis had shorter progression-free survival (hazard ratio 4.5 after adjusting for tumor type) and shorter overall survival, although the overall survival association weakened after adjustment.

Top: a swimmer plot of 28 patients over days since diagnosis, colored by diagnosis, with plasma and CSF samples marked as ctDNA detected, not detected or indeterminate, and with surgery, chemotherapy and radiation indicated. Bottom left: bar chart of patients with at least one ctDNA-positive sample: medulloblastoma 9 of 11, DMG 3 of 7, high-grade glioma 2 of 5, ependymoma 1 of 3, pilocytic astrocytoma 0 of 3. Bottom right: Kaplan-Meier curves in which patients with ctDNA-positive plasma after diagnosis (n = 7) have shorter progression-free survival (log-rank P = 0.044) and overall survival (P = 0.036) than patients without (n = 19).
Sampling timelines for all 28 patients (top), ctDNA detection by tumor type (bottom left), and progression-free and overall survival by plasma ctDNA status after diagnosis (bottom right).

Two patients show what serial sampling can add. In a patient with a thalamic high-grade glioma, tumor variants were absent from plasma during radiation and targeted therapy. They appeared as the tumor progressed on MRI, together with a rise in total cfDNA yield.

One patient with high-grade glioma. A timeline shows one CSF and nine plasma draws from diagnosis to day 696. Four axial MRI scans show progression. The tumor copy number profile shows gains on chromosomes 1, 2, 7 and 9. A heatmap shows an EGFR missense variant and a C12orf43 3' UTR variant, present in the tumor, absent from early plasma samples and detected in plasma from day 323 onward. A line chart shows plasma cfDNA yield stable near 5 to 9 ng, then rising to about 22 ng at day 575.
Plasma ctDNA and total cfDNA yield rose with radiographic progression in a patient with high-grade glioma.

In a patient with diffuse midline glioma, CSF taken shortly after diagnosis carried the tumor’s mutations and its copy number profile. Plasma stayed negative until the fragment size signal rose late in the disease course, when leptomeningeal spread was visible on spinal MRI.

One patient with diffuse midline glioma. A timeline shows one CSF and four plasma draws. A heatmap shows seven tumor variants, including H3-3A K28M, TP53 and PIK3CA, detected in CSF at day 12 and absent from all plasma samples. A bar chart shows the short to long fragment ratio: 1.37 in CSF, 0.25 to 0.27 in the first three plasma samples, and 0.59 in plasma at day 209, above the cutoff of 0.381. Three sagittal spinal MRI scans show leptomeningeal metastases at four and six months. Copy number profiles of the tumor and the CSF closely match.
CSF reproduced the tumor's mutations and copy number profile at diagnosis. The plasma fragment ratio crossed the cutoff at day 209, when leptomeningeal dissemination was visible on MRI.

Together: combining mutations, copy number and fragment size detects more ctDNA than any one of them, CSF remains the most sensitive fluid, and plasma ctDNA after diagnosis identifies children at higher risk of progression. This is a pilot cohort from a single center, and we defined the thresholds within it. Our next step is prospective validation in larger, multi-institutional cohorts.

Huge credit to Shea and Sharvari for leading this work, and to Ross Mangum, who co-supervised it from Phoenix Children’s. We thank Adriana Marshall, Ben Posorske, Alyssa Hernandez, Lisa Keller, Shari Kyman, Kazuhiro Sabet, Patricia Cornejo, Tania Contente-Cuomo, Jonathan Keats, Peter Chiarelli, Jason Hauptman and Michael Berens, and above all the patients and families who took part. This research was supported by the Ivy Foundation and Students Supporting Brain Tumor Research, and sequencing was performed in the TGen Collaborative Sequencing Center.

Read more: the paper · the preprint on bioRxiv

Social media link: BlueSky · Twitter