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.

Advances in Molecular Understanding of Polycythemia Vera, Essential Thrombocythemia, and Primary Myelofibrosis: Towards Precision Medicine

Myeloproliferative neoplasms (MPNs), including Polycythemia Vera (PV), Essential Thrombocythemia (ET), and Primary Myelofibrosis (PMF), are characterized by the clonal proliferation of hematopoietic stem cells leading to an overproduction of hematopoietic cells. The last two decades have seen significant advances in our understanding of the molecular pathogenesis of these diseases, with the discovery of key mutations in the JAK2, CALR, and MPL genes being pivotal. This review provides a comprehensive update on the molecular landscape of PV, ET, and PMF, highlighting the diagnostic, prognostic, and therapeutic implications of these genetic findings. We delve into the challenges of diagnosing and treating patients with prognostic mutations, clonal evolution, and the impact of emerging technologies like next-generation sequencing and single-cell genomics on the field. The future of MPN management lies in leveraging these molecular insights to develop personalized treatment strategies, aiming for precision medicine that optimizes outcomes for patients. This article synthesizes current knowledge on molecular diagnostics in MPNs, underscoring the critical role of genetic profiling in enhancing patient care and pointing towards future research directions that promise to further refine our approach to these complex disorders.

Stochastic modelling of single-cell gene expression adaptation reveals non-genomic contribution to evolution of tumor subclones

Cancer progression is an evolutionary process driven by the selection of cells adapted to gain growth advantage. We present the first formal study on the adaptation of gene expression in subclonal evolution. We model evolutionary changes in gene expression as stochastic Ornstein-Uhlenbeck processes, jointly leveraging the evolutionary history of subclones and single-cell expression data. Applying our model to sublines derived from single cells of a mouse melanoma revealed that sublines with distinct phenotypes are underlined by different patterns of gene expression adaptation, indicating non-genetic mechanisms of cancer evolution. Interestingly, sublines previously observed to be resistant to anti-CTLA-4 treatment showed adaptive expression of genes related to invasion and non-canonical Wnt signaling, whereas sublines that responded to treatment showed adaptive expression of genes related to proliferation and canonical Wnt signaling. Our results suggest that clonal phenotypes emerge as the result of specific adaptivity patterns of gene expression.

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.

Differential Competitive Growth of Transgenic Subclones of Neuroblastoma Cells Expressing Different Levels of Cathepsin D Co-Cultured in 2D and 3D in Response to EGF: Implications in Tumor Heterogeneity and Metastasis

Neuroblastoma (NB) is an embryonal tumor arising from the sympathetic central nervous system. The epidermal growth factor (EGF) plays a role in NB growth and metastatic behavior. Recently, we have demonstrated that cathepsin D (CD) contrasts EGF-induced NB cell growth in 2D by downregulating EGFR/MAPK signaling. Aggressive NB is highly metastatic to the bone and the brain. In the metastatic process, adherent cells detach to form clusters of suspended cells that adhere once they reach the metastatic site and form secondary colonies. Whether CD is involved in the survival of metastatic NB clones is not known. Therefore, in this study, we addressed how CD differentially affects cell growth in suspension versus the adherent condition. To mimic tumor heterogeneity, we co-cultured transgenic clones silenced for or overexpressing CD. We compared the growth kinetics of such mixed clones in 2D and 3D models in response to EGF, and we found that the Over CD clone had an advantage for growth in suspension, while the CD knocked-down clone was favored for the adherent growth in 2D. Interestingly, on switching from 3D to 2D culture conditions, the expression of E-cadherin and of N-cadherin increased in the KD-CD and Over CD clones, respectively. The fact that CD plays a dual role in cancer cell growth in 2D and 3D conditions indicates that during clonal evolution, subclones expressing different level of CD may arise, which confers survival and growth advantages depending on the metastatic step. By searching the TCGA database, we found up to 38 miRNAs capable of downregulating CD. Interestingly, these miRNAs are associated with biological processes controlling cell adhesion and cell migration. The present findings support the view that during NB growth on a substrate or when spreading as floating neurospheres, CD expression is epigenetically modulated to confer survival advantage. Thus, epigenetic targeting of CD could represent an additional strategy to prevent NB metastases.

Clonal origin and genomic diversity in Lynch syndrome-associated endometrial cancer with multiple synchronous tumors: Identification of the pathogenicity of MLH1 p.L582H

Lynch syndrome-associated endometrial cancer patients often present multiple synchronous tumors and this assessment can affect treatment strategies. We present a case of a 27-year-old woman with tumors in the uterine corpus, cervix, and ovaries who was diagnosed with endometrial cancer and exhibited cervical invasion and ovarian metastasis. Her family history suggested Lynch syndrome, and genetic testing identified a variant of uncertain significance, MLH1 p.L582H. We conducted immunohistochemical staining, microsatellite instability analysis, and Sanger sequencing for Lynch syndrome-associated cancers in three generations of the family and identified consistent MLH1 loss. Whole-exome sequencing for the corpus, cervical, and ovarian tumors of the proband identified a copy-neutral loss of heterozygosity (LOH) occurring at the MLH1 position in all tumors. This indicated that the germline variant and the copy-neutral LOH led to biallelic loss of MLH1 and was the cause of cancer initiation. All tumors shared a portion of somatic mutations with high mutant allele frequencies, suggesting a common clonal origin. There were no mutations shared only between the cervix and ovary samples. The profiles of mutant allele frequencies shared between the corpus and cervix or ovary indicated that two different subclones originating from the corpus independently metastasized to the cervix or ovary. Additionally, all tumors presented unique mutations in endometrial cancer-associated genes such as ARID1A and PIK3CA. In conclusion, we demonstrated clonal origin and genomic diversity in a Lynch syndrome-associated endometrial cancer, suggesting the importance of evaluating multiple sites in Lynch syndrome patients with synchronous tumors.

Pathogenic Variants Associated with Epigenetic Control and the NOTCH Pathway Are Frequent in Classic Hodgkin Lymphoma

Classic Hodgkin lymphoma (cHL) constitutes a B-cell neoplasm derived from germinal center lymphocytes. Despite high cure rates (80-90%) obtained with the current multiagent protocols, a significant proportion of cHL patients experience recurrences, characterized by a lower sensitivity to second-line treatments. The genomic background of chemorefractory cHL is still poorly understood, limiting personalized treatment strategies based on molecular features. In this study, using a targeted next-generation sequencing (NGS) panel specifically designed for cHL research, we compared chemosensitive and chemorefractory diagnostic tissue samples of cHL patients. Furthermore, we longitudinally examined paired diagnosis-relapsesamples of chemorefractory cHL in order to define patterns of dynamic evolution and clonal selection. Pathogenic variants in NOTCH1 and NOTCH2 genes frequently arise in cHL. Mutations in genes associated with epigenetic regulation (CREBBP and EP300) are particularly frequent in relapsed/refractory cHL. The appearance of novel clones characterized by mutations previously not identified at diagnosis is a common feature in cHL cases showing chemoresistance to frontline treatments. Our results expand current molecular and pathogenic knowledge of cHL and support the performance of molecular studies in cHL prior to the initiation of first-line therapies.

From primary myelofibrosis to chronic myeloid leukemia, BCR::ABL1+ B-Lymphoblastic leukemia, and back to primary myelofibrosis: An illustration of dynamic clonal evolution

The simultaneous detection of BCR::ABL1 and JAK2 V617F was rarely reported and their clonal relationship and dynamic clonal shift were not characterized. Here, we described a unique case with the initial presentation as JAK2 V617F+ primary myelofibrosis, followed by the emergence of BCR::ABL1+ chronic myeloid leukemia. The patient then developed BCR::ABL1+ B-lymphoblastic leukemia. Treatment for B-lymphoblastic leukemia prompted a regression to the state of primary myelofibrosis. In light of these observations, we proposed a clonal evolution model for this case.

