In vivo single-cell CRISPR uncovers distinct TNF programmes in tumour evolution.

The tumour evolution model posits that malignant transformation is preceded by randomly distributed driver mutations in cancer genes, which cause clonal expansions in phenotypically normal tissues. Although clonal expansions can remodel entire tissues1-3, the mechanisms that result in only a small number of clones transforming into malignant tumours remain unknown. Here we develop an in vivo single-cell CRISPR strategy to systematically investigate tissue-wide clonal dynamics of the 150 most frequently mutated squamous cell carcinoma genes. We couple ultrasound-guided in utero lentiviral microinjections, single-cell RNA sequencing and guide capture to longitudinally monitor clonal expansions and document their underlying gene programmes at single-cell transcriptomic resolution. We uncover a tumour necrosis factor (TNF) signalling module, which is dependent on TNF receptor 1 and involving macrophages, that acts as a generalizable driver of clonal expansions in epithelial tissues. Conversely, during tumorigenesis, the TNF signalling module is downregulated. Instead, we identify a subpopulation of invasive cancer cells that switch to an autocrine TNF gene programme associated with epithelial-mesenchymal transition. Finally, we provide in vivo evidence that the autocrine TNF gene programme is sufficient to mediate invasive properties and show that the TNF signature correlates with shorter overall survival of patients with squamous cell carcinoma. Collectively, our study demonstrates the power of applying in vivo single-cell CRISPR screening to mammalian tissues, unveils distinct TNF programmes in tumour evolution and highlights the importance of understanding the relationship between clonal expansions in epithelia and tumorigenesis.

Clonal evolution and relapse in early-stage follicular lymphoma - a tree with many branches†.

Follicular lymphoma (FL) is an indolent B-cell neoplasm characterised by multistep evolution from premalignant precursor cells carrying the hallmark t(14;18) translocation in the majority of cases. In a new article in The Journal of Pathology, samples of relapsed early-stage FL - primary manifestation and relapse with or without transformation - initially treated with radiotherapy only, were studied for clonal relationships and evolution. Using somatic mutations and the rearranged immunoglobulin sequences as markers, the majority of paired lymphoma samples showed so-called branched evolution from a common, possibly premalignant progenitor cell, with both shared and private mutations. In addition, clonally unrelated cases were identified. This and previous studies with similar findings clearly document that relapse or transformation of FL in many instances not necessarily represents a linear progression of disease due to acquisition of additional mutations and therapy resistance, but rather new outgrowths derived from a pool of clonally related, long-lived, and low proliferating precursor cells, or even unrelated second neoplasms. © 2024 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

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.

Single-cell transcriptome analysis reveals subtype-specific clonal evolution and microenvironmental changes in liver metastasis of pancreatic adenocarcinoma and their clinical implications.

BACKGROUND: Intratumoral heterogeneity (ITH) and tumor microenvironment (TME) of pancreatic ductal adenocarcinoma (PDAC) play important roles in tumor evolution and patient outcomes. However, the precise characterization of diverse cell populations and their crosstalk associated with PDAC progression and metastasis is still challenging. METHODS: We performed single-cell RNA sequencing (scRNA-seq) of treatment-naïve primary PDAC samples with and without paired liver metastasis samples to understand the interplay between ITH and TME in the PDAC evolution and its clinical associations. RESULTS: scRNA-seq analysis revealed that even a small proportion (22%) of basal-like malignant ductal cells could lead to poor chemotherapy response and patient survival and that epithelial-mesenchymal transition programs were largely subtype-specific. The clonal homogeneity significantly increased with more prevalent and pronounced copy number gains of oncogenes, such as KRAS and ETV1, and losses of tumor suppressor genes, such as SMAD2 and MAP2K4, along PDAC progression and metastasis. Moreover, diverse immune cell populations, including naïve SELLhi regulatory T cells (Tregs) and activated TIGIThi Tregs, contributed to shaping immunosuppressive TMEs of PDAC through cellular interactions with malignant ductal cells in PDAC evolution. Importantly, the proportion of basal-like ductal cells negatively correlated with that of immunoreactive cell populations, such as cytotoxic T cells, but positively correlated with that of immunosuppressive cell populations, such as Tregs. CONCLUSION: We uncover that the proportion of basal-like subtype is a key determinant for chemotherapy response and patient outcome, and that PDAC clonally evolves with subtype-specific dosage changes of cancer-associated genes by forming immunosuppressive microenvironments in its progression and metastasis.

