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.     

DNMT3A mutations in acute myeloid leukemia

New studies reveal that 20% of individuals with acute myeloid leukemia harbor somatic mutations in DNMT3A (encoding DNA methyltransferase 3A). Although these leukemias have some gene expression and DNA methylation changes, a direct link between mutant DNMT3A, epigenetic changes and pathogenesis remains to be established.     

Exome sequencing identifies recurrent mutations of the splicing factor SF3B1 gene in chronic lymphocytic leukemia

Here we perform whole-exome sequencing of samples from 105 individuals with chronic lymphocytic leukemia (CLL), the most frequent leukemia in adults in Western countries. We found 1,246 somatic mutations potentially affecting gene function and identified 78 genes with predicted functional alterations in more than one tumor sample. Among these genes, SF3B1, encoding a subunit of the spliceosomal U2 small nuclear ribonucleoprotein (snRNP), is somatically mutated in 9.7% of affected individuals. Further analysis in 279 individuals with CLL showed that SF3B1 mutations were associated with faster disease progression and poor overall survival. This work provides the first comprehensive catalog of somatic mutations in CLL with relevant clinical correlates and defines a large set of new genes that may drive the development of this common form of leukemia. The results reinforce the idea that targeting several well-known genetic pathways, including mRNA splicing, could be useful in the treatment of CLL and other malignancies.     

Exome sequencing of gastric adenocarcinoma identifies recurrent somatic mutations in cell adhesion and chromatin remodeling genes

Gastric cancer is a major cause of global cancer mortality. We surveyed the spectrum of somatic alterations in gastric cancer by sequencing the exomes of 15 gastric adenocarcinomas and their matched normal DNAs. Frequently mutated genes in the adenocarcinomas included TP53 (11/15 tumors), PIK3CA (3/15) and ARID1A (3/15). Cell adhesion was the most enriched biological pathway among the frequently mutated genes. A prevalence screening confirmed mutations in FAT4, a cadherin family gene, in 5% of gastric cancers (6/110) and FAT4 genomic deletions in 4% (3/83) of gastric tumors. Frequent mutations in chromatin remodeling genes (ARID1A, MLL3 and MLL) also occurred in 47% of the gastric cancers. We detected ARID1A mutations in 8% of tumors (9/110), which were associated with concurrent PIK3CA mutations and microsatellite instability. In functional assays, we observed both FAT4 and ARID1A to exert tumor-suppressor activity. Somatic inactivation of FAT4 and ARID1A may thus be key tumorigenic events in a subset of gastric cancers.     

Exome sequencing identifies recurrent somatic MAP2K1 and MAP2K2 mutations in melanoma

We performed exome sequencing to detect somatic mutations in protein-coding regions in seven melanoma cell lines and donor-matched germline cells. All melanoma samples had high numbers of somatic mutations, which showed the hallmark of UV-induced DNA repair. Such a hallmark was absent in tumor sample-specific mutations in two metastases derived from the same individual. Two melanomas with non-canonical BRAF mutations harbored gain-of-function MAP2K1 and MAP2K2 (MEK1 and MEK2, respectively) mutations, resulting in constitutive ERK phosphorylation and higher resistance to MEK inhibitors. Screening a larger cohort of individuals with melanoma revealed the presence of recurring somatic MAP2K1 and MAP2K2 mutations, which occurred at an overall frequency of 8%. Furthermore, missense and nonsense somatic mutations were frequently found in three candidate melanoma genes, FAT4, LRP1B and DSC1.     

Exome sequencing identifies MAX mutations as a cause of hereditary pheochromocytoma

Hereditary pheochromocytoma (PCC) is often caused by germline mutations in one of nine susceptibility genes described to date, but there are familial cases without mutations in these known genes. We sequenced the exomes of three unrelated individuals with hereditary PCC (cases) and identified mutations in MAX, the MYC associated factor X gene. Absence of MAX protein in the tumors and loss of heterozygosity caused by uniparental disomy supported the involvement of MAX alterations in the disease. A follow-up study of a selected series of 59 cases with PCC identified five additional MAX mutations and suggested an association with malignant outcome and preferential paternal transmission of MAX mutations. The involvement of the MYC-MAX-MXD1 network in the development and progression of neural crest cell tumors is further supported by the lack of functional MAX in rat PCC (PC12) cells and by the amplification of MYCN in neuroblastoma and suggests that loss of MAX function is correlated with metastatic potential.     

