Characterizing genetic intra-tumor heterogeneity across 2,658 human cancer genomes.

Intra-tumor heterogeneity (ITH) is a mechanism of therapeutic resistance and therefore an important clinical challenge. However, the extent, origin, and drivers of ITH across cancer types are poorly understood. To address this, we extensively characterize ITH across whole-genome sequences of 2,658 cancer samples spanning 38 cancer types. Nearly all informative samples (95.1%) contain evidence of distinct subclonal expansions with frequent branching relationships between subclones. We observe positive selection of subclonal driver mutations across most cancer types and identify cancer type-specific subclonal patterns of driver gene mutations, fusions, structural variants, and copy number alterations as well as dynamic changes in mutational processes between subclonal expansions. Our results underline the importance of ITH and its drivers in tumor evolution and provide a pan-cancer resource of comprehensively annotated subclonal events from whole-genome sequencing data.

The evolutionary history of 2,658 cancers.

Cancer develops through a process of somatic evolution1,2. Sequencing data from a single biopsy represent a snapshot of this process that can reveal the timing of specific genomic aberrations and the changing influence of mutational processes3. Here, by whole-genome sequencing analysis of 2,658 cancers as part of the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA)4, we reconstruct the life history and evolution of mutational processes and driver mutation sequences of 38 types of cancer. Early oncogenesis is characterized by mutations in a constrained set of driver genes, and specific copy number gains, such as trisomy 7 in glioblastoma and isochromosome 17q in medulloblastoma. The mutational spectrum changes significantly throughout tumour evolution in 40% of samples. A nearly fourfold diversification of driver genes and increased genomic instability are features of later stages. Copy number alterations often occur in mitotic crises, and lead to simultaneous gains of chromosomal segments. Timing analyses suggest that driver mutations often precede diagnosis by many years, if not decades. Together, these results determine the evolutionary trajectories of cancer, and highlight opportunities for early cancer detection.

Growth dynamics in naturally progressing chronic lymphocytic leukaemia.

How the genomic features of a patient's cancer relate to individual disease kinetics remains poorly understood. Here we used the indolent growth dynamics of chronic lymphocytic leukaemia (CLL) to analyse the growth rates and corresponding genomic patterns of leukaemia cells from 107 patients with CLL, spanning decades-long disease courses. We found that CLL commonly demonstrates not only exponential expansion but also logistic growth, which is sigmoidal and reaches a certain steady-state level. Each growth pattern was associated with marked differences in genetic composition, the pace of disease progression and the extent of clonal evolution. In a subset of patients, whose serial samples underwent next-generation sequencing, we found that dynamic changes in the disease course of CLL were shaped by the genetic events that were already present in the early slow-growing stages. Finally, by analysing the growth rates of subclones compared with their parental clones, we quantified the growth advantage conferred by putative CLL drivers in vivo.

Insights into Chemically Fueled Supramolecular Polymers.

Supramolecular polymerization can be controlled in space and time by chemical fuels. A nonassembled monomer is activated by the fuel and subsequently self-assembles into a polymer. Deactivation of the molecule either in solution or inside the polymer leads to disassembly. Whereas biology has already mastered this approach, fully artificial examples have only appeared in the past decade. Here, we map the available literature examples into four distinct regimes depending on their activation/deactivation rates and the equivalents of deactivating fuel. We present increasingly complex mathematical models, first considering only the chemical activation/deactivation rates (i.e., transient activation) and later including the full details of the isodesmic or cooperative supramolecular processes (i.e., transient self-assembly). We finish by showing that sustained oscillations are possible in chemically fueled cooperative supramolecular polymerization and provide mechanistic insights. We hope our models encourage the quantification of activation, deactivation, assembly, and disassembly kinetics in future studies.

Activation of the MAPK pathway mediates resistance to PI3K inhibitors in chronic lymphocytic leukemia.

Inhibitors of Bruton tyrosine kinase (BTK) and phosphatidylinositol 3-kinase δ (PI3Kδ) that target the B-cell receptor (BCR) signaling pathway have revolutionized the treatment of chronic lymphocytic leukemia (CLL). Mutations associated with resistance to BTK inhibitors have been identified, but limited data are available on mechanisms of resistance to PI3Kδ inhibitors. Here we present findings from longitudinal whole-exome sequencing of cells from patients with multiply relapsed CLL (N = 28) enrolled in trials of PI3K inhibitors. The nonresponder subgroup was characterized by baseline activating mutations in MAP2K1, BRAF, and KRAS genes in 60% of patients. PI3Kδ inhibition failed to inhibit ERK phosphorylation (pERK) in nonresponder CLL cells with and without mutations, whereas treatment with a MEK inhibitor rescued ERK inhibition. Overexpression of MAP2K1 mutants in vitro led to increased basal and inducible pERK and resistance to idelalisib. These data demonstrate that MAPK/ERK activation plays a key role in resistance to PI3Kδ inhibitors in CLL and provide a rationale for therapy with a combination of PI3Kδ and ERK inhibitors.

