Volume 13 Supplement 4
Domain landscapes of somatic mutations in cancer
© Nehrt et al; licensee BioMed Central Ltd. 2012
Published: 18 June 2012
Large-scale tumor sequencing projects are now underway to identify genetic mutations that drive tumor initiation and development. Most studies take a gene-based approach to identifying driver mutations, highlighting genes mutated in a large percentage of tumor samples as those likely to contain driver mutations. However, this gene-based approach usually does not consider the position of the mutation within the gene or the functional context the position of the mutation provides. Here we introduce a novel method for mapping mutations to distinct protein domains, not just individual genes, in which they occur, thus providing the functional context for how the mutation contributes to disease. Furthermore, aggregating mutations from all genes containing a specific protein domain enables the identification of mutations that are rare at the gene level, but that occur frequently within the specified domain. These highly mutated domains potentially reveal disruptions of protein function necessary for cancer development.
We mapped somatic mutations from the protein coding regions of 100 colon adenocarcinoma tumor samples to the genes and protein domains in which they occurred, and constructed topographical maps to depict the “mutational landscapes” of gene and domain mutation frequencies. We found significant mutation frequency in a number of genes previously known to be somatically mutated in colon cancer patients including APC, TP53 and KRAS. In addition, we found significant mutation frequency within specific domains located in these genes, as well as within other domains contained in genes having low mutation frequencies. These domain “peaks” were enriched with functions important to cancer development including kinase activity, DNA binding and repair, and signal transduction.
Using our method to create the domain landscapes of mutations in colon cancer, we were able to identify somatic mutations with high potential to drive cancer development. Interestingly, the majority of the genes involved have a low mutation frequency. Therefore, themethod shows good potential for identifying rare driver mutations in current, large-scale tumor sequencing projects. In addition, mapping mutations to specific domains provides the necessary functional context for understanding how the mutations contribute to the disease, and may reveal novel or more refined gene and domain target regions for drug development.
The advent of high-throughput, whole-genome DNA sequencing has enabled the evaluation of normal and tumor tissue samples from hundreds of patients in a single study, revealing both germline and somatic mutations with potential involvement in cancer susceptibility, initiation and development. However, distinguishing the handful of somatic mutations expected to initiate and maintain tumor growth, so-called driver mutations, from mutations that play no role in cancer development, passenger mutations, is still a major hurdle to fully understanding the mechanisms of the disease and to the design of more effective treatments. Most current, state-of-the-art studies take a gene-centric approach to the problem [1–6], identifying potential driver mutations as those that occur in genes mutated in a high percentage of the tumor samples. A pathway analysis typically follows to add functional context to the mutated genes. For example, mutations in the APC gene have been shown to be highly prevalent in colorectal tumors [1, 2, 7]. Unfortunately, this approach is limited to a small subset of genes and inherently disregards gene mutations occurring in a low percentage of tumor samples. Identifying rare mutations at the gene level, those that do not recur in the same gene in many patients, with high functional relevance to the oncogenic process is extremely difficult using current gene-centric approaches. This is indeed one of the most crucial problems in the fight against cancer today (http://provocativequestions.nci.nih.gov).
Furthermore, gene-centric approaches to classifying driver and passenger mutations make no distinction between mutations in different sites on the gene, disregarding important information about the functional context of the site of the mutation. A recent study by Vidal’s team demonstrated the potential of gene-centric approaches to mischaracterize mutations . The authors showed that changes causing a complete knockout of a protein (node removal) are often phenotypically distinct from mutations that disrupt specific regions of the proteins thereby eliminating any interaction(s) in which the protein participates (edgetic perturbations). In particular, the authors emphasized the importance of taking into account the modularity of proteins when studying mutation-phenotype relationships, showing several examples where mutations in the same protein but in different protein domains, which are protein regions conserved within and across species , produce distinct disease phenotypes. This result also demonstrates how pathway analyses of mutated genes can potentially provide an incomplete picture of the functional implications of mutations at the gene level. Distinct interactions for a protein in the pathway can either be preserved or disrupted depending on whether or not mutations affect the specific domain mediating the interaction.
In this study, we introduce a new approach for the analysis of cancer somatic mutations based on the study of these mutations at the protein domain level. We argue that since protein domains define the structural and functional units of the proteins, mapping mutations not only to the genes in which they occur, but also to individual protein domains, adds functional information critical to the accurate assessment of the impact of the mutations. Analysis of the positions of individual oncogenic mutations discovered from several independent studies has revealed significant clustering of these mutations at specific positions within the catalytic domain of several protein kinases [1, 11, 12]. Other studies specifically designed to show the clustering of mutations at distinct positions within protein domains have revealed numerous other domain positions highly mutated across a variety of disease types [12, 13]. In a functional analysis of candidate colon cancer genes identified by Sjöblom et al. in 2006, significant enrichments of proteins containing the MH1 and MH2 domains were found . However, as later noted by Chittenden et al., whereas 70% of the mutations in MH2-containing proteins fell within the MH2 domain, the enrichment of the MH1 domain turned out to be misleading as further analysis revealed that none of mutations in the MH1-containing proteins occurred inside of the MH1 domain itself . Similarly, many other distinct domains have been shown to be significantly enriched in cancer-associated genes including kinase domains and domains involved in transcriptional regulation and DNA maintenance and repair . Domain enrichment analysis, however, is commonly performed after a significant set of genes has been identified, and does not consider whether mutations in the genes actually occur inside the enriched domains. This approach can result in misleading assumptions about domain associations to cancer.
Our approach for identifying mutations relevant to cancer development specifically maps somatic mutations to the individual domains in which they occur, resulting in a more accurate measure of enrichment of mutations at the domain level. By mapping mutations to individual domains, our method avoids potentially misleading conclusions from gene-based domain enrichment analyses and provides an inherent functional explanation for how the mutations contribute to disease. We performed an exome-wide study of somatic mutations, including single nucleotide variants (SNVs) and short insertions and deletions (indels), from 100 colon adenocarcinoma patients obtained from The Cancer Genome Atlas (TCGA) project . To implement our approach, mutations in each domain are aggregated from all human proteins containing the domain. As a result we aggregate a wide range of human proteins, from those from the same protein families sharing high sequence similarity to highly dissimilar proteins sharing only the domain in which the mutations are present. By doing so, our domain-centric method can also reveal novel gene candidates for involvement in cancer development, through the identification of a highly mutated domain shared with other genes known to be significant in cancer.
We also introduce new terminology to describe the distribution of somatic mutations at the protein domain level: the “domain mutational landscape” consisting of a topographic representation of mutation frequencies within individual protein domains from whole-genome, cancer sequencing studies, and “domain peaks” defined as protein domains mutated at high-frequency in tumor genomes of the same or different tumor types. Sjöblom and collaborators performed the first large-scale analysis of breast and colorectal cancer mutations . Using a small sample of individual tumors, they identified 191 candidate genes (CAN genes) significantly mutated in breast and colorectal tumor samples. Wood et al. later followed up the study using a larger set of transcripts and revised statistics for identifying significantly mutated genes . The authors identified an additional 89 CAN genes, and postulated that the genomic landscape of cancer is composed of a few commonly mutated gene “mountains”, including APC, KRAS, TP53 and others for colorectal cancer, but is dominated by a larger number of infrequently mutated gene “hills”. Here, we compare the gene-based mutational landscape of a much larger set of colon tumor samples to the landscape revealed by Wood et al., and show significant clustering of mutations in many of the previously identified CAN genes.
In addition, we also describe the domain mutational landscape of colon cancer, and demonstrate how this landscape reveals major properties that cannot be revealed by gene-based landscapes. We show how highly mutated domain peaks can be missed by gene-centric methods when the individual genes containing the domains are not mutated at high frequencies. We also show how focusing only on significantly mutated genes can miss instances where mutations occurring within a shared domain are actually the more relevant functional contributors to the cancer. Therefore, due to its ability to identify and functionally characterize somatic mutations with high potential to drive cancer development, we expect our novel domain-centric method to become an integral tool for the analysis of data from future large-scale cancer sequencing studies.
Mutation counts for colon cancer
Total nonsynonymous SNVs
Total frameshift insertions
Total nonframeshift insertions
Total frameshift deletions
Total nonframeshift deletions
Total stop-loss SNVs
Total stop-gain SNVs
Mutations in domain regions
Average mutations per patient
216 (± 552)
Number of mutations per patient
Here, we show results of mapping mutations not only to individual genes, but also to the specific protein domains in which they occurred. We also constructed the domain-based mutational landscape for colon cancer from a set of 100 tumor samples.
Gene mutation landscapes
Selected genes highly mutated in colon cancer tumor genomes
We adapted the local false discovery rate analysis from Efron et al.  to identify genes and domains with significant mutation frequency. We expected these regions to contain driver mutations under the assumption that non-functional, or passenger, mutations would be uniformly distributed throughout the genome. We also normalized the mutation counts by the representative protein length for genes and by the cumulative domain length for domains, to control for the assumption that longer regions should contain more mutations. This also ensures that the domain landscape is not biased towards more frequently occurring domains. Using a local false discovery rate threshold of 0.1, we identified 154 genes with significant, length normalized mutation frequencies in colon cancer tumor samples (see Additional file 2), and 151 such genes in the breast cancer set (see Additional file 3). The top normalized mutation frequencies occurred in KRAS, OR8U1 and TAS2R43(Table 2), with significant mutation frequency for other genes with well-known relevance to colon cancer including TP53, APC and BRAF . In addition, we identified significant mutation frequencies in six CAN genes previously identified in the Sjöblom and Wood studies: APC, KRAS, TP53, FBXW7, SMAD4 and GRID1. Of the top five ranked CAN genes, four ranked in the top 20 highest mutation frequency genes in our study, and only PIK3CA did not achieve significance despite having nine mutations. However, the PI3K_p85B domain located within the PIK3CA gene was found to have significant mutation frequency. Two other CAN genes, SMAD2 and SMAD3, also did not have significant mutation frequency at gene level, but did have significant mutation frequency within the MH2 domain contained within each gene.
Domain mutation landscapes
As we did for the mutation frequencies for individual genes, we constructed the domain mutational landscape maps by plotting the domain mutation frequencies onto two-dimensional maps where each domain, not gene, is represented as a square of arbitrary size in a grid with coordinates in the x-y axes. The domain peak heights correspond to the mutation frequencies for individual domains containing mutations in the colon (Figure 1B) and breast cancer (Figure 1D) sets, respectively. The count of mutations in each domain was normalized by the cumulative length of all occurrences of the domain from all previously identified representative proteins in the genome. Again using a local false discovery rate threshold of 0.1, we identified 45 domains with significant, length normalized mutation frequencies (see Additional file 4) in the colon cancer set, and 41 such domains in the breast cancer set (see Additional file 5). We found the domain mutation landscape for colon cancer to be dominated by a few peaks corresponding to the P53, APC_crr and CENP-B_N domains, but also to contain a much larger number of smaller domain hills.
Selected domains highly mutated in colon cancer tumor genomes
Cumulative Domain Length
Genes (Number of Mutations)
CTTNBP2NL (1), CTTNBP2 (2), FILIP1 (5), FILIP1L (1)
SMAD4 (9), SMAD9 (2), SMAD2 (1), SMAD3 (1), GARS (1)
KRAS (28), NRAS (7), RAB27B (2), RAB11B (2), RABL3 (2)
+ 46 genes with 36 additional mutations
MSH4 (2), MSH6 (2), MSH5(1), MSH2 (1)
TP53 (27), TP63 (1)
TAS2R43 (36), TAS2R31 (31), TAS2R30 (25), TAS2R19 (21)
+ 19 genes with 27 additional mutations
WFDC8 (3), SLPI (1), WFDC5 (1), KAL1 (1)
Comparison of colon and breast cancer landscapes
GO term enrichment analysis
To determine functions overrepresented in our sets of significant gene and domain peaks in comparison to genes and domains not reaching significance, we first obtained all available Gene Ontology term (GO term) annotations for all genes and domains containing somatic mutations in our colon and breast cancer sets. A subsequent enrichment analysis of GO terms annotated to gene peaks in our colon cancer landscape revealed an overrepresentation of genes annotated with GO terms related to signal transduction, kinase activity, DNA damage response and the regulation of apoptosis. A list of the top overrepresented GO terms from the enrichment analysis of gene peaks for colon and breast cancer can be found in Additional file 1 – Table S4 and Table S5 respectively. A similar analysis of GO terms annotated to colon cancer domain peaks revealed an overrepresentation of domains annotated with GO terms related to DNA binding, DNA repair, enzyme regulator activity and other cancer-related functional terms such as signal transduction and kinase regulator activity. The complete lists of overrepresented GO terms from the enrichment analysis of domain peaks for colon and breast cancer can be found in Additional file 1 Table_S4 and Table_S5, respectively.
Previous studies by Sjöblom et al. and Wood et al. identified significant clustering of mutations in the “genomic landscapes” of human breast and colorectal cancers. Despite the need of larger samples to reach more accurate conclusions [1, 19, 20], these early studies demonstrated the potential of genome-wide studies to capture decades of research into the association of individual genes to cancer in one study. However, due to the rarity of mutations in the gene hills, the authors concluded that these less frequently mutated genes might be better studied within their pathway contexts to elucidate their functional roles in cancer. Today, it is still a major challenge for genome-wide studies of somatic mutations in cancers to identify rare somatic mutations, those gene mutations occurring in a low percentage of tumor samples, that still contribute to cancer initiation and progression.
By mapping mutations not only to the genes, but also to the individual domains they occurred in, we were able to construct the mutational landscapes for both genes and domains for 100 colon cancer patients (Figure 1A and 1B). We also constructed the gene and domain mutational landscapes for 522 breast cancer patients (Figure 1C and 1D) for comparison to another cancer type. Mapping the mutations to specific domains had the advantage of adding the critical functional context necessary for explaining how the mutations potentially contribute to disease. While a relatively small number of significantly mutated domains were shared in both the colon and breast cancer patients, the method also shows the potential of the domain landscape to find commonalities between different cancers at the functional level that might not be apparent at the gene level. Construction of the domain landscape also revealed many properties that are not apparent from traditional gene-based analyses by examining the individual contributions of mutations from distinct genes that fall within a shared domain. These properties include expected instances where a highly mutated gene contained a highly mutated domain, but also unexpected instances where a shared domain is highly mutated, but the individual genes are not, or even where after the removal of mutations from highly mutated genes, some genes still contain mutations within the shared domain. Examination of the domain landscape also revealed instances where all the genes contributed mutations relatively equally to the domain, and where only one or two genes contributed the majority of mutations. We also found instances where highly mutated domains are shared by genes in the same family, and by genes from different families.
Comparison of our gene-based landscape for colon cancer to the landscape constructed by Wood et al. revealed similar topographies: a few highly mutated gene mountains along with a much larger number of still significantly mutated gene hills. There was a relatively small overlap between the 154 genes identified by our study and the 140 CAN genes; only six genes were found to be significant in both studies. As noted, the two studies also used different tumor samples and different statistical models to determine significant mutation frequencies. Yet, despite these differences, four of the top five colorectal CAN genes (APC, KRAS, TP53 and FBXW7) ranked in the top twenty genes with the highest normalized mutation frequency, and the fifth top CAN gene (PIK3CA) was identified to have a significantly mutated domain. We also identified seven genes with significant mutation frequency from the Cancer Gene Census  list known to have somatic mutations in colorectal cancers including the top five CAN genes, NRAS and BRAF. A GO term enrichment analysis of all 154 significantly mutated genes in our study identified enrichment in many biological processes and molecular functions known to be disrupted in cancer development including signal transduction, regulation of apoptosis, regulation of cell proliferation and DNA damage response.
Our analysis of the gene landscape resulted in the re-identification of genes with known cancer association and confirmation on enrichment of genes involved in processes critical to cancer development, which validates that our method can identify significantly mutated genes relevant to cancer, and also provides evidence that the method can be applied to other specified regions within the genome, including domain regions. The main focus of this study, however, was the construction of the domain mutational landscape for colon cancer and its comparison with the gene-based mutational landscape. In total, we identified 45 domains with significant mutation frequency in the colon tumor samples. Again, the landscape was characterized by mountains and hills, similar to that of the gene-based landscape, with the highest peaks in the P53, APC_crr and CENP-B_N domains. The CENP-B_N domain, a known DNA-binding domain , receives mutations from the TIGD7 and JRKL genes. Although TIGD7 and JRKL are both homologs of the Jrk “jerky” gene associated with epilepsy in mice , they do not have known relevance in cancer development. The peaks for P53 and APC_crr were not surprising due to the well-known tumor suppressing functions of the genes containing the domains, TP53 and APC, respectively. However, mapping mutations to the individual domains illustrates the value of our domain-centric method to provide the essential functional context to explain the role the mutations in cancer development. The GO term enrichment analysis for significantly mutated domains confirmed enrichment of significantly mutated domains with functions important to cancer development including kinase activity, DNA binding and repair, and signal transduction.
Our study of the domain landscape of cancer mutations also highlights the relevance of considering the modularity of the proteins when studying somatic mutations. Is the whole protein responsible for the disruption that promotes tumor growth, or are only some of the functional units of the proteins relevant? For instance, the P53 domain, also known as the P53 DNA-binding domain, contains over 90% of the known TP53 mutations , even though the P53 DNA-binding domain covers approximately half of the P53 protein (193 of 393 amino acids). In our study, 27 of the 31 mutations in the P53 protein occurred within the P53 DNA-binding domain. Mutation within the domain has been shown to have multiple detrimental effects including reduced DNA binding affinity, protein misfolding, protein instability and loss of ability to oligomerize (reviewed in ). The APC gene contains seven repeats of the APC_crr domain that bind to the Arm domains of the beta-catenin protein in addition to thirteen other distinct domains . Truncating mutations mainly within the region of the protein containing the second and third repeats of the APC_crr domain, also referred to as the “mutation cluster region”, are known to eliminate APC’s ability to bind and down-regulate beta-catenin, critically impairing its function as a tumor suppressor gene in the Wnt signalling pathway [26, 27].
Together, these examples demonstrate both the advantages and a potential drawback for our domain-based approach. While the traditional, gene centric view of mutation does not consider the location of mutations within the PIK3CA gene, our domain-centric approach captures the functional modularity of protein domains and enables us to reveal specific domains critical to the cancer development process. Our approach also identifies domains with significant mutation frequency that might be missed by approaches based on mutation prevalence, as illustrated by the identification of significant mutation frequency in the PI3K_p85B domain in breast cancer patients. Yet, the power of our approach is derived from aggregating mutations from all genes containing a particular domain, therefore currently restricting our method to identifying significant mutation frequency inside domain regions. More work will be needed to extend the scope of our approach to other regions of the genome.
Comparison of the gene and domain landscapes also enabled us to identify a small number of domains, seven in total, which retained mutations even after the removal of mutations contributed from significantly mutated genes. The WAP domain in particular retained a significant number of mutations aggregated from the WFDC5, SLPI and KAL1 genes even after the removal of mutations from the significantly mutated WFDC8 gene. The WAP, whey acidic protein-type, domain contains four disulfide bonds at its core, characteristic of genes with protease inhibitor activity . WFDC8 has no known association to cancer, however, WFDC5 has been shown to be upregulated in genes undergoing P53 induced apoptosis , and SLPI has been shown to promote malignancy in a lung cancer cell line due to its protease inhibitor function . In addition, mutations in KAL1 are responsible for Kallmann syndrome . Therefore, because of the known cancer and disease relevance of mutations in the WAP domain of other genes, the presence of mutations in the WAP domain of WFDC8 encourage further study of the role of WFDC8 in colon cancer development.
The examples discussed above, in which significant domain peaks correspond to at least one significant gene peak only constitute 14 of the 45 significantly mutated domains from the colon tumor set. The other 31 domain peaks correspond to genes without significant mutation frequencies which are undetected in the gene landscape. Because these domains do not occur in significantly mutated genes, they would likely not be found by traditional, gene-centric studies, but may reveal the disruption of potentially critical functional mechanisms within the cancer tissues. One of these peaks corresponds to the cortactin-binding protein-2 domain, CortBP2, that was mutated in four genes, CTTNBP2NL (1 mutation), CTTNBP2 (2 mutations), FILIP1 (5 mutations), and FILIP1L (1 mutation). Interestingly, FILIP1L is a highly conserved protein known to inhibit proliferation and migration and increase apoptosis in endothelial cells . This anti-angiogenic protein acts as a tumor suppressor and its loss of function has been implicated in ovarian cancer, head and neck squamous cell carcinoma and oligodendrogliomas [38, 39]. While the mutation frequency for the FILIP1L gene was not significant in our study, CortBP2 ranked in the top 75 domains with the highest mutation frequency, suggesting a novel role in colon cancer development for FILIP1L and the other genes containing mutations in the CortBP2 domain. As with any in silico analysis, however, the identification of domains and genes with suspected roles in cancer development can only generate new hypotheses that must ultimately be experimentally validated.
New methods are critically needed to distinguish mutations that drive tumor initiation and development from the millions of variants being identified in current large-scale tumor sequencing projects. Our novel, domain mutational landscape approach for identifying potential driver mutations in significantly mutated domains reveals many properties that traditional gene landscapes cannot reveal while also adding the functional context necessary for understanding how individual mutations contribute to cancer development. We also compared the mutational landscapes for breast and colon cancer, demonstrating the potential for the domain landscape to identify functional similarities among different cancer types. Determining which mutations are most important for tumorigenesis will shed new light on the selective pressures experienced during the process and will ultimately provide a new set of gene and domain targets for drug development.
Cancer mutation datasets
Controlled access, whole-exome mutation data for 100 colon adenocarcinoma patients and 522 breast invasive carcinoma patients were downloaded from the TCGA Data Portal (http://tcga-data.nci.nih.gov/tcga/) using the mutation files from the ucsc.edu_COAD.IlluminaGA_DNASeq.Level_126.96.36.199 and ucsc.edu_COAD.SOLiD_DNASeq.Level_188.8.131.52 directories downloaded on September 13th, 2011 and the genome.wustl.edu_BRCA.IlluminaGA_DNASeq.Level_184.108.40.206 directory downloaded on September 21st, 2011 respectively. The individual patient records listed both somatic and germline SNVs and short insertions and deletions (indels) in addition to the tumor and matched normal tissue genotypes for each mutation. The patient records were filtered to remove germline mutations and mutations that did not pass the quality control filter. A union set of somatic SNVs and indels was created for each cancer type from the individual patient records by identifying all somatic mutations present in at least one patient. The genotypes of all patients were then examined for each mutation in the union set to count the number of times the mutation occurred somatically in the patient population. The ANNOVAR program was used to filter out mutations present in dbSNP (release 130) in order to remove mutations known to be polymorphic, therefore unlikely to be cancer driver mutations . ANNOVAR was also used to annotate the mutations within protein coding regions with their associated effects. SNVs were classified as either causing (nonsynonymous SNVs) or not causing (synonymous SNVs) amino acid changes, or either causing the gain (stop-gain) or loss (stop-loss) of a stop codon. Insertions and deletions were classified to either cause (frameshift) or not cause (nonframeshift) a shift in the reading frame. Synonymous SNVs were filtered out of the union set of mutations as they were assumed to be unlikely to affect cancer development. Somatic mutation counts after annotation and filtering are provided for colon (Table 1) and breast cancer (see Additional file 1 – Table S1).
Protein and protein domain datasets
A set of 33,963 human proteins from the RefSeq database  was downloaded via NCBI’s E-utilities service. Multiple sequence alignments for all Pfam  protein domains were downloaded from the Conserved Domains Database (CDD) , and hidden Markov models (HMMs) were built for the domains using the hmmerbuild tool (HMMER version 2.3.2)  using default parameters with the global option. HMMER's hmmpfam tool was then used to search for complete Pfam domains in all RefSeq proteins using an E-value cutoff of 0.001. Of the 11,912 Pfam domains, 4,265 mapped to at least one Refseq protein.
CAN gene datasets
The sets of breast and colorectal cancer candidate genes were downloaded from the Supporting Online Materials section of the Wood et al. study, tables S4A and S4B.
Mapping mutations to individual proteins and domains
All somatic, protein coding mutatons were mapped to individual proteins using ANNOVAR’s exonic variant annotations via the RefSeq transcript accession numbers corresponding to our RefSeq protein accessions . Mutations were mapped to specific protein domains using the hmmpfam alignment output. The method for mapping mutations to their domain positions was previously described for our Domain Mapping of Disease Mutations (DMDM) database (http://bioinf.umbc.edu/DMDM) . We defined the “cancer gene set” as the set of all genes containing at least one somatic mutation after filtering of synonymous SNVs and known polymorphisms. For each gene in the cancer gene set, the longest RefSeq protein isoform was then identified, defining the “cancer protein set”. Only mutations in the longest isoform were considered to avoid replicating individual mutations in the DNA across multiple proteins. Insertions and deletions were mapped to proteins and domains using the starting position of the mutation.
Calculation of normalized mutation frequency
To determine genes and domains frequently mutated in the patient populations, we first obtained the count of somatic mutations from all patients falling within the protein coding region of each gene, defined in the protein cancer set, and within each protein domain. Because longer regions of DNA are generally expected to accumulate more mutations than shorter regions, we then normalized the gene mutation counts by dividing each count by the length of the gene’s corresponding protein in the protein cancer set. Domain mutation counts were normalized by dividing by the cumulative length of all occurrences of the domain within the protein cancer set.
Calculation of significantly mutated genes and domains
We used aheuristic cut-off of 150 amino acids for the minimum protein or domain length to be included in our analysis. We estimated the null distribution using the “locfdr” package from R and applied these statistics to identify all protein and domains with a local false discovery rate < 0.1.
GO term enrichment analysis
Gene Ontology terms (GO terms)  for all human genes were downloaded from the BioMart portal , and were mapped to RefSeq proteins by their corresponding gene symbols. GO terms for individual Pfam domains were obtained from the “pfam2go” file on the Gene Ontology website. The annotations contained in the pfam2go file are derived from mappings of Pfam domains to InterPro domains , which are manually annotated with GO terms. Many domains still have unknown or unclear function, therefore we were only able to obtain GO term annotations for approximately 40% of the domains containing at least one somatic mutation. To account for differences in the specificity of GO term annotations, each annotated protein and domain was subsequently assigned all GO terms from all possible paths from the root of the ontology (not including the root term itself) to the annotated GO terms using the full GO ontology in OBO v1.2 format. Individual GO terms from the biological process and molecular function ontologies were tested for enrichment in proteins and domains with significant mutation frequencies using the “calculateStatistic” function of the Text::NSP::Measures::2D::Fisher::right Perl module. For each GO term, the counts for the R function corresponded to a 2 x 2 contingency table based on the counts of proteins or domains either assigned the GO term, or not assigned the GO term, and either having a significant or insignificant mutation frequency. Only proteins and domains containing at least one somatic mutation were considered.
List of abbreviations
This work was supported by the National Institutes of Health (NIH) 1K22CA143148 to M.G.K (PI). Nathan Nehrt is funded by the Research Participation Program administered by ORISE through an interagency agreement between the DOE and FDA. The mention of commercial products herein is not to be construed as either an actual or implied endorsement of such products by the Department of Health and Human Services. This is a contribution of the Food and Drug Administration and is not subject to copyright. The results published here are in whole or part based upon data generated by The Cancer Genome Atlas pilot project established by the NCI and NHGRI. Information about TCGA and the investigators and institutions who constitute the TCGA research network can be found at http://cancergenome.nih.gov/.
This article has been published as part of BMC Genomics Volume 13 Supplement 4, 2012: SNP-SIG 2011: Identification and annotation of SNPs in the context of structure, function and disease. The full contents of the supplement are available online at http://www.biomedcentral.com/bmcgenomics/supplements/13/S4.
- Sjoblom T, Jones S, Wood LD, Parsons DW, Lin J, Barber TD, Mandelker D, Leary RJ, Ptak J, Silliman N: The consensus coding sequences of human breast and colorectal cancers. Science. 2006, New York, NY, 314 (5797): 268-274. 10.1126/science.1133427.Google Scholar
- Wood LD, Parsons DW, Jones S, Lin J, Sjoblom T, Leary RJ, Shen D, Boca SM, Barber T, Ptak J: The genomic landscapes of human breast and colorectal cancers. Science. 2007, New York, NY, 318 (5853): 1108-1113. 10.1126/science.1145720.Google Scholar
- Ding L, Getz G, Wheeler DA, Mardis ER, McLellan MD, Cibulskis K, Sougnez C, Greulich H, Muzny DM, Morgan MB, et al: Somatic mutations affect key pathways in lung adenocarcinoma. Nature. 2008, 455 (7216): 1069-1075. 10.1038/nature07423.PubMed CentralView ArticlePubMedGoogle Scholar
- Parsons DW, Li M, Zhang X, Jones S, Leary RJ, Lin JC, Boca SM, Carter H, Samayoa J, Bettegowda C: The genetic landscape of the childhood cancer medulloblastoma. Science. 331 (6016): 435-439.Google Scholar
- Berger MF, Lawrence MS, Demichelis F, Drier Y, Cibulskis K, Sivachenko AY, Sboner A, Esgueva R, Pflueger D, Sougnez C: The genomic complexity of primary human prostate cancer. Nature. 470 (7333): 214-220.Google Scholar
- Jones S, Zhang X, Parsons DW, Lin JC, Leary RJ, Angenendt P, Mankoo P, Carter H, Kamiyama H, Jimeno A, et al: Core signaling pathways in human pancreatic cancers revealed by global genomic analyses. Science. 2008, 321 (5897): 1801-1806. 10.1126/science.1164368.PubMed CentralView ArticlePubMedGoogle Scholar
- Powell SM, Zilz N, Beazer-Barclay Y, Bryan TM, Hamilton SR, Thibodeau SN, Vogelstein B, Kinzler KW: APC mutations occur early during colorectal tumorigenesis. Nature. 1992, 359 (6392): 235-237. 10.1038/359235a0.View ArticlePubMedGoogle Scholar
- Zhong Q, Simonis N, Li QR, Charloteaux B, Heuze F, Klitgord N, Tam S, Yu H, Venkatesan K, Mou D, et al: Edgetic perturbation models of human inherited disorders. Mol Syst Biol. 2009, 5: 321-PubMed CentralView ArticlePubMedGoogle Scholar
- Marchler-Bauer A, Anderson JB, Derbyshire MK, DeWeese-Scott C, Gonzales NR, Gwadz M, Hao L, He S, Hurwitz DI, Jackson JD, et al: CDD: a conserved domain database for interactive domain family analysis. Nucleic acids research. 2007, 35 (Database issue): D237-240.PubMed CentralView ArticlePubMedGoogle Scholar
- Davies H, Hunter C, Smith R, Stephens P, Greenman C, Bignell G, Teague J, Butler A, Edkins S, Stevens C, et al: Somatic mutations of the protein kinase gene family in human lung cancer. Cancer Res. 2005, 65 (17): 7591-7595.PubMedGoogle Scholar
- Dixit A, Yi L, Gowthaman R, Torkamani A, Schork NJ, Verkhivker GM: Sequence and structure signatures of cancer mutation hotspots in protein kinases. PLoS One. 2009, 4 (10): e7485-10.1371/journal.pone.0007485.PubMed CentralView ArticlePubMedGoogle Scholar
- Yue P, Forrest WF, Kaminker JS, Lohr S, Zhang Z, Cavet G: Inferring the functional effects of mutation through clusters of mutations in homologous proteins. Human mutation. 2010, 31 (3): 264-271. 10.1002/humu.21194.View ArticlePubMedGoogle Scholar
- Peterson TA, Adadey A, Santana-Cruz I, Sun Y, Winder A, Kann MG: DMDM: domain mapping of disease mutations. Bioinformatics. 2010, Oxford, England, 26 (19): 2458-2459. 10.1093/bioinformatics/btq447.Google Scholar
- Chittenden TW, Howe EA, Culhane AC, Sultana R, Taylor JM, Holmes C, Quackenbush J: Functional classification analysis of somatically mutated genes in human breast and colorectal cancers. Genomics. 2008, 91 (6): 508-511. 10.1016/j.ygeno.2008.03.002.PubMed CentralView ArticlePubMedGoogle Scholar
- Futreal PA, Coin L, Marshall M, Down T, Hubbard T, Wooster R, Rahman N, Stratton MR: A census of human cancer genes. Nat Rev Cancer. 2004, 4 (3): 177-183. 10.1038/nrc1299.PubMed CentralView ArticlePubMedGoogle Scholar
- Collins FS, Barker AD: Mapping the cancer genome. Pinpointing the genes involved in cancer will help chart a new course across the complex landscape of human malignancies. Sci Am. 2007, 296 (3): 50-57. 10.1038/scientificamerican0307-50.View ArticlePubMedGoogle Scholar
- Finn RD, Tate J, Mistry J, Coggill PC, Sammut SJ, Hotz HR, Ceric G, Forslund K, Eddy SR, Sonnhammer EL, et al: The Pfam protein families database. Nucleic acids research. 2008, 36 (Database issue): D281-288.PubMed CentralPubMedGoogle Scholar
- Forrest WF, Cavet G: Comment on "The consensus coding sequences of human breast and colorectal cancers". Science. 2007, New York, NY, 317 (5844): 1500-author reply 1500Google Scholar
- Getz G, Hofling H, Mesirov JP, Golub TR, Meyerson M, Tibshirani R, Lander ES: Comment on "The consensus coding sequences of human breast and colorectal cancers". Science. 2007, New York, NY, 317 (5844): 1500-Google Scholar
- Rubin AF, Green P: Comment on "The consensus coding sequences of human breast and colorectal cancers". Science. 2007, New York, NY, 317 (5844): 1500-Google Scholar
- Efron B, Tibshirani R: Empirical bayes methods and false discovery rates for microarrays. Genetic epidemiology. 2002, 23 (1): 70-86. 10.1002/gepi.1124.View ArticlePubMedGoogle Scholar
- Iwahara J, Kigawa T, Kitagawa K, Masumoto H, Okazaki T, Yokoyama S: A helix-turn-helix structure unit in human centromere protein B (CENP-B). EMBO J. 1998, 17 (3): 827-837. 10.1093/emboj/17.3.827.PubMed CentralView ArticlePubMedGoogle Scholar
- Dou T, Gu S, Zhou Z, Ji C, Zeng L, Ye X, Xu J, Ying K, Xie Y, Mao Y: Isolation and characterization of a Jerky and JRK/JH8 like gene, tigger transposable element derived 7, TIGD7. Biochem Genet. 2004, 42 (7-8): 279-285.View ArticlePubMedGoogle Scholar
- Petitjean A, Mathe E, Kato S, Ishioka C, Tavtigian SV, Hainaut P, Olivier M: Impact of mutant p53 functional properties on TP53 mutation patterns and tumor phenotype: lessons from recent developments in the IARC TP53 database. Hum Mutat. 2007, 28 (6): 622-629. 10.1002/humu.20495.View ArticlePubMedGoogle Scholar
- Lubin DJ, Butler JS, Loh SN: Folding of tetrameric p53: oligomerization and tumorigenic mutations induce misfolding and loss of function. J Mol Biol. 395 (4): 705-716.Google Scholar
- Munemitsu S, Albert I, Souza B, Rubinfeld B, Polakis P: Regulation of intracellular beta-catenin levels by the adenomatous polyposis coli (APC) tumor-suppressor protein. Proc Natl Acad Sci U S A. 1995, 92 (7): 3046-3050. 10.1073/pnas.92.7.3046.PubMed CentralView ArticlePubMedGoogle Scholar
- Rubinfeld B, Albert I, Porfiri E, Munemitsu S, Polakis P: Loss of beta-catenin regulation by the APC tumor suppressor protein correlates with loss of structure due to common somatic mutations of the gene. Cancer Res. 1997, 57 (20): 4624-4630.PubMedGoogle Scholar
- Tanaka H, Yoshida M, Tanimura H, Fujii T, Sakata K, Tachibana Y, Ohwada J, Ebiike H, Kuramoto S, Morita K: The selective class I PI3K inhibitor CH5132799 targets human cancers harboring oncogenic PIK3CA mutations. Clin Cancer Res. 17 (10): 3272-3281.Google Scholar
- Carpenter CL, Duckworth BC, Auger KR, Cohen B, Schaffhausen BS, Cantley LC: Purification and characterization of phosphoinositide 3-kinase from rat liver. J Biol Chem. 1990, 265 (32): 19704-19711.PubMedGoogle Scholar
- Wymann MP, Pirola L: Structure and function of phosphoinositide 3-kinases. Biochim Biophys Acta. 1998, 1436 (1-2): 127-150. 10.1016/S0005-2760(98)00139-8.View ArticlePubMedGoogle Scholar
- Karakas B, Bachman KE, Park BH: Mutation of the PIK3CA oncogene in human cancers. Br J Cancer. 2006, 94 (4): 455-459. 10.1038/sj.bjc.6602970.PubMed CentralView ArticlePubMedGoogle Scholar
- Pirola L, Zvelebil MJ, Bulgarelli-Leva G, Van Obberghen E, Waterfield MD, Wymann MP: Activation loop sequences confer substrate specificity to phosphoinositide 3-kinase alpha (PI3Kalpha ). Functions of lipid kinase-deficient PI3Kalpha in signaling. J Biol Chem. 2001, 276 (24): 21544-21554. 10.1074/jbc.M011330200.View ArticlePubMedGoogle Scholar
- Simpson KJ, Nicholas KR: The comparative biology of whey proteins. J Mammary Gland Biol Neoplasia. 2002, 7 (3): 313-326. 10.1023/A:1022856801175.View ArticlePubMedGoogle Scholar
- Horikoshi N, Cong J, Kley N, Shenk T: Isolation of differentially expressed cDNAs from p53-dependent apoptotic cells: activation of the human homologue of the Drosophila peroxidasin gene. Biochem Biophys Res Commun. 1999, 261 (3): 864-869. 10.1006/bbrc.1999.1123.View ArticlePubMedGoogle Scholar
- Devoogdt N, Hassanzadeh Ghassabeh G, Zhang J, Brys L, De Baetselier P, Revets H: Secretory leukocyte protease inhibitor promotes the tumorigenic and metastatic potential of cancer cells. Proc Natl Acad Sci U S A. 2003, 100 (10): 5778-5782. 10.1073/pnas.1037154100.PubMed CentralView ArticlePubMedGoogle Scholar
- Hardelin JP, Levilliers J, del Castillo I, Cohen-Salmon M, Legouis R, Blanchard S, Compain S, Bouloux P, Kirk J, Moraine C, et al: X chromosome-linked Kallmann syndrome: stop mutations validate the candidate gene. Proc Natl Acad Sci U S A. 1992, 89 (17): 8190-8194. 10.1073/pnas.89.17.8190.PubMed CentralView ArticlePubMedGoogle Scholar
- Kwon M, Hanna E, Lorang D, He M, Quick JS, Adem A, Stevenson C, Chung JY, Hewitt SM, Zudaire E, et al: Functional characterization of filamin a interacting protein 1-like, a novel candidate for antivascular cancer therapy. Cancer Res. 2008, 68 (18): 7332-7341. 10.1158/0008-5472.CAN-08-1087.PubMed CentralView ArticlePubMedGoogle Scholar
- Burton ER, Gaffar A, Lee SJ, Adeshuko F, Whitney KD, Chung JY, Hewitt SM, Huang GS, Goldberg GL, Libutti SK, et al: Downregulation of Filamin A interacting protein 1-like is associated with promoter methylation and induces an invasive phenotype in ovarian cancer. Molecular cancer research : MCR. 2011, 9 (8): 1126-1138. 10.1158/1541-7786.MCR-11-0162.PubMed CentralView ArticlePubMedGoogle Scholar
- Ducray F, Idbaih A, de Reynies A, Bieche I, Thillet J, Mokhtari K, Lair S, Marie Y, Paris S, Vidaud M, et al: Anaplastic oligodendrogliomas with 1p19q codeletion have a proneural gene expression profile. Molecular cancer. 2008, 7: 41-10.1186/1476-4598-7-41.PubMed CentralView ArticlePubMedGoogle Scholar
- Wang K, Li M, Hakonarson H: ANNOVAR: functional annotation of genetic variants from high-throughput sequencing data. Nucleic Acids Res. 38 (16): e164-Google Scholar
- Pruitt KD, Tatusova T, Maglott DR: NCBI reference sequences (RefSeq): a curated non-redundant sequence database of genomes, transcripts and proteins. Nucleic acids research. 2007, 35 (Database issue): D61-65.PubMed CentralView ArticlePubMedGoogle Scholar
- Eddy SR: Hidden Markov models. Curr Opin Struct Biol. 1996, 6 (3): 361-365. 10.1016/S0959-440X(96)80056-X.View ArticlePubMedGoogle Scholar
- Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, et al: Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nature genetics. 2000, 25 (1): 25-29. 10.1038/75556.PubMed CentralView ArticlePubMedGoogle Scholar
- Kasprzyk A: BioMart: driving a paradigm change in biological data management. Database. 2011, Oxford, bar049-Google Scholar
- Hunter S, Apweiler R, Attwood TK, Bairoch A, Bateman A, Binns D, Bork P, Das U, Daugherty L, Duquenne L, et al: InterPro: the integrative protein signature database. Nucleic Acids Res. 2009, 37 (Database issue): D211-215.PubMed CentralView ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.