Exome sequencing in sporadic autism spectrum disorders identifies severe de novo mutations

Evidence for the etiology of autism spectrum disorders (ASDs) has consistently pointed to a strong genetic component complicated by substantial locus heterogeneity. We sequenced the exomes of 20 individuals with sporadic ASD (cases) and their parents, reasoning that these families would be enriched for de novo mutations of major effect. We identified 21 de novo mutations, 11 of which were protein altering. Protein-altering mutations were significantly enriched for changes at highly conserved residues. We identified potentially causative de novo events in 4 out of 20 probands, particularly among more severely affected individuals, in FOXP1, GRIN2B, SCN1A and LAMC3. In the FOXP1 mutation carrier, we also observed a rare inherited CNTNAP2 missense variant, and we provide functional support for a multi-hit model for disease risk. Our results show that trio-based exome sequencing is a powerful approach for identifying new candidate genes for ASDs and suggest that de novo mutations may contribute substantially to the genetic etiology of ASDs.     

Exome sequencing identifies truncating mutations in PRRT2 that cause paroxysmal kinesigenic dyskinesia

Paroxysmal kinesigenic dyskinesia is the most common type of paroxysmal movement disorder and is often misdiagnosed clinically as epilepsy. Using whole-exome sequencing followed by Sanger sequencing, we identified three truncating mutations within PRRT2 (NM_145239.2) in eight Han Chinese families with histories of paroxysmal kinesigenic dyskinesia: c.514_517delTCTG (p.Ser172Argfs*3) in one family, c.649dupC (p.Arg217Profs*8) in six families and c.972delA (p.Val325Serfs*12) in one family. These truncating mutations co-segregated exactly with the disease in these families and were not observed in 1,000 control subjects of matched ancestry. PRRT2 is a newly discovered gene consisting of four exons encoding the proline-rich transmembrane protein 2, which encompasses 340 amino acids and contains two predicted transmembrane domains. PRRT2 is highly expressed in the developing nervous system, and a truncating mutation alters the subcellular localization of the PRRT2 protein. The function of PRRT2 and its role in paroxysmal kinesigenic dyskinesia should be further investigated.     

Exome sequencing identifies ACAD9 mutations as a cause of complex I deficiency

An isolated defect of respiratory chain complex I activity is a frequent biochemical abnormality in mitochondrial disorders. Despite intensive investigation in recent years, in most instances, the molecular basis underpinning complex I defects remains unknown. We report whole-exome sequencing of a single individual with severe, isolated complex I deficiency. This analysis, followed by filtering with a prioritization of mitochondrial proteins, led us to identify compound heterozygous mutations in ACAD9, which encodes a poorly understood member of the mitochondrial acyl-CoA dehydrogenase protein family. We demonstrated the pathogenic role of the ACAD9 variants by the correction of the complex I defect on expression of the wildtype ACAD9 protein in fibroblasts derived from affected individuals. ACAD9 screening of 120 additional complex I-defective index cases led us to identify two additional unrelated cases and a total of five pathogenic ACAD9 alleles.     

Exome sequencing identifies NMNAT1 mutations as a cause of Leber congenital amaurosis

Leber congenital amaurosis (LCA) is an autosomal recessive retinal dystrophy that manifests with genetic heterogeneity. We sequenced the exome of an individual with LCA and identified nonsense (c.507G>A, p.Trp169*) and missense (c.769G>A, p.Glu257Lys) mutations in NMNAT1, which encodes an enzyme in the nicotinamide adenine dinucleotide (NAD) biosynthesis pathway implicated in protection against axonal degeneration. We also found NMNAT1 mutations in ten other individuals with LCA, all of whom carry the p.Glu257Lys variant.     

Exome sequencing identifies GRIN2A as frequently mutated in melanoma

The incidence of melanoma is increasing more than any other cancer, and knowledge of its genetic alterations is limited. To systematically analyze such alterations, we performed whole-exome sequencing of 14 matched normal and metastatic tumor DNAs. Using stringent criteria, we identified 68 genes that appeared to be somatically mutated at elevated frequency, many of which are not known to be genetically altered in tumors. Most importantly, we discovered that TRRAP harbored a recurrent mutation that clustered in one position (p. Ser722Phe) in 6 out of 167 affected individuals (～4%), as well as a previously unidentified gene, GRIN2A, which was mutated in 33% of melanoma samples. The nature, pattern and functional evaluation of the TRRAP recurrent mutation suggest that TRRAP functions as an oncogene. Our study provides, to our knowledge, the most comprehensive map of genetic alterations in melanoma to date and suggests that the glutamate signaling pathway is involved in this disease.     

