Heterogeneity of genomic evolution and mutational profiles in multiple myeloma

Impact of Genes Highly Correlated with MMSET Myeloma on the Survival of Non-MMSET Myeloma Patients

PURPOSE

The poor prognosis of multiple myeloma with t(4;14) is driven by the fusion of genes encoding multiple myeloma SET domain (MMSET) and immunoglobulin heavy chain. Specific genes affected by MMSET and their clinical implications in non-MMSET myeloma remain undetermined.

EXPERIMENTAL DESIGN

We obtained gene expression profiles of 1,032 newly diagnosed myeloma patients enrolled in Total Therapy 2, Total Therapy 3, Myeloma IX, and HOVON65-GMMGHD4 trials and 156 patients from Multiple Myeloma Resource Collection. Probes that correlated most with MMSET myeloma were selected on the basis of a multivariable linear regression and Bonferroni correction and refined on the basis of the strength of association with survival in non-MMSET patients.

RESULTS

Ten MMSET-like probes were associated with poor survival in non-MMSET myeloma. Non-MMSET myeloma patients in the highest quartile of the 10-gene signature (MMSET-like myeloma) had 5-year overall survival similar to that of MMSET myeloma [highest quartile vs. lowest quartile HR = 2.0; 95% confidence interval (CI), 1.5-2.8 in MMSET-like myeloma; HR = 2.3; 95% CI, 1.6-3.3 in MMSET myeloma]. Analyses of MMSET-like gene signature suggested the involvement of p53 and MYC pathways.

CONCLUSIONS

MMSET-like gene signature captures a subset of high-risk myeloma patients underrepresented by conventional risk stratification platforms and defines a distinct biologic subtype. Clin Cancer Res; 22(16); 4039-44. ©2016 AACR.

The kinome 'at large' in cancer

Over the past decade, rapid advances in genomics, proteomics and functional genomics technologies that enable in-depth interrogation of cancer genomes and proteomes and high-throughput analysis of gene function have enabled characterization of the kinome 'at large' in human cancers, providing crucial insights into how members of the protein kinase superfamily are dysregulated in malignancy, the context-dependent functional role of specific kinases in cancer and how kinome remodelling modulates sensitivity to anticancer drugs. The power of these complementary approaches, and the insights gained from them, form the basis of this Analysis article.

The KDM3A-KLF2-IRF4 axis maintains myeloma cell survival

KDM3A is implicated in tumorigenesis; however, its biological role in multiple myeloma (MM) has not been elucidated. Here we identify KDM3A–KLF2–IRF4 axis dependence in MM. Knockdown of KDM3A is toxic to MM cells in vitro and in vivo. KDM3A maintains expression of KLF2 and IRF4 through H3K9 demethylation, and knockdown of KLF2 triggers apoptosis. Moreover, KLF2 directly activates IRF4 and IRF4 reciprocally upregulates KLF2, forming a positive autoregulatory circuit. The interaction of MM cells with bone marrow milieu mediates survival of MM cells. Importantly, silencing of KDM3A, KLF2 or IRF4 both decreases MM cell adhesion to bone marrow stromal cells and reduces MM cell homing to the bone marrow, in association with decreased ITGB7 expression in MAF-translocated MM cell lines. Our results indicate that the KDM3A–KLF2–IRF4 pathway plays an essential role in MM cell survival and homing to the bone marrow, and therefore represents a therapeutic target.

Several histone modifiers have been implicated in the survival of multiple myeloma cells. Here, the authors reveal a role for the histone demethylase KDM3A in the survival of this haematologic cancer, and show that mechanistically KDM3A removes H3K9 methylation from the promoters of KLF2 and IRF4, genes essential for myeloma cell survival.

Targeting the Bone Marrow Microenvironment

Unprecedented advances in multiple myeloma (MM) therapy during the last 15 years are predominantly based on our increasing understanding of the pathophysiologic role of the bone marrow (BM) microenvironment. Indeed, new treatment paradigms, which incorporate thalidomide, immunomodulatory drugs (IMiDs), and proteasome inhibitors, target the tumor cell as well as its BM microenvironment. Ongoing translational research aims to understand in more detail how disordered BM-niche functions contribute to MM pathogenesis and to identify additional derived targeting agents. One of the most exciting advances in the field of MM treatment is the emergence of immune therapies including elotuzumab, daratumumab, the immune checkpoint inhibitors, Bispecific T-cell engagers (BiTes), and Chimeric antigen receptor (CAR)-T cells. This chapter will review our knowledge on the pathophysiology of the BM microenvironment and discuss derived novel agents that hold promise to further improve outcome in MM.

Integrated genomic characterization of IDH1-mutant glioma malignant progression

Gliomas represent approximately 30% of all central nervous system tumors and 80% of malignant brain tumors. To understand the molecular mechanisms underlying the malignant progression of low-grade gliomas with mutations in IDH1 (encoding isocitrate dehydrogenase 1), we studied paired tumor samples from 41 patients, comparing higher-grade, progressed samples to their lower-grade counterparts. Integrated genomic analyses, including whole-exome sequencing and copy number, gene expression and DNA methylation profiling, demonstrated nonlinear clonal expansion of the original tumors and identified oncogenic pathways driving progression. These include activation of the MYC and RTK-RAS-PI3K pathways and upregulation of the FOXM1- and E2F2-mediated cell cycle transitions, as well as epigenetic silencing of developmental transcription factor genes bound by Polycomb repressive complex 2 in human embryonic stem cells. Our results not only provide mechanistic insight into the genetic and epigenetic mechanisms driving glioma progression but also identify inhibition of the bromodomain and extraterminal (BET) family as a potential therapeutic approach.

