Search

Lung cancer is the leading cause of cancer-associated mortality worldwide 1 . Here we analysed 1,644 tumour regions sampled at surgery or during follow-up from the first 421 patients with non-small cell lung cancer prospectively enrolled into the TRACERx study. This project aims to decipher lung cancer evolution and address the primary study endpoint: determining the relationship between intratumour heterogeneity and clinical outcome. In lung adenocarcinoma, mutations in 22 out of 40 common cancer genes were under significant subclonal selection, including classical tumour initiators such as TP53 and KRAS. We defined evolutionary dependencies between drivers, mutational processes and whole genome doubling (WGD) events. Despite patients having a history of smoking, 8% of lung adenocarcinomas lacked evidence of tobacco-induced mutagenesis. These tumours also had similar detection rates for EGFR mutations and for RET, ROS1, ALK and MET oncogenic isoforms compared with tumours in never-smokers, which suggests that they have a similar aetiology and pathogenesis. Large subclonal expansions were associated with positive subclonal selection. Patients with tumours harbouring recent subclonal expansions, on the terminus of a phylogenetic branch, had significantly shorter disease-free survival. Subclonal WGD was detected in 19% of tumours, and 10% of tumours harboured multiple subclonal WGDs in parallel. Subclonal, but not truncal, WGD was associated with shorter disease-free survival. Copy number heterogeneity was associated with extrathoracic relapse within 1 year after surgery. These data demonstrate the importance of clonal expansion, WGD and copy number instability in determining the timing and patterns of relapse in non-small cell lung cancer and provide a comprehensive clinical cancer evolutionary data resource., (© 2023. The Author(s).)

Intra-tumor heterogeneity (ITH) is a pan-cancer predictor of survival, with high ITH being correlated to a dismal prognosis. The level of ITH is, hence, a clinically relevant characteristic of a malignancy. ITH of karyotypes is driven by chromosomal instability (CIN). However, not all new karyotypes generated by CIN are viable or competitive, which limits the amount of ITH. Here, we review the cellular processes and ecological properties that determine karyotype ITH. We propose a framework to understand karyotype ITH, in which cells with new karyotypes emerge through CIN, are selected by cell intrinsic and cell extrinsic selective pressures, and propagate through a cancer in competition with other malignant cells. We further discuss how CIN modulates the cell phenotype and immune microenvironment, and the implications this has for the subsequent selection of karyotypes. Together, we aim to provide a comprehensive overview of the biological processes that shape the level of karyotype heterogeneity.

Background: Endocrine therapy forms the backbone of adjuvant treatment for oestrogen-receptor-positive (ER+) breast cancer. However, it remains unclear whether adjuvant treatment improves survival rates in low-risk patients. Low intra-tumour heterogeneity (ITH) has been shown to confer low risk for recurrent disease. Here, it is studied if chromosomal copy-number ITH (CNH) can identify low-risk ER+, lymph-node-negative breast cancer patients who do not benefit from adjuvant endocrine therapy., Methods: Lymph-node-negative ER+ patients from the observational METABRIC dataset were retrospectively analysed (n = 708). CNH was determined from a single bulk copy-number measurement for each patient. Survival rates were compared between patients that did or did not receive adjuvant endocrine therapy for CNH-low, middle and high groups with Cox proportional-hazards models, using propensity-score weights to correct for confounders., Results: Adjuvant endocrine therapy improved the relapse-free survival (RFS) for CNH-high patients treatment (HR = 0.55), but not for CNH-low patients treatment (HR = 0.88). For CNH-low patients adjuvant endocrine therapy was associated with impaired OS (HR = 1.62)., Conclusions: This retrospective study of lymph-node-negative, ER+ breast cancer finds that patients identified as low risk using CNH do not benefit from adjuvant endocrine therapy., (© 2022. The Author(s), under exclusive licence to Springer Nature Limited.)