Clonal expansion in normal tissues.

Cancer originates from a single ancestral cell that acquires a driver mutation, which confers a growth or survival advantage, followed by the acquisition of additional driver mutations by descendant cells. Recently, it has become evident that somatic cell mutations accumulate in normal tissues with aging and exposure to environmental factors, such as alcohol, smoking, and UV rays, increases the mutation rate. Clones harboring driver mutations expand with age, leading to tissue remodeling. Lineage analysis of myeloproliferative neoplasms and der(1;16)-positive breast cancer revealed that driver mutations were acquired early in our lives and that the development of cancer takes decades, unveiling the previously unknown early process of cancer development. Evidence that clonal hematopoiesis affects various diseases, including nonneoplastic diseases, highlights the potential role of the identification and functional analysis of mutated clones in unraveling unknown pathologies. In this review, we summarize the recent updates on clonal expansion in normal tissues and the natural history of cancer revealed through lineage analysis of noncancerous and cancerous tissues.

Co-occurrence of Erdheim-Chester disease and clonally evolving acute myeloid leukemia with FLT3-ITD and PTPN11 mutations.

Erdheim-Chester disease (ECD) is a rare histiocytosis that tends to co-exist with other myeloid malignancies. Here, we use genetic and transcriptomic sequencing to delineate a case of co-occurring BRAFV600E-mutated ECD and acute myeloid leukemia (AML), followed by AML remission and relapse. The AML relapse involved the extinction of clones with KMT2A-AFDN and FLT3-ITD, and the predominance of PTPN11-mutated subclones with distinct transcriptomic features. This case report has highlighted the screening for other myeloid malignancies at the diagnosis of ECD and the clinical significance of PTPN11-mutated AML subclones that require meticulous monitoring.

Relationship between clonal evolution and drug resistance in bladder cancer: A genomic research review.

Bladder cancer stands as a prevalent global malignancy, exhibiting notable sex-based variations in both incidence and prognosis. Despite substantial strides in therapeutic approaches, the formidable challenge of drug resistance persists. The genomic landscape of bladder cancer, characterized by intricate clonal heterogeneity, emerges as a pivotal determinant in fostering this resistance. Clonal evolution, encapsulating the dynamic transformations within subpopulations of tumor cells over time, is implicated in the emergence of drug-resistant traits. Within this review, we illuminate contemporary insights into the role of clonal evolution in bladder cancer, elucidating its influence as a driver in tumor initiation, disease progression, and the formidable obstacle of therapy resistance.

Protocol to study the inheritance and propagation of non-genetically encoded states using barcode decay lineage tracing.

Here, we present a protocol to perform barcode decay lineage tracing followed by single-cell transcriptome analysis (BdLT-Seq). We describe steps for BdLT-Seq experimental design, building barcoded episome reporters, performing episome transfection, and barcode retrieval. We then describe procedures for sequencing library construction while providing options for sample multiplexing and data analysis. This BdLT-Seq technique enables the assessment of clonal evolution in a directional manner while preserving isogeneity, thus allowing the comparison of non-genetic molecular features between isogenic cell lineages. For complete details on the use and execution of this protocol, please refer to Shlyakhtina et al. (2023).1.

Case report: Cytopenias in VEXAS syndrome - a WHO 2022 based approach in a single-center cohort.

VEXAS syndrome is an acquired autoinflammatory disease characterized in most cases by cytopenias and macrocytic anemia. Dyshematopoiesis is a frequent finding in chronic inflammatory conditions and therefore, cytopenias are not easily classified in VEXAS patients. Here we report a series of 7 patients affected by VEXAS associated cytopenias, treated at our center. The use of NGS, together with morphological assays, integrated with the WHO 2022 criteria, allowed to identify three subsets of VEXAS associated cytopenias: ICUS (idiopathic cytopenia of uncertain significance), CCUS (clonal cytopenia of uncertain significance) at high risk of clonal evolution, and MDS. This approach could help to better understand the nature of VEXAS associated cytopenias and to guide the use of specific targeted treatments in order to achieve long lasting responses.

Hematopoietic Clonal Evolution Goes Spatial.

SUMMARY: The spatial distribution of cells carrying clonal hematopoiesis mutations in the bone marrow and the potential role of interactions with the microenvironment are largely unknown. This study takes clonal evolution to the spatial level by describing a novel technique examining the spatial location of mutated clones in the bone marrow and the first evidence that mutated hematopoietic clones are spatially constrained and have heterogenous locations within millimeters of distance. See related article by Young et al., p. 153 (10).

