Novel mutations target distinct subgroups of medulloblastoma

Medulloblastoma is a malignant childhood brain tumour comprising four discrete subgroups. Here, to identify mutations that drive medulloblastoma, we sequenced the entire genomes of 37 tumours and matched normal blood. One-hundred and thirty-six genes harbouring somatic mutations in this discovery set were sequenced in an additional 56 medulloblastomas. Recurrent mutations were detected in 41 genes not yet implicated in medulloblastoma; several target distinct components of the epigenetic machinery in different disease subgroups, such as regulators of H3K27 and H3K4 trimethylation in subgroups 3 and 4 (for example, KDM6A and ZMYM3), and CTNNB1-associated chromatin re-modellers in WNT-subgroup tumours (for example, SMARCA4 and CREBBP). Modelling of mutations in mouse lower rhombic lip progenitors that generate WNT-subgroup tumours identified genes that maintain this cell lineage (DDX3X), as well as mutated genes that initiate (CDH1) or cooperate (PIK3CA) in tumorigenesis. These data provide important new insights into the pathogenesis of medulloblastoma subgroups and highlight targets for therapeutic development.     

Subtypes of medulloblastoma have distinct developmental origins

Medulloblastoma encompasses a collection of clinically and molecularly diverse tumour subtypes that together comprise the most common malignant childhood brain tumour. These tumours are thought to arise within the cerebellum, with approximately 25% originating from granule neuron precursor cells (GNPCs) after aberrant activation of the Sonic Hedgehog pathway (hereafter, SHH subtype). The pathological processes that drive heterogeneity among the other medulloblastoma subtypes are not known, hindering the development of much needed new therapies. Here we provide evidence that a discrete subtype of medulloblastoma that contains activating mutations in the WNT pathway effector CTNNB1 (hereafter, WNT subtype) arises outside the cerebellum from cells of the dorsal brainstem. We found that genes marking human WNT-subtype medulloblastomas are more frequently expressed in the lower rhombic lip (LRL) and embryonic dorsal brainstem than in the upper rhombic lip (URL) and developing cerebellum. Magnetic resonance imaging (MRI) and intra-operative reports showed that human WNT-subtype tumours infiltrate the dorsal brainstem, whereas SHH-subtype tumours are located within the cerebellar hemispheres. Activating mutations in Ctnnb1 had little impact on progenitor cell populations in the cerebellum, but caused the abnormal accumulation of cells on the embryonic dorsal brainstem which included aberrantly proliferating Zic1(+) precursor cells. These lesions persisted in all mutant adult mice; moreover, in 15% of cases in which Tp53 was concurrently deleted, they progressed to form medulloblastomas that recapitulated the anatomy and gene expression profiles of human WNT-subtype medulloblastoma. We provide the first evidence, to our knowledge, that subtypes of medulloblastoma have distinct cellular origins. Our data provide an explanation for the marked molecular and clinical differences between SHH- and WNT-subtype medulloblastomas and have profound implications for future research and treatment of this important childhood cancer.     

Dissecting the genomic complexity underlying medulloblastoma

Medulloblastoma is an aggressively growing tumour, arising in the cerebellum or medulla/brain stem. It is the most common malignant brain tumour in children, and shows tremendous biological and clinical heterogeneity. Despite recent treatment advances, approximately 40% of children experience tumour recurrence, and 30% will die from their disease. Those who survive often have a significantly reduced quality of life. Four tumour subgroups with distinct clinical, biological and genetic profiles are currently identified. WNT tumours, showing activated wingless pathway signalling, carry a favourable prognosis under current treatment regimens. SHH tumours show hedgehog pathway activation, and have an intermediate prognosis. Group 3 and 4 tumours are molecularly less well characterized, and also present the greatest clinical challenges. The full repertoire of genetic events driving this distinction, however, remains unclear. Here we describe an integrative deep-sequencing analysis of 125 tumour-normal pairs, conducted as part of the International Cancer Genome Consortium (ICGC) PedBrain Tumor Project. Tetraploidy was identified as a frequent early event in Group 3 and 4 tumours, and a positive correlation between patient age and mutation rate was observed. Several recurrent mutations were identified, both in known medulloblastoma-related genes (CTNNB1, PTCH1, MLL2, SMARCA4) and in genes not previously linked to this tumour (DDX3X, CTDNEP1, KDM6A, TBR1), often in subgroup-specific patterns. RNA sequencing confirmed these alterations, and revealed the expression of what are, to our knowledge, the first medulloblastoma fusion genes identified. Chromatin modifiers were frequently altered across all subgroups. These findings enhance our understanding of the genomic complexity and heterogeneity underlying medulloblastoma, and provide several potential targets for new therapeutics, especially for Group 3 and 4 patients.     

