Evolutionary trajectories of small cell lung cancer under therapy.

The evolutionary processes that underlie the marked sensitivity of small cell lung cancer (SCLC) to chemotherapy and rapid relapse are unknown1-3. Here we determined tumour phylogenies at diagnosis and throughout chemotherapy and immunotherapy by multiregion sequencing of 160 tumours from 65 patients. Treatment-naive SCLC exhibited clonal homogeneity at distinct tumour sites, whereas first-line platinum-based chemotherapy led to a burst in genomic intratumour heterogeneity and spatial clonal diversity. We observed branched evolution and a shift to ancestral clones underlying tumour relapse. Effective radio- or immunotherapy induced a re-expansion of founder clones with acquired genomic damage from first-line chemotherapy. Whereas TP53 and RB1 alterations were exclusively part of the common ancestor, MYC family amplifications were frequently not constituents of the founder clone. At relapse, emerging subclonal mutations affected key genes associated with SCLC biology, and tumours harbouring clonal CREBBP/EP300 alterations underwent genome duplications. Gene-damaging TP53 alterations and co-alterations of TP53 missense mutations with TP73, CREBBP/EP300 or FMN2 were significantly associated with shorter disease relapse following chemotherapy. In summary, we uncover key processes of the genomic evolution of SCLC under therapy, identify the common ancestor as the source of clonal diversity at relapse and show central genomic patterns associated with sensitivity and resistance to chemotherapy.

Distinct Genetically Determined Origins of Myd88/BCL2-Driven Aggressive Lymphoma Rationalize Targeted Therapeutic Intervention Strategies.

Genomic profiling revealed the identity of at least 5 subtypes of diffuse large B-cell lymphoma (DLBCL), including the MCD/C5 cluster characterized by aberrations in MYD88, BCL2, PRDM1, and/or SPIB. We generated mouse models harboring B cell-specific Prdm1 or Spib aberrations on the background of oncogenic Myd88 and Bcl2 lesions. We deployed whole-exome sequencing, transcriptome, flow-cytometry, and mass cytometry analyses to demonstrate that Prdm1- or Spib-altered lymphomas display molecular features consistent with prememory B cells and light-zone B cells, whereas lymphomas lacking these alterations were enriched for late light-zone and plasmablast-associated gene sets. Consistent with the phenotypic evidence for increased B cell receptor signaling activity in Prdm1-altered lymphomas, we demonstrate that combined BTK/BCL2 inhibition displays therapeutic activity in mice and in five of six relapsed/refractory DLBCL patients. Moreover, Prdm1-altered lymphomas were immunogenic upon transplantation into immuno-competent hosts, displayed an actionable PD-L1 surface expression, and were sensitive to antimurine-CD19-CAR-T cell therapy, in vivo. SIGNIFICANCE: Relapsed/refractory DLBCL remains a major medical challenge, and most of these patients succumb to their disease. Here, we generated mouse models, faithfully recapitulating the biology of MYD88-driven human DLBCL. These models revealed robust preclinical activity of combined BTK/BCL2 inhibition. We confirmed activity of this regimen in pretreated non-GCB-DLBCL patients. See related commentary by Leveille et al., p. 8. This article is highlighted in the In This Issue feature, p. 1.

Clonal dynamics of BRAF-driven drug resistance in EGFR-mutant lung cancer.

Activation of MAPK signaling via BRAF mutations may limit the activity of EGFR inhibitors in EGFR-mutant lung cancer patients. However, the impact of BRAF mutations on the selection and fitness of emerging resistant clones during anti-EGFR therapy remains elusive. We tracked the evolution of subclonal mutations by whole-exome sequencing and performed clonal analyses of individual metastases during therapy. Complementary functional analyses of polyclonal EGFR-mutant cell pools showed a dose-dependent enrichment of BRAFV600E and a loss of EGFR inhibitor susceptibility. The clones remain stable and become vulnerable to combined EGFR, RAF, and MEK inhibition. Moreover, only osimertinib/trametinib combination treatment, but not monotherapy with either of these drugs, leads to robust tumor shrinkage in EGFR-driven xenograft models harboring BRAFV600E mutations. These data provide insights into the dynamics of clonal evolution of EGFR-mutant tumors and the therapeutic implications of BRAF co-mutations that may facilitate the development of treatment strategies to improve the prognosis of these patients.

CaMuS: simultaneous fitting and de novo imputation of cancer mutational signature.