A significant proportion of colorectal cancer (CRC) patients develop peritoneal metastases (PM) in the course of their disease. PMs are associated with a poor quality of life, significant morbidity and dismal disease outcome. To improve care for this patient group, a better understanding of the molecular characteristics of CRC-PM is required. Here we present a comprehensive molecular characterization of a cohort of 52 patients. This reveals that CRC-PM represent a distinct CRC molecular subtype, CMS4, but can be further divided in three separate categories, each presenting with unique features. We uncover that the CMS4-associated structural protein Moesin plays a key role in peritoneal dissemination. Finally, we define specific evolutionary features of CRC-PM which indicate that polyclonal metastatic seeding underlies these lesions. Together our results suggest that CRC-PM should be perceived as a distinct disease entity., (© 2022. The Author(s).)

Overall, the prognosis of patients suffering from stage II colon cancer is relatively favorable. However, a proportion of patients develop a recurrence following surgery. Clinical and histopathological properties that identify high-risk patients are of limited value and better biomarkers are urgently required. In a recent issue of The Journal of Pathology, Lahoz et al proposed that copy-number-based biomarkers could be employed for patient stratification. The authors studied copy-number alterations (CNAs) at the genomic scale by measuring the total CNA load (the aberrant genome fraction), and at a smaller scale by identifying common arm- or cytoband-level alterations. Both the overall CNA load and specific chromosomal regions were associated with an increased risk of recurrence. Most interestingly, it was demonstrated that copy-number intratumor heterogeneity, as defined by subclonal CNAs, is associated with poor disease outcome. This study demonstrates that structural genomic aberrations are promising biomarkers for patient stratification in early colon cancer. © 2022 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland., (© 2022 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.)

The tissue dynamics that govern maintenance and regeneration of the pancreas remain largely unknown. In particular, the presence and nature of a cellular hierarchy remains a topic of debate. Previous lineage tracing strategies in the pancreas relied on specific marker genes for clonal labeling, which left other populations untested and failed to account for potential widespread phenotypical plasticity. Here we employed a tracing system that depends on replication-induced clonal marks. We found that, in homeostasis, steady acinar replacement events characterize tissue dynamics, to which all acinar cells have an equal ability to contribute. Similarly, regeneration following pancreatitis was best characterized by an acinar self-replication model because no evidence of a cellular hierarchy was detected. In particular, rapid regeneration in the pancreas was found to be driven by an accelerated rate of acinar fission-like events. These results provide a comprehensive and quantitative model of cell dynamics in the exocrine pancreas., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2021 The Authors. Published by Elsevier Inc. All rights reserved.)

Effective treatments for pancreatic ductal adenocarcinoma (PDAC) are lacking, and targeted agents have demonstrated limited efficacy. It has been speculated that a rare population of cancer stem cells (CSCs) drives growth, therapy resistance, and rapid metastatic progression in PDAC. These CSCs demonstrate high clonogenicity in vitro and tumorigenic potential in vivo. However, their relevance in established PDAC tissue has not been determined. Here, we use marker-independent stochastic clonal labeling, combined with quantitative modeling of tumor expansion, to uncover PDAC tissue growth dynamics. We find that in contrast to the CSC model, all PDAC cells display clonogenic potential in situ. Furthermore, the proximity to activated cancer-associated fibroblasts determines tumor cell clonogenicity. This means that the microenvironment is dominant in defining the clonogenic activity of PDAC cells. Indeed, manipulating the stroma by Hedgehog pathway inhibition alters the tumor growth mode, revealing that tumor-stroma crosstalk shapes tumor growth dynamics and clonal architecture., Competing Interests: Declaration of interests M.F.B. has received research funding from Celgene and acted as a consultant to Servier. L.V. has received speaker and consultancy fees from Genentech. The remaining authors declare no competing interests in the context of this publication., (Copyright © 2021 The Author(s). Published by Elsevier Inc. All rights reserved.)

We recently introduced a method to derive intra-tumor heterogeneity (ITH) from a single copy number measurement. This method stratifies patients for survival and could potentially help to identify low and high-risk patients with clinical relevance., Competing Interests: E.v.D., L.V. and D.M.M. are listed as inventors in a pending patent application (NL82151) filed by Oncode Institute on behalf of the Academisch Medisch Centrum, covering the content of the paper. T.v.d.B. declares no conflict of interest., (© 2021 The Author(s). Published with license by Taylor & Francis Group, LLC.)