Longitudinal genetically detectable minimal residual disease by fluorescence in situ hybridization confers a poor prognosis in myeloma

Background

Deeper depth of response (DpR) after induction therapy, especially gain of negative minimal residual disease (MRD), has been linked to prolonged survival in multiple myeloma (MM). However, flow-MRD examination focuses on the numbers but not on the biological characteristics of residual plasma cells (PCs).

Objectives

To explore whether the genetic features of residual tumor cells affect the survival time of patients with MM.

Design

A retrospective cohort study.

Methods

We investigated the clonality of cytogenetic abnormalities (CAs) of the residual PCs using interphase fluorescence in situ hybridization (iFISH) in the National Longitudinal Cohort of Hematological Diseases in China (NCT04645199). Here, a longitudinal cohort of 269 patients with patient-paired diagnostic and post-induction iFISH results was analyzed.

Results

Persistent CAs after induction therapy were detected in about half of the patients (118/269, 43%), and patients with undetectable CAs showed significantly improved survival compared with those with genetically detectable MRD [median progression-free survival (mPFS): 59.7 versus 35.7 months, p < 0.001; median overall survival (mOS): 97.1 versus 68.8 months, p = 0.011]. In addition, different patterns of therapy-induced clonal evolution were observed by comparing the clonal structure of residual PCs with paired baseline samples. Patients who maintained at a high risk during follow-up had the worst survival (mPFS: 30.5 months; mOS: 54.4 months), while those who returned to lower risk or had iFISH- at both time points had the best survival (mPFS: 62.0 months, mOS: not reached).

Conclusion

These findings highlighted the prognostic value of genetic testing in residual tumor cells, which may provide a deep understanding of clonal evolution and guide clinical therapeutic strategies.

Cutaneous Manifestations of Myeloid Neoplasms Exhibit Broad and Divergent Morphologic and Immunophenotypic Features but Share Ancestral Clonal Mutations With Bone Marrow

In this study, we performed a comprehensive molecular analysis of paired skin and peripheral blood/bone marrow (BM) samples from 17 patients with cutaneous myeloid or cutaneous histiocytic-dendritic neoplasms. The cutaneous manifestations included 10 patients with cutaneous acute myeloid leukemia (c-AML), 2 patients with full or partial Langerhans cell differentiation, 2 patients with blastic plasmacytoid dendritic cell neoplasms (BPDCN), 1 patient with both Langerhans cell differentiation and BPDCN, and 2 patients with full or partial indeterminate dendritic cell differentiation. Seven of the 10 c-AML patients (70%) exhibited concurrent or subsequent marrow involvement by acute myeloid leukemia, with all 7 cases (100%) demonstrating shared clonal mutations in both the skin and BM. However, clonal relatedness was documented in one additional case that never had any BM involvement. Nevertheless, NPM1 mutations were identified in 7 of the 10 (70%) of these c-AML cases while one had KMT2A rearrangement and one showed inv(16). All 3 patients (100%) with Langerhans cell neoplasms, 2 patients with BPDCN (100%), and one of the 2 patients (50%) with other cutaneous dendritic cell neoplasms also demonstrated shared mutations between the skin and concurrent or subsequent myeloid neoplasms. Both BM and c-AML shared identical founding drivers, with a predominance of NPM1, DNMT3A, and translocations associated with monocytic differentiation, with common cutaneous-only mutations involving genes in the signal transduction and epigenetic pathways. Cutaneous histiocytic-dendritic neoplasms shared founding drivers in ASXL1, TET2, and/or SRSF2. However, in the Langerhans cell histiocytosis or histiocytic sarcoma cases, there exist recurrent secondary RAS pathway hits, whereas cutaneous BPDCN cases exhibit copy number or structural variants. These results enrich and broaden our understanding of clonally related cutaneous manifestations of myeloid neoplasms and further illuminate the highly diverse spectrum of morphologic and immunophenotypic features they exhibit.

Integrative NGS testing reveals clonal dynamics of adverse genomic defects contributing to a natural progression in treatment-naïve CLL patients

Large-scale next-generation sequencing (NGS) studies revealed extensive genetic heterogeneity, driving a highly variable clinical course of chronic lymphocytic leukaemia (CLL). The evolution of subclonal populations contributes to diverse therapy responses and disease refractoriness. Besides, the dynamics and impact of subpopulations before therapy initiation are not well understood. We examined changes in genomic defects in serial samples of 100 untreated CLL patients, spanning from indolent to aggressive disease. A comprehensive NGS panel LYNX, which provides targeted mutational analysis and genome-wide chromosomal defect assessment, was employed. We observed dynamic changes in the composition and/or proportion of genomic aberrations in most patients (62%). Clonal evolution of gene variants prevailed over the chromosomal alterations. Unsupervised clustering based on aberration dynamics revealed four groups of patients with different clinical behaviour. An adverse cluster was associated with fast progression and early therapy need, characterized by the expansion of TP53 defects, ATM mutations, and 18p- alongside dynamic SF3B1 mutations. Our results show that clonal evolution is active even without therapy pressure and that repeated genetic testing can be clinically relevant during long-term patient monitoring. Moreover, integrative NGS testing contributes to the consolidated evaluation of results and accurate assessment of individual patient prognosis.

SETD2 mutations do not contribute to clonal fitness in response to chemotherapy in childhood B cell acute lymphoblastic leukemia

Mutations in genes encoding epigenetic regulators are commonly observed at relapse in B cell acute lymphoblastic leukemia (B-ALL). Loss-of-function mutations in SETD2, an H3K36 methyltransferase, have been observed in B-ALL and other cancers. Previous studies on mutated SETD2 in solid tumors and acute myelogenous leukemia support a role in promoting resistance to DNA damaging agents. We did not observe chemoresistance, an impaired DNA damage response, nor increased mutation frequency in response to thiopurines using CRISPR-mediated knockout in wild-type B-ALL cell lines. Likewise, restoration of SETD2 in cell lines with hemizygous mutations did not increase sensitivity. SETD2 mutations affected the chromatin landscape and transcriptional output that was unique to each cell line. Collectively our data does not support a role for SETD2 mutations in driving clonal evolution and relapse in B-ALL, which is consistent with the lack of enrichment of SETD2 mutations at relapse in most studies.

Exploring the intratumoral heterogeneity of DNA ploidy in prostate cancer

BACKGROUND

Prostate cancer is morphologically and molecularly heterogeneous. Genomic heterogeneity might be mirrored by variability in DNA ploidy. Aneuploidy is a hallmark of genomic instability and associated with tumor aggressiveness. Little attention has been paid to the biological significance of the diploid tumor cell population that often coexists with aneuploid populations. Here, we investigated the role of DNA ploidy in tumor heterogeneity and clonal evolution.

METHODS

Three radical prostatectomy specimens with intratumoral heterogeneity based on nuclear features on H&E were selected. DNA content of each subpopulation was determined by DNA image cytometry and silver in situ hybridization (SISH). Genomic evolution was inferred from array comparative genomic hybridization (aCGH). Additionally, immunohistochemistry was used to examine the stemness-associated marker ALDH1A1.

RESULTS

Nuclear morphology reliably predicted DNA ploidy status in all three cases. In one case, aCGH analysis revealed several shared deletions and one amplification in both the diploid and the aneuploid population, suggesting that these populations could be related. In the other two cases, a statement about relatedness was not possible. Furthermore, ALDH1A1 was expressed in 2/3 cases and exclusively observed in their diploid populations.

CONCLUSIONS

