Visualizing tumor evolution with the fishplot package for R

A patient-derived T cell lymphoma biorepository uncovers pathogenetic mechanisms and host-related therapeutic vulnerabilities

Peripheral T cell lymphomas (PTCLs) comprise heterogeneous malignancies with limited therapeutic options. To uncover targetable vulnerabilities, we generate a collection of PTCL patient-derived tumor xenografts (PDXs) retaining histomorphology and molecular donor-tumor features over serial xenografting. PDX demonstrates remarkable heterogeneity, complex intratumor architecture, and stepwise trajectories mimicking primary evolutions. Combining functional transcriptional stratification and multiparametric imaging, we identify four distinct PTCL microenvironment subtypes with prognostic value. Mechanistically, we discover a subset of PTCLs expressing Epstein-Barr virus-specific T cell receptors and uncover the capacity of cancer-associated fibroblasts of counteracting treatments. PDXs' pre-clinical testing captures individual vulnerabilities, mirrors donor patients' clinical responses, and defines effective patient-tailored treatments. Ultimately, we assess the efficacy of CD5KO- and CD30- Chimeric Antigen Receptor T Cells (CD5KO-CART and CD30_CART, respectively), demonstrating their therapeutic potential and the synergistic role of immune checkpoint inhibitors for PTCL treatment. This repository represents a resource for discovering and validating intrinsic and extrinsic factors and improving the selection of drugs/combinations and immune-based therapies.

Evolutionary Pressures Shape Undifferentiated Pleomorphic Sarcoma Development and Radiotherapy Response

Radiotherapy is an integral component in the treatment of many types of cancer, with approximately half of patients with cancer receiving radiotherapy. Systemic therapy applies pressure that can select for resistant tumor subpopulations, underscoring the importance of understanding how radiation impacts tumor evolution to improve treatment outcomes. We integrated temporal genomic profiling of 120 spatially distinct tumor regions from 20 patients with undifferentiated pleomorphic sarcomas (UPS), longitudinal circulating tumor DNA analysis, and evolutionary biology computational pipelines to study UPS evolution during tumorigenesis and in response to radiotherapy. Most unirradiated UPSs displayed initial linear evolution, followed by subsequent branching evolution with distinct mutational processes during early and late development. Metrics of genetic divergence between regions provided evidence of strong selection pressures during UPS development that further increased during radiotherapy. Subclone abundance changed after radiotherapy with subclone contraction tied to alterations in calcium signaling, and inhibiting calcium transporters radiosensitized sarcoma cells. Finally, circulating tumor DNA analysis accurately measured subclone abundance and enabled noninvasive monitoring of subclonal changes. These results demonstrate that radiation exerts selective pressures on UPSs and suggest that targeting radioresistant subclonal populations could improve outcomes after radiotherapy. Significance: Radiotherapy mediates tumor evolution by leading to the expansion of resistant subclonal cancer cell populations, indicating that developing approaches to target resistant subclones will be crucial to improve radiotherapy response.

Jellyfish: integrative visualization of spatio-temporal tumor evolution and clonal dynamics

SUMMARY

Spatial and temporal intra-tumor heterogeneity drives tumor evolution and therapy resistance. Existing visualization tools often fail to capture both dimensions simultaneously. To address this, we developed Jellyfish, a tool that integrates phylogenetic and sample trees into a single plot, providing a holistic view of tumor evolution and capturing both spatial and temporal evolution. Available as a JavaScript library and R package, Jellyfish generates interactive visualizations from tumor phylogeny and clonal composition data. We demonstrate its ability to visualize complex subclonal dynamics using data from ovarian high-grade serous carcinoma.

AVAILABILITY AND IMPLEMENTATION

Jellyfish is freely available with MIT license at https://github.com/HautaniemiLab/jellyfish (JavaScript library) and https://github.com/HautaniemiLab/jellyfisher (R package).

Clonal Evolution of Myeloid Malignancies Treated With Microtransplantation: A Single-Centre Experience

Microtransplantation is a cellular therapy used in acute myeloid leukaemia and myelodysplastic syndromes as a maintenance therapy in patients ineligible for a regular allogeneic stem cell transplantation. We performed a monocentric retrospective study of acute myeloid leukaemia, myelodysplastic syndromes, and chronic myelomonocytic leukaemia patients who underwent microtransplantations at Nice University Hospital. We analysed the evolution of the disease mutational status after microtransplantation. We report 18 patients who underwent microtransplantation courses, with a total of 47 microtransplantations performed between February 2020 and June 2022. We observed long-term remissions even in high-risk patients. Founder mutations persisted throughout the follow-up, whereas it was more variable for other nonfounder mutations, with most of the nonfounder mutations variant allele frequency decreasing over time. Safety data were reassuring; no graft versus host disease was recorded, and cytokine release syndromes were manageable. Relapses or progressions were associated with the emergence or increase of a TP53 mutated clone. Microtransplantation is a promising therapy for patients ineligible for regular allogeneic stem transplantation. Further larger and randomised studies are required to establish its place as a maintenance therapy in myeloid malignancies.

Effects of local mutations in quadratic iterations

We introduce mutations in the process of discrete iterations of complex quadratic maps in the family fc(z)=z2+c. More specifically, we consider a "correct" function fc1 acting on the complex plane. A "mutation" fc0 is a different ("erroneous") map acting on a locus of given radius r around a mutation focal point ξ∗. The effect of the mutation is interpolated radially to eventually recover the original map fc1 when reaching an outer radius R. We call the resulting map a "mutated" map. In the theoretical framework of mutated iterations, we study how a mutation affects the temporal evolution of the system and the asymptotic behavior of its orbits. We use the prisoner set of the system to quantify simultaneously the long-term behavior of the entire space under mutated maps. We analyze how the position, timing, and size of the mutation can alter the system's long-term evolution (as encoded in the topology of its prisoner set). The framework is then discussed as a metaphoric model for studying the impact of copying errors in natural replication systems.

Evolution of genome and immunogenome in esophageal squamous cell carcinomas driven by neoadjuvant chemoradiotherapy

Neoadjuvant chemoradiotherapy (NCRT) followed by surgery is a standard treatment for locally advanced esophageal squamous cell carcinomas (ESCCs). However, the evolution of genome and immunogenome in ESCCs driven by NCRT remains incompletely elucidated. We performed whole-exome sequencing of 51 ESCC tumors collected before and after NCRT, 36 of which were subjected to transcriptome sequencing. Clonal analysis identified clonal extinction in 13 ESCC patients wherein all pre-NCRT clones disappeared after NCRT, and clonal persistence in 9 patients wherein clones endured following NCRT. The clone-persistent patients showed higher pre-NCRT genomic intratumoral heterogeneity and worse prognosis than the clone-extinct ones. In contrast to the clone-extinct patients, the clone-persistent patients demonstrated a high proportion of subclonal neoantigens within pre-treatment specimens. Transcriptome analysis revealed increased immune infiltrations and up-regulated immune-related pathways after NCRT, especially in the clone-extinct patients. The number of T cell receptor-neoantigen interactions was higher in the clone-extinct patients than in the clone-persistent ones. The decrease in T cell repertoire evenness positively correlated to the decreased number of clonal neoantigens after NCRT, especially in the clone-extinct patients. In conclusion, we identified two prognosis-related clonal dynamic modes driven by NCRT in ESCCs. This study extended our knowledge of the ESCC genome and immunogenome evolutions driven by NCRT.

Tumor architecture and emergence of strong genetic alterations are bottlenecks for clonal evolution in primary prostate cancer