The landscape of cancer genes and mutational processes in breast cancer

All cancers carry somatic mutations in their genomes. A subset, known as driver mutations, confer clonal selective advantage on cancer cells and are causally implicated in oncogenesis, and the remainder are passenger mutations. The driver mutations and mutational processes operative in breast cancer have not yet been comprehensively explored. Here we examine the genomes of 100 tumours for somatic copy number changes and mutations in the coding exons of protein-coding genes. The number of somatic mutations varied markedly between individual tumours. We found strong correlations between mutation number, age at which cancer was diagnosed and cancer histological grade, and observed multiple mutational signatures, including one present in about ten per cent of tumours characterized by numerous mutations of cytosine at TpC dinucleotides. Driver mutations were identified in several new cancer genes including AKT2, ARID1B, CASP8, CDKN1B, MAP3K1, MAP3K13, NCOR1, SMARCD1 and TBX3. Among the 100 tumours, we found driver mutations in at least 40 cancer genes and 73 different combinations of mutated cancer genes. The results highlight the substantial genetic diversity underlying this common disease.     

Somatic mutations affect key pathways in lung adenocarcinoma

Determining the genetic basis of cancer requires comprehensive analyses of large collections of histopathologically well-classified primary tumours. Here we report the results of a collaborative study to discover somatic mutations in 188 human lung adenocarcinomas. DNA sequencing of 623 genes with known or potential relationships to cancer revealed more than 1,000 somatic mutations across the samples. Our analysis identified 26 genes that are mutated at significantly high frequencies and thus are probably involved in carcinogenesis. The frequently mutated genes include tyrosine kinases, among them the EGFR homologue ERBB4; multiple ephrin receptor genes, notably EPHA3; vascular endothelial growth factor receptor KDR; and NTRK genes. These data provide evidence of somatic mutations in primary lung adenocarcinoma for several tumour suppressor genes involved in other cancers--including NF1, APC, RB1 and ATM--and for sequence changes in PTPRD as well as the frequently deleted gene LRP1B. The observed mutational profiles correlate with clinical features, smoking status and DNA repair defects. These results are reinforced by data integration including single nucleotide polymorphism array and gene expression array. Our findings shed further light on several important signalling pathways involved in lung adenocarcinoma, and suggest new molecular targets for treatment.     

The deubiquitinase USP9X suppresses pancreatic ductal adenocarcinoma

Pancreatic ductal adenocarcinoma (PDA) remains a lethal malignancy despite much progress concerning its molecular characterization. PDA tumours harbour four signature somatic mutations in addition to numerous lower frequency genetic events of uncertain significance. Here we use Sleeping Beauty (SB) transposon-mediated insertional mutagenesis in a mouse model of pancreatic ductal preneoplasia to identify genes that cooperate with oncogenic Kras(G12D) to accelerate tumorigenesis and promote progression. Our screen revealed new candidate genes for PDA and confirmed the importance of many genes and pathways previously implicated in human PDA. The most commonly mutated gene was the X-linked deubiquitinase Usp9x, which was inactivated in over 50% of the tumours. Although previous work had attributed a pro-survival role to USP9X in human neoplasia, we found instead that loss of Usp9x enhances transformation and protects pancreatic cancer cells from anoikis. Clinically, low USP9X protein and messenger RNA expression in PDA correlates with poor survival after surgery, and USP9X levels are inversely associated with metastatic burden in advanced disease. Furthermore, chromatin modulation with trichostatin A or 5-aza-2'-deoxycytidine elevates USP9X expression in human PDA cell lines, indicating a clinical approach for certain patients. The conditional deletion of Usp9x cooperated with Kras(G12D) to accelerate pancreatic tumorigenesis in mice, validating their genetic interaction. We propose that USP9X is a major tumour suppressor gene with prognostic and therapeutic relevance in PDA.     