Exome sequencing identifies ACSF3 as a cause of combined malonic and methylmalonic aciduria

We used exome sequencing to identify the genetic basis of combined malonic and methylmalonic aciduria (CMAMMA). We sequenced the exome of an individual with CMAMMA and followed up with sequencing of eight additional affected individuals (cases). This included one individual who was identified and diagnosed by searching an exome database. We identify mutations in ACSF3, encoding a putative methylmalonyl-CoA and malonyl-CoA synthetase as a cause of CMAMMA. We also examined a canine model of CMAMMA, which showed pathogenic mutations in a predicted ACSF3 ortholog. ACSF3 mutant alleles occur with a minor allele frequency of 0.0058 in ～1,000 control individuals, predicting a CMAMMA population incidence of ～1:30,000. ACSF3 deficiency is the first human disorder identified as caused by mutations in a gene encoding a member of the acyl-CoA synthetase family, a diverse group of evolutionarily conserved proteins, and may emerge as one of the more common human metabolic disorders.     

Exome sequencing identifies frequent mutation of ARID1A in molecular subtypes of gastric cancer

Gastric cancer is a heterogeneous disease with multiple environmental etiologies and alternative pathways of carcinogenesis. Beyond mutations in TP53, alterations in other genes or pathways account for only small subsets of the disease. We performed exome sequencing of 22 gastric cancer samples and identified previously unreported mutated genes and pathway alterations; in particular, we found genes involved in chromatin modification to be commonly mutated. A downstream validation study confirmed frequent inactivating mutations or protein deficiency of ARID1A, which encodes a member of the SWI-SNF chromatin remodeling family, in 83% of gastric cancers with microsatellite instability (MSI), 73% of those with Epstein-Barr virus (EBV) infection and 11% of those that were not infected with EBV and microsatellite stable (MSS). The mutation spectrum for ARID1A differs between molecular subtypes of gastric cancer, and mutation prevalence is negatively associated with mutations in TP53. Clinically, ARID1A alterations were associated with better prognosis in a stage-independent manner. These results reveal the genomic landscape, and highlight the importance of chromatin remodeling, in the molecular taxonomy of gastric cancer.     

A screen of the complete protein kinase gene family identifies diverse patterns of somatic mutations in human breast cancer

We examined the coding sequence of 518 protein kinases, approximately 1.3 Mb of DNA per sample, in 25 breast cancers. In many tumors, we detected no somatic mutations. But a few had numerous somatic mutations with distinctive patterns indicative of either a mutator phenotype or a past exposure.     

Exome sequencing of extreme phenotypes identifies DCTN4 as a modifier of chronic Pseudomonas aeruginosa infection in cystic fibrosis

Exome sequencing has become a powerful and effective strategy for the discovery of genes underlying Mendelian disorders. However, use of exome sequencing to identify variants associated with complex traits has been more challenging, partly because the sample sizes needed for adequate power may be very large. One strategy to increase efficiency is to sequence individuals who are at both ends of a phenotype distribution (those with extreme phenotypes). Because the frequency of alleles that contribute to the trait are enriched in one or both phenotype extremes, a modest sample size can potentially be used to identify novel candidate genes and/or alleles. As part of the National Heart, Lung, and Blood Institute (NHLBI) Exome Sequencing Project (ESP), we used an extreme phenotype study design to discover that variants in DCTN4, encoding a dynactin protein, are associated with time to first P. aeruginosa airway infection, chronic P. aeruginosa infection and mucoid P. aeruginosa in individuals with cystic fibrosis.     

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.     

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.     

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.     

Somatic mutations in PTPN11 in juvenile myelomonocytic leukemia,myelodysplastic syndromes and acute myeloid leukemia

We report here that individuals with Noonan syndrome and juvenile myelomonocytic leukemia (JMML) have germline mutations in PTPN11 and that somatic mutations in PTPN11 account for 34% of non-syndromic JMML. Furthermore, we found mutations in PTPN11 in a small percentage of individuals with myelodysplastic syndrome (MDS) and de novo acute myeloid leukemia (AML). Functional analyses documented that the two most common mutations in PTPN11 associated with JMML caused a gain of function.