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.

Disentangling oncogenic amplicons in esophageal adenocarcinoma.

Esophageal adenocarcinoma is a prominent example of cancer characterized by frequent amplifications in oncogenes. However, the mechanisms leading to amplicons that involve breakage-fusion-bridge cycles and extrachromosomal DNA are poorly understood. Here, we use 710 esophageal adenocarcinoma cases with matched samples and patient-derived organoids to disentangle complex amplicons and their associated mechanisms. Short-read sequencing identifies ERBB2, MYC, MDM2, and HMGA2 as the most frequent oncogenes amplified in extrachromosomal DNAs. We resolve complex extrachromosomal DNA and breakage-fusion-bridge cycles amplicons by integrating of de-novo assemblies and DNA methylation in nine long-read sequenced cases. Complex amplicons shared between precancerous biopsy and late-stage tumor, an enrichment of putative enhancer elements and mobile element insertions are potential drivers of complex amplicons' origin. We find that patient-derived organoids recapitulate extrachromosomal DNA observed in the primary tumors and single-cell DNA sequencing capture extrachromosomal DNA-driven clonal dynamics across passages. Prospectively, long-read and single-cell DNA sequencing technologies can lead to better prediction of clonal evolution in esophageal adenocarcinoma.

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).

Single-Cell DNA Sequencing Reveals an Evolutionary Pattern of CHIP in Transplant Eligible Multiple Myeloma Patients.

Clonal hematopoiesis of indeterminate potential (CHIP) refers to the phenomenon where a hematopoietic stem cell acquires fitness-increasing mutation(s), resulting in its clonal expansion. CHIP is frequently observed in multiple myeloma (MM) patients, and it is associated with a worse outcome. High-throughput amplicon-based single-cell DNA sequencing was performed on circulating CD34+ cells collected from twelve MM patients before autologous stem cell transplantation (ASCT). Moreover, in four MM patients, longitudinal samples either before or post-ASCT were collected. Single-cell sequencing and data analysis were assessed using the MissionBio Tapestri® platform, with a targeted panel of 20 leukemia-associated genes. We detected CHIP pathogenic mutations in 6/12 patients (50%) at the time of transplant. The most frequently mutated genes were TET2, EZH2, KIT, DNMT3A, and ASXL1. In two patients, we observed co-occurring mutations involving an epigenetic modifier (i.e., DNMT3A) and/or a gene involved in splicing machinery (i.e., SF3B1) and/or a tyrosine kinase receptor (i.e., KIT) in the same clone. Longitudinal analysis of paired samples revealed a positive selection of mutant high-fitness clones over time, regardless of their affinity with a major or minor sub-clone. Copy number analysis of the panel of all genes did not show any numerical alterations present in stem cell compartment. Moreover, we observed a tendency of CHIP-positive patients to achieve a suboptimal response to therapy compared to those without. A sub-clone dynamic of high-fitness mutations over time was confirmed.

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.

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.

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.

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.

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.

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.

Clonal hematopoiesis of indeterminate potential in persons with HIV.

BACKGROUND: Clonal hematopoiesis of indeterminate potential (CHIP) has been associated with older age, inflammation and with risk of coronary artery disease (CAD). We aimed to characterize the burden of CHIP, and to explore the association between CHIP, inflammatory markers, and CAD in older persons with HIV (PWH). METHODS: From the Copenhagen Comorbidity in HIV Infection (COCOMO) study, we included 190 individuals older than 55 years of age. We defined CHIP as variant allele fraction at least 2%. CAD was categorized according to the most severe coronary artery lesion on coronary computed tomography (CT) angiography as no coronary atherosclerosis; any atherosclerosis defined as at least 1% stenosis and obstructive CAD defined as at least 50% stenosis. RESULTS: In the entire population (median age 66 years, 87% men), we identified a total of 62 mutations distributed among 49 (26%) participants. The three most mutated genes were DNMT3A , TET2 , and ASXL1 , accounting for 49, 25, and 16% of mutations, respectively. Age and sex were the only variables associated with CHIP. IL-1ß, IL-1Ra, IL-2, IL-6, IL-10, soluble CD14, soluble CD163 and TNF-α were not associated with CHIP, and CHIP was not associated with any atherosclerosis or with obstructive CAD in adjusted analyses. CONCLUSION: In older, well treated, Scandinavian PWH, more than one in four had at least one CHIP mutation. We did not find evidence of an association between CHIP and inflammatory markers or between CHIP and CAD. CHIP is an unlikely underlying mechanism to explain the association between inflammation and CAD in treated HIV disease.

Mutations, Bottlenecks, and Clonal Sweeps: How Environmental Carcinogens and Genomic Changes Shape Clonal Evolution during Tumor Progression.