Whole-genome sequencing identifies recurrent mutations in chronic lymphocytic leukaemia

Chronic lymphocytic leukaemia (CLL), the most frequent leukaemia in adults in Western countries, is a heterogeneous disease with variable clinical presentation and evolution. Two major molecular subtypes can be distinguished, characterized respectively by a high or low number of somatic hypermutations in the variable region of immunoglobulin genes. The molecular changes leading to the pathogenesis of the disease are still poorly understood. Here we performed whole-genome sequencing of four cases of CLL and identified 46 somatic mutations that potentially affect gene function. Further analysis of these mutations in 363 patients with CLL identified four genes that are recurrently mutated: notch 1 (NOTCH1), exportin 1 (XPO1), myeloid differentiation primary response gene 88 (MYD88) and kelch-like 6 (KLHL6). Mutations in MYD88 and KLHL6 are predominant in cases of CLL with mutated immunoglobulin genes, whereas NOTCH1 and XPO1 mutations are mainly detected in patients with unmutated immunoglobulins. The patterns of somatic mutation, supported by functional and clinical analyses, strongly indicate that the recurrent NOTCH1, MYD88 and XPO1 mutations are oncogenic changes that contribute to the clinical evolution of the disease. To our knowledge, this is the first comprehensive analysis of CLL combining whole-genome sequencing with clinical characteristics and clinical outcomes. It highlights the usefulness of this approach for the identification of clinically relevant mutations in cancer.     

Whole-genome analysis informs breast cancer response to aromatase inhibition

To correlate the variable clinical features of oestrogen-receptor-positive breast cancer with somatic alterations, we studied pretreatment tumour biopsies accrued from patients in two studies of neoadjuvant aromatase inhibitor therapy by massively parallel sequencing and analysis. Eighteen significantly mutated genes were identified, including five genes (RUNX1, CBFB, MYH9, MLL3 and SF3B1) previously linked to haematopoietic disorders. Mutant MAP3K1 was associated with luminal A status, low-grade histology and low proliferation rates, whereas mutant TP53 was associated with the opposite pattern. Moreover, mutant GATA3 correlated with suppression of proliferation upon aromatase inhibitor treatment. Pathway analysis demonstrated that mutations in MAP2K4, a MAP3K1 substrate, produced similar perturbations as MAP3K1 loss. Distinct phenotypes in oestrogen-receptor-positive breast cancer are associated with specific patterns of somatic mutations that map into cellular pathways linked to tumour biology, but most recurrent mutations are relatively infrequent. Prospective clinical trials based on these findings will require comprehensive genome sequencing.     

Oncogenically active MYD88 mutations in human lymphoma

The activated B-cell-like (ABC) subtype of diffuse large B-cell lymphoma (DLBCL) remains the least curable form of this malignancy despite recent advances in therapy. Constitutive nuclear factor (NF)-κB and JAK kinase signalling promotes malignant cell survival in these lymphomas, but the genetic basis for this signalling is incompletely understood. Here we describe the dependence of ABC DLBCLs on MYD88, an adaptor protein that mediates toll and interleukin (IL)-1 receptor signalling, and the discovery of highly recurrent oncogenic mutations affecting MYD88 in ABC DLBCL tumours. RNA interference screening revealed that MYD88 and the associated kinases IRAK1 and IRAK4 are essential for ABC DLBCL survival. High-throughput RNA resequencing uncovered MYD88 mutations in ABC DLBCL lines. Notably, 29% of ABC DLBCL tumours harboured the same amino acid substitution, L265P, in the MYD88 Toll/IL-1 receptor (TIR) domain at an evolutionarily invariant residue in its hydrophobic core. This mutation was rare or absent in other DLBCL subtypes and Burkitt's lymphoma, but was observed in 9% of mucosa-associated lymphoid tissue lymphomas. At a lower frequency, additional mutations were observed in the MYD88 TIR domain, occurring in both the ABC and germinal centre B-cell-like (GCB) DLBCL subtypes. Survival of ABC DLBCL cells bearing the L265P mutation was sustained by the mutant but not the wild-type MYD88 isoform, demonstrating that L265P is a gain-of-function driver mutation. The L265P mutant promoted cell survival by spontaneously assembling a protein complex containing IRAK1 and IRAK4, leading to IRAK4 kinase activity, IRAK1 phosphorylation, NF-κB signalling, JAK kinase activation of STAT3, and secretion of IL-6, IL-10 and interferon-β. Hence, the MYD88 signalling pathway is integral to the pathogenesis of ABC DLBCL, supporting the development of inhibitors of IRAK4 kinase and other components of this pathway for the treatment of tumours bearing oncogenic MYD88 mutations.     

