Skip to main content
Download PDF

Signs of positive selection of somatic mutations in human cancers detected by EST sequence analysis

BMC Cancer volume 6, Article number: 36 (2006) Cite this article

Abstract

Background

Carcinogenesis typically involves multiple somatic mutations in caretaker (DNA repair) and gatekeeper (tumor suppressors and oncogenes) genes. Analysis of mutation spectra of the tumor suppressor that is most commonly mutated in human cancers, p53, unexpectedly suggested that somatic evolution of the p53 gene during tumorigenesis is dominated by positive selection for gain of function. This conclusion is supported by accumulating experimental evidence of evolution of new functions of p53 in tumors. These findings prompted a genome-wide analysis of possible positive selection during tumor evolution.

Methods

A comprehensive analysis of probable somatic mutations in the sequences of Expressed Sequence Tags (ESTs) from malignant tumors and normal tissues was performed in order to access the prevalence of positive selection in cancer evolution. For each EST, the numbers of synonymous and non-synonymous substitutions were calculated. In order to identify genes with a signature of positive selection in cancers, these numbers were compared to: i) expected numbers and ii) the numbers for the respective genes in the ESTs from normal tissues.

Results

We identified 112 genes with a signature of positive selection in cancers, i.e., a significantly elevated ratio of non-synonymous to synonymous substitutions, in tumors as compared to 37 such genes in an approximately equal-sized EST collection from normal tissues. A substantial fraction of the tumor-specific positive-selection candidates have experimentally demonstrated or strongly predicted links to cancer.

Conclusion

The results of EST analysis should be interpreted with extreme caution given the noise introduced by sequencing errors and undetected polymorphisms. Furthermore, an inherent limitation of EST analysis is that multiple mutations amenable to statistical analysis can be detected only in relatively highly expressed genes. Nevertheless, the present results suggest that positive selection might affect a substantial number of genes during tumorigenic somatic evolution.

Peer Review reports

Background

It is well established that most cancers are triggered by somatic or, less commonly, germline mutations in caretaker and gatekeeper genes [1–6]. The caretakers are broadly defined DNA repair genes that are responsible for maintenance of genome stability. Mutations in the caretaker genes, which are considered to be typical tumor suppressors, compromise genome stability and, more specifically, increase the probability of mutation in the gatekeepers which include both tumor suppressor genes and oncogenes [3, 7]. Tumor suppressors are genes that control cell proliferation, in particular, by causing cell death in response to DNA damage; accordingly, mutational inactivation of tumor suppressors may cause transformation. In contrast, oncogenes are genes that, when mutated, acquire new functions promoting cell proliferation and, eventually, transformation [4].

Since the pioneering work of Theodore Boveri in the beginning of the 20th century[8], tumorigenesis often has been viewed as a somatic version of Darwinian evolution [9–12]. This perspective implies positive selection of mutations that are beneficial from the standpoint of an individual cell, i.e., mutations that promote cell proliferation such as those activating the tumorigenic potential of oncogenes and those inactivating tumor suppressors. In the context of modern evolutionary synthesis, it is equally obvious that tumor evolution should involve substantial purifying selection against mutations impairing proliferation. Although the Darwinian view of tumorigenesis seems to be increasingly gaining foothold, the interplay of selective forces acting on mutations in specific genes is not understood in detail.

Altogether, mutations in more than 200 human genes have been implicated in cancer [13]. Currently, inactivation of tumor suppressors is considered to be the main driving force of tumorigenesis. The most prominent and best studied tumor suppressor is p53, a multifunctional transactivator of transcription and regulator of cell proliferation, programmed cell death, and repair [14–16]. The p53 gene is mutated in nearly 60% of human tumors. Many independent studies have shown that, in addition to its tumor suppressor properties, p53 may also behave as an oncogene[17]. Specifically, gain of new biochemical (e.g., transactivation of transcription of genes that are not affected by wild-type p53) and biological (e.g,, stimulation of cell proliferation) functions resulting from p53 mutations has been demonstrated [18–22]. Compelling evidence of p53 gain-of-function during tumorigenesis has been provided by recent reports on mouse models of Li-Fraumeni syndrome (LFS), a familial cancer predisposition syndrome caused by germline p53 mutations. These studies revealed substantial changes in the tumor spectra of mice carrying common p53 mutations, indicating that gain-of-function by p53 is important for tumorigenesis[23, 24].

The conclusion that gain-of-function in p53 mutants is important for tumorigenesis is strongly supported by the results of bioinformatic analysis of the mutation spectra of the p53 gene [25, 26]. These studies yielded three lines of evidence compatible with biologically relevant gain-of-function in p53 mutants in tumors:

i) somatic mutations of p53 detected in various cancers showed a highly significant excess of non-synonymous over synonymous substitutions, which is the signature of positive selection[27], ii) amino acid replacements caused by cancer-associated mutations clustered within evolutionarily conserved, functionally important regions of p53, and iii) mutational hotspots, the sites of frequent mutation which are subject to particularly strong positive selection, differed depending on the type of tumor, which suggests acquisition of distinct new functions by p53 in different tumors.

These observations prompted us to ask whether positive selection could also be detected in somatic mutants of other cancer-related genes in tumors. Genes evolving under positive selection during cancer progression could be viewed as candidate new oncogenes. To delineate the repertoire of such genes, we performed a genome-wide search for positive selection during cancer evolution by comparing the sequences of Expressed Sequence Tags (EST[28]) from tumors to the corresponding genomic sequences. The rationale of this analysis is to detect somatic mutations in ESTs and identify genes that show a significant excess of non-synonymous over synonymous substitutions in tumors. In principle, EST libraries provide ample material for analyzing somatic mutations in tumors and normal tissues. The problem with this approach is that differences between EST sequences and the sequences of the respective reference genes from the human genome may be caused by a variety of reasons other than somatic mutation including sequencing errors, incorrect assignment of an EST to a reference gene, and single-nucleotide polymorphisms (SNPs).

Several recent, large scale studies employed EST collections for detecting cancer-associated SNPs and cancer-specific alternative splice forms. In particular, Xu and Lee identified 316 human splice variant forms with a statistically significant cancer association; the structures of the most abundant of these were supported by sequences of the corresponding mRNAs isolated from tumors [29]. Another, larger-scale study by Gupta et al. reported 1120 tumor-specific splice isoforms with a high rate of validation by mRNA sequencing. However, when mRNA analysis was performed, the tissue specificity of many of these transcripts, particularly, those of low abundance, could not be confirmed[30]. A study by Brentani et al. took a different approach by using ESTs to identify SNPs in a predefined set of cancer-related genes; this resulted in the identification of 237 previously known and 505 new SNPs in these genes[31]. A comprehensive analysis by Qiu and coworkers involved cross-mapping of the EST database (dbEST) and the database SNP (dbSNP), yielding a statistically significant association with tumors for 4865 SNPs[32].

These studies emphasize the potential of EST analysis for detecting genomic and expression features associated with cancer. However, they are not particularly informative in terms of uncovering potential causative roles of individual genes in tumorigenesis. We were interested in mining dbEST for somatic mutations that could be positively selected in cancers, which would make the respective genes candidate oncogenes. The inherent problem of such analysis is distinguishing somatic mutations from sequencing errors and SNPs. However, the latter two sources of sequence variation are not expected to produce a signature of positive selection. Indeed, whatever biases are prevalent among sequencing errors, they would not effect the non-synonymous to synonymous substitutions. The issue with SNPs, obviously, is more complex. However, most if not all human SNPs appear to be either selectively neutral or slightly deleterious and do not show signs of frequent positive selection[33, 34]. Accordingly, the signature of positive selection, namely, an elevated non-synonymous/synonymous substitution ratio [27, 35], is expected to be detectable among somatic mutations even in the presence of some contamination by sequencing errors and SNPs.

With this premise, we partitioned the EST sequence libraries available through the dbEST database (NCBI, NIH, Bethesda) into those originating from tumors (hereinafter cancer ESTs) and those from normal tissues (normal ESTs), and identified genes with a significant excess of non-synonymous substitutions in each of the two sets. The results suggest that positive selection is more pronounced in somatic evolution of tumors than it is in normal tissues. Many genes with a signature of positive selection in tumors have established or strongly predicted links to cancer.

Results

Signatures of purifying and positive selection in cancer ests

