Frequent somatic mutations in MAP3K5 and MAP3K9 in metastatic melanoma identified by exome sequencing

We sequenced eight melanoma exomes to identify new somatic mutations in metastatic melanoma. Focusing on the mitogen-activated protein (MAP) kinase kinase kinase (MAP3K) family, we found that 24% of melanoma cell lines have mutations in the protein-coding regions of either MAP3K5 or MAP3K9. Structural modeling predicted that mutations in the kinase domain may affect the activity and regulation of these protein kinases. The position of the mutations and the loss of heterozygosity of MAP3K5 and MAP3K9 in 85% and 67% of melanoma samples, respectively, together suggest that the mutations are likely to be inactivating. In in vitro kinase assays, MAP3K5 I780F and MAP3K9 W333* variants had reduced kinase activity. Overexpression of MAP3K5 or MAP3K9 mutants in HEK293T cells reduced the phosphorylation of downstream MAP kinases. Attenuation of MAP3K9 function in melanoma cells using siRNA led to increased cell viability after temozolomide treatment, suggesting that decreased MAP3K pathway activity can lead to chemoresistance in melanoma.     

Expression profiling of medulloblastoma: PDGFRA and the RAS/MAPK pathway as therapeutic targets for metastatic disease.

Little is known about the genetic regulation of medulloblastoma dissemination, but metastatic medulloblastoma is highly associated with poor outcome. We obtained expression profiles of 23 primary medulloblastomas clinically designated as either metastatic (M+) or non-metastatic (M0) and identified 85 genes whose expression differed significantly between classes. Using a class prediction algorithm based on these genes and a leave-one-out approach, we assigned sample class to these tumors (M+ or M0) with 72% accuracy and to four additional independent tumors with 100% accuracy. We also assigned the metastatic medulloblastoma cell line Daoy to the metastatic class. Notably, platelet-derived growth factor receptor alpha (PDGFRA) and members of the downstream RAS/mitogen-activated protein kinase (MAPK) signal transduction pathway are upregulated in M+ tumors. Immunohistochemical validation on an independent set of tumors shows significant overexpression of PDGFRA in M+ tumors compared to M0 tumors. Using in vitro assays, we show that platelet-derived growth factor alpha (PDGFA) enhances medulloblastoma migration and increases downstream MAP2K1 (MEK1), MAP2K2 (MEK2), MAPK1 (p42 MAPK) and MAPK3 (p44 MAPK) phosphorylation in a dose-dependent manner. Neutralizing antibodies to PDGFRA blocks MAP2K1, MAP2K2 and MAPK1/3 phosphorylation, whereas U0126, a highly specific inhibitor of MAP2K1 and MAP2K2, also blocks MAPK1/3. Both inhibit migration and prevent PDGFA-stimulated migration. These results provide the first insight into the genetic regulation of medulloblastoma metastasis and are the first to suggest a role for PDGFRA and the RAS/MAPK signaling pathway in medulloblastoma metastasis. Inhibitors of PDGFRA and RAS proteins should therefore be considered for investigation as possible novel therapeutic strategies against medulloblastoma.     

Germline BAP1 mutations predispose to malignant mesothelioma

Because only a small fraction of asbestos-exposed individuals develop malignant mesothelioma, and because mesothelioma clustering is observed in some families, we searched for genetic predisposing factors. We discovered germline mutations in the gene encoding BRCA1 associated protein-1 (BAP1) in two families with a high incidence of mesothelioma, and we observed somatic alterations affecting BAP1 in familial mesotheliomas, indicating biallelic inactivation. In addition to mesothelioma, some BAP1 mutation carriers developed uveal melanoma. We also found germline BAP1 mutations in 2 of 26 sporadic mesotheliomas; both individuals with mutant BAP1 were previously diagnosed with uveal melanoma. We also observed somatic truncating BAP1 mutations and aberrant BAP1 expression in sporadic mesotheliomas without germline mutations. These results identify a BAP1-related cancer syndrome that is characterized by mesothelioma and uveal melanoma. We hypothesize that other cancers may also be involved and that mesothelioma predominates upon asbestos exposure. These findings will help to identify individuals at high risk of mesothelioma who could be targeted for early intervention.     

Haploinsufficiency of ATP1A2 encoding the Na+/K+ pump alpha2 subunit associated with familial hemiplegic migraine type 2