The identification of the mutational processes operating in tumour cells has implications for cancer diagnosis and therapy. These processes leave mutational patterns on the cancer genomes, which are referred to as mutational signatures. Recently, 81 mutational signatures have been inferred using computational algorithms on sequencing data of 23,879 samples. However, these published signatures may not always offer a comprehensive view on the biological processes underlying tumour types that are not included or underrepresented in the reference studies. To circumvent this problem, we designed CaMuS (Cancer Mutational Signatures) to construct de novo signatures while simultaneously fitting publicly available mutational signatures. Furthermore, we propose to estimate signature similarity by comparing probability distributions using the Hellinger distance. We applied CaMuS to infer signatures of mutational processes in poorly studied cancer types. We used whole genome sequencing data of 56 neuroblastoma, thus providing evidence for the versatility of CaMuS. Using simulated data, we compared the performance of CaMuS to sigfit, a recently developed algorithm with comparable inference functionalities. CaMuS and sigfit reconstructed the simulated datasets with similar accuracy; however two main features may argue for CaMuS over sigfit: (i) superior computational performance and (ii) a reliable parameter selection method to avoid spurious signatures.

Integrative Analysis of Pleomorphic Dermal Sarcomas Reveals Fibroblastic Differentiation and Susceptibility to Immunotherapy.

PURPOSE: Pleomorphic dermal sarcoma (PDS) is a rare malignant cutaneous tumor with an unknown cell of origin. Locally defined tumors can be treated by curative excisions, whereas advanced stages of the disease are difficult to treat, using standard regimens. EXPERIMENTAL DESIGN: We performed whole-exome sequencing on a cohort of 28 individuals and corresponding transcriptomic analysis on 21 patients, as well as quantitative IHC image analysis on 27 patients. RESULTS: PDS exhibits a universally high mutational load (42.7 mutations/mega base) with an inflamed, immunogenic tumor microenvironment. Three cases of PDS showed response to immune checkpoint blockade. Local mutation rate variation together with mRNA expression data demonstrate that PDS form a distinct entity, with PDGFRB as a lineage marker. In addition, we found that PDS is of mesenchymal, fibroblastic differentiation. CONCLUSIONS: PDS is of fibroblastic differentiation and exhibits a strong susceptibility to immunotherapy, including a high mutational burden and an inflamed tumor microenvironment.

Dissecting intratumour heterogeneity of nodal B-cell lymphomas at the transcriptional, genetic and drug-response levels.

Tumour heterogeneity encompasses both the malignant cells and their microenvironment. While heterogeneity between individual patients is known to affect the efficacy of cancer therapy, most personalized treatment approaches do not account for intratumour heterogeneity. We addressed this issue by studying the heterogeneity of nodal B-cell lymphomas by single-cell RNA-sequencing and transcriptome-informed flow cytometry. We identified transcriptionally distinct malignant subpopulations and compared their drug-response and genomic profiles. Malignant subpopulations from the same patient responded strikingly differently to anti-cancer drugs ex vivo, which recapitulated subpopulation-specific drug sensitivity during in vivo treatment. Infiltrating T cells represented the majority of non-malignant cells, whose gene-expression signatures were similar across all donors, whereas the frequencies of T-cell subsets varied significantly between the donors. Our data provide insights into the heterogeneity of nodal B-cell lymphomas and highlight the relevance of intratumour heterogeneity for personalized cancer therapy.

K-ras Mutation Subtypes in NSCLC and Associated Co-occuring Mutations in Other Oncogenic Pathways.

INTRODUCTION: Although KRAS mutations in NSCLC have been considered mutually exclusive driver mutations for a long time, there is now growing evidence that KRAS-mutated NSCLC represents a genetically heterogeneous subgroup. We sought to determine genetic heterogeneity with respect to cancer-related co-mutations and their correlation with different KRAS mutation subtypes. METHODS: Diagnostic samples from 4507 patients with NSCLC were analyzed by next-generation sequencing by using a panel of 14 genes and, in a subset of patients, fluorescence in situ hybridization. Next-generation sequencing with an extended panel of 14 additional genes was performed in 101 patients. Molecular data were correlated with clinical data. Whole-exome sequencing was performed in two patients. RESULTS: We identified 1078 patients with KRAS mutations, of whom 53.5% had at least one additional mutation. Different KRAS mutation subtypes showed different patterns of co-occurring mutations. Besides mutations in tumor protein p53 gene (TP53) (39.4%), serine/threonine kinase 11 gene (STK11) (19.8%), kelch like ECH associated protein 1 gene (KEAP1) (12.9%), and ATM serine/threonine kinase gene (ATM) (11.9%), as well as MNNG HOS Transforming gene (MET) amplifications (15.4%) and erb-b2 receptor tyrosine kinase 2 gene (ERBB2) amplifications (13.8%, exclusively in G12C), we found rare co-occurrence of targetable mutations in EGFR (1.2%) and BRAF (1.2%). Whole-exome sequencing of two patients with co-occurring phosphatidylinositol-4,5-bisphosphate 3-kinase catalytic subunit alpha gene (PIK3CA) mutation revealed clonality of mutated KRAS in one patient and subclonality in the second, suggesting different evolutionary backgrounds. CONCLUSION: KRAS-mutated NSCLC represents a genetically heterogeneous subgroup with a high frequency of co-occurring mutations in cancer-associated pathways, partly associated with distinct KRAS mutation subtypes. This diversity might have implications for understanding the variability of treatment outcome in KRAS-mutated NSCLC and for future trial design.