The ESTs from both tumors and normal cells showed a much lower ratio of non-synonymous to synonymous substitutions than expected under the model of neutral evolution (Table 1), indicating that most of these substitution were genuine mutations subject to purifying selection rather than sequencing errors or neutral SNPs. This notion was reinforced by the substantial, highly significant deficit of nonsense mutations in both EST collections compared to the neutral expectation (Table 1). Remarkably, however, the excess of synonymous over non-synonymous substitutions was less pronounced in cancer ESTs than in normal ESTs. In other words, cancer ESTs showed a significantly greater non-synonymous/synonymous substitution ratio than the normal ESTs (Table 1). (P = 3.7 × 10-32 by Fisher's two-tail exact test) This observation suggests that, compared to the somatic evolution of normal cells, somatic evolution of cancers is characterized by relaxed purifying selection and/or substantial positive selection in some genes. The latter possibility was of special interest because positive selection of somatic mutations in cancers might imply that the change in function of the respective genes was relevant for tumorigenesis and could lead to prediction of previously undetected oncogenes. Therefore we systematically screened cancer and normal EST sequences for indications of positive selection by counting synonymous and non-synonymous substitutions after controlling for sequence quality and subtracting the known SNPs (see Materials and Methods for details).

Table 1 Putative somatic mutations in cancer and normal ESTs.
Full size table

This screening identified 112 genes with a significant excess of non-synonymous over synonymous substitutions compared to the random expectation in tumors and 37 such genes among normal ESTs (Tables 2 and 3; see Additional file 1 ). The difference between the fractions of genes with a significant excess of non-synonymous substitutions in cancer and normal ESTs was highly statistically significant (P < 10-7 by the Fisher's exact test). Furthermore, a comparison of the mutation spectra in cancer-specific ESTs in the 112 concatenated genes with the mutation spectra of the same genes in ESTs from normal tissues revealed highly significant differences (P < 10-5)[36]. Many of the mutated ESTs contain so-called mutational hotspots which, for the purpose of this study, were operationally defined as sites with three or more mutations (Tables 2 and 3; see Additional file 1 ). Specifically, the 112 genes with excess of non-synonymous substitutions in cancer ESTs contained 341 hotspots, whereas the same genes contained 206 hotspots in ESTs from normal tissues. Only 63 hotspots were represented in both sets of ESTs, indicating the presence of a large number of cancer-specific hotspots.

Excess of non-synonymous over synonymous substitutions is considered to be a signature of positive selection [27, 35]. For 51 of the 112 genes with such a signature in cancer ESTs, there was also a statistically significant excess of non-synonymous substitutions in a direct comparison with the normal ESTs derived from the same genes (Table 2; see Additional file 1 and Methods for details). Notably, these genes did not seem to have an excess of nonsense mutations (Table 2; see Additional file 1) which suggests that they are not subject to strong selection for loss of function and that at least some of the positively selected amino acid replacements might be associated with gain-of-function.

The methodology employed here dictates that the list of genes with cancer-specific positive selection (CASPS) is dominated by genes that are highly expressed, particularly, in tumors; typically, statistical significance of the positive selection signature could be demonstrated only for genes for which numerous ESTs were available (Table 2; see Additional file 1). Remarkably, however, the list included 30 moderately or even weakly expressed genes that had no synonymous substitutions but had from 4 to 22 non-synonymous substitutions (Table 2; see Additional file 1).

Cancer connections of CASPS genes and biological implications

Objectively assessing the relevance of a particular gene to tumorigenesis is no easy task. Numerous genes are linked to one or another aspect of cell proliferation, and the expression of many others is perturbed in tumors, which does not necessarily point to an actual role in tumorigenesis. With all these caveats in mind, we nevertheless collated the available data on biological properties of the genes with an apparent signature of positive selection in cancer-derived and normal ESTs and examined their established and potential connections to tumorigenesis. This examination indicated that 42 of the 112 CASPS genes had definitive, experimentally supported connections to tumorigenesis (labeled 'yes' in Table 2; see Additional file 1 ), and for 21 more genes, indirect but strongly suggestive evidence of cancer connections was available (labeled 'likely'). The CASPS genes include those for several proteins involved in DNA repair, programmed cell death, and various forms of signal transduction, among which ubiquitin signaling was particularly prominently represented. All these proteins are directly linked to the control of cell proliferation. Many other CASPS genes do not have such well-defined roles but are substantially overexpressed or amplified in certain types of tumors, which is compatible with involvement in tumorigenesis (Table 2; see Additional file 1).

Figures 1, 2, 3 show the distributions of putative somatic mutations in the sequences of three CASPS genes. These distributions illustrate the complexity of mutational patterns, with distinct spectra seen in cancer and normal ESTs, and a variety of cancer-specific and normal-specific hotspots (Figs. 1,2). Note, however, the presence of a SNP in the same position in both cancer and normal ESTs (Fig. 3)

Figure 1
figure1

Distribution of non-synonymous substitutions in cancer and normal ESTs. from the protein tyrosine phosphatase IV A2 gene. The substitutions are shown in the form of amino acid replacements. Identical letters in the same position correspond to substitutions from different EST libraries. The alignments were constructed using the MACAW program [44]. Species name abbreviations: CANFA, Canis familiaris, GALGA, Gallus gallus, HOMSA, Homo sapiens, RATNO, Rattus norvegicus, XENTR, Xenopus tropicalis. TETNI, Tetraodon nigrovirides.

Full size image
Figure 2
figure2

Distribution of non-synonymous substitutions in cancer and normal ESTs. from the microsomal glutathione S transferase gene. The designations are as in Fig. 1. Additional species name abbreviations: ANOGA, Anopheles gambiae, DANRE, Danio rerio, DROME, Drosophila melanogaster, MUSMU, Mus musculus, STRPU, Strongylocentrotus purpuratus, XENLA, Xenopus laevis.

Full size image
Figure 3
figure3

Distribution of non-synonymous substitutions in cancer and normal ESTs. from the claudin 7 gene. The designations are as in Fig. 1. The site with multiple valine substitutions corresponds to a SNP. Additional species name abbreviations: HALRO, Halocynthia roretzi, TAKRU, Takifugu rubripes.

Full size image

The CASPS genes included only one well-characterized oncogene, the ret protooncogene (Table 2; see Additional file 1). We compared the list of CASPS genes to the comprehensive list of cancer-related set reported in a recent census [13]. In our analysis of EST sequences, mutation spectra with more than 5 mutations were detected for 83 of the 249 cancer-related genes, and 5 of these belonged to the CASPS list. The probability of observing 5 or more genes from a list of 249 among the 112 CASPS genes is ~0.025, suggesting a weak but non-random connection between the CASPS genes detected here and previously characterized cancer-related genes.

Discussion

The interpretation of the findings on CASPS genes described here requires extreme caution. Although filters were applied to separate somatic mutations from sequencing errors and SNPs (see Methods for details), it is impossible to guarantee that the final list is free of these irrelevant sources of variation. Furthermore, taking into consideration the number of analyzed ESTs, identification of 112 genes with apparent signs of positive selection is, in itself, not particularly surprising. The strongest indication we obtained that some of the CASPS genes are likely to be associated with tumorigenesis is the significant excess of genes with the positive selection signature among cancer ESTs compared to the ESTs from normal tissues (112 against 37). Based on this ratio and assuming that the apparent signature of positive selection in normal ESTs represents the background noise, it should be expected that ~70% of the CASPS genes are, indeed, subject to positive selection during the somatic evolution of tumors. Additionally, the evidence seems convincing for those genes that, individually, showed a significant difference in the non-synonymous to synonymous substitution ratio between cancer and normal ESTs (Table 2; see Additional file 1). From a different perspective, however, it is not certain that somatic mutations in normal tissues are not selected for. Furthermore, it cannot be ruled out that some of the genes that seem to evolve under positive selection in normal tissues are associated with the development of precancerous conditions.

Assuming that there is, indeed, a signal of tumor-specific positive selection in our list of CASPS genes, these are likely to be the tip of the proverbial iceberg of genes that evolve under this regimen in various cancers. Although the current EST database is large and represents most of human genes, it is far from being satisfactory for the purpose of analysis of somatic evolution. In the present study, we had no choice but to lump together ESTs from all types of cancers because the amount of variation in individual tumor types was insufficient for statistical analysis. Furthermore, as already indicated, this analysis is capable of detecting selection only for relatively highly expressed genes. Many genes on our CASPS list and more genes that did not make it contained only several non-synonymous substitutions with no synonymous substitutions. Obviously, the statistical power of the present analysis was insufficient to identify positive selection in these genes.

It is expected that, once the EST or complete cDNA data becomes sufficient for separate analysis of tumors of different origins or, ideally, different cell types and tumor progression stages from individual patients, approaches similar to those employed in this work will provide a wealth of information on somatic evolution of the cancer genome. Establishing ancestor-descendant relationships within individuals will allow one to arrive to definitive conclusions regarding the selection forces in action during tumorigenesis.

Conclusion

With all due caution, those genes in the CASPS list that met both criteria (significant excess of synonymous mutations and significant difference between cancer and normal ESTs) could be interesting candidates for a detailed analysis aimed at characterization of new oncogenes or genes with other, still poorly understood roles in tumorigenesis. Furthermore, the results of this work emphasize the value of massive EST and mRNA sequencing from various tumor types (or, ideally, from individual tumors) for identifying the complete catalog of genes with a causal role in tumorigenesis.

Methods

Data and sequence comparisons

The 25801 non-redundant (all identical sequences were removed) coding sequences (CDS) of human genes from the human genome draft build 35, the April 2004 freeze, obtained at the NCBI ftp server [37] were used as reference sequences to be compared with the EST sequences. The EST sequences were from the dbEST release of August, 2004 [38]. EST Library information was extracted and loaded into a mysql database [39]. Each library was manually curated and assigned to either cancer – related (1413 entries) or normal tissue (1370 entries) bins.

The CDS set was searched against dbEST using the BLASTN program with the default parameters[40]. Unigene Build #173 [41] was used to assign ESTs to a particular locus. Overall, 1844057 Unigene EST hits were identified (Table 1).

The database of single-nucleotide polymorphisms (dbSNP) in the fasta format was downloaded from [42] as of May, 2004 (build 121). Altogether, 20573 non-redundant SNPs were identified in the analyzed set of the CDS to the dbSNP consortium specifications [43] by performing a MEGABLAST search of CDS set against dbSNP (command line: megablast -U T -F m -J F -X 180 -r 10 -q -20 -P 1000 -R T -W 28). All alignments containing gaps were discarded. Using a custom PERL script, the MEGABLAST report was parsed for reliable SNPs by ensuring the identity of the RefSeq sequence with the sequence in dbSNP.

Identification of probable somatic mutations in EST sequences

To ensure that all analyzed substitutions came from high-quality sequence, a single-nucleotide substitution in an EST sequence was considered a probable somatic mutation if it was flanked, from each side, with 15 nucleotide stretches of perfect identity between the EST sequence and the reference sequence (CDS) and, in addition, a 50-nucleotide stretch with at most 3 mismatches flanking the identical 15 mers on each side. Substitutions that coincided with SNPs from dbSNP and redundant substitutions from the same EST library were discarded. The latter, highly conservative filter was applied to eliminate possible additional, relatively rare SNPs, which are not reported in the current release of dbSNP, and to ensure clonality of all analyzed mutations. The effects of these filters on the analyzed mutations spectra are shown in Table 4. Statistical significance of differences between substitution spectra was determined using a modified χ2 test [36].

Table 4 The effects of filtering for SNPs and redundant substitutions on the substitution spectra
Full size table

Statistical significance of the differences between ratios of non-synonymous to synonymous substitutions

For each EST, the numbers of synonymous and non-synonymous and substitutions was calculated. These numbers was compared to: i) expected numbers and ii) the numbers for the respective genes in the ESTs from normal (cancer) tissues. The expected values were calculated using a Monte-Carlo random permutation procedure, which was repeated 1000 times for each mutation spectrum. Each step involved random shuffling of transitions/transversions along the appropriate nucleotide sites in the CDS, e.g., for a mutation A->G, the acceptable sites for permutations were those that contained A. Binomial 1-tailed test was used for assessing the statistical significance of non-synonymous vs synonymous substitutions bias in the form:

P ( N o b s , n ) = ∑ k = N o b s n n ! k ! ( n − k ) ! p k q n − k , MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaacqWGqbaucqGGOaakcqWGobGtdaWgaaWcbaGaem4Ba8MaemOyaiMaem4CamhabeaakiabcYcaSiabd6gaUjabcMcaPiabg2da9maaqahabaWaaSaaaeaacqWGUbGBcqGGHaqiaeaacqWGRbWAcqGGHaqicqGGOaakcqWGUbGBcqGHsislcqWGRbWAcqGGPaqkcqGGHaqiaaaaleaacqWGRbWAcqGH9aqpcqWGobGtcqWGVbWBcqWGIbGycqWGZbWCaeaacqWGUbGBa0GaeyyeIuoakiabdchaWnaaCaaaleqabaGaem4AaSgaaOGaemyCae3aaWbaaSqabeaacqWGUbGBcqGHsislcqWGRbWAaaGccqGGSaalaaa@5775@

Where n = N obs +S obs is the total number of non-synonymous and synonymous substitutions observed,

p = N exp n MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaacqWGWbaCcqGH9aqpdaWcaaqaaiabd6eaonaaBaaaleaacyGGLbqzcqGG4baEcqGGWbaCaeqaaaGcbaGaemOBa4gaaaaa@361F@
(1)

