Early evolutionary branching across spatial domains predisposes to clonal replacement under chemotherapy in neuroblastoma.

Neuroblastoma (NB) is one of the most lethal childhood cancers due to its propensity to become treatment resistant. By spatial mapping of subclone geographies before and after chemotherapy across 89 tumor regions from 12 NBs, we find that densely packed territories of closely related subclones present at diagnosis are replaced under effective treatment by islands of distantly related survivor subclones, originating from a different most recent ancestor compared to lineages dominating before treatment. Conversely, in tumors that progressed under treatment, ancestors of subclones dominating later in disease are present already at diagnosis. Chemotherapy treated xenografts and cell culture models replicate these two contrasting scenarios and show branching evolution to be a constant feature of proliferating NB cells. Phylogenies based on whole genome sequencing of 505 individual NB cells indicate that a rich repertoire of parallel subclones emerges already with the first oncogenic mutations and lays the foundation for clonal replacement under treatment.

A Gradual Transition Toward Anaplasia in Wilms Tumor Through Tolerance to Genetic Damage.

Patients with Wilms tumor (WT) in general have excellent survival, but the prognosis of patients belonging to the subgroup of WT with diffuse anaplasia (DA) is poor due to frequent resistance to chemotherapy. We hypothesized that DA WT cells might undergo changes, such as acquiring a persistent tolerance to DNA damage and copy number aberrations (CNAs), which could eventually lead to their resistance to chemotherapy treatment. Tissue sections from chemotherapy-treated DA WTs (n = 12) were compared with chemotherapy-treated nonanaplastic WTs (n = 15) in a tissue microarray system, enabling analysis of 769 tumor regions. All regions were scored for anaplastic features and immunohistochemistry was used to quantify p53 expression, proliferation index (Ki67), and DNA double-strand breaks (γH2AX). CNAs were assessed by array-based genotyping and TP53 mutations using targeted sequencing. Proliferation index and the frequency of DNA double-strand breaks (γH2AX dot expression) increased with higher anaplasia scores. Almost all (95.6%) areas with full-scale anaplasia had TP53 mutations or loss of heterozygosity, along with an increased amount of CNAs. Interestingly, areas with wild-type TP53 with loss of heterozygosity and only one feature of anaplasia (anaplasia score 1) also had significantly higher proliferation indices, more DNA double-strand breaks, and more CNAs than regions without any anaplastic features (score 0); such areas may be preanaplastic cell populations under selective pressure for TP53 mutations. In conclusion, we suggest that chemoresistance of DA WTs may be partly explained by a high proliferative capability of anaplastic cells, which also have a high burden of double-stranded DNA breaks and CNAs, and that there is a gradual emergence of anaplasia in WT.

Resolving the Pathogenesis of Anaplastic Wilms Tumors through Spatial Mapping of Cancer Cell Evolution.

PURPOSE: While patients with intermediate-risk (IR) Wilms tumors now have an overall survival (OS) rate of almost 90%, those affected by high-stage tumors with diffuse anaplasia have an OS of only around 50%. We here identify key events in the pathogenesis of diffuse anaplasia by mapping cancer cell evolution over anatomic space in Wilms tumors. EXPERIMENTAL DESIGN: We spatially mapped subclonal landscapes in a retrospective cohort of 20 Wilms tumors using high-resolution copy-number profiling and TP53 mutation analysis followed by clonal deconvolution and phylogenetic reconstruction. Tumor whole-mount sections (WMS) were utilized to characterize the distribution of subclones across anatomically distinct tumor compartments. RESULTS: Compared with non-diffuse anaplasia Wilms tumors, tumors with diffuse anaplasia showed a significantly higher number of genetically distinct tumor cell subpopulations and more complex phylogenetic trees, including high levels of phylogenetic species richness, divergence, and irregularity. All regions with classical anaplasia showed TP53 alterations. TP53 mutations were frequently followed by saltatory evolution and parallel loss of the remaining wild-type (WT) allele in different regions. Morphologic features of anaplasia increased with copy-number aberration (CNA) burden and regressive features. Compartments demarcated by fibrous septae or necrosis/regression were frequently (73%) associated with the emergence of new clonal CNAs, although clonal sweeps were rare within these compartments. CONCLUSIONS: Wilms tumors with diffuse anaplasia display significantly more complex phylogenies compared with non-diffuse anaplasia Wilms tumors, including features of saltatory and parallel evolution. The subclonal landscape of individual tumors was constrained by anatomic compartments, which should be considered when sampling tissue for precision diagnostics.