Prostate cancer (PCA) exhibits high levels of intratumoral heterogeneity. In this study, we developed a mathematical model to study the growth and genetic evolution of PCA. We explored the possible evolutionary patterns and demonstrated that tumor architecture represents a major bottleneck for divergent clonal evolution. Early consecutive acquisition of strong genetic alterations serves as a proxy for the formation of aggressive tumors. A limited number of clonal hierarchy patterns were identified. A biopsy study of synthetic tumors shows complex spatial intermixing of clones and delineates the importance of biopsy extent. Deep whole-exome multiregional next-generation DNA sequencing of the primary tumors from five patients was performed to validate the results, supporting our main findings from mathematical modeling. In conclusion, our model provides qualitatively realistic predictions of PCA genomic evolution, closely aligned with the evidence available from patient samples. We share the code of the model for further studies. A record of this paper's transparent peer review process is included in the supplemental information.

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.

Unveiling the dynamics and molecular landscape of a rare chronic lymphocytic leukemia subpopulation driving refractoriness: insights from single-cell RNA sequencing

Early identification of resistant cancer cells is currently a major challenge, as their expansion leads to refractoriness. To capture the dynamics of these cells, we made a comprehensive analysis of disease progression and treatment response in a chronic lymphocytic leukemia (CLL) patient using a combination of single-cell and bulk genomic methods. At diagnosis, the patient presented with unfavorable genetic markers, including notch receptor 1 (NOTCH1) mutation and loss(11q). The initial and subsequent treatment lines did not lead to a durable response and the patient developed refractory disease. Refractory CLL cells featured substantial dysregulation in B-cell phenotypic markers such as human leukocyte antigen (HLA) genes, immunoglobulin (IG) genes, CD19 molecule (CD19), membrane spanning 4-domains A1 (MS4A1; previously known as CD20), CD79a molecule (CD79A) and paired box 5 (PAX5), indicating B-cell de-differentiation and disease transformation. We described the clonal evolution and characterized in detail two cell populations that emerged during the refractory disease phase, differing in the presence of high genomic complexity. In addition, we successfully tracked the cells with high genomic complexity back to the time before treatment, where they formed a rare subpopulation. We have confirmed that single-cell RNA sequencing enables the characterization of refractory cells and the monitoring of their development over time.

Endocrine persistence in ER+ breast cancer is accompanied by metabolic vulnerability in oxidative phosphorylation

Despite adjuvant treatment with endocrine therapies, estrogen receptor-positive (ER+) breast cancers recur in a significant proportion of patients. Recurrences are attributable to clinically undetectable endocrine-tolerant persister cancer cells that retain tumor-forming potential. Therefore, strategies targeting such persister cells may prevent recurrent disease. Using CRISPR-Cas9 genome-wide knockout screening in ER+ breast cancer cells, we identified a survival mechanism involving metabolic reprogramming with reliance upon mitochondrial respiration in endocrine-tolerant persister cells. Quantitative proteomic profiling showed reduced levels of glycolytic proteins in persisters. Metabolic tracing of glucose revealed an energy-depleted state in persisters where oxidative phosphorylation was required to generate ATP. A phase II clinical trial was conducted to evaluate changes in mitochondrial markers in primary ER+/HER2-breast tumors induced by neoadjuvant endocrine therapy ( NCT04568616 ). In an analysis of tumor specimens from 32 patients, tumors exhibiting residual cell proliferation after aromatase inhibitor-induced estrogen deprivation with letrozole showed increased mitochondrial content. Genetic profiling and barcode lineage tracing showed that endocrine-tolerant persistence occurred stochastically without genetic predisposition. Mice bearing cell line- and patient-derived xenografts were used to measure the anti-tumor effects of mitochondrial complex I inhibition in the context of endocrine therapy. Pharmacological inhibition of complex I suppressed the tumor-forming potential of persisters and synergized with the anti-estrogen fulvestrant to induce regression of patient-derived xenografts. These findings indicate that mitochondrial metabolism is essential in endocrine-tolerant persister ER+ breast cancer cells and warrant the development of treatment strategies to leverage this vulnerability in the context of endocrine-sensitive disease.

Statement of Significance

Endocrine-tolerant persister cancer cells that survive endocrine therapy can cause recurrent disease. Persister cells exhibit increased energetic dependence upon mitochondria for survival and tumor re-growth potential.

Exploring evolutionary trajectories in ovarian cancer patients by longitudinal analysis of ctDNA

OBJECTIVES

We analysed whether temporal heterogeneity of ctDNA encodes evolutionary patterns in ovarian cancer.

METHODS

Targeted sequencing of 275 cancer-associated genes was performed in a primary tumor biopsy and in ctDNA of six longitudinal plasma samples from 15 patients, using the Illumina platform.

RESULTS

While there was low overall concordance between the mutational spectrum of the primary tumor biopsies vs. ctDNA, TP53 variants were the most commonly shared somatic alterations. Up to three variant clusters were detected in each tumor biopsy, likely representing predominant clones of the primary tumor, most of them harbouring a TP53 variant. By tracing these clusters in ctDNA, we propose that liquid biopsy may allow to assess the contribution of ancestral clones of the tumor to relapsed abdominal masses, revealing two evolutionary patterns. In pattern#1, clusters detected in the primary tumor biopsy were likely relapse seeding clones, as they contributed a major share to ctDNA at relapse. In pattern#2, similar clusters were present in tumors and ctDNA; however, they were entirely cleared from liquid biopsy after chemotherapy and were undetectable at relapse. ctDNA private variants were present among both patterns, with some of them mirroring subclonal expansions after chemotherapy.

CONCLUSIONS

We demonstrate that tracing the temporal heterogeneity of ctDNA, even below exome scale resolution, deciphers evolutionary trajectories in ovarian cancer. Furthermore, we describe two evolutionary patterns that may help to identify relapse seeding clones for targeted therapy.

SJPedPanel: A Pan-Cancer Gene Panel for Childhood Malignancies to Enhance Cancer Monitoring and Early Detection

PURPOSE

The purpose of the study was to design a pan-cancer gene panel for childhood malignancies and validate it using clinically characterized patient samples.

EXPERIMENTAL DESIGN

In addition to 5,275 coding exons, SJPedPanel also covers 297 introns for fusions/structural variations and 7,590 polymorphic sites for copy-number alterations. Capture uniformity and limit of detection are determined by targeted sequencing of cell lines using dilution experiment. We validate its coverage by in silico analysis of an established real-time clinical genomics (RTCG) cohort of 253 patients. We further validate its performance by targeted resequencing of 113 patient samples from the RTCG cohort. We demonstrate its power in analyzing low tumor burden specimens using morphologic remission and monitoring samples.

RESULTS

Among the 485 pathogenic variants reported in RTCG cohort, SJPedPanel covered 86% of variants, including 82% of 90 rearrangements responsible for fusion oncoproteins. In our targeted resequencing cohort, 91% of 389 pathogenic variants are detected. The gene panel enabled us to detect ∼95% of variants at allele fraction (AF) 0.5%, whereas the detection rate is ∼80% at AF 0.2%. The panel detected low-frequency driver alterations from morphologic leukemia remission samples and relapse-enriched alterations from monitoring samples, demonstrating its power for cancer monitoring and early detection.

CONCLUSIONS

SJPedPanel enables the cost-effective detection of clinically relevant genetic alterations including rearrangements responsible for subtype-defining fusions by targeted sequencing of ∼0.15% of human genome for childhood malignancies. It will enhance the analysis of specimens with low tumor burdens for cancer monitoring and early detection.

Comparative analysis of the mutational landscape and evolutionary patterns of pancreatic ductal adenocarcinoma metastases in the liver or peritoneum

Background

Pancreatic ductal adenocarcinoma (PDAC) often presents with liver or peritoneal metastases at diagnosis. Despite similar treatment approaches, patient outcomes vary between these metastatic sites. To improve targeted therapies for metastatic PDAC, a comprehensive analysis of the genetic profiles and evolutionary patterns at these distinct metastatic locations is essential.

Methods

We performed whole exome sequencing on 44 tissue samples from 27 PDAC patients, including primary tumours and matched liver or peritoneal metastases. We analysed somatic mutation profiles, signatures, and affected pathways for each group, and examined clonal evolution using subclonal architectures and phylogenetic trees.