Emerging and extensive clonal evolution in the pancreas.

Pancreatic cancer is one of the most lethal malignancies, yet much remains to be learned regarding how its precursors develop. In a recent Nature publication, Braxton and Kiemen et al. found that the normal, adult pancreas harbors hundreds to thousands of pancreatic cancer precursors evolving by a variety of routes.

Modeling stress-induced responses: plasticity in continuous state space and gradual clonal evolution.

Mathematical models of cancer and bacterial evolution have generally stemmed from a gene-centric framework, assuming clonal evolution via acquisition of resistance-conferring mutations and selection of their corresponding subpopulations. More recently, the role of phenotypic plasticity has been recognized and models accounting for phenotypic switching between discrete cell states (e.g., epithelial and mesenchymal) have been developed. However, seldom do models incorporate both plasticity and mutationally driven resistance, particularly when the state space is continuous and resistance evolves in a continuous fashion. In this paper, we develop a framework to model plastic and mutational mechanisms of acquiring resistance in a continuous gradual fashion. We use this framework to examine ways in which cancer and bacterial populations can respond to stress and consider implications for therapeutic strategies. Although we primarily discuss our framework in the context of cancer and bacteria, it applies broadly to any system capable of evolving via plasticity and genetic evolution.

Mutation order in acute myeloid leukemia identifies uncommon patterns of evolution and illuminates phenotypic heterogeneity.

Acute myeloid leukemia (AML) has a poor prognosis and a heterogeneous mutation landscape. Although common mutations are well-studied, little research has characterized how the sequence of mutations relates to clinical features. Using published, single-cell DNA sequencing data from three institutions, we compared clonal evolution patterns in AML to patient characteristics, disease phenotype, and outcomes. Mutation trees, which represent the order of select mutations, were created for 207 patients from targeted panel sequencing data using 1 639 162 cells, 823 mutations, and 275 samples. In 224 distinct orderings of mutated genes, mutations related to DNA methylation typically preceded those related to cell signaling, but signaling-first cases did occur, and had higher peripheral cell counts, increased signaling mutation homozygosity, and younger patient age. Serial sample analysis suggested that NPM1 and DNA methylation mutations provide an advantage to signaling mutations in AML. Interestingly, WT1 mutation evolution shared features with signaling mutations, such as WT1-early being proliferative and occurring in younger individuals, trends that remained in multivariable regression. Some mutation orderings had a worse prognosis, but this was mediated by unfavorable mutations, not mutation order. These findings add a dimension to the mutation landscape of AML, identifying uncommon patterns of leukemogenesis and shedding light on heterogeneous phenotypes.

Single-cell multiomics reveals the interplay of clonal evolution and cellular plasticity in hepatoblastoma.

Hepatoblastomas (HB) display heterogeneous cellular phenotypes that influence the clinical outcome, but the underlying mechanisms are poorly understood. Here, we use a single-cell multiomic strategy to unravel the molecular determinants of this plasticity. We identify a continuum of HB cell states between hepatocytic (scH), liver progenitor (scLP) and mesenchymal (scM) differentiation poles, with an intermediate scH/LP population bordering scLP and scH areas in spatial transcriptomics. Chromatin accessibility landscapes reveal the gene regulatory networks of each differentiation pole, and the sequence of transcription factor activations underlying cell state transitions. Single-cell mapping of somatic alterations reveals the clonal architecture of each tumor, showing that each genetic subclone displays its own range of cellular plasticity across differentiation states. The most scLP subclones, overexpressing stem cell and DNA repair genes, proliferate faster after neo-adjuvant chemotherapy. These results highlight how the interplay of clonal evolution and epigenetic plasticity shapes the potential of HB subclones to respond to chemotherapy.

Impact of Clonal Heterogeneity in Multiple Myeloma.

Multiple myeloma is characterized by a highly heterogeneous disease distribution within the bone marrow-containing skeletal system. In this review, we introduce the molecular mechanisms underlying clonal heterogeneity and the spatio-temporal evolution of myeloma. We discuss the clinical impact of clonal heterogeneity, which is thought to be one of the biggest obstacles to overcome therapy resistance and to achieve cure.

Genetic intratumour heterogeneity and clonal evolution in extramammary Paget's disease.