Macrophage infiltration promotes regrowth in MYCN-amplified neuroblastoma after chemotherapy.

Despite aggressive treatment, the 5-year event-free survival rate for children with high-risk neuroblastoma is <50%. While most high-risk neuroblastoma patients initially respond to treatment, often with complete clinical remission, many eventually relapse with therapy-resistant tumors. Novel therapeutic alternatives that prevent the recurrence of therapy-resistant tumors are urgently needed. To understand the adaptation of neuroblastoma under therapy, we analyzed the transcriptomic landscape in 46 clinical tumor samples collected before (PRE) or after (POST) treatment from 22 neuroblastoma patients. RNA sequencing revealed that many of the top-upregulated biological processes in POST MYCN amplified (MNA+) tumors compared to PRE MNA+ tumors were immune-related, and there was a significant increase in numerous genes associated with macrophages. The infiltration of macrophages was corroborated by immunohistochemistry and spatial digital protein profiling. Moreover, POST MNA+ tumor cells were more immunogenic compared to PRE MNA+ tumor cells. To find support for the macrophage-induced outgrowth of certain subpopulations of immunogenic tumor cells following treatment, we examined the genetic landscape in multiple clinical PRE and POST tumor samples from nine neuroblastoma patients revealing a significant correlation between an increased amount of copy number aberrations (CNA) and macrophage infiltration in POST MNA+ tumor samples. Using an in vivo neuroblastoma patient-derived xenograft (PDX) chemotherapy model, we further show that inhibition of macrophage recruitment with anti-CSF1R treatment prevents the regrowth of MNA+ tumors following chemotherapy. Taken together, our work supports a therapeutic strategy for fighting the relapse of MNA+ neuroblastoma by targeting the immune microenvironment.

The immune cell atlas of human neuroblastoma.

Understanding the complete immune cell composition of human neuroblastoma (NB) is crucial for the development of immunotherapeutics. Here, we perform single-cell RNA sequencing (scRNA-seq) on 19 human NB samples coupled with multiplex immunohistochemistry, survival analysis, and comparison with normal fetal adrenal gland data. We provide a comprehensive immune cell landscape and characterize cell-state changes from normal tissue to NB. Our analysis reveals 27 immune cell subtypes, including distinct subpopulations of myeloid, NK, B, and T cells. Several different cell types demonstrate a survival benefit. In contrast to adult cancers and previous NB studies, we show an increase in inflammatory monocyte cell state when contrasting normal and tumor tissue, while no differences in cytotoxicity and exhaustion score for T cells, nor in Treg activity, are observed. Our receptor-ligand interaction analysis reveals a highly complex interactive network of the NB microenvironment from which we highlight several interactions that we suggest for future therapeutic studies.

Branching Copy-Number Evolution and Parallel Immune Profiles across the Regional Tumor Space of Resected Pancreatic Cancer.