Results

KRAS mutations remained the predominant driver alteration, with a prevalence of 89 % across all tumours. Notably, we observed site-specific differences in mutation frequencies, with KRAS alterations detected in 77.8 % (7/9) of peritoneal metastases and 87.5 % (7/8) of liver metastases. TP53 mutations exhibited a similar pattern, occurring in 55.6 % (5/9) of peritoneal and 37.5 % (3/8) of liver metastases. Intriguingly, we identified site-specific alterations in DNA repair pathway genes, including ATM and BRCA1, with distinct mutational profiles in liver versus peritoneal metastases. Furthermore, liver metastases demonstrated a significantly higher tumor mutational burden (TMB) compared to peritoneal metastases (median [IQR]: 2.14 [1.77-2.71] vs. 1.29 [1.21-1.69] mutations/Mb; P = 0.048).

Conclusions

In conclusion, metastasis of pancreatic cancer may be influenced by variables other than KRAS mutations, such as TP53. PDAC peritoneal and liver metastases may differ in potential therapeutic biomarkers. Further inquiry is needed on the biological mechanisms underlying metastasis and the treatment of diverse metastases.

Somatic mutation phasing and haplotype extension using linked-reads in multiple myeloma

Somatic mutation phasing informs our understanding of cancer-related events, like driver mutations. We generated linked-read whole genome sequencing data for 23 samples across disease stages from 14 multiple myeloma (MM) patients and systematically assigned somatic mutations to haplotypes using linked-reads. Here, we report the reconstructed cancer haplotypes and phase blocks from several MM samples and show how phase block length can be extended by integrating samples from the same individual. We also uncover phasing information in genes frequently mutated in MM, including DIS3, HIST1H1E, KRAS, NRAS, and TP53, phasing 79.4% of 20,705 high-confidence somatic mutations. In some cases, this enabled us to interpret clonal evolution models at higher resolution using pairs of phased somatic mutations. For example, our analysis of one patient suggested that two NRAS hotspot mutations occurred on the same haplotype but were independent events in different subclones. Given sufficient tumor purity and data quality, our framework illustrates how haplotype-aware analysis of somatic mutations in cancer can be beneficial for some cancer cases.

Sequential Analysis of cfDNA Reveals Clonal Evolution in Patients with Neuroblastoma Receiving ALK-Targeted Therapy

PURPOSE

The study of cell-free DNA (cfDNA) enables sequential analysis of tumor cell-specific genetic alterations in patients with neuroblastoma.

EXPERIMENTAL DESIGN

Eighteen patients with relapsing neuroblastoma having received lorlatinib, a third-generation ALK inhibitor, were identified (SACHA national registry and/or in the institution). cfDNA was analyzed at relapse for nine patients and sequentially for five patients (blood/bone marrow plasma) by performing whole-genome sequencing library construction followed by ALK-targeted ddPCR of the hotspot mutations [F1174L, R1275Q, and I1170N; variant allele fraction (VAF) detection limit 0.1%] and whole-exome sequencing (WES) to evaluate disease burden and clonal evolution, following comparison with tumor/germline WES.

RESULTS

Overall response rate to lorlatinib was 33% (CI, 13%-59%), with response observed in 6/10 cases without versus 0/8 cases with MYCN amplification (MNA). ALK VAFs correlated with the overall clinical disease status, with a VAF < 0.1% in clinical remission, versus higher VAFs (>30%) at progression. Importantly, sequential ALK ddPCR detected relapse earlier than clinical imaging. cfDNA WES revealed new SNVs, not seen in the primary tumor, in all instances of disease progression after lorlatinib treatment, indicating clonal evolution, including alterations in genes linked to tumor aggressivity (TP53) or novel targets (EGFR). Gene pathway analysis revealed an enrichment for genes targeting cell differentiation in emerging clones, and cell adhesion in persistent clones. Evidence of clonal hematopoiesis could be observed in follow-up samples.

CONCLUSIONS

We demonstrate the clinical utility of combining ALK cfDNA ddPCR for disease monitoring and cfDNA WES for the study of clonal evolution and resistance mechanisms in patients with neuroblastoma receiving ALK-targeted therapy.

Spatial intra-tumour heterogeneity and treatment-induced genomic evolution in oesophageal adenocarcinoma: implications for prognosis and therapy

BACKGROUND

Oesophageal adenocarcinoma (OAC) is a highly heterogeneous cancer with poor survival. Standard curative treatment is chemotherapy with or without radiotherapy followed by oesophagectomy. Genomic heterogeneity is a feature of OAC and has been linked to treatment resistance.

METHODS

Whole-genome sequencing data from 59 treatment-naïve and 18 post-treatment samples from 29 OAC patients was analysed. Twenty-seven of these were enrolled in the DOCTOR trial, sponsored by the Australasian Gastro-Intestinal Trials Group. Two biopsies from each treatment-naïve tumour were assessed to define 'shared' (between both samples) and 'private' (present in one sample) mutations.

RESULTS

Mutational signatures SBS2/13 (APOBEC) and SBS3 (BRCA) were almost exclusively detected in private mutation populations of treatment-naïve tumours. Patients presenting these signatures had significantly worse disease specific survival. Furthermore, mutational signatures associated with platinum-based chemotherapy treatment as well as high platinum enrichment scores were only detected in post-treatment samples. Additionally, clones with high putative neoantigen binding scores were detected in some treatment-naïve samples suggesting immunoediting of clones.

CONCLUSIONS

This study demonstrates the high intra-tumour heterogeneity in OAC, as well as indicators for treatment-induced changes during tumour evolution. Intra-tumour heterogeneity remains a problem for successful treatment strategies in OAC.

Genomic determinants of response and resistance to inotuzumab ozogamicin in B-cell ALL

ABSTRACT

Inotuzumab ozogamicin (InO) is an antibody-drug conjugate that delivers calicheamicin to CD22-expressing cells. In a retrospective cohort of InO-treated patients with B-cell acute lymphoblastic leukemia, we sought to understand the genomic determinants of the response and resistance to InO. Pre- and post-InO-treated patient samples were analyzed by whole genome, exome, and/or transcriptome sequencing. Acquired CD22 mutations were observed in 11% (3/27) of post-InO-relapsed tumor samples, but not in refractory samples (0/16). There were multiple CD22 mutations per sample and the mechanisms of CD22 escape included epitope loss (protein truncation and destabilization) and epitope alteration. Two CD22 mutant cases were post-InO hyper-mutators resulting from error-prone DNA damage repair (nonhomologous/alternative end-joining repair, or mismatch repair deficiency), suggesting that hypermutation drove escape from CD22-directed therapy. CD22-mutant relapses occurred after InO and subsequent hematopoietic stem cell transplantation (HSCT), suggesting that InO eliminated the predominant clones, leaving subclones with acquired CD22 mutations that conferred resistance to InO and subsequently expanded. Acquired loss-of-function mutations in TP53, ATM, and CDKN2A were observed, consistent with a compromise of the G1/S DNA damage checkpoint as a mechanism for evading InO-induced apoptosis. Genome-wide CRISPR/Cas9 screening of cell lines identified DNTT (terminal deoxynucleotidyl transferase) loss as a marker of InO resistance. In conclusion, genetic alterations modulating CD22 expression and DNA damage response influence InO efficacy. Our findings highlight the importance of defining the basis of CD22 escape and eradication of residual disease before HSCT. The identified mechanisms of escape from CD22-targeted therapy extend beyond antigen loss and provide opportunities to improve therapeutic approaches and overcome resistance. These trials were registered at www.ClinicalTrials.gov as NCT01134575, NCT01371630, and NCT03441061.

LinG3D: visualizing the spatio-temporal dynamics of clonal evolution

BACKGROUND