Magneto-capillary particle dynamics at curved interfaces: inference and criticism of dynamical models.

Time-varying fields drive the motion of magnetic particles adsorbed on liquid drops due to interfacial constraints that couple magnetic torques to capillary forces. Such magneto-capillary particle dynamics and the associated fluid flows are potentially useful for propelling drop motion, mixing drop contents, and enhancing mass transfer between phases. The design of such functions benefits from the development and validation of predictive models. Here, we apply methods of Bayesian data analysis to identify and validate a dynamical model that accurately predicts the field-driven motion of a magnetic particle adsorbed at the interface of a spherical droplet. Building on previous work, we consider candidate models that describe particle tilting at the interface, field-dependent contributions to the magnetic moment, gravitational forces, and their combinations. The analysis of each candidate is informed by particle tracking data for a magnetic Janus sphere moving in a precessing field at different frequencies and angles. We infer the uncertain parameters of each model, criticize their ability to describe and predict experimental data, and select the most probable candidate, which accounts for gravitational forces and the tilting of the Janus sphere at the interface. We show how this favored model can predict complex particle trajectories with micron-level accuracy across the range of driving fields considered. We discuss how knowledge of this "best" model can be used to design experiments that inform accurate parameter estimates or achieve desired particle trajectories.

Programmable topotaxis of magnetic rollers in time-varying fields.

We describe how spatially uniform, time-periodic magnetic fields can be designed to power and direct the migration of ferromagnetic spheres up (or down) local gradients in the topography of a solid substrate. Our results are based on a dynamical model that considers the time-varying magnetic torques on the particle and its motion through the fluid at low Reynolds number. We use both analytical theory and numerical simulation to design magnetic fields that maximize the migration velocity up (or down) an inclined plane. We show how "topotaxis" of spherical particles relies on differences in the hydrodynamic resistance to rotation about axes parallel and perpendicular to the plane. Importantly, the designed fields can drive multiple independent particles to move simultaneously in different directions as determined by gradients in their respective environments. Experiments on ferromagnetic spheres provide evidence for topotactic motions up inclined substrates. The ability to program the autonomous navigation of driven particles within anisotropic environments is relevant to the design of colloidal robots.

Automating Bayesian inference and design to quantify acoustic particle levitation.

Self-propulsion of micro- and nanoparticles powered by ultrasound provides an attractive strategy for the remote manipulation of colloidal matter using biocompatible energy inputs. Quantitative understanding of particle motion and its dependence on size, shape, and composition requires accurate characterization of the acoustic field, which depends sensitively on the experimental setup. Here, we show how automated experiments based on Bayesian inference and design can accurately and efficiently characterize the acoustic field within resonant chambers used to propel acoustic nanomotors. Repeated cycles of observation, inference, and design (OID) are guided by a physical model that describes the rate at which levitating particles approach the nodal plane. Using video microscopy, we observe the relaxation of tracer particles to this plane following the application of the acoustic field. We use sequential Monte Carlo methods to infer model parameters such as the amplitude and frequency of the resonant chamber while accounting for particle-level measurement noise and population-level heterogeneity in the field. Guided by simulated outcomes, we select the optimal design for the next experiment as to maximize the information gain in the relevant parameters. We show how this iterative process serves to discriminate between competing hypotheses and efficiently converges to accurate parameter estimates using only few automated experiments. We discuss the need for model criticism to ensure the validity of the guiding model throughout automated cycles of observation, inference, and design. This work demonstrates how Bayesian methods can learn the parameters of nonlinear, hierarchical models used to describe video microscopy data of active colloids.

Mechanisms of acquired resistance to rapalogs in metastatic renal cell carcinoma.

The mechanistic target of rapamycin (mTOR) is an established therapeutic target in renal cell carcinoma (RCC). Mechanisms of secondary resistance to rapalog therapy in RCC have not been studied previously. We identified six patients with metastatic RCC who initially responded to mTOR inhibitor therapy and then progressed, and had pre-treatment and post-treatment tumor samples available for analysis. We performed deep whole exome sequencing on the paired tumor samples and a blood sample. Sequence data was analyzed using Mutect, CapSeg, Absolute, and Phylogic to identify mutations, copy number changes, and their changes over time. We also performed in vitro functional assays on PBRM1 in RCC cell lines. Five patients had clear cell and one had chromophobe RCC. 434 somatic mutations in 416 genes were identified in the 12 tumor samples. 201 (46%) of mutations were clonal in both samples while 129 (30%) were acquired in the post-treatment samples. Tumor heterogeneity or sampling issues are likely to account for some mutations that were acquired in the post-treatment samples. Three samples had mutations in TSC1; one in PTEN; and none in MTOR. PBRM1 was the only gene in which mutations were acquired in more than one post-treatment sample. We examined the effect of PBRM1 loss in multiple RCC cell lines, and could not identify any effect on rapalog sensitivity in in vitro culture assays. We conclude that mTOR pathway gene mutations did not contribute to rapalog resistance development in these six patients with advanced RCC. Furthermore, mechanisms of resistance to rapalogs in RCC remain unclear and our results suggest that PBRM1 loss may contribute to sensitivity through complex transcriptional effects.