Pancreatic ductal adenocarcinoma (PDAC) remains a highly lethal disease. The only option for curative treatment is resection of the tumor followed by standard adjuvant chemotherapy. Yet, early relapse due to chemoresistance is almost inevitable. Herein, we delineated the genetic intratumor heterogeneity in resected PDAC, with the aim to identify evolutionary patterns that may be associated with overall survival (OS) following treatment with curative intent. Potential relationships with the adjacent immune microenvironment were also examined. The genetic and immune landscapes of the regional tumor space were analyzed in nine patients with resected PDAC. Targeted deep sequencing and genome wide SNP array were followed by clonal deconvolution and phylogenetic analysis. A mathematical complexity score was developed to calculate the network extent of each phylogeny. Spatial variation in abundancy and tumor nest infiltration of immune cells was analyzed by double IHC staining. Copy-number heterogeneity was denoted as the major contributing factor to the branching architectures of the produced phylogenetic trees. Increased tree complexity was significantly inversely associated with OS, and larger regional maximum aberrations (higher treetops) were associated with increased PD-L1 expression on tumor cells. Contrastingly, an FREM1 gene amplification, found in one patient, coincided with a particularly vigorous immune response. Findings from this limited case series suggest that complex evolutionary patterns may be associated with a shorter survival in surgically treated patients with PDAC. Some hypothesis-generating associations with the surrounding immune microenvironment were also detected. IMPLICATIONS: Evolutionary copy-number patterns may be associated with survival in patients with resected PDAC.

DEVOLUTION-A method for phylogenetic reconstruction of aneuploid cancers based on multiregional genotyping data.

Phylogenetic reconstruction of cancer cell populations remains challenging. There is a particular lack of tools that deconvolve clones based on copy number aberration analyses of multiple tumor biopsies separated in time and space from the same patient. This has hampered investigations of tumors rich in aneuploidy but few point mutations, as in many childhood cancers and high-risk adult cancer. Here, we present DEVOLUTION, an algorithm for subclonal deconvolution followed by phylogenetic reconstruction from bulk genotyping data. It integrates copy number and sequencing information across multiple tumor regions throughout the inference process, provided that the mutated clone fraction for each mutation is known. We validate DEVOLUTION on data from 56 pediatric tumors comprising 253 tumor biopsies and show a robust performance on simulations of bulk genotyping data. We also benchmark DEVOLUTION to similar bioinformatic tools using an external dataset. DEVOLUTION holds the potential to facilitate insights into the development, progression, and response to treatment, particularly in tumors with high burden of chromosomal copy number alterations.

Tracing the evolution of aneuploid cancers by multiregional sequencing with CRUST.

Clonal deconvolution of mutational landscapes is crucial to understand the evolutionary dynamics of cancer. Two limiting factors for clonal deconvolution that have remained unresolved are variation in purity and chromosomal copy number across different samples of the same tumor. We developed a semi-supervised algorithm that tracks variant calls through multi-sample spatiotemporal tumor data. While normalizing allele frequencies based on purity, it also adjusts for copy number changes at clonal deconvolution. Absent à priori copy number data, it renders in silico copy number estimations from bulk sequences. Using published and simulated tumor sequences, we reliably segregated clonal/subclonal variants even at a low sequencing depth (~50×). Given at least one pure tumor sample (>70% purity), we could normalize and deconvolve paired samples down to a purity of 40%. This renders a reliable clonal reconstruction well adapted to multi-regionally sampled solid tumors, which are often aneuploid and contaminated by non-cancer cells.

Immune checkpoint inhibitors in Wilms' tumor and Neuroblastoma: What now?

BACKGROUND: Therapeutic activation of tumor-infiltrating lymphocytes using monoclonal antibodies targeting PD1 or PD-L1 (immune checkpoint inhibitors-ICIs) has revolutionized treatment of specific solid tumors in adult cancer patients, and much hope has been placed on a similar effect in relapsed or refractory solid pediatric tumors. Recent clinical trials have disappointingly shown an almost nonexistent response rate, while case reports have demonstrated that some pediatric patients do achieve durable responses when treated with this type of drug. AIM: To elucidate this paradox, we mapped the landscape of expressed neoantigens as well as the levels of immune cell infiltration in the two most common extracranial solid pediatric tumors: Wilms tumor and neuroblastoma using state-of-the-art in silico analysis of a large cohort of patients with these tumors. METHODS: By integration of whole exome sequencing and RNA-sequencing, we mapped the landscape of neoantigens in the TARGET cohorts for these diagnoses and correlated these findings with known genetic prognostic markers. RESULTS: Our analysis shows that these tumors typically have much lower levels of expressed neoantigens than commonly seen in adult cancers, but we also identify subgroups with significantly higher levels of neoantigens. For neuroblastomas, the cases with higher levels of neoantigens were confined to the group without MYCN-amplification and for Wilms tumor restricted to the TP53-mutated cases. Furthermore, we demonstrate that neuroblastomas have an unexpectedly high level of CD8+ tumor-infiltrating lymphocytes, even when compared to adult tumor types where ICI is an approved treatment. CONCLUSION: These results could be important to consider when designing future clinical trials of ICI treatment in pediatric patients with relapsed or refractory solid tumors.