Cancers are spatially heterogenous, thus their clonal evolution, especially following anti-cancer treatments, depends on where the mutated cells are located within the tumor tissue. For example, cells exposed to different concentrations of drugs, such as cells located near the vessels in contrast to those residing far from the vasculature, can undergo a different evolutionary path. However, classical representations of cell lineage trees do not account for this spatial component of emerging cancer clones. Here, we propose routines to trace spatial and temporal clonal evolution in computer simulations of the tumor evolution models.

RESULTS

The LinG3D (Lineage Graphs in 3D) is an open-source collection of routines (in MATLAB, Python, and R) that enables spatio-temporal visualization of clonal evolution in a two-dimensional tumor slice from computer simulations of the tumor evolution models. These routines draw traces of tumor clones in both time and space, and may include a projection of a selected microenvironmental factor, such as the drug or oxygen distribution within the tumor, if such a microenvironmental factor is used in the tumor evolution model. The utility of LinG3D has been demonstrated through examples of simulated tumors with different number of clones and, additionally, in experimental colony growth assay.

CONCLUSIONS

This routine package extends the classical lineage trees, that show cellular clone relationships in time, by adding the space component to show the locations of cellular clones within the 2D tumor tissue patch from computer simulations of tumor evolution models.

CASi: A framework for cross-timepoint analysis of single-cell RNA sequencing data

Single-cell RNA sequencing (scRNA-seq) technology has been widely used to study the differences in gene expression at the single cell level, providing insights into the research of cell development, differentiation, and functional heterogeneity. Various pipelines and workflows of scRNA-seq analysis have been developed but few considered multi-timepoint data specifically. In this study, we develop CASi, a comprehensive framework for analyzing multiple timepoints' scRNA-seq data, which provides users with: (1) cross-timepoint cell annotation, (2) detection of potentially novel cell types emerged over time, (3) visualization of cell population evolution, and (4) identification of temporal differentially expressed genes (tDEGs). Through comprehensive simulation studies and applications to a real multi-timepoint single cell dataset, we demonstrate the robust and favorable performance of the proposal versus existing methods serving similar purposes.

Human ASXL1-Mutant Hematopoiesis Is Driven by a Truncated Protein Associated with Aberrant Deubiquitination of H2AK119

Mutations in additional sex combs like 1 (ASXL1) confer poor prognosis both in myeloid malignancies and in premalignant clonal hematopoiesis (CH). However, the mechanisms by which these mutations contribute to disease initiation remain unresolved, and mutation-specific targeting has remained elusive. To address this, we developed a human disease model that recapitulates the disease trajectory from ASXL1-mutant CH to lethal myeloid malignancy. We demonstrate that mutations in ASXL1 lead to the expression of a functional, truncated protein and determine that truncated ASXL1 leads to global redistribution of the repressive chromatin mark H2AK119Ub, increased transposase-accessible chromatin, and activation of both myeloid and stem cell gene-expression programs. Finally, we demonstrate that H2AK119Ub levels are tied to truncated ASXL1 expression levels and leverage this observation to demonstrate that inhibition of the PRC1 complex might be an ASXL1-mutant-specific therapeutic vulnerability in both premalignant CH and myeloid malignancy.

SIGNIFICANCE

Mutant ASXL1 is a common driver of CH and myeloid malignancy. Using primary human HSPCs, we determine that truncated ASXL1 leads to redistribution of H2AK119Ub and may affect therapeutic vulnerability to PRC1 inhibition.

Treatment with the apoptosis inhibitor Asunercept reduces clone sizes in patients with lower risk Myelodysplastic Neoplasms

In low-risk Myelodysplastic Neoplasms (MDS), increased activity of apoptosis-promoting factors such as tumor necrosis factor (TNFα) and pro-apoptotic Fas ligand (CD95L) have been described as possible pathomechanisms leading to impaired erythropoiesis. Asunercept (APG101) is a novel therapeutic fusion protein blocking CD95, which has previously shown partial efficacy in reducing transfusion requirement in a clinical phase I trial for low-risk MDS patients (NCT01736436; 2012-11-26). In the current study we aimed to evaluate the effect of Asunercept therapy on the clonal bone marrow composition to identify potential biomarkers to predict response. Bone marrow samples of n = 12 low-risk MDS patients from the above referenced clinical trial were analyzed by serial deep whole exome sequencing in a total of n = 58 time points. We could distinguish a mean of 3.5 molecularly defined subclones per patient (range 2-6). We observed a molecular response defined as reductions of dominant clone sizes by a variant allele frequency (VAF) decrease of at least 10% (mean 20%, range: 10.5-39.2%) in dependency of Asunercept treatment in 9 of 12 (75%) patients. Most of this decline in clonal populations was observed after completion of 12 weeks treatment. Particularly early and pronounced reductions of clone sizes were found in subclones driven by mutations in genes involved in regulation of methylation (n = 1 DNMT3A, n = 1 IDH2, n = 1 TET2). Our results suggest that APG101 could be efficacious in reducing clone sizes of mutated hematopoietic cells in the bone marrow of Myelodysplastic Neoplasms, which warrants further investigation.

Visualizing the Spatio-Temporal Dynamics of Clonal Evolution with LinG3D software

Cancer clonal evolution, especially following anti-cancer treatments, depends on the locations of the mutated cells within the tumor tissue. Cells near the vessels, exposed to higher concentrations of drugs, will undergo a different evolutionary path than cells residing far from the vasculature in the areas of lower drug levels. However, classical representations of cell lineage trees do not account for this spatial component of emerging cancer clones. Here, we propose the LinG3D (Lineage Graphs in 3D) algorithms to trace clonal evolution in space and time. These are an open-source collection of routines (in MATLAB, Python, and R) that enables spatio-temporal visualization of clonal evolution in a two-dimensional tumor slice from computer simulations of the tumor evolution models. These routines draw traces of tumor clones in both time and space, with an option to include a projection of a selected microenvironmental factor, such as the drug or oxygen distribution within the tumor. The utility of LinG3D has been demonstrated through examples of simulated tumors with different number of clones and, additionally, in experimental colony growth assay. This routine package extends the classical lineage trees, that show cellular clone relationships in time, by adding the space component to show the locations of cellular clones within the 2D tumor tissue patch from computer simulations of tumor evolution models.

SJPedPanel: A pan-cancer gene panel for childhood malignancies

Background

Large scale genomics projects have identified driver alterations for most childhood cancers that provide reliable biomarkers for clinical diagnosis and disease monitoring using targeted sequencing. However, there is lack of a comprehensive panel that matches the list of known driver genes. Here we fill this gap by developing SJPedPanel for childhood cancers.

Results

SJPedPanel covers 5,275 coding exons of 357 driver genes, 297 introns frequently involved in rearrangements that generate fusion oncoproteins, commonly amplified/deleted regions (e.g., MYCN for neuroblastoma, CDKN2A and PAX5 for B-/T-ALL, and SMARCB1 for AT/RT), and 7,590 polymorphism sites for interrogating tumors with aneuploidy, such as hyperdiploid and hypodiploid B-ALL or 17q gain neuroblastoma. We used driver alterations reported from an established real-time clinical genomics cohort (n=253) to validate this gene panel. Among the 485 pathogenic variants reported, our panel covered 417 variants (86%). For 90 rearrangements responsible for oncogenic fusions, our panel covered 74 events (82%). We re-sequenced 113 previously characterized clinical specimens at an average depth of 2,500X using SJPedPanel and recovered 354 (91%) of the 389 reported pathogenic variants. We then investigated the power of this panel in detecting mutations from specimens with low tumor purity (as low as 0.1%) using cell line-based dilution experiments and discovered that this gene panel enabled us to detect ∼80% variants with allele fraction of 0.2%, while the detection rate decreases to ∼50% when the allele fraction is 0.1%. We finally demonstrate its utility in disease monitoring on clinical specimens collected from AML patients in morphologic remission.

Conclusions