In this proof-of-concept study, we demonstrate the feasibility to predict the DNA ploidy status of distinct populations within one tumor by H&E morphology. Future studies are needed to further investigate the clonal relationship between the diploid and the aneuploid subpopulation and test the hypothesis that the aneuploid population is derived from the diploid one. Finally, our analyses pointed to an enrichment of the stemness-associated marker ALDH1A1 in diploid populations, which warrants further investigation in future studies.

Overview of Tumor Heterogeneity in High-Grade Serous Ovarian Cancers

Ovarian cancers encompass a group of neoplasms originating from germinal tissues and exhibiting distinct clinical, pathological, and molecular features. Among these, epithelial ovarian cancers (EOCs) are the most prevalent, comprising five distinct tumor histotypes. Notably, high-grade serous ovarian cancers (HGSOCs) represent the majority, accounting for over 70% of EOC cases. Due to their silent and asymptomatic behavior, HGSOCs are generally diagnosed in advanced stages with an evolved and complex genomic state, characterized by high intratumor heterogeneity (ITH) due to chromosomal instability that distinguishes HGSOCs. Histologically, these cancers exhibit significant morphological diversity both within and between tumors. The histologic patterns associated with solid, endometrioid, and transitional (SET) and classic subtypes of HGSOCs offer prognostic insights and may indicate specific molecular profiles. The evolution of HGSOC from primary to metastasis is typically characterized by clonal ITH, involving shared or divergent mutations in neoplastic sub-clones within primary and metastatic sites. Disease progression and therapy resistance are also influenced by non-clonal ITH, related to interactions with the tumor microenvironment and further genomic changes. Notably, significant alterations occur in nonmalignant cells, including cancer-associated fibroblast and immune cells, during tumor progression. This review provides an overview of the complex nature of HGSOC, encompassing its various aspects of intratumor heterogeneity, histological patterns, and its dynamic evolution during progression and therapy resistance.

Temporal and structural patterns of hepatitis B virus integrations in hepatocellular carcinoma

Chronic infection of hepatitis B virus (HBV) is the major cause of hepatocellular carcinoma (HCC). Notably, 90% of HBV-positive HCC cases exhibit detectable HBV integrations, hinting at the potential early entanglement of these viral integrations in tumorigenesis and their subsequent oncogenic implications. Nevertheless, the precise chronology of integration events during HCC tumorigenesis, alongside their sequential structural patterns, has remained elusive thus far. In this study, we applied whole-genome sequencing to multiple biopsies extracted from six HBV-positive HCC cases. Through this approach, we identified point mutations and viral integrations, offering a blueprint for the intricate tumor phylogeny of these samples. The emergent narrative paints a rich tapestry of diverse evolutionary trajectories characterizing the analyzed tumors. We uncovered oncogenic integration events in some samples that appear to happen before and during the initiation stage of tumor development based on their locations in reconstituted trajectories. Furthermore, we conducted additional long-read sequencing of selected samples and unveiled integration-bridged chromosome rearrangements and tandem repeats of the HBV sequence within integrations. In summary, this study revealed premalignant oncogenic and sequential complex integrations and highlighted the contributions of HBV integrations to HCC development and genome instability.

Reuniting philosophy and science to advance cancer research

Cancers rely on multiple, heterogeneous processes at different scales, pertaining to many biomedical fields. Therefore, understanding cancer is necessarily an interdisciplinary task that requires placing specialised experimental and clinical research into a broader conceptual, theoretical, and methodological framework. Without such a framework, oncology will collect piecemeal results, with scant dialogue between the different scientific communities studying cancer. We argue that one important way forward in service of a more successful dialogue is through greater integration of applied sciences (experimental and clinical) with conceptual and theoretical approaches, informed by philosophical methods. By way of illustration, we explore six central themes: (i) the role of mutations in cancer; (ii) the clonal evolution of cancer cells; (iii) the relationship between cancer and multicellularity; (iv) the tumour microenvironment; (v) the immune system; and (vi) stem cells. In each case, we examine open questions in the scientific literature through a philosophical methodology and show the benefit of such a synergy for the scientific and medical understanding of cancer.

Molecular and immunological mechanisms of clonal evolution in multiple myeloma

Multiple myeloma (MM) is a hematologic malignancy characterized by the proliferation of clonal plasma cells in the bone marrow (BM). It is known that early genetic mutations in post-germinal center B/plasma cells are the cause of myelomagenesis. The acquisition of additional chromosomal abnormalities and distinct mutations further promote the outgrowth of malignant plasma cell populations that are resistant to conventional treatments, finally resulting in relapsed and therapy-refractory terminal stages of MM. In addition, myeloma cells are supported by autocrine signaling pathways and the tumor microenvironment (TME), which consists of diverse cell types such as stromal cells, immune cells, and components of the extracellular matrix. The TME provides essential signals and stimuli that induce proliferation and/or prevent apoptosis. In particular, the molecular pathways by which MM cells interact with the TME are crucial for the development of MM. To generate successful therapies and prevent MM recurrence, a thorough understanding of the molecular mechanisms that drive MM progression and therapy resistance is essential. In this review, we summarize key mechanisms that promote myelomagenesis and drive the clonal expansion in the course of MM progression such as autocrine signaling cascades, as well as direct and indirect interactions between the TME and malignant plasma cells. In addition, we highlight drug-resistance mechanisms and emerging therapies that are currently tested in clinical trials to overcome therapy-refractory MM stages.

Clonal evolution and hierarchy in myeloid malignancies

Myeloid malignancies, a group of hematopoietic disorders that includes acute myeloid leukemia (AML), myelodysplastic syndromes (MDS), and myeloproliferative neoplasms (MPNs), are caused by the accumulation of genetic and epigenetic changes in hematopoietic stem and progenitor cells (HSPCs) over time. Despite the relatively low number of genomic drivers compared with other forms of cancer, the process by which these changes shape the genomic architecture of myeloid malignancies remains elusive. Recent advancements in clonal hematopoiesis research and the use of cutting-edge single cell technologies have shed new light on the developmental process of myeloid malignancies. In this review, we delve into the intricacies of clonal evolution in myeloid malignancies and its implications for the development of new diagnostic and therapeutic approaches.

Divergent evolution of metachronous follicular lymphoma and extranodal marginal zone lymphoma of mucosa-associated lymphoid tissue from a common precursor

The translocation t(14;18)(q32:q21)/IGH::BCL2 occurs at the pre-B stage of B-cell development in the bone marrow and is insufficient for malignant transformation, although it leads to the formation of in situ follicular B-cell neoplasia (ISFN). Despite that, the translocation is the genetic hallmark of follicular lymphoma (FL), it occurs infrequently in metachronous/synchronous lymphomas, including extranodal marginal zone lymphoma of mucosa-associated lymphoid tissue (EMZL), mantle cell lymphoma, and Hodgkin's lymphoma. In each of these scenarios, the two lymphomas often appear to be clonally related by analyses of IGH::BCL2 and/or rearranged IG genes. However, it remains largely unknown whether one lymphoma originates from the other or they develop independently. We studied five cases of metachronous EMZL and FL. In four cases, the two lymphomas were clonally related, as shown by identical IGH::BCL2 and/or rearranged IG genes or shared mutations. There were common and unique mutations between the paired EMZL and FL, indicating that they developed independently from a common premalignant cell population, harbouring IGH::BCL2 in three cases. Furthermore, case 1 presented with three metachronous FLs, and all of them originated from a common precursor cell population via divergent evolution. Our findings highlight the multi-malignant potential of IGH::BCL2-positive B-cells. © 2023 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

Unravelling the Complexity of Colorectal Cancer: Heterogeneity, Clonal Evolution, and Clinical Implications