Tumor Evolutionary Principles: How Intratumor Heterogeneity Influences Cancer Treatment and Outcome

Recent studies have shown that intratumor heterogeneity contributes to drug resistance in advanced disease. Intratumor heterogeneity may foster the selection of a resistant subclone, sometimes detectable prior to treatment. Next-generation sequencing is enabling the phylogenetic reconstruction of a cancer's life history and has revealed different modes of cancer evolution. These studies have shown that cancer evolution is not always stochastic and has certain constraints. Consideration of cancer evolution may enable the better design of clinical trials and cancer therapeutics. In this review, we summarize the different modes of cancer evolution and how this might impact clinical outcomes. Furthermore, we will discuss several therapeutic strategies for managing emergent intratumor heterogeneity.

Genomic Aberrations in Multiple Myeloma

Multiple myeloma (MM) is a genetically complex disease. The past few years have seen an evolution in cancer research with the emergence of next-generation sequencing (NGS), enabling high throughput sequencing of tumors-including whole exome, whole genome, RNA, and single-cell sequencing as well as genome-wide association study (GWAS). A few inherited variants have been described, counting for some cases of familial disease. Hierarchically, primary events in MM can be divided into hyperdiploid (HDR) and nonhyperdiploid subtypes. HRD tumors are characterized by trisomy of chromosomes 3, 5, 7, 9, 11, 15, 19, and/or 21. Non-HRD tumors harbor IGH translocations, mainly t(4;14), t(6;14), t(11;14), t(14;16), and t(14;20). Secondary events participate to the tumor progression and consist in secondary translocation involving MYC, copy number variations (CNV) and somatic mutations (such as mutations in KRAS, NRAS, BRAF, P53). Moreover, the dissection of clonal heterogeneity helps to understand the evolution of the disease. The following review provides a comprehensive review of the genomic landscape in MM.

Osteoclasts control reactivation of dormant myeloma cells by remodelling the endosteal niche

Multiple myeloma is largely incurable, despite development of therapies that target myeloma cell-intrinsic pathways. Disease relapse is thought to originate from dormant myeloma cells, localized in specialized niches, which resist therapy and repopulate the tumour. However, little is known about the niche, and how it exerts cell-extrinsic control over myeloma cell dormancy and reactivation. In this study, we track individual myeloma cells by intravital imaging as they colonize the endosteal niche, enter a dormant state and subsequently become activated to form colonies. We demonstrate that dormancy is a reversible state that is switched ‘on' by engagement with bone-lining cells or osteoblasts, and switched ‘off' by osteoclasts remodelling the endosteal niche. Dormant myeloma cells are resistant to chemotherapy that targets dividing cells. The demonstration that the endosteal niche is pivotal in controlling myeloma cell dormancy highlights the potential for targeting cell-extrinsic mechanisms to overcome cell-intrinsic drug resistance and prevent disease relapse.

Therapy resistant dormant myeloma cells contribute to disease relapse. Here, the authors use intravital microscopy to track the location of these cells and demonstrate that they hone to the endosteal niche within the bone.

Clock-like mutational processes in human somatic cells

During the course of a lifetime, somatic cells acquire mutations. Different mutational processes may contribute to the mutations accumulated in a cell, with each imprinting a mutational signature on the cell's genome. Some processes generate mutations throughout life at a constant rate in all individuals, and the number of mutations in a cell attributable to these processes will be proportional to the chronological age of the person. Using mutations from 10,250 cancer genomes across 36 cancer types, we investigated clock-like mutational processes that have been operating in normal human cells. Two mutational signatures show clock-like properties. Both exhibit different mutation rates in different tissues. However, their mutation rates are not correlated, indicating that the underlying processes are subject to different biological influences. For one signature, the rate of cell division may influence its mutation rate. This study provides the first survey of clock-like mutational processes operating in human somatic cells.

Implications of Tumor Clonal Heterogeneity in the Era of Next-Generation Sequencing

Recent whole-genome sequencing (WGS) studies have demonstrated that tumors typically comprise a founding clone and multiple subclones (i.e., clonal heterogeneity is common). The possible combination of mutations in each tumor clone is enormous, making each tumor genetically unique. Clonal heterogeneity likely has a role in cancer progression, relapse, metastasis, and chemoresistance due to functional differences in genetically unique subclones. In current clinical practice, gene mutations are only classified as being present or absent, ignoring the clonal complexity of cancers. In this review, we address how tumor clonality is measured using next-generation sequencing (NGS) data, highlight that clonal heterogeneity is common across multiple tumor types, and discuss the potential clinical implications of tumor clonal heterogeneity.

Mutational Spectrum, Copy Number Changes, and Outcome: Results of a Sequencing Study of Patients With Newly Diagnosed Myeloma

PURPOSE

At the molecular level, myeloma is characterized by copy number abnormalities and recurrent translocations into the immunoglobulin heavy chain locus. Novel methods, such as massively parallel sequencing, have begun to describe the pattern of tumor-acquired mutations, but their clinical relevance has yet to be established.

METHODS

We performed whole-exome sequencing for 463 patients who presented with myeloma and were enrolled onto the National Cancer Research Institute Myeloma XI trial, for whom complete molecular cytogenetic and clinical outcome data were available.

RESULTS