Survival rates of cancer patients vary widely within and between malignancies. While genetic aberrations are at the root of all cancers, individual genomic features cannot explain these distinct disease outcomes. In contrast, intra-tumour heterogeneity (ITH) has the potential to elucidate pan-cancer survival rates and the biology that drives cancer prognosis. Unfortunately, a comprehensive and effective framework to measure ITH across cancers is missing. Here, we introduce a scalable measure of chromosomal copy number heterogeneity (CNH) that predicts patient survival across cancers. We show that the level of ITH can be derived from a single-sample copy number profile. Using gene-expression data and live cell imaging we demonstrate that ongoing chromosomal instability underlies the observed heterogeneity. Analysing 11,534 primary cancer samples from 37 different malignancies, we find that copy number heterogeneity can be accurately deduced and predicts cancer survival across tissues of origin and stages of disease. Our results provide a unifying molecular explanation for the different survival rates observed between cancer types.

Importance: Quality improvement programs for colorectal cancer surgery have been introduced with benchmarking based on quality indicators, such as mortality. Detailed (pre)operative characteristics may offer relevant information for proper case-mix correction., Objective: To investigate the added value of machine learning to predict quality indicators for colorectal cancer surgery and identify previously unrecognized predictors of 30-day mortality based on a large, nationwide colorectal cancer registry that collected extensive data on comorbidities., Design, Setting, and Participants: All patients who underwent resection for primary colorectal cancer registered in the Dutch ColoRectal Audit between January 1, 2011, and December 31, 2016, were included. Multiple machine learning models (multivariable logistic regression, elastic net regression, support vector machine, random forest, and gradient boosting) were made to predict quality indicators. Model performance was compared with conventionally used scores. Risk factors were identified by logistic regression analyses and Shapley additive explanations (ie, SHAP values). Statistical analysis was performed between March 1 and September 30, 2020., Main Outcomes and Measures: The primary outcome of this cohort study was 30-day mortality. Prediction models were trained on a training set by performing 5-fold cross-validation, and outcomes were measured by the area under the receiver operating characteristic curve on the test set. Machine learning was further used to identify risk factors, measured by odds ratios and SHAP values., Results: This cohort study included 62 501 records, most patients were male (35 116 [56.2%]), were aged 61 to 80 years (41 560 [66.5%]), and had an American Society of Anesthesiology score of II (35 679 [57.1%]). A 30-day mortality rate of 2.7% (n = 1693) was found. The area under the curve of the best machine learning model for 30-day mortality (0.82; 95% CI, 0.79-0.85) was significantly higher than the American Society of Anesthesiology score (0.74; 95% CI, 0.71-0.77; P < .001), Charlson Comorbidity Index (0.66; 95% CI, 0.63-0.70; P < .001), and preoperative score to predict postoperative mortality (0.73; 95% CI, 0.70-0.77; P < .001). Hypertension, myocardial infarction, chronic obstructive pulmonary disease, and asthma were comorbidities with a high risk for increased mortality. Machine learning identified specific risk factors for a complicated course, intensive care unit admission, prolonged hospital stay, and readmission. Laparoscopic surgery was associated with a decreased risk for all adverse outcomes., Conclusions and Relevance: This study found that machine learning methods outperformed conventional scores to predict 30-day mortality after colorectal cancer surgery, identified specific patient groups at risk for adverse outcomes, and provided directions to optimize benchmarking in clinical audits.

Breast implant-associated anaplastic large cell lymphoma (BIA-ALCL) is a very rare type of T-cell lymphoma that is uniquely caused by a single environmental stimulus. Here, we present a comprehensive genetic analysis of a relatively large series of BIA-ALCL (n = 29), for which genome-wide chromosomal copy number aberrations (CNAs) and mutational profiles for a subset (n = 7) were determined. For comparison, CNAs for anaplastic lymphoma kinase (ALK)- nodal anaplastic large cell lymphomas (ALCLs; n = 24) were obtained. CNAs were detected in 94% of BIA-ALCLs, with losses at chromosome 20q13.13 in 66% of the samples. Loss of 20q13.13 is characteristic of BIA-ALCL compared with other classes of ALCL, such as primary cutaneous ALCL and systemic type ALK+ and ALK- ALCL. Mutational patterns confirm that the interleukin-6-JAK1-STAT3 pathway is deregulated. Although this is commonly observed across various types of T-cell lymphomas, the extent of deregulation is significantly higher in BIA-ALCL, as indicated by phosphorylated STAT3 immunohistochemistry. The characteristic loss of chromosome 20 in BIA-ALCL provides further justification to recognize BIA-ALCL as a separate disease entity. Moreover, CNA analysis may serve as a parameter for future diagnostic assays for women with breast implants to distinguish seroma caused by BIA-ALCL from other causes of seroma accumulation, such as infection or trauma., (© 2020 by The American Society of Hematology.)