Colorectal cancer (CRC) is a global health concern and a leading cause of death worldwide. The disease's course and response to treatment are significantly influenced by its heterogeneity, both within a single lesion and between primary and metastatic sites. Biomarkers, such as mutations in KRAS, NRAS, and BRAF, provide valuable guidance for treatment decisions in patients with metastatic CRC. While high concordance exists between mutational status in primary and metastatic lesions, some heterogeneity may be present. Circulating tumor DNA (ctDNA) analysis has proven invaluable in identifying genetic heterogeneity and predicting prognosis in RAS-mutated metastatic CRC patients. Tumor heterogeneity can arise from genetic and non-genetic factors, affecting tumor development and response to therapy. To comprehend and address clonal evolution and intratumoral heterogeneity, comprehensive genomic studies employing techniques such as next-generation sequencing and computational analysis are essential. Liquid biopsy, notably through analysis of ctDNA, enables real-time clonal evolution and treatment response monitoring. However, challenges remain in standardizing procedures and accurately characterizing tumor subpopulations. Various models elucidate the origin of CRC heterogeneity, highlighting the intricate molecular pathways involved. This review focuses on intrapatient cancer heterogeneity and genetic clonal evolution in metastatic CRC, with an emphasis on clinical applications.

The Contribution of Multiplexing Single Cell RNA Sequencing in Acute Myeloid Leukemia

Decades ago, the treatment for acute myeloid leukemia relied on cytarabine and anthracycline. However, advancements in medical research have introduced targeted therapies, initially employing monoclonal antibodies such as ant-CD52 and anti-CD123, and subsequently utilizing specific inhibitors that target molecular mutations like anti-IDH1, IDH2, or FLT3. The challenge lies in determining the role of these therapeutic options, considering the inherent tumor heterogeneity associated with leukemia diagnosis and the clonal drift that this type of tumor can undergo. Targeted drugs necessitate an examination of various therapeutic targets at the individual cell level rather than assessing the entire population. It is crucial to differentiate between the prognostic value and therapeutic potential of a specific molecular target, depending on whether it is found in a terminally differentiated cell with limited proliferative potential or a stem cell with robust capabilities for both proliferation and self-renewal. However, this cell-by-cell analysis is accompanied by several challenges. Firstly, the scientific aspect poses difficulties in comparing different single cell analysis experiments despite efforts to standardize the results through various techniques. Secondly, there are practical obstacles as each individual cell experiment incurs significant financial costs and consumes a substantial amount of time. A viable solution lies in the ability to process multiple samples simultaneously, which is a distinctive feature of the cell hashing technique. In this study, we demonstrate the applicability of the cell hashing technique for analyzing acute myeloid leukemia cells. By comparing it to standard single cell analysis, we establish a strong correlation in various parameters such as quality control, gene expression, and the analysis of leukemic blast markers in patients. Consequently, this technique holds the potential to become an integral part of the biological assessment of acute myeloid leukemia, contributing to the personalized and optimized management of the disease, particularly in the context of employing targeted therapies.

Evolutionary consequences of loss of sexual reproduction on male-related traits in parthenogenetic lineages of the pea aphid

Transition from sexual reproduction to parthenogenesis constitutes a major life-history change with deep evolutionary consequences for sex-related traits, which are expected to decay. The pea aphid Acyrthosiphon pisum shows intraspecific reproductive polymorphism, with cold-resistant cyclically parthenogenetic (CP) lineages that alternate sexual and asexual generations and cold-sensitive obligately parthenogenetic (OP) lineages that produce only asexual females but still males. Here, the genotyping of 219 pea aphid lineages collected in cold-winter and mild-winter regions revealed contrasting population structures. Samples from cold-winter regions consisted mostly of distinct multilocus genotypes (MLGs) usually represented by a single sample (101 different MLGs for 111 samples) and were all phenotyped as CP. In contrast, fewer MLGs were found in mild-winter regions (28 MLGs for 108 samples), all but one being OP. Since the males produced by OP lineages are unlikely to pass on their genes (sexual females being rare in mild-winter regions), we tested the hypothesis that their traits could degenerate due to lack of selection by comparing male production and male reproductive success between OP and CP lineages. Male production was indeed reduced in OP lineages, but a less clear pattern was observed for male reproductive success: females mated with OP males laid fewer eggs (fertilized or not) but OP and CP males fertilized the same proportion of eggs. These differences may stem from the type of selective forces: male production may be counter-selected whereas male performances may evolve under the slower process of relaxed selection. The overall effective reproductive capacity of OP males could result from recent sex loss in OP lineages or underestimated reproductive opportunities.

Clonal evolution in aplastic anemia: failed tumor surveillance or maladaptive recovery?

Clonal evolution to secondary paroxysmal nocturnal hemoglobinuria (PNH) or myeloid neoplasia (MN) represents one of the long-term complications of patients with aplastic anemia (AA). The recent evidence in the field of immunology and the application of next-generation sequencing have shed light on the molecular underpinnings of these clonal complications, revealing clinical and molecular risk factors as well as potential immunological players. Particularly, whether MN evolution represents a failed tumor surveillance or a maladaptive recovery is still a matter of controversy in the field of bone marrow failure syndromes. However, recent studies have explored the precise dynamics of the immune-molecular forces governing such processes over time, generating knowledge useful for potential early therapeutic strategies. In this review, we will discuss the immune pathophysiology of AA and the emergence of clonal hematopoiesis with regard to the adaptive and maladaptive mechanisms at the basis of secondary evolution trajectories operating under the immune pressure.

Case Report: Longitudinal monitoring of clonal evolution by circulating tumor DNA for resistance to anti-EGFR antibody in a case of metastatic colorectal cancer

Background

Treatment with anti-EGFR antibody has been shown to prolong survival in patients with RAS wild-type metastatic colorectal cancer (mCRC). However, even patients who initially respond to anti-EGFR antibody therapy, almost without exception, develop resistance to the therapy and then fail to respond. Secondary mutations in the mitogen-activated protein (MAPK) signaling pathway (mainly in NRAS and BRAF) have been implicated in anti-EGFR resistance. However, the process by which resistant clones develop during therapy has not been elucidated, and considerable intrapatient and interpatient heterogeneity exists. Circulating tumor DNA (ctDNA) testing has recently allowed the noninvasive detection of heterogeneous molecular alterations that underlie the evolution of resistance to anti-EGFR. In this report, we describe our observation of genomic alterations in KRAS and NRAS in a patient with acquired resistance to anti-EGFR antibody drugs by tracking clonal evolution using serial ctDNA anaylsis.

Case presentation

A 54-year-old woman was initially diagnosed with sigmoid colon cancer with multiple liver metastases. After receiving first-line mFOLFOX + cetuximab, second-line FOLFIRI + ramucirumab, third-line trifluridine/tipiracil + bevacizumab, fourth-line regorafenib, and fifth-line CAPOX + bevacizumab, she was rechallenged with CPT-11 + cetuximab. The best response to anti-EGFR rechallenge therapy was a partial response. RAS in the ctDNA was assessed during treatment. The RAS status changed from wild type to mutant type, back to wild type, and again to mutant type (NRAS/KRAS codon 61) during the course of treatment.

Conclusion

In this report, tracking of ctDNA allowed us to describe clonal evolution in a case in which we observed genomic alterations in KRAS and NRAS in a patient who acquired resistance to anti-EGFR antibody drugs during treatment. It is reasonable to consider repeat molecular interrogation during progression in patients with mCRC by using ctDNA analysis, which could help to identify patients who may benefit from a rechallenge strategy.

Genetic, Phenotypic, and Clinical Heterogeneity of NPM1-Mutant Acute Myeloid Leukemias

