Identification of Constrained Cancer Driver Genes Based on Mutation Timing

Cancer drivers are genomic alterations that provide cells containing them with a selective advantage over their local competitors, whereas neutral passengers do not change the somatic fitness of cells. Cancer-driving mutations are usually discriminated from passenger mutations by their higher degree of recurrence in tumor samples. However, there is increasing evidence that many additional driver mutations may exist that occur at very low frequencies among tumors. This observation has prompted alternative methods for driver detection, including finding groups of mutually exclusive mutations and incorporating prior biological knowledge about gene function or network structure. Dependencies among drivers due to epistatic interactions can also result in low mutation frequencies, but this effect has been ignored in driver detection so far. Here, we present a new computational approach for identifying genomic alterations that occur at low frequencies because they depend on other events. Unlike passengers, these constrained mutations display punctuated patterns of occurrence in time. We test this driver–passenger discrimination approach based on mutation timing in extensive simulation studies, and we apply it to cross-sectional copy number alteration (CNA) data from ovarian cancer, CNA and single-nucleotide variant (SNV) data from breast tumors and SNV data from colorectal cancer. Among the top ranked predicted drivers, we find low-frequency genes that have already been shown to be involved in carcinogenesis, as well as many new candidate drivers. The mutation timing approach is orthogonal and complementary to existing driver prediction methods. It will help identifying from cancer genome data the alterations that drive tumor progression.     

Cancer‐associated mutations are preferentially distributed in protein kinase functional sites

Protein kinases are a superfamily involved in many crucial cellular processes, including signal transmission and regulation of cell cycle. As a consequence of this role, kinases have been reported to be associated with many types of cancer and are considered as potential therapeutic targets. We analyzed the distribution of pathogenic somatic point mutations (drivers) in the protein kinase superfamily with respect to their location in the protein, such as in structural, evolutionary, and functionally relevant regions. We find these driver mutations are more clearly associated with key protein features than other somatic mutations (passengers) that have not been directly linked to tumor progression. This observation fits well with the expected implication of the alterations in protein kinase function in cancer pathogenicity. To explain the relevance of the detected association of cancer driver mutations at the molecular level in the human kinome, we compare these with genetically inherited mutations (SNPs). We find that the subset of nonsynonymous SNPs that are associated to disease, but sufficiently mild to the point of being widespread in the population, tend to avoid those key protein regions, where they could be more detrimental for protein function. This tendency contrasts with the one detected for cancer associated‐driver‐mutations, which seems to be more directly implicated in the alteration of protein function. The detailed analysis of protein kinase groups and a number of relevant examples, confirm the relation between cancer associated‐driver‐mutations and key regions for protein kinase structure and function. Proteins 2009. © 2009 Wiley‐Liss, Inc.     

Structural and Functional Impact of Cancer-Related Missense Somatic Mutations

A number of large-scale cancer somatic genome sequencing projects are now identifying genetic alterations in cancers. Evaluation of the effects of these mutations is essential for understanding their contribution to tumorigenesis. We have used SNPs3D, a software suite originally developed for analyzing nonsynonymous germ-line variants, to identify single-base mutations with a high impact on protein structure and function. Two machine learning methods are used: one identifying mutations that destabilize protein three-dimensional structure and the other utilizing sequence conservation and detecting all types of effects on in vivo protein function. Incorporation of detailed structure information into the analysis allows detailed interpretation of the functional effects of mutations in specific cases.  Data from a set of breast and colorectal tumors were analyzed. In known cancer genes, mutations approaching 100% of mutations are found to impact protein function, supporting the view that these methods are appropriate for identifying driver mutations. Overall, 50-60% of all somatic missense mutations are predicted to have a high impact on structural stability or to more generally affect the function of the corresponding proteins. This value is similar to the fraction of all possible missense mutations that have a high impact and is much higher than the corresponding one for human population single-nucleotide polymorphisms, at about 30%. The majority of mutations in tumor suppressors destabilize protein structure, while mutations in oncogenes operate in more varied ways, including destabilization of less active conformational states. The set of high-impact mutations encompasses the possible drivers.     

Statistical method on nonrandom clustering with application to somatic mutations in cancer

Background  

Human cancer is caused by the accumulation of tumor-specific mutations in oncogenes and tumor suppressors that confer a selective growth advantage to cells. As a consequence of genomic instability and high levels of proliferation, many passenger mutations that do not contribute to the cancer phenotype arise alongside mutations that drive oncogenesis. While several approaches have been developed to separate driver mutations from passengers, few approaches can specifically identify activating driver mutations in oncogenes, which are more amenable for pharmacological intervention.     