All tumour cells in a patient have shared and non-shared genetic alterations, and the diversity of mutations is described as intratumoural heterogeneity (ITH). Multiregion sequencing is a genome sequencing analytical technique used for multiple, spatially-separated biopsy tissues that may further our understanding of ITH and tumour evolution. Although genetic mutations in extramammary Paget's disease (EMPD) have recently been detected by next-generation sequencing analysis, there have been no reports of ITH based on multiregion sequencing in EMPD. Thus, we clarified the landscape of ITH and tumour evolution in EMPD. We performed whole-exome sequencing on 35 tissues (30 tumour tissues and five normal skin samples as a paired control), collected from five patients with EMPD. The rate of private mutations was significantly higher than that of ubiquitous and shared mutations. Ubiquitous mutations were not present in driver genes, and most driver genes exhibited private and shared mutations. The most frequent base substitution was C>T in almost all lesions, and most mutational signatures corresponded to signature 1, 2, 3, and 8. The types of proposed aetiology in most lesions were based on age and AID/APOBEC family and BRCA1/BRCA2 mutations. Evolutionary trees were characterized by short trunks and long branches due to the extremely high ratio of private mutations. In contrast, pathogenic factors, such as base substitutions, mutational signatures, and proposed aetiology, were shared. Tumour evolution in EMPD appears to be characterized by a high level of genetic ITH with shared background factors.

Joint analysis of mutational and transcriptional landscapes in human cancer reveals key perturbations during cancer evolution.

BACKGROUND: Tumors are able to acquire new capabilities, including traits such as drug resistance and metastasis that are associated with unfavorable clinical outcomes. Single-cell technologies have made it possible to study both mutational and transcriptomic profiles, but as most studies have been conducted on model systems, little is known about cancer evolution in human patients. Hence, a better understanding of cancer evolution could have important implications for treatment strategies. RESULTS: Here, we analyze cancer evolution and clonal selection by jointly considering mutational and transcriptomic profiles of single cells acquired from tumor biopsies from 49 lung cancer samples and 51 samples with chronic myeloid leukemia. Comparing the two profiles, we find that each clone is associated with a preferred transcriptional state. For metastasis and drug resistance, we find that the number of mutations affecting related genes increases as the clone evolves, while changes in gene expression profiles are limited. Surprisingly, we find that mutations affecting ligand-receptor interactions with the tumor microenvironment frequently emerge as clones acquire drug resistance. CONCLUSIONS: Our results show that lung cancer and chronic myeloid leukemia maintain a high clonal and transcriptional diversity, and we find little evidence in favor of clonal sweeps. This suggests that for these cancers selection based solely on growth rate is unlikely to be the dominating driving force during cancer evolution.

Chromoplexy Is a Frequent Early Clonal Event in EWSR1-Rearranged Round Cell Sarcomas That Can Be Detected Using Clinically Validated Targeted Sequencing Panels.

Chromoplexy is a phenomenon defined by large-scale chromosomal chained rearrangements. A previous study observed chromoplectic events in a subset of Ewing sarcomas (ES), which was linked to an increased relapse rate. Chromoplexy analysis could potentially facilitate patient risk stratification, particularly if it could be detected with clinically applied targeted next-generation sequencing (NGS) panels. Using DELLY, a structural variant (SV) calling algorithm that is part of the MSK-IMPACT pipeline, we characterized the spectrum of SVs in EWSR1-fused round cell sarcomas, including 173 ES and 104 desmoplastic small round cell tumors (DSRCT), to detect chromoplexy and evaluate its association with clinical and genomic features. Chromoplectic events were detected in 31% of the ES cases and 19% of the DSRCT cases. EWSR1 involvement accounted for 76% to 93% of these events, being rearranged with diverse noncanonical gene partners across the genome, involving mainly translocations but also intrachromosomal deletions and inversions. A major breakpoint cluster was located on EWSR1 exons 8-13. In a subset of cases, the SVs disrupted adjacent loci, forming deletion bridges. Longitudinal sequencing and breakpoint allele fraction analysis showed that chromoplexy is an early event that remains detectable throughout disease progression and likely develops simultaneously with the driver fusion. The presence of chromoplexy was validated in an external ES patient cohort with whole exome sequencing. Chromoplexy was significantly more likely to be present in cases that were metastatic at presentation. Together, this study identifies chromoplexy as a frequent genomic alteration in diverse EWSR1-rearranged tumors that can be captured by targeted NGS panels. SIGNIFICANCE: Chromoplexy is detectable using targeted NGS in a substantial portion of EWSR1-rearranged round cell sarcomas as an early and persistent clonal event, expanding the genomic complexity of fusion-associated sarcomas.