The current classification of acute myeloid leukemia (AML) relies largely on genomic alterations. AML with mutated nucleophosmin 1 (NPM1-mut) is the largest of the genetically defined groups, involving about 30% of adult AMLs and is currently recognized as a distinct entity in the actual AML classifications. NPM1-mut AML usually occurs in de novo AML and is associated predominantly with a normal karyotype and relatively favorable prognosis. However, NPM1-mut AMLs are genetically, transcriptionally, and phenotypically heterogeneous. Furthermore, NPM1-mut is a clinically heterogenous group. Recent studies have in part clarified the consistent heterogeneities of these AMLs and have strongly supported the need for an additional stratification aiming to improve the therapeutic response of the different subgroups of NPM1-mut AML patients.

Maintenance and expansion of genetic and trait variation following domestication in a clonal crop

Clonal propagation enables favourable crop genotypes to be rapidly selected and multiplied. However, the absence of sexual propagation can lead to low genetic diversity and accumulation of deleterious mutations, which may eventually render crops less resilient to pathogens or environmental change. To better understand this trade-off, we characterize the domestication and contemporary genetic diversity of Enset (Ensete ventricosum), an indigenous African relative of bananas (Musa) and a principal starch staple for 20 million Ethiopians. Wild enset reproduction occurs strictly by sexual outcrossing, but for cultivation, it is propagated clonally and associated with diversification and specialization into hundreds of named landraces. We applied tGBS sequencing to generate genome-wide genotypes for 192 accessions from across enset's cultivated distribution, and surveyed 1340 farmers on enset agronomic traits. Overall, reduced heterozygosity in the domesticated lineage was consistent with a domestication bottleneck that retained 37% of wild diversity. However, an excess of putatively deleterious missense mutations at low frequency present as heterozygotes suggested an accumulation of mutational load in clonal domesticated lineages. Our evidence indicates that the major domesticated lineages initially arose through historic sexual recombination associated with a domestication bottleneck, followed by the amplification of favourable genotypes through an extended period of clonal propagation. Among domesticated lineages, we found a significant phylogenetic signal for multiple farmer-identified food, nutrition and disease resistance traits and little evidence of contemporary recombination. The development of future-climate adapted genotypes may require crop breeding, but outcrossing risks exposing deleterious alleles as homozygotes. This trade-off may partly explain the ubiquity and persistence of clonal propagation over recent centuries of comparative climate stability.

Novel NUP98::ASH1L Gene Fusion in Acute Myeloid Leukemia Detected by Optical Genome Mapping

Optical genome mapping (OGM) recently has demonstrated the potential to improve genetic diagnostics in acute myeloid leukemia (AML). In this study, OGM was utilized as a tool for the detection of genome-wide structural variants and disease monitoring. A previously unrecognized NUP98::ASH1L fusion was detected in an adult patient with secondary AML. OGM identified the fusion of NUP98 to Absent, Small, or Homeotic-Like Histone Lysine Methyltransferase (ASH1L) as result of a complex structural rearrangement between chromosomes 1 and 11. A pipeline for the measurement of rare structural variants (Rare Variant Pipeline, Bionano Genomics, San Diego, CA, USA) was used for detection. As NUP98 and other fusions are relevant for disease classification, this demonstrates the necessity for methods such as OGM for cytogenetic diagnostics in AML. Furthermore, other structural variants showed discordant variant allele frequencies at different time points over the course of the disease and treatment pressure, indicating clonal evolution. These results support OGM to be a valuable tool for primary diagnostics in AML as well as longitudinal testing for disease monitoring and deepening our understanding of genetically heterogenous diseases.

An unusual ectopic thymoma clonal evolution analysis: A case report

Thymomas and thymic carcinomas are rare and primary tumors of the mediastinum which is derived from the thymic epithelium. Thymomas are the most common primary anterior mediastinal tumor, while ectopic thymomas are rarer. Mutational profiles of ectopic thymomas may help expand our understanding of the occurrence and treatment options of these tumors. In this report, we sought to elucidate the mutational profiles of two ectopic thymoma nodules to gain deeper understanding of the molecular genetic information of this rare tumor and to provide guidance treatment options. We presented a case of 62-year-old male patient with a postoperative pathological diagnosis of type A mediastinal thymoma and ectopic pulmonary thymoma. After mediastinal lesion resection and thoracoscopic lung wedge resection, the mediastinal thymoma was completely removed, and the patient recovered from the surgery and no recurrence was found by examination until now. Whole exome sequencing was performed on both mediastinal thymoma and ectopic pulmonary thymoma tissue samples of the patient and clonal evolution analysis were further conducted to analyze the genetic characteristics. We identified eight gene mutations that were co-mutated in both lesions. Consistent with a previous exome sequencing analysis of thymic epithelial tumor, HRAS was also observed in both mediastinal lesion and lung lesion tissues. We also evaluated the intratumor heterogeneity of non-silent mutations. The results showed that the mediastinal lesion tissue has higher degree of heterogeneity and the lung lesion tissue has relatively low amount of variant heterogeneity in the detected variants. Through pathology and genomics sequencing detection, we initially revealed the genetic differences between mediastinal thymoma and ectopic thymoma, and clonal evolution analysis showed that these two lesions originated from multi-ancestral regions.

Molecular-defined clonal evolution in patients with classical myeloproliferative neoplasms

Classical myeloproliferative neoplasms (MPNs) are characterized by distinct clinical phenotypes. The discovery of driver mutations in JAK2, CALR and MPL genes provided new insights into their pathogenesis. Next-generation sequencing (NGS) identified additional somatic mutations, most frequently in epigenetic modulator genes. In this study, a cohort of 95 MPN patients was genetically characterized using targeted NGS. Clonal hierarchies of detected mutations were subsequently analysed using colony forming progenitor assays derived from single cells to study mutation acquisition. Further, the hierarchy of mutations within distinct cell lineages was evaluated. NGS revealed mutations in three epigenetic modulator genes (TET2, DNMT3A, ASXL1) as most common co-mutations to the classical driver mutations. JAK2^V617F as well as DNMT3A and TET2 mutations were detected as primary events in disease formation and most cases presented with a linear mutation pattern. Mutations appear mostly in the myeloid lineages but can also appear in lymphoid subpopulations. In one case with a double mutant MPL gene, mutations exclusively appeared in the monocyte lineage. Overall, this study confirms the mutational heterogeneity of classical MPNs and highlights the role of JAK2^V617F and epigenetic modifier genes as early events in hematologic disease formation.

Evolutionary landscape of clonal hematopoiesis in 3,359 individuals from the general population

Knowledge about evolution of clonal hematopoiesis, which may drive malignant progression, is crucial for clinical decision-making. We investigated the landscape of clonal evolution by error-corrected sequencing on 7,045 sequential samples from 3,359 individuals in the prospective population-based Lifelines cohort, with a special focus on cytosis and cytopenia. Spliceosome (SRSF2/U2AF1/SF3B1) and JAK2 mutated clones show highest growth rates over a median 3.6-year period, while clone sizes for DNMT3A and TP53 increase only marginally, independent of cytosis or cytopenia. Nevertheless, large differences are observed between individuals carrying the same mutation, indicative of modulation by non-mutation-related factors. Clonal expansion is not dependent on classical cancer risk factors (e.g., smoking). Risk for incident myeloid malignancy diagnosis is highest for JAK2, spliceosome, or TP53 mutations and absent for DNMT3A, and it is mostly preceded by cytosis or cytopenia. The results provide important insight into high-risk evolutionary patterns to guide monitoring of "CHIP" and "CCUS."

Clinical and genomic characterization of mutational signatures across human cancers