Design and analysis issues in genome-wide somatic mutation studies of cancer

The availability of the human genome sequence and progress in sequencing and bioinformatic technologies have enabled genome-wide investigation of somatic mutations in human cancers. This article briefly reviews challenges arising in the statistical analysis of mutational data of this kind. A first challenge is that of designing studies that efficiently allocate sequencing resources. We show that this can be addressed by two-stage designs and demonstrate via simulations that even relatively small studies can produce lists of candidate cancer genes that are highly informative for future research efforts. A second challenge is to distinguish mutated genes that are selected for by cancer (drivers) from mutated genes that have no role in the development of cancer and simply happened to mutate (passengers). We suggest that this question is best approached as a classification problem and discuss some of the difficulties of more traditional testing-based approaches. A third challenge is to identify biologic processes affected by the driver genes. This can be pursued by gene set analyses. These can reliably identify functional groups and pathways that are enriched for mutated genes even when the individual genes involved in those pathways or sets are not mutated at sufficient frequencies to provide conclusive evidence as drivers.     

Distinguishing driver and passenger mutations in an evolutionary history categorized by interference

In many biological scenarios, from the development of drug resistance in pathogens to the progression of healthy cells toward cancer, quantifying the selection acting on observed mutations is a central question. One difficulty in answering this question is the complexity of the background upon which mutations can arise, with multiple potential interactions between genetic loci. We here present a method for discerning selection from a population history that accounts for interference between mutations. Given sequences sampled from multiple time points in the history of a population, we infer selection at each locus by maximizing a likelihood function derived from a multilocus evolution model. We apply the method to the question of distinguishing between loci where new mutations are under positive selection (drivers) and loci that emit neutral mutations (passengers) in a Wright-Fisher model of evolution. Relative to an otherwise equivalent method in which the genetic background of mutations was ignored, our method inferred selection coefficients more accurately for both driver mutations evolving under clonal interference and passenger mutations reaching fixation in the population through genetic drift or hitchhiking. In a population history recorded by 750 sets of sequences of 100 individuals taken at intervals of 100 generations, a set of 50 loci were divided into drivers and passengers with a mean accuracy of >0.95 across a range of numbers of driver loci. The potential application of our model, either in full or in part, to a range of biological systems, is discussed.     

Cancer driver mutations: predictions and reality

Cancer cells accumulate many genetic alterations throughout their lifetime, but only a few of them drive cancer progression, termed driver mutations. Driver mutations may vary between cancer types and patients, can remain latent for a long time and become drivers at particular cancer stages, or may drive oncogenesis only in conjunction with other mutations. The high mutational, biochemical, and histological tumor heterogeneity makes driver mutation identification very challenging. In this review we summarize recent efforts to identify driver mutations in cancer and annotate their effects. We underline the success of computational methods to predict driver mutations in finding novel cancer biomarkers, including in circulating tumor DNA (ctDNA). We also report on the boundaries of their applicability in clinical research.     

Missense mutations in cancer suppressor gene TP53 are colocalized with exonic splicing enhancers (ESEs)

Mutation databases can be viewed as footprints of functional organization of a gene and thus can be used to infer its functional organization. We studied the association of exonic splicing enhancers (ESEs) with missense mutations in the tumor suppressor gene TP53 using the International Agency for Research on Cancer (IARC) mutation database. The goals of the study were: (i) to verify the hypothesis that deleterious missense mutations are colocalized with ESEs; (ii) to identify potentially functional ESE sites in the open reading frame (ORF) of the TP53. If some sequence functions as a splicing enhancer, then nucleotide substitutions in the site will disturb splicing, abrogate p53 function, and cause an increased susceptibility to cancer. Therefore, among cancers showing p53 mutations, more missense mutations are expected within functional ESE sites as compared to non-functional ESE motifs. Using several statistical tests, we found that missense mutations in TP53 are strongly colocalized with ESEs, and that only a small fraction of ESE sites contributes to the association. There are usually one or two ESEs per exon showing a statistically significant association with missense mutations--so-called significant ESE sites. In many respects significant ESE sites are different from those that do not show association with missense mutations. We found that positions of significant ESE sites are codon-dependent--significant ESEs preferentially start from the first position of a codon, whereas non-significant ESEs show no position dependence. Significant ESEs showed a more limited set of sequences compared to non-significant ESEs. These findings suggest that there is a limited number of missense mutations that influence ESE sites and our analysis provides further insight into the types of sites that harbor exonic enhancer elements.     

