Multiple myeloma is an incurable malignancy of plasma cells, and its pathogenesis is poorly understood. Here we report the massively parallel sequencing of 38 tumor genomes and their comparison to matched normal DNAs. Several new and unexpected oncogenic mechanisms were suggested by the pattern of somatic mutation across the dataset. These include the mutation of genes involved in protein translation (seen in nearly half of the patients), genes involved in histone methylation, and genes involved in blood coagulation. In addition, a broader than anticipated role of NF-κB signaling was suggested by mutations in 11 members of the NF-κB pathway. Of potential immediate clinical relevance, activating mutations of the kinase BRAF were observed in 4% of patients, suggesting the evaluation of BRAF inhibitors in multiple myeloma clinical trials. These results indicate that cancer genome sequencing of large collections of samples will yield new insights into cancer not anticipated by existing knowledge.

Activation of NF-kappaB has been noted in many tumor types, however only rarely has this been linked to an underlying genetic mutation. An integrated analysis of high-density oligonucleotide array CGH and gene expression profiling data from 155 multiple myeloma samples identified a promiscuous array of abnormalities contributing to the dysregulation of NF-kappaB in approximately 20% of patients. We report mutations in ten genes causing the inactivation of TRAF2, TRAF3, CYLD, cIAP1/cIAP2 and activation of NFKB1, NFKB2, CD40, LTBR, TACI, and NIK that result primarily in constitutive activation of the noncanonical NF-kappaB pathway, with the single most common abnormality being inactivation of TRAF3. These results highlight the critical importance of the NF-kappaB pathway in the pathogenesis of multiple myeloma.

SUMMARY We performed massively parallel sequencing of paired tumor/normal samples from 203 multiple myeloma (MM) patients and identified significantly mutated genes and copy number alterations, and discovered putative tumor suppressor genes by determining homozygous deletions and loss-of-heterozygosity. We observed frequent mutations in KRAS (particularly in previously treated patients), NRAS, BRAF, FAM46C, TP53 and DIS3 (particularly in non-hyperdiploid MM). Mutations were often present in subclonal populations, and multiple mutations within the same pathway (e.g. KRAS, NRAS and BRAF) were observed in the same patient. In vitro modeling predicts only partial treatment efficacy of targeting subclonal mutations, and even growth promotion of non-mutated subclones in some cases. These results emphasize the importance of heterogeneity analysis for treatment decisions.

Understanding the profile of oncogene and tumor suppressor gene mutations with their interactions and impact on the prognosis of multiple myeloma (MM) can improve the definition of disease subsets and identify pathways important in disease pathobiology. Using integrated genomics of 1273 newly diagnosed patients with MM, we identified 63 driver genes, some of which are novel, including ,, ,, and Oncogene mutations are significantly more clonal than tumor suppressor mutations, indicating they may exert a bigger selective pressure. Patients with more driver gene abnormalities are associated with worse outcomes, as are identified mechanisms of genomic instability. Oncogenic dependencies were identified between mutations in driver genes, common regions of copy number change, and primary translocation and hyperdiploidy events. These dependencies included associations with t(4;14) and mutations in, , and; t(11;14) with mutations in and; t(14;16) with mutations in ,, , and; and hyperdiploidy with gain 11q, mutations in , and rearrangements. These associations indicate that the genomic landscape of myeloma is predetermined by the primary events upon which further dependencies are built, giving rise to a nonrandom accumulation of genetic hits. Understanding these dependencies may elucidate potential evolutionary patterns and lead to better treatment regimens.