Headache attacks and autonomic dysfunctions characterize migraine, a very common, disabling disorder with a prevalence of 12% in the general population of Western countries. About 20% of individuals affected with migraine experience aura, a visual or sensory-motor neurological dysfunction that usually precedes or accompanies the headache. Although the mode of transmission is controversial, population-based and twin studies have implicated genetic factors, especially in migraine with aura. Familial hemiplegic migraine is a hereditary form of migraine characterized by aura and some hemiparesis. Here we show that mutations in the gene ATP1A2 that encodes the alpha2 subunit of the Na+/K+ pump are associated with familial hemiplegic migraine type 2 (FHM2) linked to chromosome 1q23 (OMIM 602481). Functional data indicate that the putative pathogenetic mechanism is triggered by a loss of function of a single allele of ATP1A2. This is the first report associating mutations of Na+K+ pump subunits to genetic diseases.     

Exon capture analysis of G protein-coupled receptors identifies activating mutations in GRM3 in melanoma

G protein-coupled receptors (GPCRs), the largest human gene family, are important regulators of signaling pathways. However, knowledge of their genetic alterations is limited. In this study, we used exon capture and massively parallel sequencing methods to analyze the mutational status of 734 GPCRs in melanoma. This investigation revealed that one family member, GRM3, was frequently mutated and that one of its mutations clustered within one position. Biochemical analysis of GRM3 alterations revealed that mutant GRM3 selectively regulated the phosphorylation of MEK, leading to increased anchorage-independent growth and migration. Melanoma cells expressing mutant GRM3 had reduced cell growth and cellular migration after short hairpin RNA-mediated knockdown of GRM3 or treatment with a selective MEK inhibitor, AZD-6244, which is currently being used in phase 2 clinical trials. Our study yields the most comprehensive map of genetic alterations in the GPCR gene family.     

RNA interference and inhibition of MEK-ERK signaling prevent abnormal skeletal phenotypes in a mouse model of craniosynostosis

Premature fusion of one or more of the cranial sutures (craniosynostosis) in humans causes over 100 skeletal diseases, which occur in 1 of approximately 2,500 live births. Among them is Apert syndrome, one of the most severe forms of craniosynostosis, primarily caused by missense mutations leading to amino acid changes S252W or P253R in fibroblast growth factor receptor 2 (FGFR2). Here we show that a small hairpin RNA targeting the dominant mutant form of Fgfr2 (Fgfr2(S252W)) completely prevents Apert-like syndrome in mice. Restoration of normal FGFR2 signaling is manifested by an alteration of the activity of extracellular signal-regulated kinases 1 and 2 (ERK1/2), implicating the gene encoding ERK and the genes downstream of it in disease expressivity. Furthermore, treatment of the mutant mice with U0126, an inhibitor of mitogen-activated protein (MAP) kinase kinase 1 and 2 (MEK1/2) that blocks phosphorylation and activation of ERK1/2, significantly inhibits craniosynostosis. These results illustrate a pathogenic role for ERK activation in craniosynostosis resulting from FGFR2 with the S252W substitution and introduce a new concept of small-molecule inhibitor-mediated prevention and therapy for diseases caused by gain-of-function mutations in the human genome.     

Mutations of PKD1 in ADPKD2 cysts suggest a pathogenic effect of trans-heterozygous mutations

Autosomal dominant polycystic kidney disease (ADPKD) is caused by mutations in PKD1 and PKD2. The products of these genes associate to form heteromeric complexes. Several models have been proposed to explain the mechanism of cyst formation. Here we find somatic mutations of PKD2 in 71% of ADPKD2 cysts analysed. Clonal somatic mutations of PKD1 were identified in a subset of cysts that lacked PKD2 mutations.     

Exome sequencing identifies recurrent somatic RAC1 mutations in melanoma

We characterized the mutational landscape of melanoma, the form of skin cancer with the highest mortality rate, by sequencing the exomes of 147 melanomas. Sun-exposed melanomas had markedly more ultraviolet (UV)-like C>T somatic mutations compared to sun-shielded acral, mucosal and uveal melanomas. Among the newly identified cancer genes was PPP6C, encoding a serine/threonine phosphatase, which harbored mutations that clustered in the active site in 12% of sun-exposed melanomas, exclusively in tumors with mutations in BRAF or NRAS. Notably, we identified a recurrent UV-signature, an activating mutation in RAC1 in 9.2% of sun-exposed melanomas. This activating mutation, the third most frequent in our cohort of sun-exposed melanoma after those of BRAF and NRAS, changes Pro29 to serine (RAC1(P29S)) in the highly conserved switch I domain. Crystal structures, and biochemical and functional studies of RAC1(P29S) showed that the alteration releases the conformational restraint conferred by the conserved proline, causes an increased binding of the protein to downstream effectors, and promotes melanocyte proliferation and migration. These findings raise the possibility that pharmacological inhibition of downstream effectors of RAC1 signaling could be of therapeutic benefit.     