Intraflagellar transport (IFT) in C. elegans chemosensory cilia is an example of functional coordination and cooperation of two motor proteins with distinct motility properties operating together in large groups to transport cargoes: a fast and processive homodimeric kinesin-2, OSM-3, and a slow and less processive heterotrimeric kinesin-2, kinesin-II. To study the mechanism of the collective dynamics of kinesin-II of C. elegans cilia in an in vitro system, we used Total Internal Reflection Fluorescence microscopy to image the motility of truncated, heterodimeric kinesin-II constructs at high motor densities. Using an analysis technique based on correlation of the fluorescence intensities, we extracted quantitative motor parameters, such as motor density, velocity and average run length, from the image. Our experiments and analyses show that kinesin-II motility parameters are far less affected by (self) crowding than OSM-3. Our observations are supported by numerical calculations based on the TASEP-LK model (Totally Asymmetric Simple Exclusion Process-Langmuir Kinetics). From a comparison of data and modelling of OSM-3 and kinesin-II, a general picture emerges of the collective dynamics of the kinesin motors driving IFT in C. elegans chemosensory cilia and the way the motors deal with crowding., Competing Interests: The authors have declared that no competing interests exist.

Lineage tracing is a powerful tool that can be used to uncover cell fates. Here, we describe a novel method for the quantitative analysis of clonal dynamics in grafted cancer tissues. The protocol involves the preparation and validation of cells for lineage tracing, establishment of grafts and label induction, analysis of clone-size distribution and fitting of the experimental data to a mathematical tumor growth model. In contrast to other lineage-tracing strategies, the method described here assesses stem cell functionality and infers tumor expansion dynamics independently of molecular markers such as putative cancer stem cell (CSC)-specific genes. The experimental system and analytical framework presented can be used to quantify clonal advantages that specific mutations provide, in both the absence and presence of (targeted) therapeutic agents. This protocol typically takes ~20 weeks to complete from cell line selection to inference of growth dynamics, depending on the grafted cancer growth rate.

Cancer evolution is predominantly studied by focusing on differences in the genetic characteristics of malignant cells within tumors. However, the spatiotemporal dynamics of clonal outgrowth that underlie evolutionary trajectories remain largely unresolved. Here, we sought to unravel the clonal dynamics of colorectal cancer (CRC) expansion in space and time by using a color-based clonal tracing method. This method involves lentiviral red-green-blue (RGB) marking of cell populations, which enabled us to track individual cells and their clonal outgrowth during tumor initiation and growth in a xenograft model. We found that clonal expansion largely depends on the location of a clone, as small clones reside in the center and large clones mostly drive tumor growth at the border. These dynamics are recapitulated in a computational model, which confirms that the clone position within a tumor rather than cell-intrinsic features, is crucial for clonal outgrowth. We also found that no significant clonal loss occurs during tumor growth and clonal dispersal is limited in most models. Our results imply that, in addition to molecular features of clones such as (epi-)genetic differences between cells, clone location and the geometry of tumor growth are crucial for clonal expansion. Our findings suggest that either microenvironmental signals on the tumor border or differences in physical properties within the tumor, are major contributors to explain heterogeneous clonal expansion. Thus, this study provides further insights into the dynamics of solid tumor growth and progression, as well as the origins of tumor cell heterogeneity in a relevant model system., Competing Interests: The authors declare no conflict of interest., (Copyright © 2019 the Author(s). Published by PNAS.)

By using marker-free lineage tracing in combination with quantitative analysis, we recently revealed cancer stem cell functionality in established human colon cancer is not intrinsically defined, but fully spatiotemporally regulated.