Accelerating the Design of Self-Guided Microrobots in Time-Varying Magnetic Fields.

Mobile robots combine sensory information with mechanical actuation to move autonomously through structured environments and perform specific tasks. The miniaturization of such robots to the size of living cells is actively pursued for applications in biomedicine, materials science, and environmental sustainability. Existing microrobots based on field-driven particles rely on knowledge of the particle position and the target destination to control particle motion through fluid environments. Often, however, these external control strategies are challenged by limited information and global actuation where a common field directs multiple robots with unknown positions. In this Perspective, we discuss how time-varying magnetic fields can be used to encode the self-guided behaviors of magnetic particles conditioned on local environmental cues. Programming these behaviors is framed as a design problem: we seek to identify the design variables (e.g., particle shape, magnetization, elasticity, stimuli-response) that achieve the desired performance in a given environment. We discuss strategies for accelerating the design process using automated experiments, computational models, statistical inference, and machine learning approaches. Based on the current understanding of field-driven particle dynamics and existing capabilities for particle fabrication and actuation, we argue that self-guided microrobots with potentially transformative capabilities are close at hand.

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.

Combined tumor and immune signals from genomes or transcriptomes predict outcomes of checkpoint inhibition in melanoma.

Immune checkpoint blockade (CPB) improves melanoma outcomes, but many patients still do not respond. Tumor mutational burden (TMB) and tumor-infiltrating T cells are associated with response, and integrative models improve survival prediction. However, integrating immune/tumor-intrinsic features using data from a single assay (DNA/RNA) remains underexplored. Here, we analyze whole-exome and bulk RNA sequencing of tumors from new and published cohorts of 189 and 178 patients with melanoma receiving CPB, respectively. Using DNA, we calculate T cell and B cell burdens (TCB/BCB) from rearranged TCR/Ig sequences and find that patients with TMBhigh and TCBhigh or BCBhigh have improved outcomes compared to other patients. By combining pairs of immune- and tumor-expressed genes, we identify three gene pairs associated with response and survival, which validate in independent cohorts. The top model includes lymphocyte-expressed MAP4K1 and tumor-expressed TBX3. Overall, RNA or DNA-based models combining immune and tumor measures improve predictions of melanoma CPB outcomes.

Longitudinal Single-Cell Dynamics of Chromatin Accessibility and Mitochondrial Mutations in Chronic Lymphocytic Leukemia Mirror Disease History.

While cancers evolve during disease progression and in response to therapy, temporal dynamics remain difficult to study in humans due to the lack of consistent barcodes marking individual clones in vivo. We employ mitochondrial single-cell assay for transposase-accessible chromatin with sequencing to profile 163,279 cells from 9 patients with chronic lymphocytic leukemia (CLL) collected across disease course and utilize mitochondrial DNA (mtDNA) mutations as natural genetic markers of cancer clones. We observe stable propagation of mtDNA mutations over years in the absence of strong selective pressure, indicating clonal persistence, but dramatic changes following tight bottlenecks, including disease transformation and relapse posttherapy, paralleled by acquisition of copy-number variants and changes in chromatin accessibility and gene expression. Furthermore, we link CLL subclones to distinct chromatin states, providing insight into nongenetic sources of relapse. mtDNA mutations thus mirror disease history and provide naturally occurring genetic barcodes to enable patient-specific study of cancer subclonal dynamics. SIGNIFICANCE: Single-cell multi-omic profiling of CLL reveals the utility of somatic mtDNA mutations as in vivo barcodes, which mark subclones that can evolve over time along with changes in accessible chromatin and gene expression profiles to capture dynamics of disease evolution. See related commentary by Hilton and Scott, p. 2965. This article is highlighted in the In This Issue feature, p. 2945.

Molecular features of exceptional response to neoadjuvant anti-androgen therapy in high-risk localized prostate cancer.