SJPedPanel enables the detection of clinically relevant genetic alterations including rearrangements responsible for subtype-defining fusions for childhood cancers by targeted sequencing of ∼0.15% of human genome. It will enhance the analysis of specimens with low tumor burdens for cancer monitoring and early detection.

Evolutionary mode and timing of dissemination of high-grade serous carcinomas

Dissemination within the peritoneal cavity is a main determinant of poor patient outcomes from high-grade serous carcinomas (HGSCs). The dissemination process is poorly understood from a cancer evolutionary perspective. We reconstructed the evolutionary trajectories across a median of 5 tumor sites and regions from each of 23 patients based on deep whole-exome sequencing. Polyclonal cancer origin was detected in 1 patient. Ovarian tumors had more complex subclonal architectures than other intraperitoneal tumors in each patient, which indicated that tumors developed earlier in the ovaries. Three common modes of dissemination were identified, including monoclonal or polyclonal dissemination of monophyletic (linear) or polyphyletic (branched) subclones. Mutation profiles of initial or disseminated clones varied greatly among cancers, but recurrent mutations were found in 7 cancer-critical genes, including TP53, BRCA1, BRCA2, and DNMT3A, and in the PI3K/AKT1 pathway. Disseminated clones developed late in the evolutionary trajectory models of most cancers, in particular in cancers with DNA damage repair deficiency. Polyclonal dissemination was predicted to occur predominantly as a single and rapid wave, but chemotherapy exposure was associated with higher genomic diversity of disseminated clones. In conclusion, we described three common evolutionary dissemination modes across HGSCs and proposed factors associated with dissemination diversity.

Druggable growth dependencies and tumor evolution analysis in patient-derived organoids of neuroendocrine neoplasms from multiple body sites

Neuroendocrine neoplasms (NENs) comprise well-differentiated neuroendocrine tumors (NETs) and poorly differentiated neuroendocrine carcinomas (NECs). Treatment options for patients with NENs are limited, in part due to lack of accurate models. We establish patient-derived tumor organoids (PDTOs) from pulmonary NETs and derive PDTOs from an understudied subtype of NEC, large cell neuroendocrine carcinoma (LCNEC), arising from multiple body sites. PDTOs maintain the gene expression patterns, intra-tumoral heterogeneity, and evolutionary processes of parental tumors. Through hypothesis-driven drug sensitivity analyses, we identify ASCL1 as a potential biomarker for response of LCNEC to treatment with BCL-2 inhibitors. Additionally, we discover a dependency on EGF in pulmonary NET PDTOs. Consistent with these findings, we find that, in an independent cohort, approximately 50% of pulmonary NETs express EGFR. This study identifies an actionable vulnerability for a subset of pulmonary NETs, emphasizing the utility of these PDTO models.

Genomic determinants of response and resistance to inotuzumab ozogamicin in B-cell ALL

Inotuzumab ozogamicin (InO) is an antibody-drug conjugate that delivers calicheamicin to CD22-expressing cells. In a retrospective cohort of InO treated patients with B-cell acute lymphoblastic leukemia, we sought to understand the genomic determinants of response to InO. Acquired CD22 mutations were observed in 11% (3/27) of post-InO relapsed tumor samples. There were multiple CD22 mutations per sample and the mechanisms of CD22 escape included protein truncation, protein destabilization, and epitope alteration. Hypermutation by error-prone DNA damage repair (alternative end-joining, mismatch repair deficiency) drove CD22 escape. Acquired loss-of-function mutations in TP53 , ATM and CDKN2A were observed, suggesting compromise of the G1/S DNA damage checkpoint as a mechanism of evading InO-induced apoptosis. In conclusion, genetic alterations modulating CD22 expression and DNA damage response influence InO efficacy. The escape strategies within and beyond antigen loss to CD22-targeted therapy elucidated in this study provide insights into improving therapeutic approaches and overcoming resistance.

KEY POINTS

We identified multiple mechanisms of CD22 antigen escape from inotuzumab ozogamicin, including protein truncation, protein destabilization, and epitope alteration.Hypermutation caused by error-prone DNA damage repair was a driver of CD22 mutation and escape.

VISUAL ABSTRACT

Acquired resistance to a GPRC5D-directed T-cell engager in multiple myeloma is mediated by genetic or epigenetic target inactivation

Bispecific antibodies targeting GPRC5D demonstrated promising efficacy in multiple myeloma, but acquired resistance usually occurs within a few months. Using a single-nucleus multi-omic strategy in three patients from the MYRACLE cohort (ClinicalTrials.gov registration: NCT03807128 ), we identified two resistance mechanisms, by bi-allelic genetic inactivation of GPRC5D or by long-range epigenetic silencing of its promoter and enhancer regions. Molecular profiling of target genes may help to guide the choice of immunotherapy and early detection of resistance in multiple myeloma.

Case report: Rare case of donor cell-derived T-cell acute lymphoblastic leukaemia in a female patient after receiving an allo-transplant from her male sibling

Donor-derived haematological neoplasms, in which recipients present with haematological malignancies that have evolved from transplant donor stem cells, have previously been described for myelodysplastic syndrome, myeloproliferative neoplasms, acute myeloid leukaemia and less often, leukaemias of lymphoid origin. Here we describe a rare and complex case of donor-derived T-cell acute lymphoblastic leukaemia with a relatively short disease latency of less than 4 years. Through genomic and in vitro analyses, we identified novel mutations in NOTCH1 as well as a novel activating mutation in STAT5B; the latter targetable with the clinically available drugs, venetoclax and ruxolitinib.

Human CD34+-derived complete plasmacytoid and conventional dendritic cell vaccine effectively induces antigen-specific CD8+ T cell and NK cell responses in vitro and in vivo

Allogeneic stem cell transplantation (alloSCT) can be curative for hemato-oncology patients due to effective graft-versus-tumor immunity. However, relapse remains the major cause of treatment failure, emphasizing the need for adjuvant immunotherapies. In this regard, post-transplantation dendritic cell (DC) vaccination is a highly interesting strategy to boost graft-versus-tumor responses. Previously, we developed a clinically applicable protocol for simultaneous large-scale generation of end-stage blood DC subsets from donor-derived CD34+ stem cells, including conventional type 1 and 2 DCs (cDC1s and cDC2s), and plasmacytoid DCs (pDCs). In addition, the total cultured end-product (DC-complete vaccine), also contains non-end-stage-DCs (i.e. non-DCs). In this study, we aimed to dissect the phenotypic identity of these non-DCs and their potential immune modulatory functions on the potency of cDCs and pDCs in stimulating tumor-reactive CD8+ T and NK cell responses, in order to obtain rationale for clinical translation of our DC-complete vaccine. The non-DC compartment was heterogeneous and comprised of myeloid progenitors and (immature) granulocyte- and monocyte-like cells. Importantly, non-DCs potentiated toll-like receptor-induced DC maturation, as reflected by increased expression of co-stimulatory molecules and enhanced cDC-derived IL-12 and pDC-derived IFN-α production. Additionally, antigen-specific CD8+ T cells effectively expanded upon DC-complete vaccination in vitro and in vivo. This effect was strongly augmented by non-DCs in an antigen-independent manner. Moreover, non-DCs did not impair in vitro DC-mediated NK cell activation, degranulation nor cytotoxicity. Notably, in vivo i.p. DC-complete vaccination activated i.v. injected NK cells. Together, these data demonstrate that the non-DC compartment potentiates DC-mediated activation and expansion of antigen-specific CD8+ T cells and do not impair NK cell responses in vitro and in vivo. This underscores the rationale for further clinical translation of our CD34+-derived DC-complete vaccine in hemato-oncology patients post alloSCT.

Single-Cell Evolutionary Analysis Reveals Drivers of Plasticity and Mediators of Chemoresistance in Small Cell Lung Cancer