Clonal evolution in relapsed acute myeloid leukaemia revealed by whole-genome sequencing

Most patients with acute myeloid leukaemia (AML) die from progressive disease after relapse, which is associated with clonal evolution at the cytogenetic level. To determine the mutational spectrum associated with relapse, we sequenced the primary tumour and relapse genomes from eight AML patients, and validated hundreds of somatic mutations using deep sequencing; this allowed us to define clonality and clonal evolution patterns precisely at relapse. In addition to discovering novel, recurrently mutated genes (for example, WAC, SMC3, DIS3, DDX41 and DAXX) in AML, we also found two major clonal evolution patterns during AML relapse: (1) the founding clone in the primary tumour gained mutations and evolved into the relapse clone, or (2) a subclone of the founding clone survived initial therapy, gained additional mutations and expanded at relapse. In all cases, chemotherapy failed to eradicate the founding clone. The comparison of relapse-specific versus primary tumour mutations in all eight cases revealed an increase in transversions, probably due to DNA damage caused by cytotoxic chemotherapy. These data demonstrate that AML relapse is associated with the addition of new mutations and clonal evolution, which is shaped, in part, by the chemotherapy that the patients receive to establish and maintain remissions.     

Subgroup-specific structural variation across 1,000 medulloblastoma genomes

Medulloblastoma, the most common malignant paediatric brain tumour, is currently treated with nonspecific cytotoxic therapies including surgery, whole-brain radiation, and aggressive chemotherapy. As medulloblastoma exhibits marked intertumoural heterogeneity, with at least four distinct molecular variants, previous attempts to identify targets for therapy have been underpowered because of small samples sizes. Here we report somatic copy number aberrations (SCNAs) in 1,087 unique medulloblastomas. SCNAs are common in medulloblastoma, and are predominantly subgroup-enriched. The most common region of focal copy number gain is a tandem duplication of SNCAIP, a gene associated with Parkinson's disease, which is exquisitely restricted to Group 4α. Recurrent translocations of PVT1, including PVT1-MYC and PVT1-NDRG1, that arise through chromothripsis are restricted to Group 3. Numerous targetable SCNAs, including recurrent events targeting TGF-β signalling in Group 3, and NF-κB signalling in Group 4, suggest future avenues for rational, targeted therapy.     

Mutant alpha subunits of Gi2 inhibit cyclic AMP accumulation

One or more of three Gi proteins, Gi1-3, mediates hormonal inhibition of adenylyl cyclase. Whether this inhibition is mediated by the alpha or by the beta gamma subunits of Gi proteins is unclear. Mutations inhibiting the intrinsic GTPase activity of another G protein, the stimulatory regulator of adenylyl cyclase (Gs), constitutively activate it by replacing either of two conserved amino acids in its alpha subunit (alpha s). These mutations create the gsp oncogene which is found in human pituitary and thyroid tumours. In a second group of human endocrine tumours, somatic mutations in the alpha subunit of Gi2 replace a residue cognate to one of those affected by gsp mutations. This implies that the mutations convert the alpha i2 gene into a dominantly acting oncogene, called gip2, and that the mutant alpha i2 subunits are constitutively active. We have therefore assessed cyclic AMP accumulation in cultured cells which stably or transiently express exogenous wild-type alpha i2 complementary DNA or either of two mutant alpha i2 cDNAs. The results show that putatively oncogenic mutations in alpha i2 constitutively activate the protein's ability to inhibit cAMP accumulation.     