We identified 15 significantly mutated genes: IRF4, KRAS, NRAS, MAX, HIST1H1E, RB1, EGR1, TP53, TRAF3, FAM46C, DIS3, BRAF, LTB, CYLD, and FGFR3. The mutational spectrum is dominated by mutations in the RAS (43%) and nuclear factor-κB (17%) pathways, but although they are prognostically neutral, they could be targeted therapeutically. Mutations in CCND1 and DNA repair pathway alterations (TP53, ATM, ATR, and ZNFHX4 mutations) are associated with a negative impact on survival. In contrast, those in IRF4 and EGR1 are associated with a favorable overall survival. We combined these novel mutation risk factors with the recurrent molecular adverse features and international staging system to generate an international staging system mutation score that can identify a high-risk population of patients who experience relapse and die prematurely.

CONCLUSION

We have refined our understanding of genetic events in myeloma and identified clinically relevant mutations that may be used to better stratify patients at presentation.

SnapShot: Multiple Myeloma

Multiple myeloma (MM) is a plasma cell malignancy characterized by a heterogeneous clinical presentation. Genetic abnormalities are not only key events in the origin and progression of the disease but are also useful tools for prognosis, risk stratification, and therapeutic decision making. Although still incurable, a revolution in the treatment of MM is currently ongoing, leading to a significant improvement of clinical outcome and survival. To view this SnapShot, open or download the PDF.

Understanding the role of hyperdiploidy in myeloma prognosis: which trisomies really matter?

The prognosis of multiple myeloma is mainly dependent upon chromosomal changes. The 2 major abnormalities driving poor outcome are del(17p) and t(4;14). However, the outcome of these high-risk patients is not absolutely uniform, with some patients presenting long survival. We hypothesized that these better outcomes might be related to concomitant "good-risk" chromosomal changes exploring hyperdiploidy. We analyzed a large series of 965 myeloma patients, including 168 patients with t(4;14) and 126 patients with del(17p), using high-throughput single-nucleotide polymorphism arrays after plasma cell sorting. As expected, trisomic chromosomes were highly associated. Using the LASSO model, we found that only chromosome 3, when trisomic, was associated with a longer progression-free survival and that 3 trisomies modulated overall survival (OS) in myeloma patients: trisomies 3 and 5 significantly improved OS, whereas trisomy 21 worsened OS. In patients with t(4;14), trisomies 3 and/or 5 seemed to overcome the poor prognosis. For the first time, using a specific modeling approach, we show that not all trisomies display the same prognostic impact. This finding could be important for routine assessment of prognosis in myeloma, and some high-risk patients with a traditional evaluation could in fact be standard-risk patients.

A mutational signature in gastric cancer suggests therapeutic strategies

Targeting defects in the DNA repair machinery of neoplastic cells, for example, those due to inactivating BRCA1 and/or BRCA2 mutations, has been used for developing new therapies in certain types of breast, ovarian and pancreatic cancers. Recently, a mutational signature was associated with failure of double-strand DNA break repair by homologous recombination based on its high mutational burden in samples harbouring BRCA1 or BRCA2 mutations. In pancreatic cancer, all responders to platinum therapy exhibit this mutational signature including a sample that lacked any defects in BRCA1 or BRCA2. Here, we examine 10,250 cancer genomes across 36 types of cancer and demonstrate that, in addition to breast, ovarian and pancreatic cancers, gastric cancer is another cancer type that exhibits this mutational signature. Our results suggest that 7–12% of gastric cancers have defective double-strand DNA break repair by homologous recombination and may benefit from either platinum therapy or PARP inhibitors.

Cancer genome analysis has demonstrated that some breast and ovarian tumours show reduced homologous recombination, a feature that can be therapeutically exploited. Here, Alexandrov et al. search for this mutational signature in 36 different cancer types and find that some gastric tumours also harbour this mutational spectrum.

Mechanisms and Clinical Applications of Genome Instability in Multiple Myeloma

Ongoing genomic instability represents a hallmark of multiple myeloma (MM) cells, which manifests largely as whole chromosome- or translocation-based aneuploidy. Importantly, although it supports tumorigenesis, progression and, response to treatment in MM patients, it remains one of the least understood components of malignant transformation in terms of molecular basis. Therefore these aspects make the comprehension of genomic instability a pioneering strategy for novel therapeutic and clinical speculations to use in the management of MM patients. Here we will review mechanisms mediating genomic instability in MM cells with an emphasis placed on pathogenic mutations affecting DNA recombination, replication and repair, telomere function and mitotic regulation of spindle attachment, centrosome function, and chromosomal segregation. We will discuss the mechanisms by which genetic aberrations give rise to multiple pathogenic events required for myelomagenesis and conclude with a discussion of the clinical applications of these findings in MM patients.

The Tao of myeloma

Multiple myeloma is a plasma cell malignancy in which significant advances have been observed during the last 15 years. Our understanding of the disease has been advanced through its molecular characterization. We have also seen improvements in patient care with the development of 2 new classes of active agents, proteasome inhibitors and immunomodulatory drugs (IMiDs), resulting in a significant improvement in overall survival of myeloma patients such that it can now be debated as to whether some subsets of myeloma patients can be cured. However, the advances in our understanding of myeloma biology occurred in parallel with advances in treatment as opposed to being directly informed by the research. Moreover, the molecular characterization of malignant plasma cells would not have predicted the effectiveness of these novel therapies.We hypothesize that proteasome inhibitors and IMiDs are highly active because malignant plasma cells are constrained by many of the characteristics of their normal counterparts and these novel therapies target both normal plasma cell biology and the cancer biology of myeloma. Thus, a better understanding of normal plasma cell biology will likely yield as many actionable targets as mapping the genomic landscape of this disease.

