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

Evaluation of simulation methods for tumor subclonal reconstruction

Most neoplastic tumors originate from a single cell, and their evolution can be genetically traced through lineages characterized by common alterations such as small somatic mutations (SSMs), copy number alterations (CNAs), structural variants (SVs), and aneuploidies. Due to the complexity of these alterations in most tumors and the errors introduced by sequencing protocols and calling algorithms, tumor subclonal reconstruction algorithms are necessary to recapitulate the DNA sequence composition and tumor evolution in silico. With a growing number of these algorithms available, there is a pressing need for consistent and comprehensive benchmarking, which relies on realistic tumor sequencing generated by simulation tools. Here, we examine the current simulation methods, identifying their strengths and weaknesses, and provide recommendations for their improvement. Our review also explores potential new directions for research in this area. This work aims to serve as a resource for understanding and enhancing tumor genomic simulations, contributing to the advancement of the field.

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.

CNETML: maximum likelihood inference of phylogeny from copy number profiles of multiple samples

Phylogenetic trees based on copy number profiles from multiple samples of a patient are helpful to understand cancer evolution. Here, we develop a new maximum likelihood method, CNETML, to infer phylogenies from such data. CNETML is the first program to jointly infer the tree topology, node ages, and mutation rates from total copy numbers of longitudinal samples. Our extensive simulations suggest CNETML performs well on copy numbers relative to ploidy and under slight violation of model assumptions. The application of CNETML to real data generates results consistent with previous discoveries and provides novel early copy number events for further investigation.

Evolutionary route of nasopharyngeal carcinoma metastasis and its clinical significance

It is critical to understand factors associated with nasopharyngeal carcinoma (NPC) metastasis. To track the evolutionary route of metastasis, here we perform an integrative genomic analysis of 163 matched blood and primary, regional lymph node metastasis and distant metastasis tumour samples, combined with single-cell RNA-seq on 11 samples from two patients. The mutation burden, gene mutation frequency, mutation signature, and copy number frequency are similar between metastatic tumours and primary and regional lymph node tumours. There are two distinct evolutionary routes of metastasis, including metastases evolved from regional lymph nodes (lymphatic route, 61.5%, 8/13) and from primary tumours (hematogenous route, 38.5%, 5/13). The hematogenous route is characterised by higher IFN-γ response gene expression and a higher fraction of exhausted CD8⁺ T cells. Based on a radiomics model, we find that the hematogenous group has significantly better progression-free survival and PD-1 immunotherapy response, while the lymphatic group has a better response to locoregional radiotherapy.

Clinically relevant treatment of PDX models reveals patterns of neuroblastoma chemoresistance

Chemotherapy resistance and relapses are common in high-risk neuroblastoma (NB). Here, we developed a clinically relevant in vivo treatment protocol mimicking the first-line five-chemotherapy treatment regimen of high-risk NB and applied this protocol to mice with MYCN-amplified NB patient-derived xenografts (PDXs). Genomic and transcriptomic analyses were used to reveal NB chemoresistance mechanisms. Intrinsic resistance was associated with high genetic diversity and an embryonic phenotype. Relapsed NB with acquired resistance showed a decreased adrenergic phenotype and an enhanced immature mesenchymal-like phenotype, resembling multipotent Schwann cell precursors. NBs with a favorable treatment response presented a lineage-committed adrenergic phenotype similar to normal neuroblasts. Novel integrated phenotypic gene signatures reflected treatment response and patient prognosis. NB organoids established from relapsed PDX tumors retained drug resistance, tumorigenicity, and transcriptional cell states. This work sheds light on the mechanisms of NB chemotherapy response and emphasizes the importance of transcriptional cell states in chemoresistance.

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.