A dynamic mutational landscape associated with an inter-regionally diverse immune response in malignant rhabdoid tumour.

Malignant rhabdoid tumour (MRT) is a childhood neoplasm of high malignancy characterised by biallelic mutation and/or loss of the epigenetic master regulator SMARCB1, accompanied by no or few other oncogenic drivers. In spite of their generally low mutational burden, an intratumoural T-cell response has been reported in a subset of MRTs, indicating that immune checkpoint inhibition may be considered a viable therapy option for some patients. We assess here the evolution over time and space of predicted neoantigens and indicators of immune checkpoint status in two MRT patients who progressed under treatment. Both patients showed an accumulation of novel clonal and subclonal mutations, including predicted neoantigens, in metastases compared to their inferred ancestral clones in the primary tumours. The first patient had peritoneal metastases from an MRT of the liver. Clonal deconvolution revealed polyclonal seeding from the primary tumour to a single metastatic site, followed by a local subclonal burst of mutations. The second patient had a renal MRT with multiple pulmonary metastases, each of which could be traced back to a single genetically unique founder cell, with formation of novel subclones in two metastases. Both patients showed a regionally heterogeneous landscape of predicted neoantigens and of tumour-infiltrating lymphocytes expressing CD8 and PD1. In both patients, some tumour regions fulfilled established criteria for PD-L1 positivity (> 1% of tumour cells), while others did not. This suggests that even in a tumour type like MRT, with a single driver mutation, there can be heterogeneity in neoantigen repertoire, immune response, and biomarkers for checkpoint blockade among sampled locations. This must be taken into account when assessing progressed MRT patients for checkpoint inhibition therapy. © 2020 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of Pathological Society of Great Britain and Ireland.

Extensive Clonal Branching Shapes the Evolutionary History of High-Risk Pediatric Cancers.

Darwinian evolution of tumor cells remains underexplored in childhood cancer. We here reconstruct the evolutionary histories of 56 pediatric primary tumors, including 24 neuroblastomas, 24 Wilms tumors, and 8 rhabdomyosarcomas. Whole-genome copy-number and whole-exome mutational profiling of multiple regions per tumor were performed, followed by clonal deconvolution to reconstruct a phylogenetic tree for each tumor. Overall, 88% of the tumors exhibited genetic variation among primary tumor regions. This variability typically emerged through collateral phylogenetic branching, leading to spatial variability in the distribution of more than 50% (96/173) of detected diagnostically informative genetic aberrations. Single-cell sequencing of 547 individual cancer cells from eight solid pediatric tumors confirmed branching evolution to be a fundamental underlying principle of genetic variation in all cases. Strikingly, cell-to-cell genetic diversity was almost twice as high in aggressive compared with clinically favorable tumors (median Simpson index of diversity 0.45 vs. 0.88; P = 0.029). Similarly, a comparison of multiregional sampling data from a total of 274 tumor regions showed that new phylogenetic branches emerge at a higher frequency per sample and carry a higher mutational load in high-risk than in low-risk tumors. Timelines based on spatial genetic variation showed that the mutations most influencing relapse risk occur at initiation of clonal expansion in neuroblastoma and rhabdomyosarcoma, whereas in Wilms tumor, they are late events. Thus, from an evolutionary standpoint, some high-risk childhood cancers are born bad, whereas others grow worse over time. SIGNIFICANCE: Different pediatric cancers with a high risk of relapse share a common generic pattern of extensively branching evolution of somatic mutations. GRAPHICAL ABSTRACT: http://cancerres.aacrjournals.org/content/canres/80/7/1512/F1.large.jpg.