Solid malignancies have been speculated to depend on cancer stem cells (CSCs) for expansion and relapse after therapy. Here we report on quantitative analyses of lineage tracing data from primary colon cancer xenograft tissue to assess CSC functionality in a human solid malignancy. The temporally obtained clone size distribution data support a model in which stem cell function in established cancers is not intrinsically, but is entirely spatiotemporally orchestrated. Functional stem cells that drive tumour expansion predominantly reside at the tumour edge, close to cancer-associated fibroblasts. Hence, stem cell properties change in time depending on the cell location. Furthermore, although chemotherapy enriches for cells with a CSC phenotype, in this context functional stem cell properties are also fully defined by the microenvironment. To conclude, we identified osteopontin as a key cancer-associated fibroblast-produced factor that drives in situ clonogenicity in colon cancer.

Using the totally asymmetric simple-exclusion-process and mean-field transport theory, we investigate the transport in closed random networks with simple crossing topology-two incoming, two outgoing segments, as a model for molecular motor motion along biopolymer networks. Inspired by in vitro observation of molecular motor motion, we model the motor behavior at the intersections by introducing different exit rates for the two outgoing segments. Our simulations of this simple network reveal surprisingly rich behavior of the transport current with respect to the global density and exit rate ratio. For asymmetric exit rates, we find a broad current plateau at intermediate motor densities resulting from the competition of two subnetwork populations. This current plateau leads to stabilization of transport properties within such networks.

We demonstrate that application of an increasing shear field on a glass leads to an intriguing dynamic first-order transition in analogy with equilibrium transitions. By following the particle dynamics as a function of the driving field in a colloidal glass, we identify a critical shear rate upon which the diffusion time scale of the glass exhibits a sudden discontinuity. Using a new dynamic order parameter, we show that this discontinuity is analogous to a first-order transition, in which the applied stress acts as the conjugate field on the system's dynamic evolution. These results offer new perspectives to comprehend the generic shear-banding instability of a wide range of amorphous materials.

We study condensation in one-dimensional transport models with a kinetic constraint. The kinetic constraint results in clustering of immobile vehicles; these clusters can grow to macroscopic condensates, indicating the onset of dynamic phase separation between free-flowing and arrested traffic. We investigate analytically the conditions under which this occurs and derive a necessary and sufficient criterion for phase separation. This criterion is applied to the well-known Nagel-Schreckenberg model of traffic flow to analytically investigate the existence of dynamic condensates. We find that true condensates occur only when acceleration out of jammed traffic happens in a single time step, in the limit of strong overbraking. Our predictions are further verified with simulation results on the growth of arrested clusters. These results provide analytic understanding of dynamic arrest and dynamic phase separation in one-dimensional traffic and transport models.

Dynamic arrest is a general phenomenon across a wide range of dynamic systems including glasses, traffic flow, and dynamics in cells, but the universality of dynamic arrest phenomena remains unclear. We connect the emergence of traffic jams in a simple traffic flow model directly to the dynamic slowing down in kinetically constrained models for glasses. In kinetically constrained models, the formation of glass becomes a true (singular) phase transition in the limit T→0. Similarly, using the Nagel-Schreckenberg model to simulate traffic flow, we show that the emergence of jammed traffic acquires the signature of a sharp transition in the deterministic limit p→1, corresponding to overcautious driving. We identify a true dynamic critical point marking the onset of coexistence between free flowing and jammed traffic, and demonstrate its analogy to the kinetically constrained glass models. We find diverging correlations analogous to those at a critical point of thermodynamic phase transitions.

The 275GHz electron-paramagnetic-resonance spectrometer we reported on in 2004 has been equipped with a new probe head, which contains a cavity especially designed for operation in continuous-wave mode. The sensitivity and signal stability that is achieved with this new probe head is illustrated with 275GHz continuous-wave spectra of a 1mM frozen solution of the complex Fe(III)-ethylenediamine tetra-acetic acid and of 10mM frozen solutions of the protein rubredoxin, which contains Fe(3+) in its active site, from three different organisms. The high quality of the spectra of the rubredoxins allows the determination of the zero-field-splitting parameters with an accuracy of 0.5GHz. The success of our approach results partially from the enhanced absolute sensitivity, which can be reached using a single-mode cavity. At least as important is the signal stability that we were able to achieve with the new probe head., (Copyright © 2011 Elsevier Inc. All rights reserved.)