Mutational signatures, the generic patterns of mutations, are the footprints of both endogenous and exogenous factors that have influenced cancer development. To date, dozens of mutational signatures have been discerned through computational methods. However, the etiology, mutational properties, clonality, immunology and prognostic value of mutation signatures across cancer types are poorly understood. To address this, we extensively characterized mutational signatures across 8836 cancer samples spanning 42 cancer types. We confirmed and extended clinical and genomic features associated with mutation signatures. Mutation distribution analysis showed that most mutation processes were depleted in exons and APOBEC signatures (SBS2 and SBS13), the Pol-η related signature (SBS9) and SBS40 tended to contribute clustered mutations. We observed that age-related signatures (SBS1 and SBS5) and SBS40 tended to induce mutations affecting cancer genes and subclonal drivers posted by specific signatures (eg, mismatch repair deficiency-related signature SBS44) were unlikely subjected to positive selection. We also revealed early mutation signatures (eg, UV light exposure-related signature SBS7a) and signatures (eg, reactive oxygen species-related signature SBS18) predominated in the late stage of tumorigenesis. Comprehensive association analysis of mutation processes with microenvironment revealed that APOBEC- and mismatch repair deficiency-related signatures were positively associated with immune parameters, while age-related signatures showed negative correlations. In addition, prognostic association analysis showed that many signatures were favorable (eg, SBS9) or adverse factors (eg, SBS18) of patient survival. Our findings enhance appreciation of the role of mutational signatures in tumor evolution and underline their potential in immunotherapy guidance and prognostic prediction.

Reconstructing Clonal Evolution-A Systematic Evaluation of Current Bioinformatics Approaches

The accurate reconstruction of clonal evolution, including the identification of newly developing, highly aggressive subclones, is essential for the application of precision medicine in cancer treatment. Reconstruction, aiming for correct variant clustering and clonal evolution tree reconstruction, is commonly performed by tedious manual work. While there is a plethora of tools to automatically generate reconstruction, their reliability, especially reasons for unreliability, are not systematically assessed. We developed clevRsim-an approach to simulate clonal evolution data, including single-nucleotide variants as well as (overlapping) copy number variants. From this, we generated 88 data sets and performed a systematic evaluation of the tools for the reconstruction of clonal evolution. The results indicate a major negative influence of a high number of clones on both clustering and tree reconstruction. Low coverage as well as an extreme number of time points usually leads to poor clustering results. An underlying branched independent evolution hampers correct tree reconstruction. A further major decline in performance could be observed for large deletions and duplications overlapping single-nucleotide variants. In summary, to explore the full potential of reconstructing clonal evolution, improved algorithms that can properly handle the identified limitations are greatly needed.

Integrated analysis of genotype and phenotype reveals clonal evolution and cytogenetically driven disruption of myeloid cell maturation in myelodysplastic syndromes

BACKGROUND

Myelodysplastic syndromes (MDS) are a heterogenous collection of clonal bone marrow diseases characterized by cytopenias, abnormal karyotypes, molecular abnormalities, and dysplasia by flow cytometry and/or morphology. The progression of MDS to severe cytopenias and/or overt leukemia is associated with the accumulation of additional cytogenetic abnormalities, suggesting clonal evolution. The impact of these accumulated abnormalities on myeloid maturation and the severity of the disease is poorly understood.

METHODS

Bone marrow specimens from 16 patients with cytogenetic abnormalities were flow cytometrically sorted into three myeloid populations: progenitors, immature myeloid cells, and mature myeloid cells. Fluorescence in situ hybridization analysis was performed on each to determine the distribution of chromosomal abnormalities during myeloid maturation.

RESULTS

Our findings revealed three distinct distributions of cytogenetic abnormalities across myeloid maturation, each of which corresponded to specific cytogenetic abnormalities. Group 1 had continuous distribution across all maturational stages and contained patients with a single cytogenetic aberration associated with good-to-intermediate prognosis; Group 2 had accumulation of abnormalities in immature cells and contained patients with high-risk monosomy 7; and Group 3 had abnormalities defining the founding clone equally distributed across maturational stages while subclonal abnormalities were enriched in progenitor cells and contained patients with multiple, non-monosomy 7, abnormalities with evidence of clonal evolution.

CONCLUSIONS

Our findings demonstrate that low-risk abnormalities (e.g., del(20q) and trisomy 8) occurring in the founding clone display a markedly different disease etiology, with respect to myeloid maturation, than monosomy 7 or abnormalities acquired in subclones, which result in a disruption of myeloid cell maturation in MDS.

Determination of ERG(+), EZH2, NKX3.1, and SPINK-1 subtypes to evaluate their association with clonal origin and disease progression in multifocal prostate cancer

BACKGROUND

The prognostic relevance of prostate cancer (PCa) molecular subtypes remains controversial, given the presence of multiple foci with the possibility of different subtypes in the same patient.

AIM

To determine the clonal origin of heterogeneity in PCa and its association with disease progression, SPOP, ERG(+), EZH2, NKX3.1, and SPINK-1 subtypes were analyzed.

METHODS

A total of 103 samples from 20 PCa patients were analyzed; foci of adjacent non-tumor prostate tissue, HGPIN, GL3, GL4, GL5, and LN were examined to determine the presence of the TMPRSS2-ERG fusion and ERG, EZH2, NKX3.1, and SPINK-1 expression levels, using RT-PCR. Mutations in exons 6 and 7 of the SPOP gene were determined by sequencing. The presence of subtypes and molecular patterns were identified by combining all subtypes analyzed. To establish the clonal origin of multifocal PCa, molecular concordance between different foci of the same patient was determined. Association of these subtypes with histopathological groups and time to biochemical recurrence (BCR) was assessed.

RESULTS

No mutation was found in SPOP in any sample. The ERG(+) subtype was the most frequent. The molecular pattern containing all four PCa subtypes was only detected in 3 samples (4%), all LN, but it was the most frequent (40%) in patients. Molecular discordance was the predominant status (55%) when all analyzed molecular characteristics were considered. It was possible to find all subtypes, starting as a preneoplastic lesion, and all but one LN molecular subtype were ERG(+) and NKX3.1 subtypes. Only the expression of the NKX3.1 gene was significantly different among the histopathological groups. No association was found between BCR time in patients and molecular subtypes or molecular concordance or between clinicopathological characteristics and molecular subtypes of ERG, EZH2, and SPINK-1.

CONCLUSION

The predominance of molecular discordance in prostatic foci per patient, which reflects the multifocal origin of PCa foci, highlights the importance of analyzing multiple samples to establish the prognostic and therapeutic relevance of molecular subtypes in a patient. All the subtypes analyzed here are of early onset, starting from preneoplastic lesions. NKX3.1 gene expression is the only molecular characteristic that shows a progression pattern by sample.

Association between human leukocyte antigen and immunosuppressive treatment outcomes in Chinese patients with aplastic anemia

Background

Activated cytotoxic T cells (CTLs) recognize the auto-antigens presented on hematopoietic stem/progenitor cells (HSPCs) through class I human leukocyte antigen (HLA) molecules and play an important role in the immune pathogenesis of aplastic anemia (AA). Previous reports demonstrated that HLA was related to the disease susceptibility and response to immunosuppressive therapy (IST) in AA patients. Recent studies have indicated that specific HLA allele deletions, which helped AA patients to evade CTL-driven autoimmune responses and escape from immune surveillance, may lead to high-risk clonal evolution. Therefore, HLA genotyping has a particular predictive value for the response to IST and the risk of clonal evolution. However, there are limited studies on this topic in the Chinese population.

Methods

To explore the value of HLA genotyping in Chinese patients with AA, 95 AA patients treated with IST were retrospectively investigated.

Results