Cancer Genomics in Clinical Context

Precision medicine requires appropriate application of genomics in clinical practice. In cancer, we have witnessed practice-changing examples of how genomic knowledge is translated into more tailored and effective therapies. The next opportunity is to embed cancer genomics in clinical context so that patient-centric longitudinal clinical, genomic, and molecular phenotypes can be compiled for adaptive learning between precision medicine research and clinical care with the goal of accelerating clinically-actionable discoveries. We describe here an adaptive learning platform, APOLLO™ (adaptive patient-oriented longitudinal learning and optimization) designed to integrate genomic research in the context of, but not in the path of, routine and investigational clinical care for purposes of enabling data-driven discovery across disciplines such that every patient can contribute to and potentially benefit from research discoveries.

Evolving Paradigms in the Management of Multiple Myeloma: Novel Agents and Targeted Therapies

Multiple myeloma (MM) is a clonal plasma cell disorder defined by bone marrow infiltration and osteolytic bone lesions and is the second most common hematologic malignancy after non-Hodgkin lymphoma. The landscape of MM treatment was transformed at the dawn of the twenty-first century by the introduction of novel agents including proteasome inhibitors (bortezomib) and immunomodulatory drugs (thalidomide, lenalidomide), which have prolonged the survival of MM patients. The recently revised International Myeloma Working Group diagnostic criteria for MM added validated biomarkers (clonal bone marrow plasma cell ≥60%, involved:uninvolved serum free light chain ratio ≥100, or >1 focal lesion on magnetic resonance imaging) to identify near inevitable progression to symptomatic MM requiring therapy. In addition, the definition of myeloma-defining CRAB features (hypercalcemia, renal failure, anemia, and bone lesions) has been refined based on advances in imaging and laboratory techniques since the 2003 IMWG consensus. Despite expanded treatment options, MM remains an incurable disease. Drug resistance and clonal evolution remain problematic, and novel therapeutic agents are needed. New approaches to myeloma treatment include anti-CD38 antibodies, next generation proteasome inhibitors, epigenetic modulation with histone deacetylase inhibitors, and targeting the tumor microenvironment. In this article, the diagnosis, staging, and prognostic stratification of newly diagnosed MM will be reviewed. Clinical data pertaining to the emerging targeted agents will be discussed, and a suggested framework for integration of these new therapeutic options will be provided.

Common progenitor cells in mature B-cell malignancies: implications for therapy

PURPOSE OF REVIEW

This review summarizes the recent progress in defining the patterns of genetic evolution giving rise to relapse in follicular lymphoma and multiple myeloma, and discusses their implications with respect to 'personalized medicine'.

RECENT FINDINGS

High-throughput sequencing studies have uncovered a large degree of clonal heterogeneity within tumors, and found that subclones have a variable contribution to relapse. Recent studies aimed at defining patterns of clonal evolution have revealed that serial tumors in some malignancies share their ancestry in a less evolved common progenitor cell (CPC) that bears only a subset of the mutations that are present in the fully evolved tumors that present clinically. This pattern of 'divergent evolution' means that the majority of 'actionable mutations' in tumor specimens are absent within the progenitors that give rise to relapse.

SUMMARY

Follicular lymphoma and multiple myeloma are clinically, biologically and genetically distinct mature B-cell malignancies. However, recent studies have found them to share important similarities in their patterns of genetic evolution. Tumor cells that constitute subclonal populations within these tumors, or between consecutive tumors, share their origins within a genetically less evolved common progenitor cell. This pattern of evolution indicates that current therapies are unable to eradicate these less evolved populations that are at the root of relapse. This suggests that in order to obtain the best results with modern 'targeted therapies' that are directed towards 'actionable mutations', these mutations should be considered within the context of the evolutionary stage at which mutations are acquired, not simply on a presence or absence basis.

SP140L, an Evolutionarily Recent Member of the SP100 Family, Is an Autoantigen in Primary Biliary Cirrhosis

The SP100 family members comprise a set of closely related genes on chromosome 2q37.1. The widely expressed SP100 and the leukocyte-specific proteins SP110 and SP140 have been associated with transcriptional regulation and various human diseases. Here, we have characterized the SP100 family member SP140L. The genome sequence analysis showed the formation of SP140L gene through rearrangements of the two neighboring genes, SP100 and SP140, during the evolution of higher primates. The SP140L expression is interferon-inducible with high transcript levels in B cells and other peripheral blood mononuclear cells. Subcellularly, SP140L colocalizes with SP100 and SP140 in nuclear structures that are devoid of SP110, PML, or p300 proteins. Similarly to SP100 and SP140 protein, we detected serum autoantibodies to SP140L in patients with primary biliary cirrhosis using luciferase immunoprecipitation system and immunoblotting assays. In conclusion, our results show that SP140L is phylogenetically recent member of SP100 proteins and acts as an autoantigen in primary biliary cirrhosis patients.

Intratumor heterogeneity and transcriptional profiling in glioblastoma: translational opportunities

The study of phenotypic and genetic intratumor heterogeneity in glioblastoma is attracting a lot of attention. Recent studies have demonstrated that transcriptional profiling analysis can help interpret the complexity of this disease. Previously proposed molecular classifiers have been recently challenged due to the unexpected degree of intratumor heterogeneity that has been described spatially and at single-cell level. Different computational methods have been employed to analyze this huge amount of data, but new experimental designs including multisampling from individual patients and single-cell experiments require new specific approaches. In light of these results, there is hope that integration of genetic, phenotypic and transcriptional data coupled with functional experiments might help define new therapeutic strategies and classify patients according to key pathways and molecular targets that can be further investigated to develop personalized and combinatorial treatment strategies.

