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.     

Exome sequencing supports a de novo mutational paradigm for schizophrenia

Despite its high heritability, a large fraction of individuals with schizophrenia do not have a family history of the disease (sporadic cases). Here we examined the possibility that rare de novo protein-altering mutations contribute to the genetic component of schizophrenia by sequencing the exomes of 53 sporadic cases, 22 unaffected controls and their parents. We identified 40 de novo mutations in 27 cases affecting 40 genes, including a potentially disruptive mutation in DGCR2, a gene located in the schizophrenia-predisposing 22q11.2 microdeletion region. A comparison to rare inherited variants indicated that the identified de novo mutations show a large excess of non-synonymous changes in schizophrenia cases, as well as a greater potential to affect protein structure and function. Our analyses suggest a major role for de novo mutations in schizophrenia as well as a large mutational target, which together provide a plausible explanation for the high global incidence and persistence of the disease.     

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 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 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 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 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 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.     

Exome sequencing identifies somatic mutations of DNA methyltransferase gene DNMT3A in acute monocytic leukemia

Abnormal epigenetic regulation has been implicated in oncogenesis. We report here the identification of somatic mutations by exome sequencing in acute monocytic leukemia, the M5 subtype of acute myeloid leukemia (AML-M5). We discovered mutations in DNMT3A (encoding DNA methyltransferase 3A) in 23 of 112 (20.5%) cases. The DNMT3A mutants showed reduced enzymatic activity or aberrant affinity to histone H3 in vitro. Notably, there were alterations of DNA methylation patterns and/or gene expression profiles (such as HOXB genes) in samples with DNMT3A mutations as compared with those without such changes. Leukemias with DNMT3A mutations constituted a group of poor prognosis with elderly disease onset and of promonocytic as well as monocytic predominance among AML-M5 individuals. Screening other leukemia subtypes showed Arg882 alterations in 13.6% of acute myelomonocytic leukemia (AML-M4) cases. Our work suggests a contribution of aberrant DNA methyltransferase activity to the pathogenesis of acute monocytic leukemia and provides a useful new biomarker for relevant cases.     

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 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 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.     

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.     

Whole-genome sequencing and comprehensive variant analysis of a Japanese individual using massively parallel sequencing

We report the analysis of a Japanese male using high-throughput sequencing to × 40 coverage. More than 99% of the sequence reads were mapped to the reference human genome. Using a Bayesian decision method, we identified 3,132,608 single nucleotide variations (SNVs). Comparison with six previously reported genomes revealed an excess of singleton nonsense and nonsynonymous SNVs, as well as singleton SNVs in conserved non-coding regions. We also identified 5,319 deletions smaller than 10 kb with high accuracy, in addition to copy number variations and rearrangements. De novo assembly of the unmapped sequence reads generated around 3 Mb of novel sequence, which showed high similarity to non-reference human genomes and the human herpesvirus 4 genome. Our analysis suggests that considerable variation remains undiscovered in the human genome and that whole-genome sequencing is an invaluable tool for obtaining a complete understanding of human genetic variation.     

Whole-genome sequencing of multiple Arabidopsis thaliana populations

The plant Arabidopsis thaliana occurs naturally in many different habitats throughout Eurasia. As a foundation for identifying genetic variation contributing to adaptation to diverse environments, a 1001 Genomes Project to sequence geographically diverse A. thaliana strains has been initiated. Here we present the first phase of this project, based on population-scale sequencing of 80 strains drawn from eight regions throughout the species' native range. We describe the majority of common small-scale polymorphisms as well as many larger insertions and deletions in the A. thaliana pan-genome, their effects on gene function, and the patterns of local and global linkage among these variants. The action of processes other than spontaneous mutation is identified by comparing the spectrum of mutations that have accumulated since A. thaliana diverged from its closest relative 10 million years ago with the spectrum observed in the laboratory. Recent species-wide selective sweeps are rare, and potentially deleterious mutations are more common in marginal populations.