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.

•

27 different immune subtypes form a complex neuroblastoma microenvironment

•

Inflammatory monocytes correlate with neuroblastoma disease progression

•

Active NK cells correlate with improved survival in human neuroblastoma

•

T cell exhaustion and Treg activity demonstrate no difference between risk groups

Abstract 6215: The immune cell atlas of human neuroblastoma

Our current knowledge of the different immune cells in neuroblastoma is based on in vitro and in vivo studies mainly focusing on a single cell type. Different studies have conveyed conflicting results and despite the use of anti-GD2 immunotherapy in the clinic, a comprehensive immune cell overview at the single-cell level is still missing and understanding the complete composition of immune cells neuroblastoma will be crucial for the development of novel immunotherapeutics against the disease. In this study, we performed single-cell RNA-sequencing on human neuroblastoma coupled with multiplex immunohistochemistry and survival analysis to provide a complete cellular and molecular immune cell landscape of human neuroblastoma. Further, we contrasted the neuroblastoma data with single-cell RNA-sequencing data from normal fetal adrenal gland to characterize cell-state changes from normal tissue to neuroblastoma. Our analysis revealed 27 immune cell subtypes including distinct subpopulations of myeloid, NK, B and T cells that were associated with a survival benefit. Multiple subtypes of B and NK cells, not seen in neuroblastoma before, were identified. We propose the presence of tertiary lymphoid structures and showed that Active NK cells correlated with improved survival. In contrast to adult cancers, we detected no difference in cytotoxicity and exhaustion score for cytotoxic T cells, nor Treg activity. However, we demonstrated an increase in inflammatory monocyte signature score from normal to tumor derived myeloid cells. Finally, we performed receptor-ligand interaction analysis between tumor, stroma and immune cells, where we highlight several interactions that we suggest for future studies of how to exploit immune cells as a therapeutic option in neuroblastoma. Our findings significantly broaden the understanding of the immune composition in neuroblastoma and provides a resource for the development of novel immunotherapeutics.

Citation Format: Bronte Manouk Verhoeven, Shenglin Mei, Thale K. Olsen, Karin K. Gustafsson, Anders Valind, Axel Lindström, David G. Nord, Shahrzad S. Fard, Catharina Hagerling, Peter V. Kharchenko, Per Kogner, John I. Johnsen, Ninib Baryawno. The immune cell atlas of human neuroblastoma [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 6215.

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.

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.

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.

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.

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.

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.

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.

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.

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.

Elevated tolerance to aneuploidy in cancer cells: estimating the fitness effects of chromosome number alterations by in silico modelling of somatic genome evolution

An unbalanced chromosome number (aneuploidy) is present in most malignant tumours and has been attributed to mitotic mis-segregation of chromosomes. However, recent studies have shown a relatively high rate of chromosomal mis-segregation also in non-neoplastic human cells, while the frequency of aneuploid cells remains low throughout life in most normal tissues. This implies that newly formed aneuploid cells are subject to negative selection in healthy tissues and that attenuation of this selection could contribute to aneuploidy in cancer. To test this, we modelled cellular growth as discrete time branching processes, during which chromosome gains and losses were generated and their host cells subjected to selection pressures of various magnitudes. We then assessed experimentally the frequency of chromosomal mis-segregation as well as the prevalence of aneuploid cells in human non-neoplastic cells and in cancer cells. Integrating these data into our models allowed estimation of the fitness reduction resulting from a single chromosome copy number change to an average of ≈30% in normal cells. In comparison, cancer cells showed an average fitness reduction of only 6% (p = 0.0008), indicative of aneuploidy tolerance. Simulations based on the combined presence of chromosomal mis-segregation and aneuploidy tolerance reproduced distributions of chromosome aberrations in >400 cancer cases with higher fidelity than models based on chromosomal mis-segregation alone. Reverse engineering of aneuploid cancer cell development in silico predicted that aneuploidy intolerance is a stronger limiting factor for clonal expansion of aneuploid cells than chromosomal mis-segregation rate. In conclusion, our findings indicate that not only an elevated chromosomal mis-segregation rate, but also a generalised tolerance to novel chromosomal imbalances contribute to the genomic landscape of human tumours.