Mutation-specific non-canonical pathway of PTEN as a distinct therapeutic target for glioblastoma

PTEN is one of the most frequently altered tumor suppressor genes in malignant tumors. The dominant-negative effect of PTEN alteration suggests that the aberrant function of PTEN mutation might be more disastrous than deletion, the most frequent genomic event in glioblastoma (GBM). This study aimed to understand the functional properties of various PTEN missense mutations and to investigate their clinical relevance. The genomic landscape of PTEN alteration was analyzed using the Samsung Medical Center GBM cohort and validated via The Cancer Genome Atlas dataset. Several hotspot mutations were identified, and their subcellular distributions and phenotypes were evaluated. We established a library of cancer cell lines that overexpress these mutant proteins using the U87MG and patient-derived cell models lacking functional PTEN. PTEN mutations were categorized into two major subsets: missense mutations in the phosphatase domain and truncal mutations in the C2 domain. We determined the subcellular compartmentalization of four mutant proteins (H93Y, C124S, R130Q, and R173C) from the former group and found that they had distinct localizations; those associated with invasive phenotypes (‘edge mutations’) localized to the cell periphery, while the R173C mutant localized to the nucleus. Invasive phenotypes derived from edge substitutions were unaffected by an anti-PI3K/Akt agent but were disrupted by microtubule inhibitors. PTEN mutations exhibit distinct functional properties regarding their subcellular localization. Further, some missense mutations (‘edge mutations’) in the phosphatase domain caused enhanced invasiveness associated with dysfunctional cytoskeletal assembly, thus suggesting it to be a potent therapeutic target.Subject terms: Cancer, Oncogenes     

Mapping the protein–protein and genetic interactions of cancer to guide precision medicine

Massive efforts to sequence cancer genomes have compiled an impressive catalogue of cancer mutations, revealing the recurrent exploitation of a handful of ‘hallmark cancer pathways’. However, unraveling how sets of mutated proteins in these and other pathways hijack pro-proliferative signaling networks and dictate therapeutic responsiveness remains challenging. Here, we show that cancer driver protein–protein interactions are enriched for additional cancer drivers, highlighting the power of physical interaction maps to explain known, as well as uncover new, disease-promoting pathway interrelationships. We hypothesize that by systematically mapping the protein–protein and genetic interactions in cancer—thereby creating Cancer Cell Maps—we will create resources against which to contextualize a patient’s mutations into perturbed pathways/complexes and thereby specify a matching targeted therapeutic cocktail.     

Statistical analysis of pathogenicity of somatic mutations in cancer

Recent large-scale sequencing studies have revealed that cancer genomes contain variable numbers of somatic point mutations distributed across many genes. These somatic mutations most likely include passenger mutations that are not cancer causing and pathogenic driver mutations in cancer genes. Establishing a significant presence of driver mutations in such data sets is of biological interest. Whereas current techniques from phylogeny are applicable to large data sets composed of singly mutated samples, recently exemplified with a p53 mutation database, methods for smaller data sets containing individual samples with multiple mutations need to be developed. By constructing distinct models of both the mutation process and selection pressure upon the cancer samples, exact statistical tests to examine this problem are devised. Tests to examine the significance of selection toward missense, nonsense, and splice site mutations are derived, along with tests assessing variation in selection between functional domains. Maximum-likelihood methods facilitate parameter estimation, including levels of selection pressure and minimum numbers of pathogenic mutations. These methods are illustrated with 25 breast cancers screened across the coding sequences of 518 kinase genes, revealing 90 base substitutions in 71 genes. Significant selection pressure upon truncating mutations was established. Furthermore, an estimated minimum of 29.8 mutations were pathogenic.     

Functional impact bias reveals cancer drivers

Identifying cancer driver genes and pathways among all somatic mutations detected in a cohort of tumors is a key challenge in cancer genomics. Traditionally, this is done by prioritizing genes according to the recurrence of alterations that they bear. However, this approach has some known limitations, such as the difficulty to correctly estimate the background mutation rate, and the fact that it cannot identify lowly recurrently mutated driver genes. Here we present a novel approach, Oncodrive-fm, to detect candidate cancer drivers which does not rely on recurrence. First, we hypothesized that any bias toward the accumulation of variants with high functional impact observed in a gene or group of genes may be an indication of positive selection and can thus be used to detect candidate driver genes or gene modules. Next, we developed a method to measure this bias (FM bias) and applied it to three datasets of tumor somatic variants. As a proof of concept of our hypothesis we show that most of the highly recurrent and well-known cancer genes exhibit a clear FM bias. Moreover, this novel approach avoids some known limitations of recurrence-based approaches, and can successfully identify lowly recurrent candidate cancer drivers.     