Neuroblastoma with flat genomic profile: a question of representativity?

Neuroblastoma is one of the most common paediatric malignancies. Detection of somatic genetic alterations in this tumour is instrumental for its risk stratification and treatment. On the other hand, an absence of detected chromosomal imbalances in neuroblastoma biopsies is difficult to interpret because it is unclear whether this situation truly reflects the tumour genome or if it is due to suboptimal sampling. We here present a neuroblastoma in the left adrenal of a newborn. The tumour was subjected to single-nucleotide polymorphism array analysis of five tumour regions with >80% tumour cells in histological mirror sections. This revealed no aberrations compared with a normal reference sample from the patient. Whole exome sequencing identified two single-nucleotide variants present in most tumour regions, corroborating that the tumour resulted from monoclonal expansion. Our data provide proof-of-principle that rare cases of neuroblastoma can have a normal whole genome copy number and allelic profile.

Convergent evolution of 11p allelic loss in multifocal Wilms tumors arising in WT1 mutation carriers.

Wilms tumors in patients with constitutional WT1 mutations are examples of Knudson's tumor suppressor paradigm, with somatic inactivation of the second allele occurring through 11p loss of heterozygosity. The time point of this second hit has remained unknown. We analyzed seven Wilms tumors from two patients with constitutional WT1 mutations by whole exome sequencing and genomic array. All tumors exhibited wild type WT1 loss through uniparental isodisomy. Each tumor had a unique genomic breakpoint in 11p, typically accompanied by a private activating mutation of CTNNB1. Hence, convergent evolution rather than field carcinogenesis underlies multifocal tumors in WT1 mutation carriers.

Four evolutionary trajectories underlie genetic intratumoral variation in childhood cancer.

A major challenge to personalized oncology is that driver mutations vary among cancer cells inhabiting the same tumor. Whether this reflects principally disparate patterns of Darwinian evolution in different tumor regions has remained unexplored1-5. We mapped the prevalence of genetically distinct clones over 250 regions in 54 childhood cancers. This showed that primary tumors can simultaneously follow up to four evolutionary trajectories over different anatomic areas. The most common pattern consists of subclones with very few mutations confined to a single tumor region. The second most common is a stable coexistence, over vast areas, of clones characterized by changes in chromosome numbers. This is contrasted by a third, less frequent, pattern where a clone with driver mutations or structural chromosome rearrangements emerges through a clonal sweep to dominate an anatomical region. The fourth and rarest pattern is the local emergence of a myriad of clones with TP53 inactivation. Death from disease was limited to tumors exhibiting the two last, most dynamic patterns.

The fetal thymus has a unique genomic copy number profile resulting from physiological T cell receptor gene rearrangement.

Somatic mosaicism, the presence of genetically distinct cells within an organism, has been increasingly associated with human morbidity, ranging from being a cause of rare syndromes to a risk factor for common disorders such as malignancy and cardiovascular disease. Previous studies interrogating the normal prevalence of somatic mosaicism have focused on adults. We here present an estimate of the baseline frequency of somatic mosaic copy number variation (CNV) at the time around birth, by sampling eight different organs from a total of five fetuses and newborns. Overall we find a significantly lower frequency of organ specific (i.e. mosaic) CNVs as compared to adults (p = 0.003; Mann-Whitney U-test). The rate of somatic CNV in adults has been estimated to around 2.2 CNV per organ assayed. In contrast, after stringent filtering, we found no organ-private CNVs in fetuses or newborns with exception of the thymus. This organ exhibited a specific genome profile in the form of deletions resulting from polyclonal T-cell receptor rearrangements. This implies that somatic non-immune related CNVs, if present at birth, are typically confined to very small cell populations within organs.