High-risk localized prostate cancer (HRLPC) is associated with a substantial risk of recurrence and disease mortality. Recent clinical trials have shown that intensifying anti-androgen therapies administered before prostatectomy can induce pathologic complete responses or minimal residual disease, called exceptional response, although the molecular determinants of these clinical outcomes are largely unknown. Here, we perform whole-exome and transcriptome sequencing on pre-treatment multi-regional tumor biopsies from exceptional responders (ERs) and non-responders (NRs, pathologic T3 or lymph node-positive disease) to intensive neoadjuvant anti-androgen therapies. Clonal SPOP mutation and SPOPL copy-number loss are exclusively observed in ERs, while clonal TP53 mutation and PTEN copy-number loss are exclusively observed in NRs. Transcriptional programs involving androgen signaling and TGF-ß signaling are enriched in ERs and NRs, respectively. These findings may guide prospective validation studies of these molecular features in large HRLPC clinical cohorts treated with neoadjuvant anti-androgens to improve patient stratification.

Framework for quality assessment of whole genome cancer sequences.

Bringing together cancer genomes from different projects increases power and allows the investigation of pan-cancer, molecular mechanisms. However, working with whole genomes sequenced over several years in different sequencing centres requires a framework to compare the quality of these sequences. We used the Pan-Cancer Analysis of Whole Genomes cohort as a test case to construct such a framework. This cohort contains whole cancer genomes of 2832 donors from 18 sequencing centres. We developed a non-redundant set of five quality control (QC) measurements to establish a star rating system. These QC measures reflect known differences in sequencing protocol and provide a guide to downstream analyses and allow for exclusion of samples of poor quality. We have found that this is an effective framework of quality measures. The implementation of the framework is available at: https://dockstore.org/containers/quay.io/jwerner_dkfz/pancanqc:1.2.2 .

RNA sequence analysis reveals macroscopic somatic clonal expansion across normal tissues.

How somatic mutations accumulate in normal cells is poorly understood. A comprehensive analysis of RNA sequencing data from ~6700 samples across 29 normal tissues revealed multiple somatic variants, demonstrating that macroscopic clones can be found in many normal tissues. We found that sun-exposed skin, esophagus, and lung have a higher mutation burden than other tested tissues, which suggests that environmental factors can promote somatic mosaicism. Mutation burden was associated with both age and tissue-specific cell proliferation rate, highlighting that mutations accumulate over both time and number of cell divisions. Finally, normal tissues were found to harbor mutations in known cancer genes and hotspots. This study provides a broad view of macroscopic clonal expansion in human tissues, thus serving as a foundation for associating clonal expansion with environmental factors, aging, and risk of disease.

Author Correction: Liquid versus tissue biopsy for detecting acquired resistance and tumor heterogeneity in gastrointestinal cancers.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Mitochondrial Reprogramming Underlies Resistance to BCL-2 Inhibition in Lymphoid Malignancies.

Mitochondrial apoptosis can be effectively targeted in lymphoid malignancies with the FDA-approved B cell lymphoma 2 (BCL-2) inhibitor venetoclax, but resistance to this agent is emerging. We show that venetoclax resistance in chronic lymphocytic leukemia is associated with complex clonal shifts. To identify determinants of resistance, we conducted parallel genome-scale screens of the BCL-2-driven OCI-Ly1 lymphoma cell line after venetoclax exposure along with integrated expression profiling and functional characterization of drug-resistant and engineered cell lines. We identified regulators of lymphoid transcription and cellular energy metabolism as drivers of venetoclax resistance in addition to the known involvement by BCL-2 family members, which were confirmed in patient samples. Our data support the implementation of combinatorial therapy with metabolic modulators to address venetoclax resistance.

Liquid versus tissue biopsy for detecting acquired resistance and tumor heterogeneity in gastrointestinal cancers.

During cancer therapy, tumor heterogeneity can drive the evolution of multiple tumor subclones harboring unique resistance mechanisms in an individual patient1-3. Previous case reports and small case series have suggested that liquid biopsy (specifically, cell-free DNA (cfDNA)) may better capture the heterogeneity of acquired resistance4-8. However, the effectiveness of cfDNA versus standard single-lesion tumor biopsies has not been directly compared in larger-scale prospective cohorts of patients following progression on targeted therapy. Here, in a prospective cohort of 42 patients with molecularly defined gastrointestinal cancers and acquired resistance to targeted therapy, direct comparison of postprogression cfDNA versus tumor biopsy revealed that cfDNA more frequently identified clinically relevant resistance alterations and multiple resistance mechanisms, detecting resistance alterations not found in the matched tumor biopsy in 78% of cases. Whole-exome sequencing of serial cfDNA, tumor biopsies and rapid autopsy specimens elucidated substantial geographic and evolutionary differences across lesions. Our data suggest that acquired resistance is frequently characterized by profound tumor heterogeneity, and that the emergence of multiple resistance alterations in an individual patient may represent the 'rule' rather than the 'exception'. These findings have profound therapeutic implications and highlight the potential advantages of cfDNA over tissue biopsy in the setting of acquired resistance.