The alleles HLA-B*15:18 and HLA-C*04:01 were associated with a superior long-term response to IST (P = 0.025; P = 0.027, respectively), while the allele HLA-B*40:01 indicated an inferior result (P = 0.02). The allele HLA-A*01:01 and HLA-B*54:01 were associated with high-risk clonal evolution (P = 0.032; P = 0.01, respectively), and the former had a higher frequency in very severe AA (VSAA) patients than that in severe AA (SAA) patients (12.7% vs 0%, P = 0.02). The HLA-DQ*03:03 and HLA-DR*09:01 alleles were associated with high-risk clonal evolution and poor long-term survival in patients aged ≥40 years. Such patients may be recommended for early allogeneic hematopoietic stem cell transplantation rather than the routine IST treatment.

Conclusion

HLA genotype has crucial value in predicting the outcome of IST and long-term survival in AA patients, and thus may assist an individualized treatment strategy.

Measurable Residual Disease and Clonal Evolution in Acute Myeloid Leukemia from Diagnosis to Post-Transplant Follow-Up: The Role of Next-Generation Sequencing

It has now been ascertained that acute myeloid leukemias-as in most type of cancers-are mixtures of various subclones, evolving by acquiring additional somatic mutations over the course of the disease. The complexity of leukemia clone architecture and the phenotypic and/or genotypic drifts that can occur during treatment explain why more than 50% of patients-in hematological remission-could relapse. Moreover, the complexity and heterogeneity of clone architecture represent a hindrance for monitoring measurable residual disease, as not all minimal residual disease monitoring methods are able to detect genetic mutations arising during treatment. Unlike with chemotherapy, which imparts a relatively short duration of selective pressure on acute myeloid leukemia clonal architecture, the immunological effect related to allogeneic hematopoietic stem cell transplant is prolonged over time and must be overcome for relapse to occur. This means that not all molecular abnormalities detected after transplant always imply inevitable relapse. Therefore, transplant represents a critical setting where a measurable residual disease-based strategy, performed during post-transplant follow-up by highly sensitive methods such as next-generation sequencing, could optimize and improve treatment outcome. The purpose of our review is to provide an overview of the role of next-generation sequencing in monitoring both measurable residual disease and clonal evolution in acute myeloid leukemia patients during the entire course of the disease, with special focus on the transplant phase.

Subclonal somatic copy number alterations emerge and dominate in recurrent osteosarcoma

Multiple large-scale tumor genomic profiling efforts have been undertaken in osteosarcoma, however, little is known about the spatial and temporal intratumor heterogeneity and how it may drive treatment resistance. We performed whole-genome sequencing of 37 tumor samples from eight patients with relapsed or refractory osteosarcoma. Each patient had at least one sample from a primary site and a metastatic or relapse site. We identified subclonal copy number alterations in all but one patient. We observed that in five patients, a subclonal copy number clone from the primary tumor emerged and dominated at subsequent relapses. MYC gain/amplification was enriched in the treatment-resistant clone in 6 out of 7 patients with more than one clone. Amplifications in other potential driver genes, such as CCNE1, RAD21, VEGFA, and IGF1R, were also observed in the resistant copy number clones. Our study sheds light on intratumor heterogeneity and the potential drivers of treatment resistance in osteosarcoma.

Fatal Progression of Mutated TP53-Associated Clonal Hematopoiesis following Anti-CD19 CAR-T Cell Therapy

We present the case of a 64-year-old man diagnosed with large B-cell lymphoma who relapsed twice after standard-of-care therapy. Due to persisting cytopenia, Next generation sequencing analysis was performed, revealing a small TP53-mutated clone. As a third-line therapy, the patient was treated with CAR-T cells, which resulted in complete remission. However, this treatment also led to the expansion of the TP53-mutated clone and therapy-related myelodysplasia with a complex aberrant karyotype. This case may serve as a paradigmatic example of clonal hematopoietic progression in a patient undergoing CAR-T cell therapy, especially in the context of a TP53-mutated clone.

Refining the Treatment of Pancreatic Cancer From Big Data to Improved Individual Survival

Pancreatic cancer is one of the most lethal cancers worldwide, most notably in Europe and North America. Great strides have been made in combining the most effective conventional therapies to improve survival at least in the short and medium term. The start of treatment can only be made once a diagnosis is made, which at this point, the tumor volume is already very high in the primary cancer and systemically. If caught at the earliest opportunity (in circa 20% patients) surgical resection of the primary followed by combination chemotherapy can achieve 5-year overall survival rates of 30%-50%. A delay in detection of even a few months after symptom onset will result in the tumor having only borderline resectabilty (in 20%-30% of patients), in which case the best survival is achieved by using short-course chemotherapy before tumor resection as well as adjuvant chemotherapy. Once metastases become visible (in 40%-60% of patients), cure is not possible, palliative cytotoxics only being able to prolong life by few months. Even in apparently successful therapy in resected and borderline resectable patients, the recurrence rate is very high. Considerable efforts to understand the nature of pancreatic cancer through large-scale genomics, transcriptomics, and digital profiling, combined with functional preclinical models, using genetically engineered mouse models and patient derived organoids, have identified the critical role of the tumor microenvironment in determining the nature of chemo- and immuno-resistance. This functional understanding has powered fresh and exciting approaches for the treatment of this cancer.

Insights into IGH clonal evolution in BCP-ALL: frequency, mechanisms, associations, and diagnostic implications

Introduction

The malignant transformation leading to a maturation arrest in B-cell precursor acute lymphoblastic leukemia (BCP-ALL) occurs early in B-cell development, in a pro-B or pre-B cell, when somatic recombination of variable (V), diversity (D), and joining (J) segment immunoglobulin (IG) genes and the B-cell rescue mechanism of V_H replacement might be ongoing or fully active, driving clonal evolution. In this study of newly diagnosed BCP-ALL, we sought to understand the mechanistic details of oligoclonal composition of the leukemia at diagnosis, clonal evolution during follow-up, and clonal distribution in different hematopoietic compartments.

Methods

Utilizing high-throughput sequencing assays and bespoke bioinformatics we identified BCP-ALL-derived clonally-related IGH sequences by their shared 'DNJ-stem'.

Results

We introduce the concept of 'marker DNJ-stem' to cover the entirety of, even lowly abundant, clonally-related family members. In a cohort of 280 adult patients with BCP-ALL, IGH clonal evolution at diagnosis was identified in one-third of patients. The phenomenon was linked to contemporaneous recombinant and editing activity driven by aberrant ongoing D_H/V_H-DJ_H recombination and V_H replacement, and we share insights and examples for both. Furthermore, in a subset of 167 patients with molecular subtype allocation, high prevalence and high degree of clonal evolution driven by ongoing D_H/V_H-DJ_H recombination were associated with the presence of KMT2A gene rearrangements, while V_H replacements occurred more frequently in Ph-like and DUX4 BCP-ALL. Analysis of 46 matched diagnostic bone marrow and peripheral blood samples showed a comparable clonal and clonotypic distribution in both hematopoietic compartments, but the clonotypic composition markedly changed in longitudinal follow-up analysis in select cases. Thus, finally, we present cases where the specific dynamics of clonal evolution have implications for both the initial marker identification and the MRD monitoring in follow-up samples.

Discussion

Consequently, we suggest to follow the marker DNJ-stem (capturing all family members) rather than specific clonotypes as the MRD target, as well as to follow both VDJ_H and DJ_H family members since their respective kinetics are not always parallel. Our study further highlights the intricacy, importance, and present and future challenges of IGH clonal evolution in BCP-ALL.

Evolutionary history of transformation from chronic lymphocytic leukemia to Richter syndrome