Time to redefine Myeloma

In November 2014 the International Myeloma Working Group (IMWG) revised the definition of multiple myeloma, such that asymptomatic patients with newly diagnosed multiple myeloma without any of the traditional 'CRAB' (hypercalcaemia, renal impairment, anaemia, bone disease) end organ damage criteria but with one of three new criteria would be recommended to start treatment. Previously, the standard of care for such patients was expectant management. These three new criteria are: greater than 60% clonal plasma cells on bone marrow biopsy, a serum free light chain (sFLC) ratio of >100 (the involved sFLC must be >100 mg/l) and greater than one unequivocal focal lesion on advanced imaging (low dose whole body computerized tomography, magnetic resonance imaging, (18) F fluorodeoxyglucose positron emission tomography). Although this would appear to affect a small number of patients, the impact of these changes are broad, leading to an increased use of advanced imaging, a debate around the management of patients previously diagnosed with smouldering myeloma, changed terminology and clinical trial design and an extension of the use of biomarkers. For the first time the philosophy of treatment in myeloma will change from treatment initiation only being triggered by overt end organ damage to an era where sub clinical risk factors will also be taken into account.

A functional variant that affects exon-skipping and protein expression of SP140 as genetic mechanism predisposing to multiple sclerosis

Several variants in strong linkage disequilibrium (LD) at the SP140 locus have been associated with multiple sclerosis (MS), Crohn's disease (CD) and chronic lymphocytic leukemia (CLL). To determine the causal polymorphism, we have integrated high-density data sets of expression quantitative trait loci (eQTL), using GEUVADIS RNA sequences and 1000 Genomes genotypes, with MS-risk variants of the high-density Immunochip array performed by the International Multiple Sclerosis Genetic Consortium (IMSGC). The variants most associated with MS were also correlated with a decreased expression of the full-length RNA isoform of SP140 and an increase of an isoform lacking exon 7. By exon splicing assay, we have demonstrated that the rs28445040 variant was the causal factor for skipping of exon 7. Western blots of peripheral blood mononuclear cells from MS patients showed a significant allele-dependent reduction of the SP140 protein expression. To confirm the association of this functional variant with MS and to compare it with the best-associated variant previously reported by GWAS (rs10201872), a case-control study including 4384 MS patients and 3197 controls was performed. Both variants, in strong LD (r(2) = 0.93), were found similarly associated with MS [P-values, odds ratios: 1.9E-9, OR = 1.35 (1.22-1.49) and 4.9E-10, OR = 1.37 (1.24-1.51), respectively]. In conclusion, our data uncover the causal variant for the SP140 locus and the molecular mechanism associated with MS risk. In addition, this study and others previously reported strongly suggest that this functional variant may be shared with other immune-mediated diseases as CD and CLL.

Subclonal diversification of primary breast cancer revealed by multiregion sequencing

The sequencing of cancer genomes may enable tailoring of therapeutics to the underlying biological abnormalities driving a particular patient's tumor. However, sequencing-based strategies rely heavily on representative sampling of tumors. To understand the subclonal structure of primary breast cancer, we applied whole-genome and targeted sequencing to multiple samples from each of 50 patients' tumors (303 samples in total). The extent of subclonal diversification varied among cases and followed spatial patterns. No strict temporal order was evident, with point mutations and rearrangements affecting the most common breast cancer genes, including PIK3CA, TP53, PTEN, BRCA2 and MYC, occurring early in some tumors and late in others. In 13 out of 50 cancers, potentially targetable mutations were subclonal. Landmarks of disease progression, such as resistance to chemotherapy and the acquisition of invasive or metastatic potential, arose within detectable subclones of antecedent lesions. These findings highlight the importance of including analyses of subclonal structure and tumor evolution in clinical trials of primary breast cancer.

APOBEC Enzymes: Mutagenic Fuel for Cancer Evolution and Heterogeneity

Deep sequencing technologies are revealing the complexities of cancer evolution, casting light on mutational processes fueling tumor adaptation, immune escape, and treatment resistance. Understanding mechanisms driving cancer diversity is a critical step toward developing strategies to attenuate tumor evolution and adaptation. One emerging mechanism fueling tumor diversity and subclonal evolution is genomic DNA cytosine deamination catalyzed by APOBEC3B and at least one other APOBEC family member. Deregulation of APOBEC3 enzymes causes a general mutator phenotype that manifests as diverse and heterogeneous tumor subclones. Here, we summarize knowledge of the APOBEC DNA deaminase family in cancer, and their role as driving forces for intratumor heterogeneity and a therapeutic target to limit tumor adaptation.

SIGNIFICANCE

APOBEC mutational signatures may be enriched in tumor subclones, suggesting APOBEC cytosine deaminases fuel subclonal expansions and intratumor heterogeneity. APOBEC family members might represent a new class of drug target aimed at limiting tumor evolution, adaptation, and drug resistance.

Synthetic Lethal Approaches Exploiting DNA Damage in Aggressive Myeloma