Clonal dynamics towards the development of venetoclax resistance in chronic lymphocytic leukemia.

Deciphering the evolution of cancer cells under therapeutic pressure is a crucial step to understand the mechanisms that lead to treatment resistance. To this end, we analyzed whole-exome sequencing data of eight chronic lymphocytic leukemia (CLL) patients that developed resistance upon BCL2-inhibition by venetoclax. Here, we report recurrent mutations in BTG1 (2 patients) and homozygous deletions affecting CDKN2A/B (3 patients) that developed during treatment, as well as a mutation in BRAF and a high-level focal amplification of CD274 (PD-L1) that might pinpoint molecular aberrations offering structures for further therapeutic interventions.

Genomic and Functional Fidelity of Small Cell Lung Cancer Patient-Derived Xenografts.

Small cell lung cancer (SCLC) patient-derived xenografts (PDX) can be generated from biopsies or circulating tumor cells (CTC), though scarcity of tissue and low efficiency of tumor growth have previously limited these approaches. Applying an established clinical-translational pipeline for tissue collection and an automated microfluidic platform for CTC enrichment, we generated 17 biopsy-derived PDXs and 17 CTC-derived PDXs in a 2-year timeframe, at 89% and 38% efficiency, respectively. Whole-exome sequencing showed that somatic alterations are stably maintained between patient tumors and PDXs. Early-passage PDXs maintain the genomic and transcriptional profiles of the founder PDX. In vivo treatment with etoposide and platinum (EP) in 30 PDX models demonstrated greater sensitivity in PDXs from EP-naïve patients, and resistance to EP corresponded to increased expression of a MYC gene signature. Finally, serial CTC-derived PDXs generated from an individual patient at multiple time points accurately recapitulated the evolving drug sensitivities of that patient's disease. Collectively, this work highlights the translational potential of this strategy.Significance: Effective translational research utilizing SCLC PDX models requires both efficient generation of models from patients and fidelity of those models in representing patient tumor characteristics. We present approaches for efficient generation of PDXs from both biopsies and CTCs, and demonstrate that these models capture the mutational landscape and functional features of the donor tumors. Cancer Discov; 8(5); 600-15. ©2018 AACR.This article is highlighted in the In This Issue feature, p. 517.

Resource constrained flux balance analysis predicts selective pressure on the global structure of metabolic networks.

BACKGROUND: A universal feature of metabolic networks is their hourglass or bow-tie structure on cellular level. This architecture reflects the conversion of multiple input nutrients into multiple biomass components via a small set of precursor metabolites. However, it is yet unclear to what extent this structural feature is the result of natural selection. RESULTS: We extend flux balance analysis to account for limited cellular resources. Using this model, optimal structure of metabolic networks can be calculated for different environmental conditions. We observe a significant structural reshaping of metabolic networks for a toy-network and E. coli core metabolism if we increase the share of invested resources for switching between different nutrient conditions. Here, hub nodes emerge and the optimal network structure becomes bow-tie-like as a consequence of limited cellular resource constraint. We confirm this theoretical finding by comparing the reconstructed metabolic networks of bacterial species with respect to their lifestyle. CONCLUSIONS: We show that bow-tie structure can give a system-level fitness advantage to organisms that live in highly competitive and fluctuating environments. Here, limitation of cellular resources can lead to an efficiency-flexibility tradeoff where it pays off for the organism to shorten catabolic pathways if they are frequently activated and deactivated. As a consequence, generalists that shuttle between diverse environmental conditions should have a more predominant bow-tie structure than specialists that visit just a few isomorphic habitats during their life cycle.