The transition from a single, initiated cell to a full-blown malignant tumor involves significant genomic evolution. Exposure to carcinogens-whether directly mutagenic or not-can drive progression toward malignancy, as can stochastic acquisition of cancer-promoting genetic events. Mouse models using both carcinogens and germline genetic manipulations have enabled precise inquiry into the evolutionary dynamics that take place as a tumor progresses from benign to malignant to metastatic stages. Tumor progression is characterized by changes in somatic point mutations and copy-number alterations, even though any single tumor can itself have a high or low burden of genomic alterations. Further, lineage-tracing, single-cell analyses and CRISPR barcoding have revealed the distinct clonal dynamics within benign and malignant tumors. Application of these tools in a range of mouse models can shed unique light on the patterns of clonal evolution that take place in both mouse and human tumors.

Chronic lymphocytic leukemia patient-derived xenografts recapitulate clonal evolution to Richter transformation.

Chronic lymphocytic leukemia (CLL) is a B-cell neoplasm with a heterogeneous clinical behavior. In 5-10% of patients the disease transforms into a diffuse large-B cell lymphoma known as Richter transformation (RT), which is associated with dismal prognosis. Here, we aimed to establish patient-derived xenograft (PDX) models to study the molecular features and evolution of CLL and RT. We generated two PDXs by injecting CLL (PDX12) and RT (PDX19) cells into immunocompromised NSG mice. Both PDXs were morphologically and phenotypically similar to RT. Whole-genome sequencing analysis at different time points of the PDX evolution revealed a genomic landscape similar to RT tumors from both patients and uncovered an unprecedented RT subclonal heterogeneity and clonal evolution during PDX generation. In PDX12, the transformed cells expanded from a very small subclone already present at the CLL stage. Transcriptomic analysis of PDXs showed a high oxidative phosphorylation (OXPHOS) and low B-cell receptor (BCR) signaling similar to the RT in the patients. IACS-010759, an OXPHOS inhibitor, reduced proliferation, and circumvented resistance to venetoclax. In summary, we have generated new RT-PDX models, one of them from CLL cells that mimicked the evolution of CLL to RT uncovering intrinsic features of RT cells of therapeutical value.

Clonal evolution dissection reveals that a high MSI2 level promotes chemoresistance in T-cell acute lymphoblastic leukemia.

ABSTRACT: T-cell acute lymphoblastic leukemia (T-ALL) is an aggressive cancer with resistant clonal propagation in recurrence. We performed high-throughput droplet-based 5' single-cell RNA with paired T-cell receptor (TCR) sequencing of paired diagnosis-relapse (Dx_Rel) T-ALL samples to dissect the clonal diversities. Two leukemic evolutionary patterns, "clonal shift" and "clonal drift" were unveiled. Targeted single-cell DNA sequencing of paired Dx_Rel T-ALL samples further corroborated the existence of the 2 contrasting clonal evolution patterns, revealing that dynamic transcriptional variation might cause the mutationally static clones to evolve chemotherapy resistance. Analysis of commonly enriched drifted gene signatures showed expression of the RNA-binding protein MSI2 was significantly upregulated in the persistent TCR clonotypes at relapse. Integrated in vitro and in vivo functional studies suggested that MSI2 contributed to the proliferation of T-ALL and promoted chemotherapy resistance through the posttranscriptional regulation of MYC, pinpointing MSI2 as an informative biomarker and novel therapeutic target in T-ALL.

Clonal evolution in inherited marrow failure syndromes predicts disease progression.

Progression to myelodysplastic syndromes (MDS) and acute myeloid leukemia is one of the most serious complications of the inherited bone marrow failure and MDS-predisposition syndromes. Given the lack of predictive markers, this risk can also be a source of great uncertainty and anxiety to patients and their providers alike. Recent data show that some acquired mutations may provide a window into this risk. While maladaptive mechanisms, such as monosomy 7, are associated with a high risk of leukemogenesis, mutations that offset the inherited defect (known as somatic genetic rescue) may attenuate this risk. Somatic mutations that are shared with age-acquired clonal hematopoiesis mutations also show syndrome-specific patterns that may provide additional data as to disease risk. This review focuses on recent progress in this area with an emphasis on the biological underpinnings and interpretation of these patterns for patient care decisions.

Somatic evolution of marine transmissible leukemias in the common cockle, Cerastoderma edule.

Transmissible cancers are malignant cell lineages that spread clonally between individuals. Several such cancers, termed bivalve transmissible neoplasia (BTN), induce leukemia-like disease in marine bivalves. This is the case of BTN lineages affecting the common cockle, Cerastoderma edule, which inhabits the Atlantic coasts of Europe and northwest Africa. To investigate the evolution of cockle BTN, we collected 6,854 cockles, diagnosed 390 BTN tumors, generated a reference genome and assessed genomic variation across 61 tumors. Our analyses confirmed the existence of two BTN lineages with hemocytic origins. Mitochondrial variation revealed mitochondrial capture and host co-infection events. Mutational analyses identified lineage-specific signatures, one of which likely reflects DNA alkylation. Cytogenetic and copy number analyses uncovered pervasive genomic instability, with whole-genome duplication, oncogene amplification and alkylation-repair suppression as likely drivers. Satellite DNA distributions suggested ancient clonal origins. Our study illuminates long-term cancer evolution under the sea and reveals tolerance of extreme instability in neoplastic genomes.