Small cell lung cancer (SCLC) is often a heterogeneous tumor, where dynamic regulation of key transcription factors can drive multiple populations of phenotypically different cells which contribute differentially to tumor dynamics. This tumor is characterized by a very low 2-year survival rate, high rates of metastasis, and rapid acquisition of chemoresistance. The heterogeneous nature of this tumor makes it difficult to study and to treat, as it is not clear how or when this heterogeneity arises. Here we describe temporal, single-cell analysis of SCLC to investigate tumor initiation and chemoresistance in both SCLC xenografts and an autochthonous SCLC model. We identify an early population of tumor cells with high expression of AP-1 network genes that are critical for tumor growth. Furthermore, we have identified and validated the cancer testis antigens (CTA) PAGE5 and GAGE2A as mediators of chemoresistance in human SCLC. CTAs have been successfully targeted in other tumor types and may be a promising avenue for targeted therapy in SCLC.

IMPLICATIONS

Understanding the evolutionary dynamics of SCLC can shed light on key mechanisms such as cellular plasticity, heterogeneity, and chemoresistance.

Intensive management facilitates bacterial invasion on soil microbial community

Intensive management has greatly altered natural forests, especially forests around the world are increasingly being converted into economic plantations. Soil microbiota are critical for community functions in all ecosystems, but the effects of microbial disturbance during economic plantation remain unclear. Here, we used Escherichia coli O157:H7, a model pathogenic species for bacterial invasion, to assess the invasion impacts on the soil microbial community under intensive management. The E. coli invasion was tracked for 135 days to explore the instant and legacy impacts on the resident community. Our results showed that bamboo economic plantations altered soil abiotic and biotic properties, especially increasing pH and community diversity. Higher pH in bamboo soils resulted in longer pathogen survivals than in natural hardwood soils, indicating that pathogen suppression during intensive management should arouse our attention. A longer invasion legacy effect on the resident community (P < 0.05) were found in bamboo soils underlines the need to quantify the soil resilience even when the invasion was unsuccessful. Deterministic processes drove community assembly in bamboo plantations, and this selection acted more strongly during by E. coli invasion than in hardwood soils. We also showed more associated co-occurrence patterns in bamboo plantations, suggesting more complex potential interactions within the microbial community. Apart from community structure, community functions are also strongly related to the resident species associated with invaders. These findings provide new perspectives to understand intensive management facilitates the bacterial invasion, and the impacts would leave potential risks on environmental and human health.

Organic pulses and bacterial invasion alleviated by the resilience of soil microbial community

Biogas slurry is a nutrient-rich secondary product of livestock feces digestion which is recycled as a crop plantation fertilizer and provides exogenous microbes to the soil. However, the effects of biogas slurry microbes on the soil resident community remain unknown. In this study, we examined the ecological consequences of long-term biogas slurry pulse on the soil resident community and found that it promoted crop yield and altered soil characteristics. The soil microbial ecosystem was altered as a result of organic amendments due to the exogenous input of microbes and nutrients. Nevertheless, the soil resident communities were highly resilient to long-term organic pulses, as evidenced by community diversity and composition. The two dominant bacterial species in biogas slurry were Sterolibacterium and Clostridium. Notably, the abundance of Clostridium in biogas slurry increased following long-term amendments, while other species such as GP1 and Subdivision3_genera_incertae_sedis decreased; which was consistent with the results of module-eigengene analysis. Long-term organic pulses shifted the balance of microbial community assembly from stochastic to deterministic processes. Overall, our findings indicated that organic pulses accompanied with bacterial invasion could be alleviated by the resilience of soil microbial communities, thereby emphasizing the importance of microbiota assemblage and network architecture.

Evolutionary states and trajectories characterized by distinct pathways stratify patients with ovarian high grade serous carcinoma

Ovarian high-grade serous carcinoma (HGSC) is typically diagnosed at an advanced stage, with multiple genetically heterogeneous clones existing in the tumors long before therapeutic intervention. Herein we integrate clonal composition and topology using whole-genome sequencing data from 510 samples of 148 patients with HGSC in the prospective, longitudinal, multiregion DECIDER study. Our results reveal three evolutionary states, which have distinct features in genomics, pathways, and morphological phenotypes, and significant association with treatment response. Nested pathway analysis suggests two evolutionary trajectories between the states. Experiments with five tumor organoids and three PI3K inhibitors support targeting tumors with enriched PI3K/AKT pathway with alpelisib. Heterogeneity analysis of samples from multiple anatomical sites shows that site-of-origin samples have 70% more unique clones than metastatic tumors or ascites. In conclusion, these analysis and visualization methods enable integrative tumor evolution analysis to identify patient subtypes using data from longitudinal, multiregion cohorts.

UBTF tandem duplications define a distinct subtype of adult de novo acute myeloid leukemia

Tandem duplications (TDs) of the UBTF gene have been recently described as a recurrent alteration in pediatric acute myeloid leukemia (AML). Here, by screening 1946 newly diagnosed adult AML, we found that UBTF-TDs occur in about 3% of patients aged 18-60 years, in a mutually exclusive pattern with other known AML subtype-defining alterations. The characteristics of 59 adults with UBTF-TD AML included young age (median 37 years), low bone marrow (BM) blast infiltration (median 25%), and high rates of WT1 mutations (61%), FLT3-ITDs (51%) and trisomy 8 (29%). BM morphology frequently demonstrates dysmyelopoiesis albeit modulated by the co-occurrence of FLT3-ITD. UBTF-TD patients have lower complete remission (CR) rates (57% after 1 course and 76% after 2 courses of intensive chemotherapy [ICT]) than UBTF-wild-type patients. In patients enrolled in the ALFA-0702 study (n = 614 patients including 21 with UBTF-TD AML), the 3-year disease-free survival (DFS) and overall survival of UBTF-TD patients were 42.9% (95%CI: 23.4-78.5%) and 57.1% (95%CI: 39.5-82.8%) and did not significantly differ from those of ELN 2022 intermediate/adverse risk patients. Finally, the study of paired diagnosis and relapsed/refractory AML samples suggests that WT1-mutated clones are frequently selected under ICT. This study supports the recognition of UBTF-TD AML as a new AML entity in adults.

Tracking early lung cancer metastatic dissemination in TRACERx using ctDNA

Circulating tumour DNA (ctDNA) can be used to detect and profile residual tumour cells persisting after curative intent therapy1. The study of large patient cohorts incorporating longitudinal plasma sampling and extended follow-up is required to determine the role of ctDNA as a phylogenetic biomarker of relapse in early-stage non-small-cell lung cancer (NSCLC). Here we developed ctDNA methods tracking a median of 200 mutations identified in resected NSCLC tissue across 1,069 plasma samples collected from 197 patients enrolled in the TRACERx study2. A lack of preoperative ctDNA detection distinguished biologically indolent lung adenocarcinoma with good clinical outcome. Postoperative plasma analyses were interpreted within the context of standard-of-care radiological surveillance and administration of cytotoxic adjuvant therapy. Landmark analyses of plasma samples collected within 120 days after surgery revealed ctDNA detection in 25% of patients, including 49% of all patients who experienced clinical relapse; 3 to 6 monthly ctDNA surveillance identified impending disease relapse in an additional 20% of landmark-negative patients. We developed a bioinformatic tool (ECLIPSE) for non-invasive tracking of subclonal architecture at low ctDNA levels. ECLIPSE identified patients with polyclonal metastatic dissemination, which was associated with a poor clinical outcome. By measuring subclone cancer cell fractions in preoperative plasma, we found that subclones seeding future metastases were significantly more expanded compared with non-metastatic subclones. Our findings will support (neo)adjuvant trial advances and provide insights into the process of metastatic dissemination using low-ctDNA-level liquid biopsy.

In Multiple Myeloma, High-Risk Secondary Genetic Events Observed at Relapse Are Present From Diagnosis in Tiny, Undetectable Subclonal Populations

PURPOSE