Richter syndrome (RS) arising from chronic lymphocytic leukemia (CLL) exemplifies an aggressive malignancy that develops from an indolent neoplasm. To decipher the genetics underlying this transformation, we computationally deconvoluted admixtures of CLL and RS cells from 52 patients with RS, evaluating paired CLL-RS whole-exome sequencing data. We discovered RS-specific somatic driver mutations (including IRF2BP2, SRSF1, B2M, DNMT3A and CCND3), recurrent copy-number alterations beyond del(9p21)(CDKN2A/B), whole-genome duplication and chromothripsis, which were confirmed in 45 independent RS cases and in an external set of RS whole genomes. Through unsupervised clustering, clonally related RS was largely distinct from diffuse large B cell lymphoma. We distinguished pathways that were dysregulated in RS versus CLL, and detected clonal evolution of transformation at single-cell resolution, identifying intermediate cell states. Our study defines distinct molecular subtypes of RS and highlights cell-free DNA analysis as a potential tool for early diagnosis and monitoring.

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.

Current concepts of anti-EGFR targeting in metastatic colorectal cancer

Anti-EGFR targeting is one of the key strategies in the treatment of metastatic colorectal cancer (mCRC). For almost two decades oncologists have struggled to implement EGFR antibodies in the mCRC continuum of care. Both sidedness and RAS mutational status rank high among the predictive factors for the clinical efficacy of EGFR inhibitors. A prospective phase III trial has recently confirmed that anti-EGFR targeting confers an overall survival benefit only in left sided RAS-wildtype tumors when given in first line. It is a matter of discussion if more clinical benefit can be reached by considering putative primary resistance mechanisms (e.g., HER2, BRAF, PIK3CA, etc.) at this early stage of treatment. The value of this procedure in daily routine clinical utility has not yet been clearly delineated. Re-exposure to EGFR antibodies becomes increasingly crucial in the disease journey of mCRC. Yet re- induction or re-challenge strategies have been problematic as they relied on mathematical models that described the timely decay of EGFR antibody resistant clones. The advent of liquid biopsy and the implementation of more accurate next-generation sequencing (NGS) based high throughput methods allows for tracing of EGFR resistant clones in real time. These displays the spatiotemporal heterogeneity of metastatic disease compared to the former standard radiographic assessment and re-biopsy. These techniques may move EGFR inhibition in mCRC into the area of precision medicine in order to apply EGFR antibodies with the increase or decrease of EGFR resistant clones. This review critically discusses established concepts of tackling the EGFR pathway in mCRC and provides insight into the growing field of liquid biopsy guided personalized approaches of EGFR inhibition in mCRC.

Triple-negative Thrombocythemia and Subsequent Acute Lymphoblastic Leukemia with Additional Somatic Mutations

Triple-negative essential thrombocythemia (ET) is a condition in which mutations in JAK2, CALR and MPL are all negative. Transformation to acute myeloid leukemia may occur during the course of ET, while B-acute lymphoblastic leukemia (B-ALL) is rare. We experienced a case diagnosed as B-ALL during the course of triple-negative ET. Notably, cytoreduction was required for the excessive increase in blood cells during the bone marrow recovery period after chemotherapies. Whole-exome sequencing identified 17 somatic mutations: 9 were identified in both ET and B-ALL samples, while 8 were specific to B-ALL, suggesting that these 8 might have caused the transformation.

Targeting stem cells in myelodysplastic syndromes and acute myeloid leukemia

The genetic architecture of cancer has been delineated through advances in high-throughput next-generation sequencing, where the sequential acquisition of recurrent driver mutations initially targeted towards normal cells ultimately leads to malignant transformation. Myelodysplastic syndromes (MDS) and acute myeloid leukemia (AML) are hematologic malignancies frequently initiated by mutations in the normal hematopoietic stem cell compartment leading to the establishment of leukemic stem cells. Although the genetic characterization of MDS and AML has led to identification of new therapeutic targets and development of new promising therapeutic strategies, disease progression, relapse, and treatment-related mortality remain a major challenge in MDS and AML. The selective persistence of rare leukemic stem cells following therapy-induced remission implies unique resistance mechanisms of leukemic stem cells towards conventional therapeutic strategies and that leukemic stem cells represent the cellular origin of relapse. Therefore, targeted surveillance of leukemic stem cells following therapy should, in the future, allow better prediction of relapse and disease progression, but is currently challenged by our restricted ability to distinguish leukemic stem cells from other leukemic cells and residual normal cells. To advance current and new clinical strategies for the treatment of MDS and AML, there is a need to improve our understanding and characterization of MDS and AML stem cells at the cellular, molecular, and genetic levels. Such work has already led to the identification of promising new candidate leukemic stem cell molecular targets that can now be exploited in preclinical and clinical therapeutic strategies, towards more efficient and specific elimination of leukemic stem cells.

Heterogeneous genetic and non-genetic mechanisms contribute to response and resistance to azacitidine monotherapy

Acute myeloid leukaemia is prevalent in older patients that are often ineligible for intensive chemotherapy and treatment options remain limited with azacitidine being at the forefront. Azacitidine has been used in the clinic for decades, however, we still lack a complete understanding of the mechanisms by which the drug exerts its anti-tumour effect. To gain insight into the mechanism of action, we defined the mutational profile of sequential samples of patients treated with azacitidine. We did not identify any mutations that could predict response and observed lack of a uniform pattern of clonal evolution. Focusing on responders, at remission, we observed three types of response: (1) an almost complete elimination of mutations (33%), (2) no change (17%), and (3) change with no discernible pattern (50%). Heterogeneous patterns were also observed at relapse, with no clonal evolution between remission and relapse in some patients. Lack of clonal evolution suggests that non-genetic mechanisms might be involved. Towards understanding such mechanisms, we investigated the immune microenvironment in a number of patients and we observed lack of a uniform response following therapy. We identified a higher frequency of cytotoxic T cells in responders and higher frequency of naïve helper T cells in non-responders.

Tracking Clonal Evolution of Multiple Myeloma Using Targeted Next-Generation DNA Sequencing

Clonal evolution drives treatment failure in multiple myeloma (MM). Here, we used a custom 372-gene panel to track genetic changes occurring during MM progression at different stages of the disease. A tumor-only targeted next-generation DNA sequencing was performed on 69 samples sequentially collected from 30 MM patients. The MAPK/ERK pathway was mostly affected with KRAS mutated in 47% of patients. Acquisition and loss of mutations were observed in 63% and 37% of patients, respectively. Four different patterns of mutation evolution were found: branching-, mutation acquisition-, mutation loss- and a stable mutational pathway. Better response to anti-myeloma therapy was more frequently observed in patients who followed the mutation loss-compared to the mutation acquisition pathway. More than two-thirds of patients had druggable genes mutated (including cases of heavily pre-treated disease). Only 7% of patients had a stable copy number variants profile. Consequently, a redistribution in stages according to R-ISS between the first and paired samples (R-ISS″) was seen. The higher the R-ISS″, the higher the risk of MM progression and death. We provided new insights into the genetics of MM evolution, especially in heavily pre-treated patients. Additionally, we confirmed that redefining R-ISS at MM relapse is of high clinical value.

Geographic Origin and Vertical Transmission of Leishmania infantum Parasites in Hunting Hounds, United States

Vertical transmission of leishmaniasis is common but is difficult to study against the background of pervasive vector transmission. We present genomic data from dogs in the United States infected with Leishmania infantum parasites; these infections have persisted in the apparent absence of vector transmission. We demonstrate that these parasites were introduced from the Old World separately and more recently than L. infantum from South America. The parasite population shows unusual genetics consistent with a lack of meiosis: a high level of heterozygous sites shared across all isolates and no decrease in linkage with genomic distance between variants. Our data confirm that this parasite population has been evolving with little or no sexual reproduction. This demonstration of vertical transmission has profound implications for the population genetics of Leishmania parasites. When investigating transmission in complex natural settings, considering vertical transmission alongside vector transmission is vital.