is the ratio of the number of non-synonymous substitutions derived from the Monte-Carlo procedure to the total number of substitutions, and q = 1-p.

One-tailed Fisher's exact test was used to assess the significance of the difference between the cancer and normal substitution spectra.

Abbreviations

CASPS:

Cancer-associate positive selection

EST:

Expressed sequence tag

SNP:

single-nucleotide polymorphism

References

  1. 1.

    Knudson AG: Cancer genetics. Am J Med Genet. 2002, 111 (1): 96-102. 10.1002/ajmg.10320.

    Article  PubMed  Google Scholar 

  2. 2.

    Vogelstein B, Kinzler KW: The Genetic Basis of Human Cancer. 2002, New York , McGraw Hill

    Google Scholar 

  3. 3.

    Kinzler KW, Vogelstein B: Cancer-susceptibility genes. Gatekeepers and caretakers. Nature. 1997, 386 (6627): 761, 763-10.1038/386761a0.

    Article  PubMed  Google Scholar 

  4. 4.

    Vogelstein B, Kinzler KW: Cancer genes and the pathways they control. Nat Med. 2004, 10 (8): 789-799. 10.1038/nm1087.

    CAS  Article  PubMed  Google Scholar 

  5. 5.

    Lengauer C, Kinzler KW, Vogelstein B: Genetic instabilities in human cancers. Nature. 1998, 396 (6712): 643-649. 10.1038/25292.

    CAS  Article  PubMed  Google Scholar 

  6. 6.

    Blagosklonny MV: Molecular theory of cancer. Cancer Biol Ther. 2005, 4 (6): 621-627.

    CAS  Article  PubMed  Google Scholar 

  7. 7.

    Levitt NC, Hickson ID: Caretaker tumour suppressor genes that defend genome integrity. Trends Mol Med. 2002, 8 (4): 179-186. 10.1016/S1471-4914(02)02298-0.

    CAS  Article  PubMed  Google Scholar 

  8. 8.

    Boveri T: Zur Frage der Entstehung maligner Tumoren. 1914, Jena , Gustav Fischer

    Google Scholar 

  9. 9.

    Cahill DP, Kinzler KW, Vogelstein B, Lengauer C: Genetic instability and darwinian selection in tumours. Trends Cell Biol. 1999, 9 (12): M57-60. 10.1016/S0962-8924(99)01661-X.

    CAS  Article  PubMed  Google Scholar 

  10. 10.

    Vineis P: Cancer as an evolutionary process at the cell level: an epidemiological perspective. Carcinogenesis. 2003, 24 (1): 1-6. 10.1093/carcin/24.1.1.

    CAS  Article  PubMed  Google Scholar 

  11. 11.

    Breivik J: The evolutionary origin of genetic instability in cancer development. Semin Cancer Biol. 2005, 15 (1): 51-60. 10.1016/j.semcancer.2004.09.008.

    CAS  Article  PubMed  Google Scholar 

  12. 12.

    Gatenby RA, Vincent TL: An evolutionary model of carcinogenesis. Cancer Res. 2003, 63 (19): 6212-6220.

    CAS  PubMed  Google Scholar 

  13. 13.

    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.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  14. 14.

    Hollstein M, Sidransky D, Vogelstein B, Harris CC: p53 mutations in human cancers. Science. 1991, 253 (5015): 49-53.

    CAS  Article  PubMed  Google Scholar 

  15. 15.

    Levine AJ: p53, the cellular gatekeeper for growth and division. Cell. 1997, 88 (3): 323-331. 10.1016/S0092-8674(00)81871-1.

    CAS  Article  PubMed  Google Scholar 

  16. 16.

    Harris SL, Levine AJ: The p53 pathway: positive and negative feedback loops. Oncogene. 2005, 24 (17): 2899-2908. 10.1038/sj.onc.1208615.

    CAS  Article  PubMed  Google Scholar 

  17. 17.

    Lane DP, Benchimol S: p53: oncogene or anti-oncogene?. Genes Dev. 1990, 4 (1): 1-8.

    CAS  Article  PubMed  Google Scholar 

  18. 18.

    Dittmer D, Pati S, Zambetti G, Chu S, Teresky AK, Moore M, Finlay C, Levine AJ: Gain of function mutations in p53. Nat Genet. 1993, 4 (1): 42-46. 10.1038/ng0593-42.

    CAS  Article  PubMed  Google Scholar 

  19. 19.

    Hsiao M, Low J, Dorn E, Ku D, Pattengale P, Yeargin J, Haas M: Gain-of-function mutations of the p53 gene induce lymphohematopoietic metastatic potential and tissue invasiveness. Am J Pathol. 1994, 145 (3): 702-714.

    CAS  PubMed  PubMed Central  Google Scholar 

  20. 20.

    Parant JM, Lozano G: Disrupting TP53 in mouse models of human cancers. Hum Mutat. 2003, 21 (3): 321-326. 10.1002/humu.10186.

    CAS  Article  PubMed  Google Scholar 

  21. 21.

    Pugacheva EN, Ivanov AV, Kravchenko JE, Kopnin BP, Levine AJ, Chumakov PM: Novel gain of function activity of p53 mutants: activation of the dUTPase gene expression leading to resistance to 5-fluorouracil. Oncogene. 2002, 21 (30): 4595-4600. 10.1038/sj.onc.1205704.

    CAS  Article  PubMed  Google Scholar 

  22. 22.

    Blagosklonny MV: p53 from complexity to simplicity: mutant p53 stabilization, gain-of-function, and dominant-negative effect. Faseb J. 2000, 14 (13): 1901-1907. 10.1096/fj.99-1078rev.

    CAS  Article  PubMed  Google Scholar 

  23. 23.

    Lang GA, Iwakuma T, Suh YA, Liu G, Rao VA, Parant JM, Valentin-Vega YA, Terzian T, Caldwell LC, Strong LC, El-Naggar AK, Lozano G: Gain of function of a p53 hot spot mutation in a mouse model of Li-Fraumeni syndrome. Cell. 2004, 119 (6): 861-872. 10.1016/j.cell.2004.11.006.

    CAS  Article  PubMed  Google Scholar 

  24. 24.

    Olive KP, Tuveson DA, Ruhe ZC, Yin B, Willis NA, Bronson RT, Crowley D, Jacks T: Mutant p53 gain of function in two mouse models of Li-Fraumeni syndrome. Cell. 2004, 119 (6): 847-860. 10.1016/j.cell.2004.11.004.

    CAS  Article  PubMed  Google Scholar 

  25. 25.

    Glazko GV, Koonin EV, Rogozin IB: Mutation hotspots in the p53 gene in tumors of different origin: correlation with evolutionary conservation and signs of positive selection. Biochim Biophys Acta. 2004, 1679 (2): 95-106.

    CAS  Article  PubMed  Google Scholar 

  26. 26.

    Koonin EV, Rogozin IB, Glazko GV: p53 gain-of-function: tumor biology and bioinformatics come together. Cell Cycle. 2005, 4 (5): 686-688.

    CAS  Article  PubMed  Google Scholar 

  27. 27.

    Hurst LD: The Ka/Ks ratio: diagnosing the form of sequence evolution. Trends Genet. 2002, 18 (9): 486-10.1016/S0168-9525(02)02722-1.

    Article  PubMed  Google Scholar 

  28. 28.

    Boguski MS, Lowe TM, Tolstoshev CM: dbEST--database for "expressed sequence tags". Nat Genet. 1993, 4 (4): 332-333. 10.1038/ng0893-332.

    CAS  Article  PubMed  Google Scholar 

  29. 29.

    Xu Q, Lee C: Discovery of novel splice forms and functional analysis of cancer-specific alternative splicing in human expressed sequences. Nucleic Acids Res. 2003, 31 (19): 5635-5643. 10.1093/nar/gkg786.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  30. 30.

    Gupta S, Zink D, Korn B, Vingron M, Haas SA: Strengths and weaknesses of EST-based prediction of tissue-specific alternative splicing. BMC Genomics. 2004, 5 (1): 72-10.1186/1471-2164-5-72.

    Article  PubMed  PubMed Central  Google Scholar 

  31. 31.

    Brentani H, Caballero OL, Camargo AA, da Silva AM, da Silva WAJ, Dias Neto E, Grivet M, Gruber A, Guimaraes PE, Hide W, Iseli C, Jongeneel CV, Kelso J, Nagai MA, Ojopi EP, Osorio EC, Reis EM, Riggins GJ, Simpson AJ, de Souza S, Stevenson BJ, Strausberg RL, Tajara EH, Verjovski-Almeida S, Acencio ML, Bengtson MH, Bettoni F, Bodmer WF, Briones MR, Camargo LP, Cavenee W, Cerutti JM, Coelho Andrade LE, Costa dos Santos PC, Ramos Costa MC, da Silva IT, Estecio MR, Sa Ferreira K, Furnari FB, Faria MJ, Galante PA, Guimaraes GS, Holanda AJ, Kimura ET, Leerkes MR, Lu X, Maciel RM, Martins EA, Massirer KB, Melo AS, Mestriner CA, Miracca EC, Miranda LL, Nobrega FG, Oliveira PS, Paquola AC, Pandolfi JR, Campos Pardini MI, Passetti F, Quackenbush J, Schnabel B, Sogayar MC, Souza JE, Valentini SR, Zaiats AC, Amaral EJ, Arnaldi LA, de Araujo AG, de Bessa SA, Bicknell DC, Ribeiro de Camaro ME, Carraro DM, Carrer H, Carvalho AF, Colin C, Costa F, Curcio C, Guerreiro da Silva ID, Pereira da Silva N, Dellamano M, El-Dorry H, Espreafico EM, Scattone Ferreira AJ, Ayres Ferreira C, Fortes MA, Gama AH, Giannella-Neto D, Giannella ML, Giorgi RR, Goldman GH, Goldman MH, Hackel C, Ho PL, Kimura EM, Kowalski LP, Krieger JE, Leite LC, Lopes A, Luna AM, Mackay A, Mari SK, Marques AA, Martins WK, Montagnini A, Mourao Neto M, Nascimento AL, Neville AM, Nobrega MP, O'Hare MJ, Otsuka AY, Ruas de Melo AI, Paco-Larson ML, Guimaraes Pereira G, Pereira da Silva N, Pesquero JB, Pessoa JG, Rahal P, Rainho CA, Rodrigues V, Rogatto SR, Romano CM, Romeiro JG, Rossi BM, Rusticci M, Guerra de Sa R, Sant' Anna SC, Sarmazo ML, Silva TC, Soares FA, Sonati Mde F, de Freitas Sousa J, Queiroz D, Valente V, Vettore AL, Villanova FE, Zago MA, Zalcberg H: The generation and utilization of a cancer-oriented representation of the human transcriptome by using expressed sequence tags. Proc Natl Acad Sci U S A. 2003, 100 (23): 13418-13423. 10.1073/pnas.1233632100.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  32. 32.

    Qiu P, Wang L, Kostich M, Ding W, Simon JS, Greene JR: Genome wide in silico SNP-tumor association analysis. BMC Cancer. 2004, 4: 4-10.1186/1471-2407-4-4.

    Article  PubMed  PubMed Central  Google Scholar 

  33. 33.

    Zhao Z, Fu YX, Hewett-Emmett D, Boerwinkle E: Investigating single nucleotide polymorphism (SNP) density in the human genome and its implications for molecular evolution. Gene. 2003, 312: 207-213. 10.1016/S0378-1119(03)00670-X.

    CAS  Article  PubMed  Google Scholar 

  34. 34.

    Zhang L, Li WH: Human SNPs Reveal no Evidence of Frequent Positive Selection. Mol Biol Evol. 2005

    Google Scholar 

  35. 35.

    Li WH: Molecular Evolution. 1997, Sunderland, MA , Sinauer

    Google Scholar 

  36. 36.

    Adams WT, Skopek TR: Statistical test for the comparison of samples from mutational spectra. J Mol Biol. 1987, 194 (3): 391-396. 10.1016/0022-2836(87)90669-3.

    CAS  Article  PubMed  Google Scholar 

  37. 37.

    NCBI ftp site: Human Genome [ftp://ftp.ncbi.nih.gov/genomes/H_sapiens/].

  38. 38.

    dbEST [http://www.ncbi.nlm.nih.gov/dbEST/].

  39. 39.

    MySQL [http://www.mysql.com].

  40. 40.

    Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997, 25 (17): 3389-3402. 10.1093/nar/25.17.3389.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  41. 41.

    Unigene Build #173 [ftp://ftp.ncbi.nih.gov/repository/UniGene/2004.LOG].

  42. 42.

    dbSNP [ftp://ftp.ncbi.nih.gov/snp/organisms/human_9606/ss_fasta/].

  43. 43.

    dbSNP Readme [ftp://ftp.ncbi.nih.gov/snp/00readme.txt].

  44. 44.

    Schuler GD, Altschul SF, Lipman DJ: A workbench for multiple alignment construction and analysis. Proteins. 1991, 9 (3): 180-190. 10.1002/prot.340090304.

    CAS  Article  PubMed  Google Scholar 

Pre-publication history

  1. The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2407/6/36/prepub

Download references

Acknowledgements

This work was supported by the Intramural Research Program of the National Library of Medicine at the National Institutes of Health/DHHS. FAK is an NSF Graduate Fellow. We thank Yuri Pavlov for helpful discussions.

Author information

Affiliations

  1. National Center for Biotechnology Information, National Library of Medicine, National Institutes of Health, Bethesda, MD, USA

    Vladimir N Babenko, Malay K Basu, Igor B Rogozin & Eugene V Koonin

  2. Section of Ecology, Behavior and Evolution, University of California San Diego, La Jolla, CA, USA

    Fyodor A Kondrashov

Authors
  1. Vladimir N Babenko
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Malay K Basu
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Fyodor A Kondrashov
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Igor B Rogozin
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Eugene V Koonin
    View author publications

    You can also search for this author in PubMed Google Scholar

Corresponding author

Correspondence to Eugene V Koonin.

Additional information

Competing interests

The author(s) declare that they have no competing interests.

Authors' contributions

The study was conceived and designed by EVK and IBR; VNB, MKB, FAK, and IBR performed the computational analysis of the EST mutation data; EVK analyzed the biological aspects of the candidate positively selected genes; VNB wrote the initial draft of the Methods and Results; EVK wrote the final manuscript which was read and approved by all authors.

Electronic supplementary material

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Authors’ original file for figure 3

Rights and permissions

Reprints and Permissions

About this article

Cite this article

Babenko, V.N., Basu, M.K., Kondrashov, F.A. et al. Signs of positive selection of somatic mutations in human cancers detected by EST sequence analysis. BMC Cancer 6, 36 (2006). https://doi.org/10.1186/1471-2407-6-36

Download citation

Keywords

  • Somatic Mutation
  • Sequencing Error
  • Synonymous Substitution
  • Mutation Spectrum
  • Significant Excess

BMC Cancer

ISSN: 1471-2407

Contact us