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

Neuroblastoma displays important clinical and genetic heterogeneity, with emergence of new mutations at tumor progression. 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. 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 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. .