Ongoing DNA damage is a common feature of epithelial cancers. Here, we show that tumor cells derived from multiple myeloma, a disease of clonal plasma cells, demonstrate DNA-replicative stress, leading to DNA damage. We identified a poor-prognosis subset of multiple myeloma with extensive chromosomal instability and replicative stress, which rely on ATR to compensate for DNA-replicative stress; conversely, silencing of ATR or treatment with a specific ATR inhibitor triggers multiple myeloma cell apoptosis. We show that oncogenes, such as MYC, induce DNA damage in multiple myeloma cells not only by increased replicative stress, but also via increased oxidative stress, and that reactive oxygen species-inducer piperlongumine triggers further DNA damage and apoptosis. Importantly, ATR inhibition combined with piperlongumine triggers synergistic multiple myeloma cytotoxicity. This synthetic lethal approach, enhancing oxidative stress while concomitantly blocking replicative stress response, provides a novel combination targeted therapy to address an unmet medical need in this subset of multiple myeloma.

SIGNIFICANCE

Multiple myeloma remains an incurable disease. We have identified a subset of multiple myeloma patients with poor prognosis, whose tumors present chromosomal instability, replicative and oxidative stress, and DNA damage. We define a synthetic lethal approach enhancing oxidative stress while targeting replicative stress response, inducing tumor cell apoptosis in this patient subset. Cancer Discov; 5(9); 972-87. ©2015 AACR.This article is highlighted in the In This Issue feature, p. 893.

Single-molecule analysis reveals widespread structural variation in multiple myeloma

Multiple myeloma (MM), a malignancy of plasma cells, is characterized by widespread genomic heterogeneity and, consequently, differences in disease progression and drug response. Although recent large-scale sequencing studies have greatly improved our understanding of MM genomes, our knowledge about genomic structural variation in MM is attenuated due to the limitations of commonly used sequencing approaches. In this study, we present the application of optical mapping, a single-molecule, whole-genome analysis system, to discover new structural variants in a primary MM genome. Through our analysis, we have identified and characterized widespread structural variation in this tumor genome. Additionally, we describe our efforts toward comprehensive characterization of genome structure and variation by integrating our findings from optical mapping with those from DNA sequencing-based genomic analysis. Finally, by studying this MM genome at two time points during tumor progression, we have demonstrated an increase in mutational burden with tumor progression at all length scales of variation.

Massive parallel IGHV gene sequencing reveals a germinal center pathway in origins of human multiple myeloma

Human multiple myeloma (MM) is characterized by accumulation of malignant terminally differentiated plasma cells (PCs) in the bone marrow (BM), raising the question when during maturation neoplastic transformation begins. Immunoglobulin IGHV genes carry imprints of clonal tumor history, delineating somatic hypermutation (SHM) events that generally occur in the germinal center (GC). Here, we examine MM-derived IGHV genes using massive parallel deep sequencing, comparing them with profiles in normal BM PCs. In 4/4 presentation IgG MM, monoclonal tumor-derived IGHV sequences revealed significant evidence for intraclonal variation (ICV) in mutation patterns. IGHV sequences of 2/2 normal PC IgG populations revealed dominant oligoclonal expansions, each expansion also displaying mutational ICV. Clonal expansions in MM and in normal BM PCs reveal common IGHV features. In such MM, the data fit a model of tumor origins in which neoplastic transformation is initiated in a GC B-cell committed to terminal differentiation but still targeted by on-going SHM. Strikingly, the data parallel IGHV clonal sequences in some monoclonal gammopathy of undetermined significance (MGUS) known to display on-going SHM imprints. Since MGUS generally precedes MM, these data suggest origins of MGUS and MM with IGHV gene mutational ICV from the same GC B-cell, arising via a distinctive pathway.

Genomic heterogeneity in multiple myeloma

Multiple myeloma (MM) is an incurable malignancy in majority of patients characterized by clonal proliferation of plasma cells. To date, treatment is established based on general conditions and age of patients. However, MM is a heterogeneous disease, featured by various subtypes and different outcomes. Thus, the understanding of MM biology is currently a major challenge to eventually cure the disease. During the last decade, karyotype studies and gene expression profiling have identified robust prognostic markers as well as a widespread genomic landscape. More recently, studies of epigenetic, transcriptional modifications and next generation sequencing have allowed characterization of critical genes and pathways, clonal heterogeneity and mutational profiles involved in myelomagenesis. Altogether, these findings constitute important tools to develop new targeted and personalized therapies in MM.

Prognostic indicators of lenalidomide for multiple myeloma: consensus and controversy

The long-term outcome of multiple myeloma (MM) has been greatly improved through new agents, one being lenalidomide (LEN). Based upon the findings of in vitro experiments, its mode of action against MM occurs through a combination of direct tumoricidal effects on myeloma cells, modulatory effects on tumor immunity and tumor microenvironment-regulatory effects. However, it has not been clearly defined whether the clinical response and long-term outcome of MM with LEN treatment truly reflect the mechanisms of action of LEN proposed by in vitro studies. To ascertain what is known and what remains to be elucidated with LEN, we review the current literature on the mode of action of LEN in association with myeloma pathophysiology, and discuss the prognostic indicators in the treatment of MM with LEN.

Glioblastoma adaptation traced through decline of an IDH1 clonal driver and macro-evolution of a double-minute chromosome

BACKGROUND

Glioblastoma (GBM) is the most common malignant brain cancer occurring in adults, and is associated with dismal outcome and few therapeutic options. GBM has been shown to predominantly disrupt three core pathways through somatic aberrations, rendering it ideal for precision medicine approaches.

METHODS

We describe a 35-year-old female patient with recurrent GBM following surgical removal of the primary tumour, adjuvant treatment with temozolomide and a 3-year disease-free period. Rapid whole-genome sequencing (WGS) of three separate tumour regions at recurrence was carried out and interpreted relative to WGS of two regions of the primary tumour.

RESULTS