Use of human tissue to assess the oncogenic activity of melanoma-associated mutations

Multiple genetic alterations occur in melanoma, a lethal skin malignancy of increasing incidence. These include mutations that activate Ras and two of its effector cascades, Raf and phosphoinositide 3-kinase (PI3K). Induction of Ras and Raf can be caused by active N-Ras and B-Raf mutants as well as by gene amplification. Activation of PI3K pathway components occurs by PTEN loss and by AKT3 amplification. Melanomas also commonly show impairment of the p16(INK4A)-CDK4-Rb and ARF-HDM2-p53 tumor suppressor pathways. CDKN2A mutations can produce p16(INK4A) and ARF protein loss. Rb bypass can also occur through activating CDK4 mutations as well as by CDK4 amplification. In addition to ARF deletion, p53 pathway disruption can result from dominant negative TP53 mutations. TERT amplification also occurs in melanoma. The extent to which these mutations can induce human melanocytic neoplasia is unknown. Here we characterize pathways sufficient to generate human melanocytic neoplasia and show that genetically altered human tissue facilitates functional analysis of mutations observed in human tumors.     

Somatic mosaic IDH1 and IDH2 mutations are associated with enchondroma and spindle cell hemangioma in Ollier disease and Maffucci syndrome

Ollier disease and Maffucci syndrome are non-hereditary skeletal disorders characterized by multiple enchondromas (Ollier disease) combined with spindle cell hemangiomas (Maffucci syndrome). We report somatic heterozygous mutations in IDH1 (c.394C>T encoding an R132C substitution and c.395G>A encoding an R132H substitution) or IDH2 (c.516G>C encoding R172S) in 87% of enchondromas (benign cartilage tumors) and in 70% of spindle cell hemangiomas (benign vascular lesions). In total, 35 of 43 (81%) subjects with Ollier disease and 10 of 13 (77%) with Maffucci syndrome carried IDH1 (98%) or IDH2 (2%) mutations in their tumors. Fourteen of 16 subjects had identical mutations in separate lesions. Immunohistochemistry to detect mutant IDH1 R132H protein suggested intraneoplastic and somatic mosaicism. IDH1 mutations in cartilage tumors were associated with hypermethylation and downregulated expression of several genes. Mutations were also found in 40% of solitary central cartilaginous tumors and in four chondrosarcoma cell lines, which will enable functional studies to assess the role of IDH1 and IDH2 mutations in tumor formation.     

Gain-of-function RAF1 mutations cause Noonan and LEOPARD syndromes with hypertrophic cardiomyopathy

Noonan and LEOPARD syndromes are developmental disorders with overlapping features, including cardiac abnormalities, short stature and facial dysmorphia. Increased RAS signaling owing to PTPN11, SOS1 and KRAS mutations causes approximately 60% of Noonan syndrome cases, and PTPN11 mutations cause 90% of LEOPARD syndrome cases. Here, we report that 18 of 231 individuals with Noonan syndrome without known mutations (corresponding to 3% of all affected individuals) and two of six individuals with LEOPARD syndrome without PTPN11 mutations have missense mutations in RAF1, which encodes a serine-threonine kinase that activates MEK1 and MEK2. Most mutations altered a motif flanking Ser259, a residue critical for autoinhibition of RAF1 through 14-3-3 binding. Of 19 subjects with a RAF1 mutation in two hotspots, 18 (or 95%) showed hypertrophic cardiomyopathy (HCM), compared with the 18% prevalence of HCM among individuals with Noonan syndrome in general. Ectopically expressed RAF1 mutants from the two HCM hotspots had increased kinase activity and enhanced ERK activation, whereas non-HCM-associated mutants were kinase impaired. Our findings further implicate increased RAS signaling in pathological cardiomyocyte hypertrophy.     

Germline mutations in BAP1 predispose to melanocytic tumors

