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.

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.

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.

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.

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.

Whole chromosome gain does not in itself confer cancer-like chromosomal instability.

Constitutional aneuploidy is typically caused by a single-event meiotic or early mitotic error. In contrast, somatic aneuploidy, found mainly in neoplastic tissue, is attributed to continuous chromosomal instability. More debated as a cause of aneuploidy is aneuploidy itself; that is, whether aneuploidy per se causes chromosomal instability, for example, in patients with inborn aneuploidy. We have addressed this issue by quantifying the level of somatic mosaicism, a proxy marker of chromosomal instability, in patients with constitutional aneuploidy by precise background-filtered dual-color FISH. In contrast to previous studies that used less precise methods, we find that constitutional trisomy, even for large chromosomes that are often trisomic in cancer, does not confer a significantly elevated rate of somatic chromosomal mosaicism in individual cases. Constitutional triploidy was associated with an increased level of somatic mosaicism, but this consisted mostly of reversion from trisomy to disomy and did not correspond to a proportionally elevated level of chromosome mis-segregation in triploids, indicating that the observed mosaicism resulted from a specific accumulation of cells with a hypotriploid chromosome number. In no case did the rate of somatic mosaicism in constitutional aneuploidy exceed that of "chromosomally stable" cancer cells. Our findings show that even though constitutional aneuploidy was in some cases associated with low-level somatic mosaicism, it was insufficient to generate the cancer-like levels expected if aneuploidy single-handedly triggered cancer-like chromosomal instability.

Confined trisomy 8 mosaicism of meiotic origin: a rare cause of aneuploidy in childhood cancer.

Whether chromosome abnormalities observed in tumor cells may in some cases reflect low-grade somatic mosaicism for anomalies present already at zygote formation, rather than acquired somatic mutations, has for long remained a speculation. We here report a patient with Wilms tumor, where constitutional somatic mosaicism of trisomy 8 was detected in a previously healthy 2 ½-year-old boy. Single Nucleotide Polymorphism (SNP) array analysis of tumor tissue revealed a complex distribution of allele frequencies for chromosome 8 that could not be explained solely by mitotic events. Combined analysis of allele frequencies, chromosome banding, and fluorescence in situ hybridization revealed that the majority of tumor cells contained four copies of chromosome 8, with three distinct haplotypes at a 2:1:1 ratio. Because the patient had not been subject to organ transplantation, these findings indicated that the tumor karyotype evolved from a cell with trisomy 8 of meiotic origin, with subsequent somatic gain of one additional chromosome copy. Haplotype analysis was consistent with trisomy 8 through nondisjunction at meiosis I. Matched normal renal tissue or peripheral blood did not contain detectable amounts of cells with trisomy 8, consistent with the complete lack of mosaic trisomy 8 syndrome features in the patient. This case provides proof of principle for the hypothesis that tumor genotypes may in rare cases reflect meiotic rather than mitotic events, also in patients lacking syndromic features. © 2014 Wiley Periodicals, Inc.

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.

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.

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.

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.

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.

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.

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.

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.