Multiple myeloma (MM) is characterized by copy number abnormalities (CNAs), some of which influence patient outcomes and are sometimes observed only at relapse(s), suggesting their acquisition during tumor evolution. However, the presence of micro-subclones may be missed in bulk analyses. Here, we use single-cell genomics to determine how often these high-risk events are missed at diagnosis and selected at relapse.

MATERIALS AND METHODS

We analyzed 81 patients with plasma cell dyscrasias using single-cell CNA sequencing. Sixty-six patients were selected at diagnosis, nine at first relapse, and six in presymptomatic stages. A total of 956 newly diagnosed patients with MM and patients with first relapse MM have been identified retrospectively with required cytogenetic data to evaluate enrichment of CNA risk events and survival impact.

RESULTS

A total of 52,176 MM cells were analyzed. Seventy-four patients (91%) had 2-16 subclones. Among these patients, 28.7% had a subclone with high-risk features (del(17p), del(1p32), and 1q gain) at diagnosis. In a patient with a subclonal 1q gain at diagnosis, we analyzed the diagnosis, postinduction, and first relapse samples, which showed a rise of the high-risk 1q gain subclone (16%, 70%, and 92%, respectively). In our clinical database, we found that the 1q gain frequency increased from 30.2% at diagnosis to 43.6% at relapse (odds ratio, 1.78; 95% CI, 1.58 to 2.00). We subsequently performed survival analyses, which showed that the progression-free and overall survival curves were superimposable between patients who had the 1q gain from diagnosis and those who seemingly acquired it at relapse. This strongly suggests that many patients had 1q gains at diagnosis in microclones that were missed by bulk analyses.

CONCLUSION

These data suggest that identifying these scarce aggressive cells may necessitate more aggressive treatment as early as diagnosis to prevent them from becoming the dominant clone.

Amp(1q) and tetraploidy are commonly acquired chromosomal abnormalities in relapsed multiple myeloma

Long-term disease control in multiple myeloma (MM) is typically an unmet medical need, and most patients experience multiple relapses. Fluorescence in situ hybridization (FISH) is the standard technique to detect chromosomal abnormalities (CAs), which are important to estimate the prognosis of MM and the allocation of risk adapted therapies. In advanced stages, the importance of CAs needs further investigation. From 148 MM patients, two or more paired samples, at least one of which was collected at relapse, were analyzed by FISH. Using targeted next-generation sequencing, we molecularly investigated samples harboring relapse-associated CAs. Sixty-one percent of the patients showed a change in the cytogenetic profile during the disease course, including 10% who acquired high-risk cytogenetics. Amp(1q) (≥4 copies of 1q21), driven by an additional increase in copy number in patients who already had 3 copies of 1q21, was the most common acquired CA with 16% affected patients. Tetraploidy, found in 10% of the samples collected at the last time-point, was unstable over the course of the disease and was associated with TP53 lesions. Our results indicate that cytogenetic progression is common in relapsed patients. The relatively high frequency of amp(1q) suggests an active role for this CA in disease progression.

Haplotype-aware analysis of somatic copy number variations from single-cell transcriptomes

Genome instability and aberrant alterations of transcriptional programs both play important roles in cancer. Single-cell RNA sequencing (scRNA-seq) has the potential to investigate both genetic and nongenetic sources of tumor heterogeneity in a single assay. Here we present a computational method, Numbat, that integrates haplotype information obtained from population-based phasing with allele and expression signals to enhance detection of copy number variations from scRNA-seq. Numbat exploits the evolutionary relationships between subclones to iteratively infer single-cell copy number profiles and tumor clonal phylogeny. Analysis of 22 tumor samples, including multiple myeloma, gastric, breast and thyroid cancers, shows that Numbat can reconstruct the tumor copy number profile and precisely identify malignant cells in the tumor microenvironment. We identify genetic subpopulations with transcriptional signatures relevant to tumor progression and therapy resistance. Numbat requires neither sample-matched DNA data nor a priori genotyping, and is applicable to a wide range of experimental settings and cancer types.

Unusual phylogenetic tree and circulating actionable ESR1 mutations in an aggressive luminal/HER2-low breast cancer: Case report

Under therapeutic pressure aggressive tumors evolve rapidly. Herein, a luminal B/HER2-low breast cancer was tracked for >3 years during a total of 6 largely unsuccessful therapy lines, from adjuvant to advanced settings. Targeted next generation sequencing (NGS) of the primary lesion, two metastases and 14 blood drawings suggested a striking, unprecedented coexistence of three evolution modes: punctuated, branched and convergent. Punctuated evolution of the trunk was supported by en bloc inheritance of a large set (19 distinct genes) of copy number alterations. Branched evolution was supported by the distribution of site-specific SNVs. Convergent evolution was characterized by a unique asynchronous expansion of three actionable (OncoKB level 3A) mutations at two consecutive ESR1 codons. Low or undetectable in all the sampled tumor tissues, ESR1 mutations expanded rapidly in blood during HER2/hormone double-blockade, and predicted life-threatening local progression at lung and liver metastatic foci. Dramatic clinical response to Fulvestrant (assigned off-label exclusively based on liquid biopsy) was associated with clearance of all 3 subclones and was in stark contrast to the poor therapeutic efficacy reported in large liquid biopsy-informed interventional trials. Altogether, deconvolution of the tumor phylogenetic tree, as shown herein, may help to customize treatment in breast cancers that rapidly develop refractoriness to multiple drugs.

Molecular Characterization of Acquired Resistance to KRASG12C-EGFR Inhibition in Colorectal Cancer

With the combination of KRASG12C and EGFR inhibitors, KRAS is becoming a druggable target in colorectal cancer. However, secondary resistance limits its efficacy. Using cell lines, patient-derived xenografts, and patient samples, we detected a heterogeneous pattern of putative resistance alterations expected primarily to prevent inhibition of ERK signaling by drugs at progression. Serial analysis of patient blood samples on treatment demonstrates that most of these alterations are detected at a low frequency except for KRASG12C amplification, a recurrent resistance mechanism that rises in step with clinical progression. Upon drug withdrawal, resistant cells with KRASG12C amplification undergo oncogene-induced senescence, and progressing patients experience a rapid fall in levels of this alteration in circulating DNA. In this new state, drug resumption is ineffective as mTOR signaling is elevated. However, our work exposes a potential therapeutic vulnerability, whereby therapies that target the senescence response may overcome acquired resistance.

SIGNIFICANCE

Clinical resistance to KRASG12C-EGFR inhibition primarily prevents suppression of ERK signaling. Most resistance mechanisms are subclonal, whereas KRASG12C amplification rises over time to drive a higher portion of resistance. This recurrent resistance mechanism leads to oncogene-induced senescence upon drug withdrawal and creates a potential vulnerability to senolytic approaches. This article is highlighted in the In This Issue feature, p. 1.

Reconstructing Phylogenetic Relationship in Bladder Cancer: A Methodological Overview

Bladder cancer (BC) expresses itself as a highly heterogeneous disease both at the histological and molecular level, often occurring as synchronous or metachronous multifocal disease with high risk of recurrence and potential to metastasize. Multiple sequencing studies focusing on both non-muscle-invasive bladder cancer (NMIBC) and muscle-invasive bladder cancer (MIBC) gave insights into the extent of both inter- and intrapatient heterogeneity, but many questions on clonal evolution in BC remain unanswered. In this review article, we provide an overview over the technical and theoretical concepts linked to reconstructing evolutionary trajectories in BC and propose a set of tools and established software for phylogenetic analysis.

Insights from a rare myeloproliferative neoplasm with coexisting BCR-ABL1 fusion gene, CALR, and TET2 mutations treated with nilotinib and ruxolitinib

Myeloproliferative neoplasms (MPNs) with concurrent BCR-ABL1 fusion gene and CALR mutation are especially rare. We report a patient with coexisting BCR-ABL1 fusion gene, CALR, and TET2 mutations who was treated with the combination of the second-generation TKI nilotinib and JAK1/JAK2 inhibitor ruxolitinib.