Common acquired melanocytic nevi are benign neoplasms that are composed of small, uniform melanocytes and are typically present as flat or slightly elevated pigmented lesions on the skin. We describe two families with a new autosomal dominant syndrome characterized by multiple, skin-colored, elevated melanocytic tumors. In contrast to common acquired nevi, the melanocytic neoplasms in affected family members ranged histopathologically from epithelioid nevi to atypical melanocytic proliferations that showed overlapping features with melanoma. Some affected individuals developed uveal or cutaneous melanomas. Segregating with this phenotype, we found inactivating germline mutations of BAP1, which encodes a ubiquitin carboxy-terminal hydrolase. The majority of melanocytic neoplasms lost the remaining wild-type allele of BAP1 by various somatic alterations. In addition, we found BAP1 mutations in a subset of sporadic melanocytic neoplasms showing histological similarities to the familial tumors. These findings suggest that loss of BAP1 is associated with a clinically and morphologically distinct type of melanocytic neoplasm.     

HRPT2, encoding parafibromin,is mutated in hyperparathyroidism-jaw tumor syndrome

We report here the identification of a gene associated with the hyperparathyroidism-jaw tumor (HPT-JT) syndrome. A single locus associated with HPT-JT (HRPT2) was previously mapped to chromosomal region 1q25-q32. We refined this region to a critical interval of 12 cM by genotyping in 26 affected kindreds. Using a positional candidate approach, we identified thirteen different heterozygous, germline, inactivating mutations in a single gene in fourteen families with HPT-JT. The proposed role of HRPT2 as a tumor suppressor was supported by mutation screening in 48 parathyroid adenomas with cystic features, which identified three somatic inactivating mutations, all located in exon 1. None of these mutations were detected in normal controls, and all were predicted to cause deficient or impaired protein function. HRPT2 is a ubiquitously expressed, evolutionarily conserved gene encoding a predicted protein of 531 amino acids, for which we propose the name parafibromin. Our findings suggest that HRPT2 is a tumor-suppressor gene, the inactivation of which is directly involved in predisposition to HPT-JT and in development of some sporadic parathyroid tumors.     

Mutations in GRIN2A and GRIN2B encoding regulatory subunits of NMDA receptors cause variable neurodevelopmental phenotypes

N-methyl-D-aspartate (NMDA) receptors mediate excitatory neurotransmission in the mammalian brain. Two glycine-binding NR1 subunits and two glutamate-binding NR2 subunits each form highly Ca2(+)-permeable cation channels which are blocked by extracellular Mg2(+) in a voltage-dependent manner. Either GRIN2B or GRIN2A, encoding the NMDA receptor subunits NR2B and NR2A, was found to be disrupted by chromosome translocation breakpoints in individuals with mental retardation and/or epilepsy. Sequencing of GRIN2B in 468 individuals with mental retardation revealed four de novo mutations: a frameshift, a missense and two splice-site mutations. In another cohort of 127 individuals with idiopathic epilepsy and/or mental retardation, we discovered a GRIN2A nonsense mutation in a three-generation family. In a girl with early-onset epileptic encephalopathy, we identified the de novo GRIN2A mutation c.1845C>A predicting the amino acid substitution p.N615K. Analysis of NR1-NR2A(N615K) (NR2A subunit with the p.N615K alteration) receptor currents revealed a loss of the Mg2(+) block and a decrease in Ca2(+) permeability. Our findings suggest that disturbances in the neuronal electrophysiological balance during development result in variable neurological phenotypes depending on which NR2 subunit of NMDA receptors is affected.     

Whole-genome sequencing of liver cancers identifies etiological influences on mutation patterns and recurrent mutations in chromatin regulators

Hepatocellular carcinoma (HCC) is the third leading cause of cancer-related death worldwide. We sequenced and analyzed the whole genomes of 27 HCCs, 25 of which were associated with hepatitis B or C virus infections, including two sets of multicentric tumors. Although no common somatic mutations were identified in the multicentric tumor pairs, their whole-genome substitution patterns were similar, suggesting that these tumors developed from independent mutations, although their shared etiological backgrounds may have strongly influenced their somatic mutation patterns. Statistical and functional analyses yielded a list of recurrently mutated genes. Multiple chromatin regulators, including ARID1A, ARID1B, ARID2, MLL and MLL3, were mutated in ～50% of the tumors. Hepatitis B virus genome integration in the TERT locus was frequently observed in a high clonal proportion. Our whole-genome sequencing analysis of HCCs identified the influence of etiological background on somatic mutation patterns and subsequent carcinogenesis, as well as recurrent mutations in chromatin regulators in HCCs.     

Genomic sequencing of colorectal adenocarcinomas identifies a recurrent VTI1A-TCF7L2 fusion