The (non)malignancy of cancerous amino acidic substitutions

The process of natural selection acts both on individual organisms within a population and on individual cells within an organism as they develop into cancer. In this work, we have taken a first step toward understanding the differences in selection pressures exerted on the human genome under these disparate circumstances. Focusing on single amino acid substitutions, we have found that cancer‐related mutations (CRMs) are frequent in evolutionarily conserved sites, whereas single amino acid polymorphisms (SAPs) tend to appear in sites having a more relaxed evolutionary pressure. Those CRMs classed as cancer driver mutations show greater enrichment for conserved sites than passenger mutations. Consistent with this, driver mutations are enriched for sites annotated as key functional residues and their neighbors, and are more likely to be located on the surface of proteins than expected by chance. Overall the pattern of CRM and polymorphism is remarkably similar, but we do see a clear signal indicative of diversifying selection for disruptive amino acid substitutions in the cancer driver mutations. The ultimate consequence of the appearance of those mutations must be advantageous for the tumor cell, leading to cell population‐growth and migration events similar to those seen in natural ecosystems. Proteins 2010. © 2009 Wiley‐Liss, Inc.     

Finding Driver Pathways in Cancer: Models and Algorithms

ABSTRACT: BACKGROUND: Cancer sequencing projects are now measuring somatic mutations in large numbers of cancer genomes. A key challenge in interpreting these data is to distinguish driver mutations, mutations important for cancer development, from passenger mutations that have accumulated in somatic cells but without functional consequences. A common approach to identify genes harboring driver mutations is a single gene test that identifies individual genes that are recurrently mutated in a significant number of cancer genomes. However, the power of this test is reduced by: (1) the necessity of estimating the background mutation rate (BMR) for each gene; (2) the mutational heterogeneity in most cancers meaning that groups of genes (e.g. pathways), rather than single genes, are the primary target of mutations. RESULTS: We investigate the problem of discovering driver pathways, groups of genes containing driver mutations, directly from cancer mutation data and without prior knowledge of pathways or other interactions between genes. We introduce two generative models of somatic mutations in cancer and study the algorithmic complexity of discovering driver pathways in both models. We show that a single gene test for driver genes is highly sensitive to the estimate of the BMR. In contrast, we show that an algorithmic approach that maximizes a straightforward measure of the mutational properties of a driver pathway successfully discovers these groups of genes without an estimate of the BMR. Moreover, this approach is also successful in the case when the observed frequencies of passenger and driver mutations are indistinguishable, a situation where single gene tests fail. CONCLUSIONS: Accurate estimation of the BMR is a challenging task. Thus, methods that do not require an estimate of the BMR, such as the ones we provide here, can give increased power for the discovery of driver genes.     

Likelihood models of somatic mutation and codon substitution in cancer genes

The role of somatic mutation in cancer is well established and several genes have been identified that are frequent targets. This has enabled large-scale screening studies of the spectrum of somatic mutations in cancers of particular organs. Cancer gene mutation databases compile the results of many studies and can provide insight into the importance of specific amino acid sequences and functional domains in cancer, as well as elucidate aspects of the mutation process. Past studies of the spectrum of cancer mutations (in particular genes) have examined overall frequencies of mutation (at specific nucleotides) and of missense, nonsense, and silent substitution (at specific codons) both in the sequence as a whole and in a specific functional domain. Existing methods ignore features of the genetic code that allow some codons to mutate to missense, or stop, codons more readily than others (i.e., by one nucleotide change, vs. two or three). A new codon-based method to estimate the relative rate of substitution (fixation of a somatic mutation in a cancer cell lineage) of nonsense vs. missense mutations in different functional domains and in different tumor tissues is presented. Models that account for several potential influences on rates of somatic mutation and substitution in cancer progenitor cells and allow biases of mutation rates for particular dinucleotide sequences (CGs and dipyrimidines), transition vs. transversion bias, and variable rates of silent substitution across functional domains (useful in detecting investigator sampling bias) are considered. Likelihood-ratio tests are used to choose among models, using cancer gene mutation data. The method is applied to analyze published data on the spectrum of p53 mutations in cancers. A novel finding is that the ratio of the probability of nonsense to missense substitution is much lower in the DNA-binding and transactivation domains (ratios near 1) than in structural domains such as the linker, tetramerization (oligomerization), and proline-rich domains (ratios exceeding 100 in some tissues), implying that the specific amino acid sequence may be less critical in structural domains (e.g., amino acid changes less often lead to cancer). The transition vs. transversion bias and effect of CpG dinucleotides on mutation rates in p53 varied greatly across cancers of different organs, likely reflecting effects of different endogenous and exogenous factors influencing mutation in specific organs.     

