|Home | About | Journals | Submit | Contact Us | Français|
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 67 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.
Approximately 69,000 individuals were diagnosed with invasive melanoma in the United States in 2010, and there were approximately 8,700 deaths from melanoma in that year1. New insights into the pathogenesis of this lethal disease are needed. Although candidate gene analyses have been powerful in identifying melanoma mutations, such those in as BRAF, cKit and ERBB4 (refs. 2–4), and whole-genome sequencing of a single melanoma has been completed5, no comprehensive analysis of this tumor type has yet been performed. We have now used whole-exome sequencing to search for additional candidate melanoma genes in an unbiased fashion.
We conducted an exome resequencing of 14 matched normal and metastatic tumor DNAs from untreated individuals with melanoma. We enriched exonic sequences using Agilent's SureSelect technology for targeted exon capture6, targeting 37 Mb of sequence from exons and their flanking regions in ~20,000 genes. We performed sequencing with the Illumina GAII platform and aligned the reads using ELAND (Illumina, Inc.) followed by cross_match (see URLs) to the reference human genome (build 36.1). On average, we generated 12 Gb of sequence per sample to a mean depth of 180× or greater to achieve exome builds with at least 90% of the exons covered by high quality genotype calls. To eliminate common germline mutations, we removed any potential somatic mutation that was observed in dbSNP130 or in the 1000 Genomes Project data. To determine which of these alterations were somatic (that is, tumor-specific), we compared these data to the matched normal tissue. From these putative alterations, we identified 5,161 potential somatic mutations in 3,568 different genes in the 14 samples sequenced.
To discriminate true mutations from the possible sequence alterations identified, we applied several steps, described in the Supplementary Note and Supplementary Table 1. These refinements gave us a 97.9% coverage rate, a 2.4% false-negative rate and a sensitivity of 81%. Furthermore, these tools removed ~18% of the alterations, leaving 4,226 putative somatic alterations for further scrutiny. A total of 3,871 alterations were heterozygous, and 355 changes were putative alterations in regions of loss of heterozygosity (Supplementary Table 2). Of these alterations, 2,813 caused amino acid changes (non-synonymous), including 2,589 missense and 175 nonsense alterations, and 49 occurred at splice sites. There were 1,387 silent (synonymous) substitutions. We observed a total of 18 small deletions and 9 insertions.
The observed somatic mutations could either be ‘driver’ mutations that have a role in melanoma neoplasia or functionally inert ‘passenger’ changes. This screen yielded a ratio of non-synonymous to synonymous changes (N:S ratio) of 2.0:1, which is not higher than the N:S ratio of 2.5:1 predicted for non-selected passenger mutations7, suggesting that most of these alterations are likely to be ‘passenger’ mutations. The number of C>T mutations was significantly greater than the numbers of other nucleotide substitutions, resulting in a high prevalence of C>T/G>A transitions (P < 0.001) (Supplementary Fig. 1). Finally, we observed a total of 116 dinucleotide substitutions; of these, 59 were CC>TT/GG>AA changes, all consistent with the previously documented ultraviolet light exposure signature8.
To search for new recurrent mutations, we looked for alterations that occurred in 2 or more of the 14 samples subjected to whole-exome sequencing. From this analysis, as expected, the BRAF p.Val600Glu alteration was captured in 7 out of the 14 samples. In addition to BRAF, we found nine more genes harboring a recurrent mutation. Seven of the new recurrent mutations were non-synonymous and two were synonymous (Table 1). Further screening of the new hotspot mutations in an additional 153 melanomas identified DCC and ZNF831 to have a third recurring mutation each. Strikingly, TRRAP, encoding the transformation/transcription domain-associated protein, contained four additional cases that harbored the identical mutation (Fig. 1 and Supplementary Fig. 2), one of which we found in the commercially available melanoma cell line A375. The recurring mutation in TRRAP was a cytosine to thymine change at position 2,165 of the transcript (uc003upp.1), leading to a serine to phenylalanine substitution at amino acid residue 722 of the protein (p.Ser722Phe). The likelihood for the occurrence of six identical mutations is approximately 5 × 10−20, suggesting that the TRRAP hotspot mutation is important.
Subsequent sequence analysis of all coding exons of TRRAP in 25 additional melanomas revealed no additional non-synonymous alterations. The positions of the mutations in TRRAP imply that they are likely to be oncogenic, as we observed no truncating mutations, and 100% of the alterations occurred in one location. In addition, the affected residue is highly conserved evolutionarily (Supplementary Fig. 3). The clustering of somatic missense mutations is similar to that observed for activating mutations, such those in as RAS (ref. 9), BRAF (ref. 2) and PIK3CA (ref. 10), suggesting that TRRAP may be a new oncogene.
TRRAP functions as part of a multiprotein coactivator complex possessing histone acetyltranferase activity that is central to the transcriptional activity of p53, c-Myc and E2F1 (refs. 11,12). To determine whether the TRRAP recurrent alteration (p.Ser722Phe) is transforming, we transiently transfected NIH 3T3 cells with empty vector, wild-type TRRAP, mutant TRRAP or oncogenic Ras with p.Gly12Val. Two weeks after transfection, mutant TRRAP transformed NIH 3T3 cells more efficiently than wild-type TRRAP (Fig. 2a). Levels of TRRAP were comparable between the transfected cells (Fig. 2b).
To assess if melanoma cells harboring endogenous TRRAP mutations are dependent on TRRAP for survival, we used short hairpin RNA (shRNA) to stably knockdown TRRAP protein levels in melanoma cells harboring either wild-type TRRAP (SK-Mel-28 and 17T) or mutant TRRAP (63T and A375). We confirmed specific targeting of TRRAP by immunoblotting of transiently transfected HEK293T cells or one of the targeted melanoma cell lines (Fig. 2c,d). We grew clones in 10% or 2.5% serum to test cell viability. Knockdown of TRRAP had minimal effect on the survival of cells expressing wild-type TRRAP but substantially increased apoptosis rates of melanoma lines carrying mutant TRRAP, particularly when the clones were grown in low serum (Fig. 2e–h and Supplementary Fig. 4). These results were corroborated by protein blot analysis, which showed increased levels of cleaved PARP, specifically in mutant cells knocked down for TRRAP as compared to vector control, as well as to TRRAP wild-type cells (Fig. 2i,j). Thus, mutant TRRAP is essential for melanoma cell survival, which is consistent with previous results showing that TRRAP knockout mice are embryonic lethal13.
To further evaluate the prevalence and spectrum of somatic mutations in genes identified in the whole-exome discovery screen, we selected genes with more mutations than expected using the observed mutation rate (binomial P < 0.05; see Online Methods) (Supplementary Table 3) and those that were mutated in more than two discovery screen samples. We identified 16 such genes and analyzed them for mutations in an additional 38 melanoma samples using PCR amplification and Sanger sequencing. We identified a total of 65 putative changes in the prevalence screen (Table 2).
Notably, we found these sixteen genes, which, except for BRAF, had not previously been identified as playing a role in melanoma, to scale up and harbor an increasing number of somatic mutations in the prevalence screen. This correlated with the initial frequency identified in the discovery screen (Table 2 and Supplementary Table 4). SIFT analysis predicted that a large fraction of the alterations would affect protein function (Supplementary Table 5).
We found GRIN2A, which was mutated in 6 of the 14 melanomas in the discovery screen, to harbor an additional 11 somatic mutations in the prevalence screen. In addition to the discovery and prevalence screen samples, we searched for GRIN2A somatic mutations in two independent validation panel sets from individuals with melanoma. Validation panel set 1, which included 39 melanomas, revealed 11 GRIN2A somatic mutations (28.2%), and validation panel set 2, which included 32 melanomas, identified 5 tumors with mutations (15.6%). Lastly, GRIN2A sequencing in 12 commercially available cell lines revealed a mutation in the 501Mel cell line. In total, we identified 34 distinct GRIN2A mutations in 135 samples (25.2%) (Supplementary Table 4). The number of C>T mutations identified in GRIN2A was significantly greater than the numbers of other nucleotide substitutions, resulting in a high prevalence of C>T/G>A transitions (P < 0.001) (Supplementary Fig. 5). Supplementary Figure 6 depicts the various genetic and analytic stages used in our study.
GRIN2A, found on chromosome 16p13.2, encodes a glutamate (N-methyl-(D)-aspartic acid (NMDA)) receptor subunit ε-1 that is part of the class of ionotropic glutamate receptors and bears the agonist binding site for glutamate14. The location of the discovered GRIN2A alterations is summarized in Figure 3. We observed two mutation clusters surrounding amino acids 371, 372 and 373 in the PBP1 iGluR NMDA NR2 domain and 1,073, 1,074 and 1,076 in the NMDAR2 C domain. In addition, we observed three recurrent alterations (p.Ser278Phe, p.GluE371Lys and p.Glu1175Lys). The affected residues within these clusters are highly conserved evolutionarily, and SIFT analysis predicts that over 51% would affect protein function (Supplementary Table 5). Furthermore, we identified five nonsense mutations.
The nature of somatic mutations in tumors may aid in classifying the identified genes as oncogenes or tumor suppressors15. Generally, oncogenes harbor adjacent missense recurrent mutations that affect one allele. Tumor suppressor genes are generally mutated throughout the gene; many of the mutations truncate the protein and generally affect both alleles. Based on this classification, the mechanism by which mutations in GRIN2A have a tumorigenic effect is unclear; however, their frequency strongly suggests that they play a major role in melanoma and are worthy of future investigation.
The comprehensive approach to acquiring data in this study provided an opportunity to investigate whether any particular new gene pathways have a role in melanoma (Online Methods). From this analysis, we found the glutamate signaling pathway to be highly significant. Glutamate is known to activate two different types of receptors: iono-tropic glutamate receptors (iGluRs) and metabotropic glutamate receptors (mGluRs). iGluRs are ligand-gated ion channels that allow cations such as calcium and potassium to pass through the plasma membrane of the cell after the binding of glutamate to the receptors. iGluRs are subdivided into three receptor types according to agonists response, one of which is N-methyl-d-aspartate (NMDA)16. Thus, GRIN2A, the most highly mutated gene in our screen, encodes a glutamate receptor subunit that binds NMDA. Another highly mutated gene identified was PLCB4, which encodes a downstream protein involved in metabotropic glutamate receptor-related signal transduction causing inositol phosphate production and protein kinase C (PKC) activation17.
Importantly, a recent genetic analysis of the G-protein coupled receptor family in melanoma identified that GRM3, a metabotropic glutamate receptor also activated by glutamate, is mutated in 16% of melanoma cases (Prickett, T.D. et al., unpublished data). In this case as well, the GRM-signaling downstream effector PLCB4 is activated. Finally, a survey of the tyrosine kinome in melanoma has pointed to ERBB4, Pyk2 and the Ephrin receptors to be highly mutated4. Interestingly, ERBB4, its ligand NRG1, Pyk2 and the Ephrins have been shown to play a crucial role in modulation of NMDA receptor signaling18–21. Thus, molecular and genetic studies have implicated crosstalk between NRG1-ERBB4, GRM3, Ephrin signaling and glutamate receptor functions. A link between the glutamate pathway and tumorigenesis has been seen in neuronal tumors, where glioma cells releasing an excess of glutamate showed more aggressive growth than parental glioma cells22. Furthermore, previous reports showing that expression of GRM1 results in melanocytes transformation23,24 also implicated glutamate signaling in melanoma. A glutamate pathway schematic summarizing these various components is depicted in Supplementary Figure 7.
Our study presents, to our knowledge, the most complete analysis of melanoma exome alterations to date, thus allowing us to (i) identify TRRAP as an unexpected target of recurring genetic alterations, (ii) reveal genes that were not previously connected with melanoma, one of which, GRIN2A, is one of the most highly mutated in melanoma to date and (iii) show that a majority of melanoma tumors had alterations in genes encoding members of the glutamate signaling pathway.
The results may enhance our biological insight of melanoma and will hopefully direct us toward new strategies to improve patient care. Some of our findings have potential therapeutic implications, as gluta-mate pathway modification has previously been shown to limit tumor growth25. Further investigation of this pathway in melanoma as well as development of such inhibitors is warranted.
Methods and any associated references are available in the online version of the paper at http://www.nature.com/naturegenetics/.
Exome sequencing identifies GRIN2A as frequently mutated in melanoma
Corresponding author: Yardena Samuels, National Human Genome Research Institute, NIH, 50 South Drive, MSC 8000, Building 50, Room 5140, Bethesda MD 20892-8000, Phone: 301-451-2628, Fax: 301-480-9864, vog.hin.liam@ysleumas
Supplementary Figure 1. Mutation spectra of single base pair substitutions in melanoma whole exome sequencing. The number of each of the six classes of base substitutions resulting in nonsynonymous changes in the whole exome screen is shown.
Supplementary Figure 2. Detection of TRRAP recurrent mutation.
Representative examples of the hotspot mutation in TRRAP. In each case, the top sequence chromatogram was obtained from normal tissue and the lower sequence chromatogram from the indicated tumors. Arrows indicate the location of missense mutations. The nucleotide and amino acid alterations are indicated below the chromatograms.
Supplementary Figure 3. Comparison of conserved Serine-722 of human TRRAP with its orthologs. The human TRRAP orthologs in 12 species were compared using indicated NCBI accession numbers by COBALT algorithm (http://www.ncbi.nlm.nih.gov/tools/cobalt/cobalt.cgi). Conserved serine at amino acid 722 in humans is shown in red and aligned with other species. Nonconserved amino acids are shown in blue.
Supplementary Figure 4. Effect of mutant TRRAP on cell growth.
Apoptosis quantification of melanoma cell lines transduced with shRNA control or shRNAs targeting TRRAP. Cells were grown in growth medium containing 10% serum for the indicated times. Apoptosis was assessed by fluorescence microscopy of Hoechst 33258-stained cells.
Supplementary Figure 5. Mutation spectra of single base pair substitutions in GRIN2A. The number of each of the six classes of base substitutions resulting in nonsynonymous changes in GRIN2A is shown.
Supplementary Figure 6. Whole exome capture and sequencing analysis of 14 melanoma samples. Schematic overview of mutation identification approach of 14 whole exomes. Discovery and Prevalence samples were obtained at The Surgery Branch, National Cancer Institute, Validation set 1 were obtained at University of Colorado Denver School of Medicine and Validation set 2 were obtained at The University of Texas MD Anderson Cancer Center.
Supplementary Figure 7. Glutamate signaling pathway mutations in melanoma. The genes that function in glutamate signaling are specified. Circled genes have somatic mutations. The number of mutations in each gene is indicated by the number adjacent to the circle.
PSD complex, post synaptic density complex. NMDAR, N-methyl-D-aspartic acid (NMDA) receptor. RTK, receptor tyrosine kinase.
Supplementary Tables 1, 2 and 3 are attached individually
Supplementary Table 4. Somatic mutations identified in the Discovery, Prevalence and Validation Screens
Supplementary Table 5. SIFT analysis of somatic mutations identified in highly mutated genes
Supplementary Table 6. Characteristics of melanoma patient samples used in this study
Supplementary Table 7. Primers used for recurrent mutation confirmation
Supplementary Table 8. Primers used for PCR and sequencing of prevalence genes
Supplementary Table 9. Primers used for TRRAP plasmid construction
Bioinformatic and experimental steps used to discriminate between true somatic mutations from the large number of possible sequence alterations identified in the whole exome screen.
We used a previously described Bayesian genotype caller, MPG to call genotypes (available at http://research.nhgri.nih.gov/software/bam2mpg/) 1. Empirical testing against Sanger-based sequencing was employed to set a score cutoff above which the alterations could be trusted. A total of 91 regions were amplified by PCR from the relevant melanoma genomic DNA sample using specific primers and directly sequenced with dye-terminator chemistry. We determined whether a mutation was somatic by examining the sequence of the gene in genomic DNA from matched normal tissue. From the 91 sequences obtained, 47 were confirmed and 44 were not. Looking for differences between matched tissue samples will select for both somatic variants, and sequencing errors, resulting in a higher non-confirmation rate. We further refined the analysis by examining the relationship between the MPG quality score and depth of coverage at each position. The quality score should rise with increasing amounts of sequence covering a given position. When the score does not rise with increasing depth of coverage, this suggests a false genotype call. As indicated in Supplementary Table 1, all the confirmed alterations correlated with a MPG score / depth of coverage ratio of 0.5 or above in both the tumor and matched normal samples, with only one exception. In contrast, 91% of the alterations that were not confirmed by Sanger sequencing, had a MPG / coverage ratio below 0.5 in either the tumor or the normal or both. Thus, of the 91 assessed alterations, using 0.5 as a cutoff, we could identify 46 somatic mutations out of 47 by conventional sequencing, indicating a 97.9% coverage rate and a 2.4% false negative rate. All our future analyses therefore relate to alterations that obtained a MPG / coverage ratio of 0.5 or above in both the tumor and normal samples. Using the 0.5 ratio cutoff, of 47 somatic substitutions that were discovered prior to this study in these samples, 38 were present in our whole exome study, giving us a sensitivity of 81%.
Tumor Tissue Collection (NCI)
A panel of pathology-confirmed metastatic melanoma tumor resections was assembled from patients referred to the Surgery Branch of the National Cancer Institute for enrollment in IRB-approved Surgery Branch clinical trials of immunotherapy. All of the patients had progressive metastatic melanoma, were 16 years or older and passed eligibility testing for the relevant protocol. Before tissue was acquired all patients signed written informed consents. The melanoma tissue panel was selected on the basis of an available pathology-confirmed, melanoma tissue culture line, paired with freshly frozen tissue from the resected melanoma metastasis and apheresiscollected, peripheral blood mononuclear cells.
Tumor Tissue Collection (MDACC)
Tumor samples and clinical data were collected under an IRB-approved protocol for patients undergoing clinically indicated surgical treatment for melanoma at the University of Texas MD Anderson Cancer Center. All of the patients provided informed consent for collection of tissues and clinical data. All specimens were reviewed by a dermatopathologist to confirm viable tumor content prior to use for molecular studies, and for standard clinicopathologic characterization of each tumor. All samples and clinical data were de-identified prior to distribution to investigators for inclusion in the study. Samples were selected for use in the study based on the availability of sufficient DNA to complete the proposed studies.
Tumor Tissue Collection (University of Colorado)
A Panel of DNA samples from melanoma patients was assembled from peripheral blood and metastatic melanoma. All patients with a confirmed diagnosis of melanoma seen in the Melanoma Research Clinic at the University of Colorado Anschutz Medical Campus who were 18 years of age or older were given an opportunity to participate. All patients signed consents for blood and tissue collection for our IRB-approved Skin Cancer Biorepository after meeting with a clinical nurse coordinator. Blood was drawn into PaxGene DNA tubes and processed according to standard protocol. Melanomas were immediately placed on ice and transported to the laboratory where portions of the tumor were selected for DNA extraction. All samples were de-identified and given a unique barcode. The scientific advisory committee of the biorepository reviewed and approved the request for samples. Panels were assembled of peripheral blood DNA and metastatic melanoma DNA from individual patients and stored at -80° C until use.
1. Teer, J.K. et al. Systematic comparison of three genomic enrichment methods for massively parallel DNA sequencing. Genome Res (2010).
We thank V. Maduro, H. Abaan, P. Cruz and J. Mullikin for generating the sequence data analyzed here. We thank V.G. Prieto for pathologic review of the biospecimens from MelCore at MD Anderson. We thank T. Wolfsberg for bioinformatics help and J. Fekecs and D. Leja for graphical assistance. This work was supported by the Intramural Research Programs of the National Human Genome Research Institute, the National Cancer Institute, National Institutes of Health, USA and The University of Texas MD Anderson Cancer Center Melanoma SPORE (P50 CA93459).
Tissue and melanoma cell lines used for the discovery and prevalence screen in this study were described previously26. Tissues used for validation set 1 were fresh frozen melanoma tumors obtained from the University of Colorado Denver Skin Cancer Biorepository, Division of Medical Oncology. Tissue was collected at University of Colorado Hospital, Anschutz Medical Campus, under that institution's Institutional Review Board protocols. DNA was isolated from enriched macrodissected tumor isolates as previously described (see URLs). Tissue processing and storage procedures have been previously described27. Tissues used for validation set 2 of melanomas were obtained from optimum cutting temperature–embedded frozen clinical specimens from the Melanoma Informatics, Tissue Resource, and Pathology Core (MelCore) at The University of Texas MD Anderson Cancer Center under that institution's Institutional Review Board approved protocols. DNA isolation from the tumor-enriched isolates has been described previously28. The clinical information associated with the melanoma tumors used in this study is provided in Supplementary Table 6.
Exome capture was performed using the SureSelect Human All Exon System (Agilent Technologies). The manufacturer's protocol for SureSelect Human All Exon System (Illumina Paired-End Sequencing Library Prep), version 1.0.1 was used, with the following modifications: Bioanalyzer steps were either performed using agarose gel or omitted. In the sample preparation step 9, samples were purified using Ampure XP beads (Agencourt/Beckman Coulter Genomics) according to the manufacturer's protocols. In step 12, samples were purified with the QIAquick MinElute kit (Qiagen Inc.). One column was used for each sample; the four 250-μl post-amplification aliquots were pooled and passed over the column in several spin steps. Samples were eluted in 12 μl buffer EB and quantitated using the Qubit dsDNA BR Assay kit (Invitrogen Corp). In the post-hybridization amplification step 2, samples were purified with AMPure XP beads as described above. Samples were then eluted in 30 μL buffer EB.
Sequencing was performed on the Illumina GAIIx platform with version 4 chemistry and version 4 flowcells according to the manufacturer's instructions. Seventy six base paired-end reads were generated.
Reads were initially aligned using ELAND (Illumina Inc.). ELAND alignments were used to place reads in bins of about five million base pairs. Unmapped reads were placed in the bin of the mate pair if the mate was mapped. Cross_match (P. Green, see URLs) was used to align the reads assigned to each bin to the corresponding ~5 Mb of genomic sequence. Cross_match alignments were converted to the SamTools bam format, and then genotypes were called using bam2mpg29 (see URLs). Bam2mpg was used to implement the most probable Genotype (MPG) algorithm and a Bayesian based method was used to determine the probability of each genotype given the data observed at that position. The quality score represents the difference of the log likelihoods of the most and second most probable genotype. The MPG was divided by the coverage at each position to calculate the MPG:coverage ratio.
To eliminate common germline mutations from consideration, alterations observed in dbSNP130 or in the 1000 Genomes Project 11_2010 data release project (see URLs) were removed. From this list of variants, we included those positions called by at least three of the four analysis methods used by the project. We further limited the list to those variants above 5% minor allele frequency. Polymorphisms were further removed by examination of the sequence of the gene in genomic DNA from matched normal tissue. Genotypes were annotated as described previously30. The types of mutation are: synonymous (changes protein coding region but not the amino acid), non-synonymous (changes protein coding region), missense variant (changes amino acid), nonsense variant, stop (introduces a stop codon), DIV-c (in-frame deletion or insertion variant in a coding region) and DIV-fs (frameshifting deletion or insertion variant in a coding region).
To determine which genes in the discovery screen were likely to be drivers, we calculated two metrics to assess the frequency and probability of mutations. The first method considered the frequency of mutations in the discovery screen, which was calculated by dividing the total number of non-synonymous mutations observed by the total number of base pairs sequenced. The second method fit a binomial distribution based on the number of base pairs sequenced for the gene of interest, the number of mutations observed and the background rate of mutation. For the background rate, we used the observed rate in the exome screen: 11.2 mutations per Mb. With this study design, our power to detect genes mutated in more than 20% of tumors was over 95%; using the two-mutation threshold, the detection probability was 80%.
Genes identified to harbor recurrent mutations were confirmed, and further screened using two primer sets listed in Supplementary Table 7 in an additional 153 melanoma samples. A subset of 16 genes found mutated in the discovery screen was selected for analysis in the prevalence screen using criteria described in the text. These genes were amplified and sequenced in an additional 38 melanoma samples using the primers listed in Supplementary Tables 7 and 8. BRAF was only evaluated for the V600E alteration. Mutational analysis, confirmation and determination of somatic status were carried out as previously described4,26. All genes had 93% coverage or above.
TRRAP expressing construct in the CβSBS vector was a kind gift from M. Cole (Norris Cotton Cancer Center Lebanon, New Hampshire, USA). The TRRAP mutation was generated as previously described26 using primers shown in Supplementary Table 9. The TRRAP mutation bearing the PCR product was inserted in the TRRAP-CβSBS plasmid. Sequence-verified mutant TRRAP was subsequently used in the various experiments. Details are available upon request.
Metastatic melanoma tumor lines were maintained as previously described26. HEK 293T cells were purchased from American Type Culture Collection (ATCC), and SK-Mel-28 and A375 cells were purchased from National Cancer Institute, Division of Cancer Treatment, Developmental Therapeutics Program, Frederick, Maryland, USA. Cell lines were maintained as previously described26.
Transfected cells or stable clones were lysed as described4. Briefly sonicated whole cell extracts in 2×SDS sample buffer were subjected to protein blotting. Primary antibodies used were TRRAP antibody (Cell Signaling, Cat. No. 3966S), PARP antibody (Cell Signaling, Cat. No. 9542) and α-tubulin antibody (Calbiochem-EMD Biosciences, Cat. No. 555627).
Constructs for stable depletion of TRRAP (Catalog# RHS4533, Construct# TRCN0000005361, TRCN0000005362) (Open Biosystems) were confirmed to knockdown TRRAP at the protein level. Lentiviral stocks were prepared as previously described4. Melanoma cell lines (63T, 17T, A375 and SK-Mel-28) were infected with control and TRRAP-specific shRNA lentiviruses. Selection and growth were done as previously described26.
Cells were collected by incubation in trypsin containing EDTA followed by centrifugation and were fixed in a solution containing 3.7% formaldehyde, 0.5% NP40 and 10 μg/ml Hoechst 33258 in PBS. Apoptotic indices were determined by visual scoring of at least 300 nuclei.
For each plasmid (CβSBS, CβSBS -TRRAP (wild type), CβSBS -TRRAP (p.Ser722Phe) and Ras p.Gly12Val), 3 μg of plasmid was transfected using the calcium phosphate precipitation method into NIH 3T3 cells cultured in T25 flasks. Twenty-four hours after transfection, 1–5% of transfected cells were seeded into T25 flasks and cultured in growth medium containing 2.5% FBS for 14–20 days. The flasks were scored under a microscope by visually counting foci.
Three sources of gene sets were selected: Gene Ontology, Kyoto Encyclopedia of Genes and Genomes, and MSigDB (see URLs), containing 10,147, 214 and 1,892 groups, respectively, when downloaded in September 2010. The sets included groups based on molecular function, cellular localization, biological processes and signaling pathways. For each of group of genes, we examined whether we assessed the respective genes in the exome analysis and whether mutations were observed. Several different statistics were then calculated to accentuate different aspects of the analysis. First, we examined the number of genes we successfully sequenced in the study compared to the number that were mutated. This proportion was then compared to the background rate of 3,308 genes mutated in the total of 16,768 genes using a hypergeometric distribution. As the purpose was for ranking, no multiple comparison correction was implemented. Second, we examined the average number of non-synonymous mutations observed for the mutated genes by dividing the total number of non-synonymous mutations into the total number of genes observed to be mutated. Third, we used a binomial calculation. For each group of genes, we determined the total number of mutations observed and the number of base pairs that were successfully sequenced. We then computed the P value as the probability of a group having at least as many mutations as were observed, given the numbers of base pairs sequenced and the background passenger frequencies, using the binomial distribution in R (see URLs). The background passenger frequencies were conservatively estimated as the total numbers of mutations observed the whole exome study divided by the total number of base pairs sequenced in the study (assuming that all of the mutations observed were passengers).
X.W., V.W., J.C.L., J.K.T., T.D.P. and Y.S. designed the study; K.S.-H., M.A.D., J.E.G., W.R., S.R. and S.A.R. collected and analyzed the melanoma samples; X.W., J.K.T., J.G., J.C.L., S.D. and the NISC Comparative Sequencing Program analyzed the genetic data; V.W. and T.D.P. performed and analyzed the functional data. All authors contributed to the final version of the paper.
URLs. Green Group (cross_match), http://www.phrap.org; Laboratory of Thomas Ried, http://www.riedlab.nci.nih.gov; bam2mpg, http://research.nhgri.nih.gov/software/bam2mpg/; 1000 Genomes Project 11_2010 data release project, ftp://ftp-trace.ncbi.nih.gov/1000genomes/ftp/release/20100804/; Gene Ontology, http://www.geneontology.org/; Kyoto Encyclopedia of Genes and Genomes, http://www.genome.jp/kegg; MSigDB, http://www.broadinstitute.org/gsea/msigdb; R, http://www.r-project.org.
Accession codes. Data from this study are deposited in dbSNP. Accession codes: BRAF, CCDS5863.1; CPT1A, CCDS8185.1; DCC, CCDS11952.1; FCRL1, CCDS1170.1; LRRN3, CCDS5754.1; NOS1, CCDS41842.1; PLCH1, CCDS33881.1; SLC17A5, CCDS4981.1; TRRAP, CCDS5659.1; ZNF831, CCDS42894.1; GRIN2A, CCDS10539.1; CCDC63, CCDS9151.1; TMEM132B, CCDS41859.1; PLCB4, CCDS13104.1; AKR1B10, CCDS5832.1; TAS2R60, CCDS5885.1; KHDRBS2, CCDS4963.1; PTPRO, CCDS8675.1; SYT4, CCDS11922.1; UGT2B10, NM_001075.4; SLC6A11, CCDS2602.1; C12orf63, CCDS9062.1; PCDHB8, CCDS4250.1.
Note: Supplementary information is available on the Nature Genetics website.
COMPETING FINANCIAL INTERESTS The authors declare no competing financial interests.