Prior studies have identified recurrent oncogenic mutations in colorectal adenocarcinoma and have surveyed exons of protein-coding genes for mutations in 11 affected individuals. Here we report whole-genome sequencing from nine individuals with colorectal cancer, including primary colorectal tumors and matched adjacent non-tumor tissues, at an average of 30.7× and 31.9× coverage, respectively. We identify an average of 75 somatic rearrangements per tumor, including complex networks of translocations between pairs of chromosomes. Eleven rearrangements encode predicted in-frame fusion proteins, including a fusion of VTI1A and TCF7L2 found in 3 out of 97 colorectal cancers. Although TCF7L2 encodes TCF4, which cooperates with β-catenin in colorectal carcinogenesis, the fusion lacks the TCF4 β-catenin-binding domain. We found a colorectal carcinoma cell line harboring the fusion gene to be dependent on VTI1A-TCF7L2 for anchorage-independent growth using RNA interference-mediated knockdown. This study shows previously unidentified levels of genomic rearrangements in colorectal carcinoma that can lead to essential gene fusions and other oncogenic events.     

Recurrent mutations in the U2AF1 splicing factor in myelodysplastic syndromes

Myelodysplastic syndromes (MDS) are hematopoietic stem cell disorders that often progress to chemotherapy-resistant secondary acute myeloid leukemia (sAML). We used whole-genome sequencing to perform an unbiased comprehensive screen to discover the somatic mutations in a sample from an individual with sAML and genotyped the loci containing these mutations in the matched MDS sample. Here we show that a missense mutation affecting the serine at codon 34 (Ser34) in U2AF1 was recurrently present in 13 out of 150 (8.7%) subjects with de novo MDS, and we found suggestive evidence of an increased risk of progression to sAML associated with this mutation. U2AF1 is a U2 auxiliary factor protein that recognizes the AG splice acceptor dinucleotide at the 3' end of introns, and the alterations in U2AF1 are located in highly conserved zinc fingers of this protein. Mutant U2AF1 promotes enhanced splicing and exon skipping in reporter assays in vitro. This previously unidentified, recurrent mutation in U2AF1 implicates altered pre-mRNA splicing as a potential mechanism for MDS pathogenesis.     

Insertional mutagenesis identifies multiple networks of cooperating genes driving intestinal tumorigenesis

The evolution of colorectal cancer suggests the involvement of many genes. To identify new drivers of intestinal cancer, we performed insertional mutagenesis using the Sleeping Beauty transposon system in mice carrying germline or somatic Apc mutations. By analyzing common insertion sites (CISs) isolated from 446 tumors, we identified many hundreds of candidate cancer drivers. Comparison to human data sets suggested that 234 CIS-targeted genes are also dysregulated in human colorectal cancers. In addition, we found 183 CIS-containing genes that are candidate Wnt targets and showed that 20 CISs-containing genes are newly discovered modifiers of canonical Wnt signaling. We also identified mutations associated with a subset of tumors containing an expanded number of Paneth cells, a hallmark of deregulated Wnt signaling, and genes associated with more severe dysplasia included those encoding members of the FGF signaling cascade. Some 70 genes had co-occurrence of CIS pairs, clustering into 38 sub-networks that may regulate tumor development.     

Integrated analysis of somatic mutations and focal copy-number changes identifies key genes and pathways in hepatocellular carcinoma

Hepatocellular carcinoma (HCC) is the most common primary liver malignancy. Here, we performed high-resolution copy-number analysis on 125 HCC tumors and whole-exome sequencing on 24 of these tumors. We identified 135 homozygous deletions and 994 somatic mutations of genes with predicted functional consequences. We found new recurrent alterations in four genes (ARID1A, RPS6KA3, NFE2L2 and IRF2) not previously described in HCC. Functional analyses showed tumor suppressor properties for IRF2, whose inactivation, exclusively found in hepatitis B virus (HBV)-related tumors, led to impaired TP53 function. In contrast, inactivation of chromatin remodelers was frequent and predominant in alcohol-related tumors. Moreover, association of mutations in specific genes (RPS6KA3-AXIN1 and NFE2L2-CTNNB1) suggested that Wnt/β-catenin signaling might cooperate in liver carcinogenesis with both oxidative stress metabolism and Ras/mitogen-activated protein kinase (MAPK) pathways. This study provides insight into the somatic mutational landscape in HCC and identifies interactions between mutations in oncogene and tumor suppressor gene mutations related to specific risk factors.