We found extensive mutational and copy-number heterogeneity within the primary tumour. We identified a TP53 mutation and two focal amplifications involving PDGFRA, KIT and CDK4, on chromosomes 4 and 12. A clonal IDH1 R132H mutation in the primary, a known GBM driver event, was detectable at only very low frequency in the recurrent tumour. After sub-clonal diversification, evidence was found for a whole-genome doubling event and a translocation between the amplified regions of PDGFRA, KIT and CDK4, encoded within a double-minute chromosome also incorporating miR26a-2. The WGS analysis uncovered progressive evolution of the double-minute chromosome converging on the KIT/PDGFRA/PI3K/mTOR axis, superseding the IDH1 mutation in dominance in a mutually exclusive manner at recurrence, consequently the patient was treated with imatinib. Despite rapid sequencing and cancer genome-guided therapy against amplified oncogenes, the disease progressed, and the patient died shortly after.

CONCLUSION

This case sheds light on the dynamic evolution of a GBM tumour, defining the origins of the lethal sub-clone, the macro-evolutionary genomic events dominating the disease at recurrence and the loss of a clonal driver. Even in the era of rapid WGS analysis, cases such as this illustrate the significant hurdles for precision medicine success.

The ribonuclease DIS3 promotes let-7 miRNA maturation by degrading the pluripotency factor LIN28B mRNA

Multiple myeloma, the second most frequent hematologic tumor after lymphomas, is an incurable cancer. Recent sequencing efforts have identified the ribonuclease DIS3 as one of the most frequently mutated genes in this disease. DIS3 represents the catalytic subunit of the exosome, a macromolecular complex central to the processing, maturation and surveillance of various RNAs. miRNAs are an evolutionarily conserved class of small noncoding RNAs, regulating gene expression at post-transcriptional level. Ribonucleases, including Drosha, Dicer and XRN2, are involved in the processing and stability of miRNAs. However, the role of DIS3 on the regulation of miRNAs remains largely unknown. Here we found that DIS3 regulates the levels of the tumor suppressor let-7 miRNAs without affecting other miRNA families. DIS3 facilitates the maturation of let-7 miRNAs by reducing in the cytoplasm the RNA stability of the pluripotency factor LIN28B, a inhibitor of let-7 processing. DIS3 inactivation, through the increase of LIN28B and the reduction of mature let-7, enhances the translation of let-7 targets such as MYC and RAS leading to enhanced tumorigenesis. Our study establishes that the ribonuclease DIS3, targeting LIN28B, sustains the maturation of let-7 miRNAs and suggests the increased translation of critical oncogenes as one of the biological outcomes of DIS3 inactivation.

The evolutionary history of lethal metastatic prostate cancer

Cancers emerge from an ongoing Darwinian evolutionary process, often leading to multiple competing subclones within a single primary tumour. This evolutionary process culminates in the formation of metastases, which is the cause of 90% of cancer-related deaths. However, despite its clinical importance, little is known about the principles governing the dissemination of cancer cells to distant organs. Although the hypothesis that each metastasis originates from a single tumour cell is generally supported, recent studies using mouse models of cancer demonstrated the existence of polyclonal seeding from and interclonal cooperation between multiple subclones. Here we sought definitive evidence for the existence of polyclonal seeding in human malignancy and to establish the clonal relationship among different metastases in the context of androgen-deprived metastatic prostate cancer. Using whole-genome sequencing, we characterized multiple metastases arising from prostate tumours in ten patients. Integrated analyses of subclonal architecture revealed the patterns of metastatic spread in unprecedented detail. Metastasis-to-metastasis spread was found to be common, either through de novo monoclonal seeding of daughter metastases or, in five cases, through the transfer of multiple tumour clones between metastatic sites. Lesions affecting tumour suppressor genes usually occur as single events, whereas mutations in genes involved in androgen receptor signalling commonly involve multiple, convergent events in different metastases. Our results elucidate in detail the complex patterns of metastatic spread and further our understanding of the development of resistance to androgen-deprivation therapy in prostate cancer.

Clonal status of actionable driver events and the timing of mutational processes in cancer evolution

Deciphering whether actionable driver mutations are found in all or a subset of tumor cells will likely be required to improve drug development and precision medicine strategies. We analyzed nine cancer types to determine the subclonal frequencies of driver events, to time mutational processes during cancer evolution, and to identify drivers of subclonal expansions. Although mutations in known driver genes typically occurred early in cancer evolution, we also identified later subclonal "actionable" mutations, including BRAF (V600E), IDH1 (R132H), PIK3CA (E545K), EGFR (L858R), and KRAS (G12D), which may compromise the efficacy of targeted therapy approaches. More than 20% of IDH1 mutations in glioblastomas, and 15% of mutations in genes in the PI3K (phosphatidylinositol 3-kinase)-AKT-mTOR (mammalian target of rapamycin) signaling axis across all tumor types were subclonal. Mutations in the RAS-MEK (mitogen-activated protein kinase kinase) signaling axis were less likely to be subclonal than mutations in genes associated with PI3K-AKT-mTOR signaling. Analysis of late mutations revealed a link between APOBEC-mediated mutagenesis and the acquisition of subclonal driver mutations and uncovered putative cancer genes involved in subclonal expansions, including CTNNA2 and ATXN1. Our results provide a pan-cancer census of driver events within the context of intratumor heterogeneity and reveal patterns of tumor evolution across cancers. The frequent presence of subclonal driver mutations suggests the need to stratify targeted therapy response according to the proportion of tumor cells in which the driver is identified.

Pathogenesis beyond the cancer clone(s) in multiple myeloma

