A large scale survey reveals that chromosomal copy-number alterations significantly affect gene modules involved in cancer initiation and progression
© Alloza et al; licensee BioMed Central Ltd. 2011
Received: 2 March 2010
Accepted: 6 May 2011
Published: 6 May 2011
Recent observations point towards the existence of a large number of neighborhoods composed of functionally-related gene modules that lie together in the genome. This local component in the distribution of the functionality across chromosomes is probably affecting the own chromosomal architecture by limiting the possibilities in which genes can be arranged and distributed across the genome. As a direct consequence of this fact it is therefore presumable that diseases such as cancer, harboring DNA copy number alterations (CNAs), will have a symptomatology strongly dependent on modules of functionally-related genes rather than on a unique "important" gene.
We carried out a systematic analysis of more than 140,000 observations of CNAs in cancers and searched by enrichments in gene functional modules associated to high frequencies of loss or gains.
The analysis of CNAs in cancers clearly demonstrates the existence of a significant pattern of loss of gene modules functionally related to cancer initiation and progression along with the amplification of modules of genes related to unspecific defense against xenobiotics (probably chemotherapeutical agents). With the extension of this analysis to an Array-CGH dataset (glioblastomas) from The Cancer Genome Atlas we demonstrate the validity of this approach to investigate the functional impact of CNAs.
The presented results indicate promising clinical and therapeutic implications. Our findings also directly point out to the necessity of adopting a function-centric, rather a gene-centric, view in the understanding of phenotypes or diseases harboring CNAs.
For many years, research into the genetic basis of many different diseases has cumulated a large corpus of knowledge on the links between particular genes and diseases. However, most diseases are multigenic and cannot be explained or characterized by a single gene but rather by a combination of genes . In a broad sense, multigenecity reflects disruptions in proteins that participate in a protein complex, in a pathway  or, in general, in any known or yet to be discovered functional unit or module.
It is widely accepted that most of the biological functionality of the cell arises from complex interactions between their molecular components that define functional entities or modules that operate as small subcellular systems . Recent evidence shows that genes mapping in close physical locations in the chromosomes tend to have a high degree of co-expression, that can often be linked to common functionality . We have recently demonstrated that chromosomal regions enriched in certain types of functions have arisen at determinate evolutionary moments and have maintained since then in the genomes along the evolution . It is a well established fact that chromosomal regions change their activity in diseases that harbor copy number alterations, such as some cancers . Moreover, examples of disease-related genes mapping together in the genome have recently been described [7, 8].
Despite these observations, the strategy followed in order to study pathologies harboring CNAs is still based on finding one or a few target genes . However, it should be expected that modules composed of functionally-related genes, rather than one or a few key genes, were affected by the chromosomal alterations.
The aim of this study is to evaluate the importance of the functional component that can be attributed to the co-localization of modules of functionally related genes in those regions of the chromosomes that are lost or gained during the curse of the disease, and the possible impact that such component can have on its symptomatology. An interesting model to test this hypothesis is cancer. Genetic instability, producing not only point mutations but mainly CNAs, has long been recognized as a fundamental feature of many cancer types . It has been suggested that only a few essential functional alterations in the physiology of the cells, the so-called hallmarks of cancer, are behind the vast catalog of cancer genotypes . Our hypothesis is that in most cases such functions will be collectively carried out by modules of functionally-related genes that lay close together in the chromosomes. Such modules could consequently be found by analyzing the genomic regions most affected by chromosomal instability.
To this end, a functional analysis of the genes located in the regions undergoing copy number alterations in different cancer types has been conducted. This analysis aims to detect enrichments in cancer-related functions significantly associated to chromosomal regions that are systematically gained or lost in cancers. The conventional way to determine the level of involvement of a gene in the disease would imply the establishment of a threshold beyond which the participation of the gene in deletions of amplifications can be declared significant. The way in which this threshold is fixed tends to be arbitrary and can compromise the conclusions of the work . To avoid the imposition of artificial thresholds, we have applied a generalization of the concept of gene set analysis (GSA), previously proposed in the context of transcriptomics  and further generalized to other genomic domains such as evolution , QTL analysis  and genotyping . In the approach used here a GSA test  is conducted to relate pre-defined sets of functionally-related genes to a variable that represents the propensity of a gene to be involved in a deletion or in an amplification. Such variable is the absolute frequency at which each gene is affected by copy number alterations (CNAs) because it is located in a lost or gained chromosomal region. This family of methods based on the analysis of functionally-related gene sets, (modules therein), are known to be more sensitive than other alternative functional enrichment methods [17–19].
Array-CGH constitute an ideal type of data, given that they provide a detailed and accurate description of the regions gained and lost and consequently a well defined characterization of the genes included in them. Unfortunately, there are still a limited number of such arrays available in the databases. For example the progenetix database , a repository of cytogenetic abnormalities in cancer, contains only 3016 samples of array-CGH at present (August 2009). The more general ArrayExpress database contains 34 experiments, only 10 of which contain the word cancer, summing up a total of only 463 samples. Both databases are highly overlapping and include data from different platforms which makes the analysis even more difficult . With the exception of initiatives such as The Cancer Genome Atlas , in which large numbers of cancers are simultaneously analyzed with Array-CGH, most experiments involve small numbers of samples. Small sample sizes of many experiments hinder the possibility of obtaining significant associations between gene modules and regions gained or lost in cancers. However, this situation will eventually change in the near future.
On the other hand, classical CGH assays provide an alternative low-resolution description of regions gained or lost. Fortunately, a large number of such observations is available in repositories such as the Mitelman database . This is a manually curated repository that relates chromosomal aberrations to tumour features, based information extracted from the literature. Here we conducted a systematic analysis of more than 140.000 observations of CNAs corresponding to losses and gains of genetic material in cancers. The results revealed a significant pattern of losses and gains of modules of functionally-related genes that will affect physiological processes associated to cancer initiation and progression. Consequently, a function-centric view, which takes into account the cooperative effects of the genes that integrate functional modules, can provide a more complete understanding of diseases that harbor CNAs.
Distribution of CNAs
Classical cytogenetic techniques are more suitable for detecting regional deletions (deletions affecting to regions of the chromosomes, accounting for 19859 out of a total of 86048 deletions reported in the Mitelman database) than for detecting regional amplifications (only 1011 out of a total of 55935) so, there is an unbalance between the number of observations in each case, which constitutes an obvious limitation in the conclusions obtained for amplifications using this type of data. Apparently this limitation did not seem to apply to whole chromosome CNAs. There was a clear significant negative correlation between the chromosome size and the number of observed whole chromosome CNA deletions. Probably this observation is reflecting the obvious general fact that the smaller chromosomes have fewer genes and their loss is potentially less damaging for the cell. On the other hand, no equivalent trend could be observed in the case of whole chromosome gains. In this case an average of ~2000 observations (SD ~1000) was recorded for all the chromosomes, except 7, 8 and 21, with 4156, 7013 and 4974 observations respectively. These relevant exceptions were due to the specific and highly recurrent presence of these chromosome aneuploidies in well defined epithelial and hematological cancers.
Particular behaviors of some chromosomes deserve to be mentioned. For example, chromosome Y displayed a large number of deletions but an abnormally low number of amplifications. An opposite trend was observed for chromosomes 8 and 21. Particularly interesting is chromosome 7 that showed an unexpectedly high number of amplifications and deletions simultaneously.
Although these observations have been extensively reported, nothing is known about the advantageous biological reasons that a particular tumor may have due to the acquisition of a chromosomal aberration along its evolution.
Functional roles found in frequently lost regions
GO terms most significantly associated to chromosomal regions frequently lost in cancers.
Number of genes
homophilic cell adhesión
3.8779 × 10-17
4.77 × 10-18
Metastasis and invasion
calcium-dependent cell-cell adhesion
2.0064 × 10-08
9.88 × 10-9
Metastasis and invasion
1.2491 × 10-05
1.46 × 10-6
sensory perception of taste
cellular component organization and biogenesis
fertilization (sensu Metazoa)
Metastasis and invasion
maintenance of fidelity during DNA-dependent DNA replication
male gamete generation
localization of cell
Metastasis and invasion
pancreatic ribonuclease activity
Metastasis and invasion?
serine-type endopeptidase inhibitor activity
sulfate porter activity
Metastasis and invasion
interferon-alpha/beta receptor binding
carboxypeptidase A activity
Metastasis and invasion?
Metastasis and invasion?
microtubule organizing center part
intermediate filament cytoskeleton
integral to plasma membrane
1.50 × 10-001
Metastasis and invasion
All the functions described have been obtained by analyzing all the cancers together. It is difficult to particularize the analysis in individual cancers given the poor resolution of the technology and the comparatively fewer cases recorded for individual cancers. Only in the case of leukemia, with 6859 cases available, was it possible to find significant terms (see Additional file 1 Table S4). Actually, a richer and more detailed picture of the functions of the modules affected by the deletions is obtained. Thus, in addition to processes related to cell adhesion, localization of cell, sulfate transport and other processes described above, more specific terms for signaling were found, such as "cell-cell signaling" (p = 5.79 × 10-03) and some descendants or "cell surface receptor linked signal transduction" (p = 5.18 × 10-03). Other interesting terms that can be associated to metastasis would be "regulation of cell migration" (p = 0.0269). Distinct modules related to phosporylation processes have also been significantly lost. Finally, processes related to the immune system, probably characteristic of this type of cancer, were also lost, such as "antigen processing and presentation" (p = 0.0159).
Additional File 2 contains the list of the genes for each significant GO term with the corresponding chromosomal locations and the number of CNAs in which they are involved.
Functional roles found in frequently amplified regions
GO terms corresponding to the "biological process" ontology significantly associated to chromosomal regions frequently amplified in cancers.
Number of genes
defense response to bacterium
1.7594 × 10 -11
2.98 × 10 -11
xenobiotic metabolic process
defense response to fungus
Chromosome-specific functional roles
While many of the observations made for functional roles lost are based on regional deletions, the observations made for functional roles gained are almost exclusively based on whole chromosome gains. This fact suggests the existence of a strong regionalization of the functional units composed by gene modules. While the concentration of functionally-related genes in neighborhoods has already been described , the existence of chromosome-specific functions, beyond the obvious ascribed to the sexual chromosomes, has not systematically been explored. Additional file 1 Table S5 shows the results obtained upon the application of a functional enrichment method  that finds significant over-representations of functional terms in each chromosome. Surprisingly only a few chromosomes -12, 13, 14, 15 and 21- are not significantly enriched in any functionally-related gene module. Some of the chromosome-specific gene modules present in frequently amplified chromosomes correspond, obviously, to functions systematically gained in cancers (e.g. defense response to bacterium or defense response to fungus in chromosome 8 or mismatch repair in chromosome 7). Although the existence of biological functions specific of one or of a few chromosomes does not seem to be a frequent occurrence (we must remember that there are several thousands of GO terms), some of these, however, are clearly over-represented in some chromosomes. This observation suggests that when a complete chromosome is lost or gained (a frequent event in cancer) the activity of whole modules of genes carrying out a particular biological role are dramatically changed.
Fine-scale analysis of CNAs of functional roles lost in glioblastoma
GO terms corresponding to the "biological process" ontology significantly associated to chromosomal regions frequently lost in the glioblastomas .
sensory perception of chemical stimulus
9.83 × 10-24
1.97 × 10-17
2.39 × 10-17
7.17 × 10-14
4.45 × 10-12
1.78 × 10-06
G-protein coupled receptor protein signaling pathway
5.44 × 10-8
1.60 × 10-5
chromatin assembly or disassembly
1.85 × 10-5
cell surface receptor linked signal transduction
establishment and/or maintenance of chromatin architecture
protein-DNA complex assembly
chromosome organization and biogenesis (sensu Eukaryota)
chromosome organization and biogenesis
lipid metabolic process
organelle organization and biogenesis
regulation of liquid surface tension
The results demonstrate the importance of having high-resolution Array-CGH data. With only 294 samples we have found well defined and highly significant terms. In the conventional CGH case we needed almost two orders of magnitude more samples to obtain similar results.
Given the poor resolution that classical cytogenetic techniques have, the analysis presented here can be considered as a coarse-grain view of the real importance of neighborhoods of modules composed by functionally-related genes. The extension of this methodology to the analysis of the data produced by Array-CGH  has demonstrated to be straightforward. Also the results are noticeably better when a reasonable number of samples are analyzed. The forthcoming deluge of data foreseeable from such high-throughput platforms will allow refining the results presented here in a near future. Moreover, more sophisticated methods of error correction in the GSA methods will result in more specific and accurate findings in the future . In this study we have used the popular FDR, that keeps false discoveries at a reasonable rate without being too conservative .
Despite its low resolution and limitations, the analysis based on cytogenetic data resulted to be sensitive enough to reveal the existence of an important number of gene modules, significantly associated to regions systematically deleted in cancers, whose functionalities (as represented by GO terms) clearly account for cancer-related biological processes and molecular roles. Similarly, regions systematically amplified by cancer cells are significantly populated by genes carrying out processes related to cell immortalization and to unspecific defense (probably against chemotherapeutical agents). When the vast catalog of cancer cell genotypes is analyzed in detail, a small number of functional alterations, also denominated hallmarks of cancer by some authors [11, 28], seems to account for all this complexity. These altered functions shared by most (if not all) cancers are: self-sufficiency in growth signals, insensitivity to antigrowth signals, evasion of apoptosis, limitless replicative potential, sustained angiogenesis, and tissue invasion and metastasis . Interestingly, losses and gains of gene modules carrying out functions (whose GO terms are compatible with the acquisition of all of them) were found. Moreover, another well-established feature of cancer, the acquisition of multi drug resistance by amplification of chromosomal regions  was also found to depend on the amplification of the corresponding gene modules (amplification of the "xenobiotic metabolic process" GO term; see table 2). Finally, the enabling characteristic for the acquisition of these functionalities, genomic instability, also considered a cancel hallmark by other authors , has been clearly detected by our analysis. The analysis of leukemia in the CGH data and the analysis of the TCGA glioblastomas data  revealed, in addition to processes related to general cancer properties, other affected processes related to particular characteristics of the cancer types.
The emerging picture from this study is a scenario in which cancer cell populations undergo a Darwinian dynamics along the distinct phases of the disease , which mainly occurs by genetic modifications caused by chromosomal instability [10, 28]. This chromosomal instability, affecting large chromosomal regions, is the substrate of a strong selective pressure operating on a reduced number of cellular functions leading to cancer . Our hypothesis, consistent with the results found, is that in a considerable number of cases such functions will be carried out by modules of genes that are located in neighborhoods along the chromosomes rather than by one or few key genes.
Regardless of the fact that some recent evidences link chromosomal physical location to function [4, 5] and that disease-related genes have been recently described as mapping in adjacent positions in the chromosomes [6–8], the putative impact that this local distribution of cellular functionality could have in the symptomatology of diseases that harbor CNAs still remains largely unexplored.
The systematic study presented here has firmly related, for the first time, molecular functional roles of gene modules involved in CNAs to biological processes related to cancer initiation and progression. Our conclusions point to the necessity of taking a more function-centric perspective to understand the functional effect of CNAs, especially in the context of cancer. Even more important is the relevance of these findings for the design of therapeutic strategies in the treatment of cancer. Also, drug discovery processes need to account for this new scenario [38, 39]. The recent observation of an unexpectedly widespread distribution and prevalence of CNA polymorphisms in the human genome  makes even more urgent the adoption of a viewpoint in which gene modules, instead of genes in isolation, are considered to be behind most phenotypes or disease symptoms. Our findings are consistent with the idea that pathways, rather than individual genes, appear to govern the course of tumorigenesis [40, 41].
Database of chromosomal alterations
The Mitelman database of Chromosome Aberrations in Cancer http://cgap.nci.nih.gov/Chromosomes/Mitelman was used as primary source of information. A total of 86048 observations (independent CNA event described in the database) of deletions (19859 of them corresponding to regional deletions affecting only to parts of the chromosomes and 66189 to whole chromosome deletions), including any type of cancer, were obtained from the database. Each individual case in the database has a description of the associated CNA in a standard format with the indication of the cytoband or cytobands affected by the deletion or the amplification; for example del(7)(p13p15). Although with a smaller coverage, amplifications were also stored in the database. A total of 55935 observations of amplifications corresponding to 1011 regional amplifications and 54924 whole chromosome amplifications were retrieved from the database.
In order to avoid possible artifacts of erroneous observations, we have only taken into account deletions or amplifications occurring at least in two cases.
Array-CGH data from TCGA
The Cancer Genome Atlas  was used as the source for Array-CGH data on cancer. A large experiment with 294 Array-CGH (Agilent 244 K) of glioblastoma was used . The data already segmented was available and was retrieved for the analysis.
Genes were mapped to their corresponding cytobands using the Ensembl database , release 45, Homo sapiens (NCBI 36).
Functional annotations and representation
The widely accepted GO functional categories  were used for the functional profiling of the genes involved in CNAs. The GO annotations for the human genes were provided by the ID-converter engine of the Babelomics package http://www.babelomics.org which contains annotations taken from the GO database . GO terms corresponding to levels 3 to 6 in the hierarchy were tested. This correspond to a total of 3549 Biological Process, 2817 Molecular Function and 394 Cellular Component GO terms
Plots of the relationships existent among GO terms can be obtained by using the GO visualization tools implemented in the Babelomics package .
Chromosomal functional enrichment
The significance in the enrichment of functional terms in chromosomes was carried out with the FatiGO  program as implemented in the Babelomics package http://www.babelomics.org. Briefly, the program builds a 2 × 2 contingency table for each functional term checked and applies an exact Fisher's test. The resulting p-values are adjusted to compensate multiple-testing effects using the popular FDR method .
Functional profiling of genes more frequently involved in CNAs
We used for this purpose a gene-set functional profiling method. This methodology, introduced in the field of microarray data analysis , comprises a family of methods that allows studying the coordinate over- or under-representation of pre-defined gene sets in the extremes of a list of genes ranked by some criteria . In particular we have used a generalized GSA test  implemented in the Babelomics package . This algorithm has the advantage of dealing directly with the ranked list, without requiring either any extra information or the original experimental values . In particular, the test seeks for significant asymmetrical distributions of GO terms across consecutive partitions of the ranked list . If the term is systematically over-represented in the upper part of the partitions (corresponding to the genes more frequently deleted) we will consider that this GO functionality is systematically deleted. The same rationale can be applied to the list ranked by frequency of amplifications. The FDR method  was used to adjust the p-values for multiple-testing effects given that a large number of GO terms are tested.
GO definitions , obtained as explained above, were used to define such gene modules composed of functionally-related genes.
Validation of the roles of biological processes in cancer
The GoPubmed , a tool that relate terms to genes using sophisticated textmining methods was used to check whether the biological processes found were linked to cancer by co-citations in the scientific literature.
This work is supported by grants BIO2008-04212 and FIS PI 08/0440 from the Spanish Ministry of Science and Innovation and PROMETEO/2010/001 from the GVA-FEDER. The CIBER de Enfermedades Raras is an initiative of the ISCIII. This work is also partly supported by a grant (RD06/0020/1019) from Red Temática de Investigación Cooperativa en Cáncer (RTICC), Instituto de Salud Carlos III (ISCIII), Spanish Ministry of Science and Innovation. EA is supported by a fellowship from the FIS of the Spanish Ministry of Health (FI06/00027).
- Vogelstein B, Lane D, Levine AJ: Surfing the p53 network. Nature. 2000, 408 (6810): 307-310. 10.1038/35042675.View ArticlePubMedGoogle Scholar
- Badano JL, Katsanis N: Beyond Mendel: an evolving view of human genetic disease transmission. Nat Rev Genet. 2002, 3 (10): 779-789.View ArticlePubMedGoogle Scholar
- Hartwell LH, Hopfield JJ, Leibler S, Murray AW: From molecular to modular cell biology. Nature. 1999, 402 (6761 Suppl): C47-52.View ArticlePubMedGoogle Scholar
- Hurst LD, Pal C, Lercher MJ: The evolutionary dynamics of eukaryotic gene order. Nat Rev Genet. 2004, 5 (4): 299-310. 10.1038/nrg1319.View ArticlePubMedGoogle Scholar
- Al-Shahrour F, Minguez P, Marques-Bonet T, Gazave E, Navarro A, Dopazo J: Selection upon genome architecture: conservation of functional neighborhoods with changing genes. PLoS Comput Biol. 2010, 6 (10): e1000953-10.1371/journal.pcbi.1000953.View ArticlePubMedPubMed CentralGoogle Scholar
- Zhou Y, Luoh SM, Zhang Y, Watanabe C, Wu TD, Ostland M, Wood WI, Zhang Z: Genome-wide identification of chromosomal regions of increased tumor expression by transcriptome analysis. Cancer Res. 2003, 63 (18): 5781-5784.PubMedGoogle Scholar
- Israeli O, Goldring-Aviram A, Rienstein S, Ben-Baruch G, Korach J, Goldman B, Friedman E: In silico chromosomal clustering of genes displaying altered expression patterns in ovarian cancer. Cancer Genet Cytogenet. 2005, 160 (1): 35-42. 10.1016/j.cancergencyto.2004.11.011.View ArticlePubMedGoogle Scholar
- Snijders AM, Schmidt BL, Fridlyand J, Dekker N, Pinkel D, Jordan RC, Albertson DG: Rare amplicons implicate frequent deregulation of cell fate specification pathways in oral squamous cell carcinoma. Oncogene. 2005, 24 (26): 4232-4242. 10.1038/sj.onc.1208601.View ArticlePubMedGoogle Scholar
- Pinkel D, Albertson DG: Array comparative genomic hybridization and its applications in cancer. Nat Genet. 2005, 37 (Suppl): S11-17.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Hanahan D, Weinberg RA: The hallmarks of cancer. Cell. 2000, 100 (1): 57-70. 10.1016/S0092-8674(00)81683-9.View ArticlePubMedGoogle Scholar
- Pan KH, Lih CJ, Cohen SN: Effects of threshold choice on biological conclusions reached during analysis of gene expression by DNA microarrays. Proc Natl Acad Sci USA. 2005, 102 (25): 8961-8965. 10.1073/pnas.0502674102.View ArticlePubMedPubMed CentralGoogle Scholar
- Mootha VK, Lindgren CM, Eriksson KF, Subramanian A, Sihag S, Lehar J, Puigserver P, Carlsson E, Ridderstrale M, Laurila E, et al: PGC-1alpha-responsive genes involved in oxidative phosphorylation are coordinately downregulated in human diabetes. Nat Genet. 2003, 34 (3): 267-273. 10.1038/ng1180.View ArticlePubMedGoogle Scholar
- Al-Shahrour F, Arbiza L, Dopazo H, Huerta-Cepas J, Minguez P, Montaner D, Dopazo J: From genes to functional classes in the study of biological systems. BMC Bioinformatics. 2007, 8: 114-10.1186/1471-2105-8-114.View ArticlePubMedPubMed CentralGoogle Scholar
- Wu C, Delano DL, Mitro N, Su SV, Janes J, McClurg P, Batalov S, Welch GL, Zhang J, Orth AP, et al: Gene set enrichment in eQTL data identifies novel annotations and pathway regulators. PLoS Genet. 2008, 4 (5): e1000070-10.1371/journal.pgen.1000070.View ArticlePubMedPubMed CentralGoogle Scholar
- Medina I, Montaner D, Bonifaci N, Pujana MA, Carbonell J, Tarraga J, Al-Shahrour F, Dopazo J: Gene set-based analysis of polymorphisms: finding pathways or biological processes associated to traits in genome-wide association studies. Nucleic Acids Res. 2009, 37: W340-344. 10.1093/nar/gkp481.View ArticlePubMedPubMed CentralGoogle Scholar
- Dopazo J: Formulating and testing hypotheses in functional genomics. Artif Intell Med. 2009, 45 (2-3): 97-107. 10.1016/j.artmed.2008.08.003.View ArticlePubMedGoogle Scholar
- Goeman JJ, Buhlmann P: Analyzing gene expression data in terms of gene sets: methodological issues. Bioinformatics. 2007, 23 (8): 980-987. 10.1093/bioinformatics/btm051.View ArticlePubMedGoogle Scholar
- Huang DW, Sherman BT, Lempicki RA: Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res. 2008, 37 (1): 1-13.View ArticlePubMed CentralGoogle Scholar
- The Progenetix database. [http://www.progenetix.net]
- Larsson O, Sandberg R: Lack of correct data format and comparability limits future integrative microarray research. Nat Biotechnol. 2006, 24 (11): 1322-1323. 10.1038/nbt1106-1322.View ArticlePubMedGoogle Scholar
- The Cancer Genome Atlas. [http://cancergenome.nih.gov/]
- Mitelman Database of Chromosome Aberrations in Cancer. [http://cgap.nci.nih.gov/Chromosomes/Mitelman]
- The GOPubMed web server. [http://www.gopubmed.org/]
- Lengauer C, Kinzler KW, Vogelstein B: Genetic instabilities in human cancers. Nature. 1998, 396 (6712): 643-649. 10.1038/25292.View ArticlePubMedGoogle Scholar
- Volpert OV, Dameron KM, Bouck N: Sequential development of an angiogenic phenotype by human fibroblasts progressing to tumorigenicity. Oncogene. 1997, 14 (12): 1495-1502. 10.1038/sj.onc.1200977.View ArticlePubMedGoogle Scholar
- Stetler-Stevenson WG: Matrix metalloproteinases in angiogenesis: a moving target for therapeutic intervention. J Clin Invest. 1999, 103 (9): 1237-1241. 10.1172/JCI6870.View ArticlePubMedPubMed CentralGoogle Scholar
- Merlo LM, Pepper JW, Reid BJ, Maley CC: Cancer as an evolutionary and ecological process. Nat Rev Cancer. 2006, 6 (12): 924-935. 10.1038/nrc2013.View ArticlePubMedGoogle Scholar
- Al-Shahrour F, Diaz-Uriarte R, Dopazo J: FatiGO: a web tool for finding significant associations of Gene Ontology terms with groups of genes. Bioinformatics. 2004, 20 (4): 578-580. 10.1093/bioinformatics/btg455.View ArticlePubMedGoogle Scholar
- The_Cancer_Genome_Atlas_Research_Network: Comprehensive genomic characterization defines human glioblastoma genes and core pathways. Nature. 2008, 455 (7216): 1061-1068. 10.1038/nature07385.View ArticleGoogle Scholar
- Prado-Lopez S, Conesa A, Arminan A, Martinez-Losa M, Escobedo-Lucea C, Gandia C, Tarazona S, Melguizo D, Blesa D, Montaner D, et al: Hypoxia Promotes Efficient Differentiation of Human Embryonic Stem Cells to Functional Endothelium. Stem Cells. 2010, 28 (3): 407-418.PubMedGoogle Scholar
- Redon R, Ishikawa S, Fitch KR, Feuk L, Perry GH, Andrews TD, Fiegler H, Shapero MH, Carson AR, Chen W, et al: Global variation in copy number in the human genome. Nature. 2006, 444 (7118): 444-454. 10.1038/nature05329.View ArticlePubMedPubMed CentralGoogle Scholar
- Gold DL, Miecznikowski JC, Liu S: Error control variability in pathway-based microarray analysis. Bioinformatics. 2009, 25 (17): 2216-2221. 10.1093/bioinformatics/btp385.View ArticlePubMedPubMed CentralGoogle Scholar
- Benjamini Y, Yekutieli D: The control of false discovery rate in multiple testing under dependency. Annals of Statistics. 2001, 29: 1165-1188. 10.1214/aos/1013699998.View ArticleGoogle Scholar
- Schimke RT: Gene amplification, drug resistance, and cancer. Cancer Res. 1984, 44 (5): 1735-1742.PubMedGoogle Scholar
- Breivik J: Don't stop for repairs in a war zone: Darwinian evolution unites genes and environment in cancer development. Proc Natl Acad Sci USA. 2001, 98 (10): 5379-5381. 10.1073/pnas.101137698.View ArticlePubMedPubMed CentralGoogle Scholar
- Duesberg P, Li R, Fabarius A, Hehlmann R: The chromosomal basis of cancer. Cell Oncol. 2005, 27 (5-6): 293-318.PubMedPubMed CentralGoogle Scholar
- Butcher EC, Berg EL, Kunkel EJ: Systems biology in drug discovery. Nat Biotechnol. 2004, 22 (10): 1253-1259. 10.1038/nbt1017.View ArticlePubMedGoogle Scholar
- Kitano H: A robustness-based approach to systems-oriented drug design. Nat Rev Drug Discov. 2007, 6 (3): 202-210. 10.1038/nrd2195.View ArticlePubMedGoogle Scholar
- Bardelli A, Velculescu VE: Mutational analysis of gene families in human cancer. Curr Opin Genet Dev. 2005, 15 (1): 5-12. 10.1016/j.gde.2004.12.009.View ArticlePubMedGoogle Scholar
- Wood LD, Parsons DW, Jones S, Lin J, Sjoblom T, Leary RJ, Shen D, Boca SM, Barber T, Ptak J, et al: The genomic landscapes of human breast and colorectal cancers. Science. 2007, 318 (5853): 1108-1113. 10.1126/science.1145720.View ArticlePubMedGoogle Scholar
- Hubbard TJ, Aken BL, Ayling S, Ballester B, Beal K, Bragin E, Brent S, Chen Y, Clapham P, Clarke L, et al: Ensembl 2009. Nucleic Acids Res. 2009, D690-697. 37 Database
- 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. Nat Genet. 2000, 25 (1): 25-29. 10.1038/75556.View ArticlePubMedPubMed CentralGoogle Scholar
- Medina I, Carbonell J, Pulido L, Madeira SC, Goetz S, Conesa A, Tarraga J, Pascual-Montano A, Nogales-Cadenas R, Santoyo J, et al: Babelomics: an integrative platform for the analysis of transcriptomics, proteomics and genomic data with advanced functional profiling. Nucleic Acids Res. 2010, 38 (Suppl): W210-213.View ArticlePubMedPubMed CentralGoogle Scholar
- Dopazo J: Functional interpretation of microarray experiments. Omics. 2006, 10 (3): 398-410. 10.1089/omi.2006.10.398.View ArticlePubMedGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1755-8794/4/37/prepub
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.