clevRvis: visualization techniques for clonal evolution

BACKGROUND

A thorough analysis of clonal evolution commonly requires integration of diverse sources of data (e.g., karyotyping, next-generation sequencing, and clinical information). Subsequent to actual reconstruction of clonal evolution, detailed analysis and interpretation of the results are essential. Often, however, only few tumor samples per patient are available. Thus, information on clonal development and therapy effect may be incomplete. Furthermore, analysis of biallelic events-considered of high relevance with respect to disease course-can commonly only be realized by time-consuming analysis of the raw results and even raw sequencing data.

RESULTS

We developed clevRvis, an R/Bioconductor package providing an extensive set of visualization techniques for clonal evolution. In addition to common approaches for visualization, clevRvis offers a unique option for allele-aware representation: plaice plots. Biallelic events may be visualized and inspected at a glance. Analyzing 4 public datasets, we show that plaice plots help to gain new insights into tumor development and investigate hypotheses on disease progression and therapy resistance. In addition to a graphical user interface, automatic phylogeny-aware color coding of the plots, and an approach to explore alternative trees, clevRvis provides 2 algorithms for fully automatic time point interpolation and therapy effect estimation. Analyzing 2 public datasets, we show that both approaches allow for valid approximation of a tumor's development in between measured time points.

CONCLUSIONS

clevRvis represents a novel option for user-friendly analysis of clonal evolution, contributing to gaining new insights into tumor development.

Post-irradiation intratumoral heterogeneity modulates response to immune checkpoint inhibition therapy in a murine melanoma model

PURPOSE

The underlying mechanism for radiation as a potentiator of immune checkpoint inhibition (ICI) is unclear. We developed a novel murine model to investigate the effects of post-irradiation intratumoral heterogeneity (ITH) on response to ICI.

EXPERIMENTAL DESIGN

Parental mouse melanoma B16F10 cells were irradiated in vitro (5Gy x 3 fractions), then an a priori determined number of resulting colonies were implanted in C57BL/6J immunocompetent mice creating syngeneic models of unirradiated (parental) and irradiated tumors with low (irradiated-L) and high (irradiated-H) ITH. Mice were treated with placebo, α-PD-L1, α-CTLA-4 or dual ICI. Murine tumors underwent whole exome sequencing (WES). Clinically correlated paired pre- and post-irradiation patient rectal adenocarcinoma samples underwent WES.

RESULTS

Irradiated-L tumors showed increased tumor mutational burden (TMB) and a sustained decrease in ITH. Irradiated-L tumors were predicted to express five neoantigens with high variant allele frequency/clonal distribution. Mice with irradiated-L and irradiated-H versus parental B16F10 tumors demonstrated longer overall survival with dual ICI. Only mice with irradiated-L tumors experienced an overall survival benefit with single agent ICI. Clinically correlated rectal adenocarcinoma samples showed similarly increased TMB and decreased ITH following irradiation.

CONCLUSIONS

Post-irradiation ITH modulates ICI response in a murine melanoma model. Irradiation may offer a mechanism to widen the therapeutic window of ICI.

SLC11A1 as a stratification indicator for immunotherapy or chemotherapy in patients with glioma

Background

Glioma is a fatal tumor originating from the brain, which accounts for most intracranial malignancies. Currently, Immunotherapy has turned into a novel and promising treatment in glioma patients. however, there are still few effective biomarkers to mirror the reaction to immunotherapy in patients with glioma. Therefore, we intended to elucidate the evaluable efficacy of SLC11A1 in glioma patients.

Methods

In this study, samples from Shanghai General Hospital and data from TCGA, GEO, CGGA datasets were used to investigate and validate the relationship between SLC11A1 and the progression of glioma. We evaluated the predictive value of SLC11A1 on the prognosis of glioma with cox regression analysis. Then the relationship between immune infiltration and SLC11A1 was also analyzed. Ultimately, we performed the prediction on the immunotherapeutic response and therapeutic drugs according to the expression of SLC11A1.

Results

Expression of SLC11A1 increased with progression and predicted unfavorable prognosis for glioma patients. The hazard ratio for SLC11A1 expression was 2.33 with 95% CI (1.92-2.58) (P < 0.001) in cox analysis. And based on expression, we found SLC11A1 stratified glioma patients into subgroups with different immune activation statuses. Moreover, we observed that patients with higher SLC11A1 levels companied with better immunotherapeutic response, while those with lower SLC11A1 levels may respond better to temozolomide.

Conclusion

This study provided evidence that SLC11A1 was a novel prognostic marker and immunotherapy response indicator for gliomas. In some cases, SLC11A1 could be an effective marker for identifying patients who might benefit from immunotherapy or chemotherapy.

Pharmacogenomic Profiling of Pediatric Acute Myeloid Leukemia to Identify Therapeutic Vulnerabilities and Inform Functional Precision Medicine

Despite the expanding portfolio of targeted therapies for adults with acute myeloid leukemia (AML), direct implementation in children is challenging due to inherent differences in underlying genetics. Here we established the pharmacologic profile of pediatric AML by screening myeloblast sensitivity to approved and investigational agents, revealing candidates of immediate clinical relevance. Drug responses ex vivo correlated with patient characteristics, exhibited age-specific alterations, and concorded with activities in xenograft models. Integration with genomic data uncovered new gene-drug associations, suggesting actionable therapeutic vulnerabilities. Transcriptome profiling further identified gene-expression signatures associated with on- and off-target drug responses. We also demonstrated the feasibility of drug screening-guided treatment for children with high-risk AML, with two evaluable cases achieving remission. Collectively, this study offers a high-dimensional gene-drug clinical data set that could be leveraged to research the unique biology of pediatric AML and sets the stage for realizing functional precision medicine for the clinical management of the disease.

SIGNIFICANCE

We conducted integrated drug and genomic profiling of patient biopsies to build the functional genomic landscape of pediatric AML. Age-specific differences in drug response and new gene-drug interactions were identified. The feasibility of functional precision medicine-guided management of children with high-risk AML was successfully demonstrated in two evaluable clinical cases. This article is highlighted in the In This Issue feature, p. 476.

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.

Lorlatinib and compound mutations in ALK+ large-cell neuroendocrine lung carcinoma: a case report

Large-cell neuroendocrine lung carcinoma (LCNEC) is a high-grade neoplasm with median survival of 1 year and limited therapeutic options. Here, we report the unusual case of a 47-yr-old female smoker with stage IV LCNEC featuring EML4-ALK variant 2 (E20:A20), wild-type TP53/RB1, and low tumor mutational burden of 3.91 mut/Mb. Despite early progression within 3 mo under crizotinib, a durable response was achieved with alectinib. Oligoprogression in the left breast 10 mo later was treated by surgery, followed by a switch to ceritinib upon multifocal progression and detection of ALK:p.V1180L in the mastectomy specimen, but without success. Another rebiopsy revealed ALK:p.L1196M, but the tumor did not respond to brigatinib or carboplatin/pemetrexed, before stabilization under lorlatinib. Diffuse progression 8 mo later with detection of ALK :p.L1196M/p.G1202R and p.L1196M/ p.D1203N evolving from the previous p.L1196M did not respond to chemoimmunotherapy, and the patient succumbed with an overall survival (OS) of 37 mo. This case illustrates the importance of molecular profiling for LCNEC regardless of smoking status, and the superiority of next-generation ALK inhibitors compared to crizotinib for ALK+ cases. Lorlatinib retained efficacy in the heavily pretreated setting, whereas its upfront use could possibly have prevented the stepwise emergence of compound ALK mutations. Furthermore, the disease course was more aggressive and OS shorter compared to the V2/TP53wt ALK+ lung adenocarcinoma, whereas crizotinib, ceritinib, and brigatinib did not confer the benefit expected according to next-generation sequencing results, which also underline the need for more potent drugs against ALK in the high-risk setting of neuroendocrine histology.