Toward an integrated map of genetic interactions in cancer cells

Cancer genomes often harbor hundreds of molecular aberrations. Such genetic variants can be drivers or passengers of tumorigenesis and create vulnerabilities for potential therapeutic exploitation. To identify genotype‐dependent vulnerabilities, forward genetic screens in different genetic backgrounds have been conducted. We devised MINGLE, a computational framework to integrate CRISPR/Cas9 screens originating from different libraries building on approaches pioneered for genetic network discovery in model organisms. We applied this method to integrate and analyze data from 85 CRISPR/Cas9 screens in human cancer cells combining functional data with information on genetic variants to explore more than 2.1 million gene‐background relationships. In addition to known dependencies, we identified new genotype‐specific vulnerabilities of cancer cells. Experimental validation of predicted vulnerabilities identified GANAB and PRKCSH as new positive regulators of Wnt/β‐catenin signaling. By clustering genes with similar genetic interaction profiles, we drew the largest genetic network in cancer cells to date. Our scalable approach highlights how diverse genetic screens can be integrated to systematically build informative maps of genetic interactions in cancer, which can grow dynamically as more data are included.     

RAS and RHO family GTPase mutations in cancer: twin sons of different mothers?

Abstract

The RAS and RHO family comprise two major branches of the RAS superfamily of small GTPases. These proteins function as regulated molecular switches and control cytoplasmic signaling networks that regulate a diversity of cellular processes, including cell proliferation and cell migration. In the early 1980s, mutationally activated RAS genes encoding KRAS, HRAS and NRAS were discovered in human cancer and now comprise the most frequently mutated oncogene family in cancer. Only recently, exome sequencing studies identified cancer-associated alterations in two RHO family GTPases, RAC1 and RHOA. RAS and RHO proteins share significant identity in their amino acid sequences, protein structure and biochemistry. Cancer-associated RAS mutant proteins harbor missense mutations that are found primarily at one of three mutational hotspots (G12, G13 and Q61) and have been identified as gain-of-function oncogenic alterations. Although these residues are conserved in RHO family proteins, the gain-of-function mutations found in RAC1 are found primarily at a distinct hotspot. Unexpectedly, the cancer-associated mutations found with RHOA are located at different hotspots than those found with RAS. Furthermore, since the RHOA mutations suggested a loss-of-function phenotype, it has been unclear whether RHOA functions as an oncogene or tumor suppressor in cancer development. Finally, whereas RAS mutations are found in a broad spectrum of cancer types, RHOA and RAC1 mutations occur in a highly restricted range of cancer types. In this review, we focus on RHOA missense mutations found in cancer and their role in driving tumorigenesis, with comparisons to cancer-associated mutations in RAC1 and RAS GTPases.     

Mutant p53: it’s not all one and the same

Mutation of the TP53 tumor suppressor gene is the most common genetic alteration in cancer, and almost 1000 alleles have been identified in human tumors. While virtually all TP53 mutations are thought to compromise wild type p53 activity, the prevalence and recurrence of missense TP53 alleles has motivated countless research studies aimed at understanding the function of the resulting mutant p53 protein. The data from these studies support three distinct, but perhaps not necessarily mutually exclusive, mechanisms for how different p53 mutants impact cancer: first, they lose the ability to execute wild type p53 functions to varying degrees; second, they act as a dominant negative (DN) inhibitor of wild type p53 tumor-suppressive programs; and third, they may gain oncogenic functions that go beyond mere p53 inactivation. Of these possibilities, the gain of function (GOF) hypothesis is the most controversial, in part due to the dizzying array of biological functions that have been attributed to different mutant p53 proteins. Herein we discuss the current state of understanding of TP53 allele variation in cancer and recent reports that both support and challenge the p53 GOF model. In these studies and others, researchers are turning to more systematic approaches to profile TP53 mutations, which may ultimately determine once and for all how different TP53 mutations act as cancer drivers and whether tumors harboring distinct mutations are phenotypically unique. From a clinical perspective, such information could lead to new therapeutic approaches targeting the effects of different TP53 alleles and/or better sub-stratification of patients harboring TP53 mutant cancers.Subject terms: Cancer genetics, Tumour-suppressor proteins