Aberrant epigenetic regulation in clear cell sarcoma of the kidney featuring distinct DNA hypermethylation and EZH2 overexpression.

The global methylation profile and the mutational status of 633 specific epigenetic regulators were analyzed in the pediatric tumor clear cell sarcoma of the kidney (CCSK). Methylation array analyses of 30 CCSKs revealed CCSK tumor DNA to be globally hypermethylated compared to Wilms tumor, normal fetal kidney, and adult kidney. The aberrant methylation pattern of CCSKs was associated with activation of genes involved in embryonic processes and with silencing of genes linked to normal kidney function. No epigenetic regulator was recurrently mutated in our cohort, but a mutation in the key epigenetic regulator EZH2 was discovered in one case. EZH2 mRNA was significantly higher in CCSK compared to Wilms tumor and normal kidney, and the EZH2 protein was strongly expressed in more than 90 % of CCSK tumor cells in 9/9 tumors analyzed. This was in striking contrast to the lack of EZH2 protein expression in Wilms tumor stromal elements, indicating that EZH2 could be explored further as a diagnostic marker and a potential drug target for CCSK.

BCOR internal tandem duplication and YWHAE-NUTM2B/E fusion are mutually exclusive events in clear cell sarcoma of the kidney.

Clear cell sarcoma of the kidney (CCSK) is the second most common pediatric renal tumor. Two recurrent genetic aberrations have been described in CCSK. One is a fusion of YWHAE and NUTM2B/E, the other is an internal tandem duplication (ITD) in the BCOR gene. Here it is shown that YWHAE-NUTM2B/E fusion and the BCOR ITD are mutually exclusive events and activated different downstream signaling systems. This has important diagnostic implications and opens up for further mechanistic studies of CCSK pathogenesis.

Intratumoral genome diversity parallels progression and predicts outcome in pediatric cancer.

Genetic differences among neoplastic cells within the same tumour have been proposed to drive cancer progression and treatment failure. Whether data on intratumoral diversity can be used to predict clinical outcome remains unclear. We here address this issue by quantifying genetic intratumoral diversity in a set of chemotherapy-treated childhood tumours. By analysis of multiple tumour samples from seven patients we demonstrate intratumoral diversity in all patients analysed after chemotherapy, typically presenting as multiple clones within a single millimetre-sized tumour sample (microdiversity). We show that microdiversity often acts as the foundation for further genome evolution in metastases. In addition, we find that microdiversity predicts poor cancer-specific survival (60%; P=0.009), independent of other risk factors, in a cohort of 44 patients with chemotherapy-treated childhood kidney cancer. Survival was 100% for patients lacking microdiversity. Thus, intratumoral genetic diversity is common in childhood cancers after chemotherapy and may be an important factor behind treatment failure.

Activation of human telomerase reverse transcriptase through gene fusion in clear cell sarcoma of the kidney.

Clear cell sarcoma of the kidney (CCSK) is a rare tumor type affecting infants and young children. Most CCSKs display few genomic aberrations, and no general underlying mechanism for tumor initiation has yet been identified, although a YWHAE-NUTM2B/NUTM2E fusion gene has been observed in a minority of cases. We performed RNA-sequencing of 22 CCSKs to investigate the presence of additional fusion transcripts. The presence of the YWHAE-NUTM2B/NUTM2E fusion was confirmed in two cases. In addition, a novel IRX2-TERT fusion transcript was identified in one case. SNP-array analyses revealed the underlying event to be an interstitial deletion in the short arm of chromosome 5 (5p15.33). TERT was dramatically upregulated under the influence of the IRX2 promoter. In line with TERT expression being driven by active IRX2 regulatory elements, we found a high expression of IRX2 in CCSKs irrespective of fusion gene status. IRX2 was also expressed in human fetal kidney - the presumed tissue of origin for CCSK. We conclude that in addition to promoter mutations and epigenetic events, TERT can also be activated in tumors via formation of fusion transcripts.