The clonal and mutational evolution spectrum of primary triple-negative breast cancers

Primary triple-negative breast cancers (TNBCs), a tumour type defined by lack of oestrogen receptor, progesterone receptor and ERBB2 gene amplification, represent approximately 16% of all breast cancers. Here we show in 104 TNBC cases that at the time of diagnosis these cancers exhibit a wide and continuous spectrum of genomic evolution, with some having only a handful of coding somatic aberrations in a few pathways, whereas others contain hundreds of coding somatic mutations. High-throughput RNA sequencing (RNA-seq) revealed that only approximately 36% of mutations are expressed. Using deep re-sequencing measurements of allelic abundance for 2,414 somatic mutations, we determine for the first time-to our knowledge-in an epithelial tumour subtype, the relative abundance of clonal frequencies among cases representative of the population. We show that TNBCs vary widely in their clonal frequencies at the time of diagnosis, with the basal subtype of TNBC showing more variation than non-basal TNBC. Although p53 (also known as TP53), PIK3CA and PTEN somatic mutations seem to be clonally dominant compared to other genes, in some tumours their clonal frequencies are incompatible with founder status. Mutations in cytoskeletal, cell shape and motility proteins occurred at lower clonal frequencies, suggesting that they occurred later during tumour progression. Taken together, our results show that understanding the biology and therapeutic responses of patients with TNBC will require the determination of individual tumour clonal genotypes.     

Somatic activation of the K-ras oncogene causes early onset lung cancer in mice

About 30% of human tumours carry ras gene mutations. Of the three genes in this family (composed of K-ras, N-ras and H-ras), K-ras is the most frequently mutated member in human tumours, including adenocarcinomas of the pancreas ( approximately 70-90% incidence), colon ( approximately 50%) and lung ( approximately 25-50%). To construct mouse tumour models involving K-ras, we used a new gene targeting procedure to create mouse strains carrying oncogenic alleles of K-ras that can be activated only on a spontaneous recombination event in the whole animal. Here we show that mice carrying these mutations were highly predisposed to a range of tumour types, predominantly early onset lung cancer. This model was further characterized by examining the effects of germline mutations in the tumour suppressor gene p53, which is known to be mutated along with K-ras in human tumours. This approach has several advantages over traditional transgenic strategies, including that it more closely recapitulates spontaneous oncogene activation as seen in human cancers.     

Mutations in the p53 gene occur in diverse human tumour types

The p53 gene has been a constant source of fascination since its discovery nearly a decade ago. Originally considered to be an oncogene, several convergent lines of research have indicated that the wild-type gene product actually functions as a tumour suppressor gene. For example, expression of the neoplastic phenotype is inhibited, rather than promoted, when rat cells are transfected with the murine wild-type p53 gene together with mutant p53 genes and/or other oncogenes. Moreover, in human tumours, the short arm of chromosome 17 is often deleted. In colorectal cancers, the smallest common region of deletion is centred at 17p13.1; this region harbours the p53 gene, and in two tumours examined in detail, the remaining (non-deleted) p53 alleles were found to contain mutations. This result was provocative because allelic deletion coupled with mutation of the remaining allele is a theoretical hallmark of tumour-suppressor genes. In the present report, we have attempted to determine the generality of this observation; that is, whether tumours with allelic deletions of chromosome 17p contain mutant p53 genes in the allele that is retained. Our results suggest that (1) most tumours with such allelic deletions contain p53 point mutations resulting in amino-acid substitutions, (2) such mutations are not confined to tumours with allelic deletion, but also occur in at least some tumours that have retained both parental 17p alleles, and (3) p53 gene mutations are clustered in four 'hot-spots' which exactly coincide with the four most highly conserved regions of the gene. These results suggest that p53 mutations play a role in the development of many common human malignancies.