Over the past 4 decades, basic research has provided crucial information regarding the cellular and molecular biology of cancer. In particular, the relevance of cancer microenvironment (including both cellular and noncellular elements) and the concept of clonal evolution and heterogeneity have emerged as important in cancer pathogenesis, immunologic escape, and resistance to therapy. Multiple myeloma (MM), a cancer of terminally differentiated plasma cells, is emblematic of the impact of cancer microenvironment and the role of clonal evolution. Although genetic and epigenetic aberrations occur in MM and evolve over time under the pressure of exogenous stimuli, they are also largely present in premalignant plasma cell dyscrasia such as monoclonal gammopathy of undetermined significance (MGUS) and smoldering multiple myeloma (SMM), suggesting that genetic mutations alone are necessary, but not sufficient, for myeloma transformation. The role of bone marrow microenvironment in mediating survival, proliferation, and resistance to therapy in myeloma is well established; and although an appealing speculation, its role in fostering the evolution of MGUS or SMM into MM is yet to be proven. In this review, we discuss MM pathogenesis with a particular emphasis on the role of bone marrow microenvironment.

Tracking the origins and drivers of subclonal metastatic expansion in prostate cancer

Tumour heterogeneity in primary prostate cancer is a well-established phenomenon. However, how the subclonal diversity of tumours changes during metastasis and progression to lethality is poorly understood. Here we reveal the precise direction of metastatic spread across four lethal prostate cancer patients using whole-genome and ultra-deep targeted sequencing of longitudinally collected primary and metastatic tumours. We find one case of metastatic spread to the surgical bed causing local recurrence, and another case of cross-metastatic site seeding combining with dynamic remoulding of subclonal mixtures in response to therapy. By ultra-deep sequencing end-stage blood, we detect both metastatic and primary tumour clones, even years after removal of the prostate. Analysis of mutations associated with metastasis reveals an enrichment of TP53 mutations, and additional sequencing of metastases from 19 patients demonstrates that acquisition of TP53 mutations is linked with the expansion of subclones with metastatic potential which we can detect in the blood.

Inferring mutational timing and reconstructing tumour evolutionary histories

Cancer evolution can be considered within a Darwinian framework. Both micro and macro-evolutionary theories can be applied to understand tumour progression and treatment failure. Owing to cancers' complexity and heterogeneity the rules of tumour evolution, such as the role of selection, remain incompletely understood. The timing of mutational events during tumour evolution presents diagnostic, prognostic and therapeutic opportunities. Here we review the current sampling and computational approaches for inferring mutational timing and the evidence from next generation sequencing-informed data on mutational timing across all tumour types. We discuss how this knowledge can be used to illuminate the genes and pathways that drive cancer initiation and relapse; and to support drug development and clinical trial design.

Cereblon binding molecules in multiple myeloma

Immunomodulation is an established treatment strategy in multiple myeloma with thalidomide and its derivatives lenalidomide and pomalidomide as its FDA approved representatives. Just recently the method of action of these cereblon binding molecules was deciphered and results from large phase 3 trials confirmed the backbone function of this drug family in various combination therapies. This review details the to-date knowledge concerning mechanism of IMiD action, clinical applications and plausible escape mechanisms in which cells may become resistant/refractory to cereblon binding molecule based treatment.

Translational implications of tumor heterogeneity

Advances in next-generation sequencing and bioinformatics have led to an unprecedented view of the cancer genome and its evolution. Genomic studies have demonstrated the complex and heterogeneous clonal landscape of tumors of different origins and the potential impact of intratumor heterogeneity on treatment response and resistance, cancer progression, and the risk of disease relapse. However, the significance of subclonal mutations, in particular mutations in driver genes, and their evolution through time and their dynamics in response to cancer therapies, is yet to be determined. The necessary tools are now available to prospectively determine whether clonal heterogeneity can be used as a biomarker of clinical outcome and to what extent subclonal somatic alterations might influence clinical outcome. Studies that use longitudinal tissue sampling, integrating both genomic and clinical data, have the potential to reveal the subclonal composition and track the evolution of tumors to address these questions and to begin to define the breadth of genetic diversity in different tumor types and its relevance to patient outcome. Such studies may provide further evidence for drug-resistance mechanisms informing combinatorial, adaptive, and tumor immune therapies placed within the context of tumor evolution.

Longitudinal analysis of 25 sequential sample-pairs using a custom multiple myeloma mutation sequencing panel (M(3)P)

Recent advances in genomic sequencing technologies now allow results from deep next-generation sequencing to be obtained within clinically meaningful timeframes, making this an attractive approach to better guide personalized treatment strategies. No multiple myeloma-specific gene panel has been established so far; we therefore designed a 47-gene-targeting gene panel, containing 39 genes known to be mutated in ≥3 % of multiple myeloma cases and eight genes in pathways therapeutically targeted in multiple myeloma (MM). We performed targeted sequencing on tumor/germline DNA of 25 MM patients in which we also had a sequential sample post treatment. Mutation analysis revealed KRAS as the most commonly mutated gene (36 % in each time point), followed by NRAS (20 and 16 %), TP53 (16 and 16 %), DIS3 (16 and 16 %), FAM46C (12 and 16 %), and SP140 (12 and 12 %). We successfully tracked clonal evolution and identified mutation acquisition and/or loss in FAM46C, FAT1, KRAS, NRAS, SPEN, PRDM1, NEB, and TP53 as well as two mutations in XBP1, a gene associated with bortezomib resistance. Thus, we present the first longitudinal analysis of a MM-specific targeted sequencing gene panel that can be used for individual tumor characterization and for tracking clonal evolution over time.