Abstract
The study of pathway disruption is key to understanding cancer biology. Advances in high throughput technologies have led to the rapid accumulation of genomic data. The explosion in available data has generated opportunities for investigation of concerted changes that disrupt biological functions, this in turns created a need for computational tools for pathway analysis. In this review, we discuss approaches to the analysis of genomic data and describe the publicly available resources for studying biological pathways.
Background
The development of cancer involves the accumulation of genetic and epigenetic alterations. Genetic events such as chromosomal rearrangements, changes in gene dosage, and sequence mutations can influence gene expression patterns, which contribute to the hallmark phenotypes of cancer1,2. The interaction between pathways and the involvement of pathways in multiple phenotypes complicate the interpretation of gene expression patterns. For example, the epidermal growth factor receptor (EGFR, HER1, ERBB-1) signaling pathway plays a role in specific phenotypes including resistance to apoptosis, increased proliferation, mitogenesis, transcription of numerous target genes, and actin reorganization, in several cancers (Fig. 1)1,3,4. In order to decipher the interaction within and between pathways, computational tools are necessary to annotate components, to identify co-regulated expression, and to identify sets of genes or pathways which are statistically over/under-represented within a dataset.

Example of EGFR-mediated signaling changes, a commonly disrupted pathway in lung cancer. The EGFR pathway could be disrupted by an increased expression of growth factor ligands. By targeting EGFR with tyrosine kinase inhibitors (TKIs) and MAb (monoclonal antibodies), EGFR activity can be eliminated. However, a downstream factor (e.g. MAPK signaling pathway) may also be activated to disrupt the pathway, thus making TKIs ineffective. Pathway data was obtained and selected from the Cancer Cell Map database and drawn using Cytoscape.
Methods for Gene Classification
A major analytical step to mine large microarray data is sample classification or identification of gene sets with characteristic biological function. Entrez Gene at the National Center for Biotechnology Information (NCBI) provides unique identifiers for genes, and is a searchable database providing gene-specific information and links to external databases, including the Gene Ontology (GO) consortium, KEGG and Reactome 5 . A limitation of Entrez Gene is that genes are searched individually, which could be time consuming. Here, we describe the Gene Ontology (GO), a structural language to annotate gene functions for batch processing, and also methods of clustering analysis. The algorithmic basis of clustering identifies a pattern associated in a data set, which could be subsequently followed by GO analysis to identify its underlying biology.
Gene Ontology Annotation
The Gene Ontology (GO) Consortium was established in 2000 to provide a controlled vocabulary for annotating homologous gene and protein sequences in different organisms6,7. GO classifies genes and gene products based on three hierarchical structures that describe a given entry's biological processes, cellular components, and molecular functions, and organizes them into a parent-child relationship 6 . Through easy on-line access (http://www.geneontology.org), the genome databases are being unified to expedite the process of retrieving information on genes and proteins based on shared biology among multiple organisms. Several software tools, including GoMiner8,9, MAPPFinder 10 , and Onto-Express11,12, have been developed to explore the GO relationships among high-throughput data. However, the biological functions of genes/proteins are often complex and annotating them into categories may oversimplify their biology. The flat-format output does not convey the richness of GO's hierarchical structure. Nevertheless, this established system of nomenclature of genes and proteins is important for the interoperability of databases, batch processing, and the future design of pathway databases.
Clustering
The biological system is integrative with tightly regulated processes, and genes with similar functions often exhibit coordinated expression patterns13–16. Transcriptional profiling studies typically aim to identify patterns of change among clinically related samples or to classify subgroups of samples15–17. Clustering of microarray data is widely used to identify groups of genes that display coordinated expression patterns performed in a supervised or unsupervised manner (Fig. 2)13,14,17–21. Unsupervised clustering is to classify data without a priori labeling of samples, whereas supervised clustering classifies data based on knowledge of samples type (e.g. cancer subtype)21–24. Clustering techniques are generally classified into two types: hierarchical and partitional25,26. Hierarchical clustering is constructed by either agglomerative (bottom-up) or divisive (top-down) approaches 25 . Agglomerative algorithms begin with separate clusters and merge them into successively larger clusters, while divisive algorithms begin with the whole dataset and divide the data into smaller clusters successively 25 . The output of agglomerative clustering is a tree of clusters called a dendrogram, in which each branch represents group of genes that have a higher order relationship (Fig. 2B)25,27. Partitional clustering directly reduces the dataset into a set of non-overlapping clusters 26 . Representative algorithms of partitional clustering include k-means clustering and self-organizing maps (SOM) 25 . k-means clustering requires the user to define k number of clusters26,28, and SOM partitions data into a two dimensional grid of clusters13,29,30. However, hierarchical clustering is more frequently used17–20,30. Detailed reviews of clustering algorithms are available and this topic will not be discussed further in this review26,31–33.

Graphical output display of heatmap, hierarchical clustering, and principal component analysis.
Dimensionality Reduction
Dimensionality reduction of data is used to minimize the number of input variables for finding coherent patterns of gene expression in an efficient manner25,34,35. Algorithms like principle component analysis (PCA) and multi-dimensional scaling (MDS) both employ this technique for classification procedures25,34,36,37. PCA visualizes multidimensional datasets by projecting data into a sub-space with 2 or 3 dimensions (Fig. 2C)34,35,37,38. The three-dimensional graphical display of MDS can be useful to portray relationships among the data points but might be complex to interpret and require subjective judgments.
Classification analysis may provide some pattern to the experimental datasets. Subsequently, the identified pattern may be further evaluated for biological interpretation using tools such as GO and/or Entrez Gene. However, the inherent limitation of pre-processed databases is subjective to the interpretation of the curator. Therefore, further validation should be considered. In a study that was conducted under the hypothesis that members in the same cluster would share related biological annotations, the majority of the clusters generated by three different clustering algorithms do not correspond well with known biology 39 . Furthermore, there is a need to improve the different clustering algorithms to enhance consistency of the results39,40. It is crucial to associate biological functions or regulatory pathways with each identified cluster of genes in order to deduce biological significance to each sample group 41 .
Construction of Pathway Database
A remarkable number of published articles have collectively yielded thousands of molecular interactions for human and for model species. The challenge is to extract these individual interactions from the literature and to comprehend the dynamics of the interlocking networks as a whole. In recent years, massive efforts have been devoted to managing, integrating, and interpreting the available scientific information in a meaningful manner (i.e. building interactomes or networks of genes and pathways)42,43. Three categories of information are essential for the construction of interactome databases: gene and protein sequences, gene and protein biological information, and molecular interaction resources (Fig. 3). The major repositories of genes and protein sequences are listed in Table 1. Examples of nucleotide sequence databases include NCBI GenBank, EMBL, and DDBJ, all of which are part of the International Nucleotide Sequence Database Collaboration to facilitate data exchange and enhance accuracy44–47. The major databases for gene and protein biological information are listed in Table 2. Gene Ontology (GO), OMIM, Entrez Gene, and Universal Protein Resource Knowledgebase (UniProtKB) are the foundation for building these hierarchical databases5,7,48,49. The main publicly available molecular interaction databases are listed in Table 3. Currently, DIP, IntAct, MINT, HPRD, and MIPS all support the Human Proteome Organization (HUPO) Proteomics Standards Initiative Molecular Interaction (PSI-MI) standard format50–55. This is a unified data standard to represent molecular interaction data in a controlled vocabulary, which facilitates data comparison, exchange, and linking queries together 51 .

Biological knowledgebases contain a myriad of specific information on each gene/protein. Sequence databases are the basis for gene and protein information. Gene and protein information is further extracted and their inter-relationships are experimentally identified, building molecular interaction databases. All of this information is the foundation of pathway databases.
Gene and protein databases.
Gene and protein information databases.
Molecular interaction database.
The wealth of biological resources can complicate the construction of pathway databases (Fig. 4). When assembling information into a pathway database, developers must be cautious to distinguish those interactions that are deduced from hypothetical situations from those that have been experimentally confirmed. Within the latter group, care must also be taken to determine whether interactions have been confirmed in a single direct experiment or a high-throughput experiment. Furthermore, the use of natural language processing (NLP) systems to automate the extraction of information from published articles and to identify relationships between gene and protein names or interactions must be reviewed for biological relevance56,57. This method is useful as a first-pass tool for mining and extracting the knowledge in the literature. However, the constantly advancing nature of research, the further refinement of biological knowledge associated with each gene or protein further refining, the incompletion of the annotation database, and the complexity of entity names in the biological domain often makes it challenging for NLP to be high-quality with huge successes.

An approach to building pathway databases. Biological knowledgebases are represented as rectangles with squared edges. Computational tools for text-mining and language control are represented as ellipses. Molecular interaction and pathway databases are represented by rectangles with rounded edges.
Descriptions of Specific Pathway Database
Pathway databases facilitate the data mining process for cancer researchers. The major pathway databases are listed in Table 4. A collection of biological pathway and network databases is summarized in Table 5, including Pathguide: The Pathway Resource List (http://www.pathguide.org) 58 . This website is updated regularly and currently about 224 biological pathway resources are accessible through the Pathguide website. Here, we focus on a subset of databases that are publicly available.
Pathway databases.
Collection of databases.
KEGG
The KEGG (Kyoto Encyclopedia of Genes and Genomes) database has been established since 1995 and has been one of the most popular knowledge databases to date 59 . The KEGG PATHWAY database consists of manually assembled pathway maps based on inspection of published literature. Pathway maps are grouped into metabolism, genetic information processing, environmental information processing, cellular processes, human diseases, and drug development. Most of the pathways associated with cancer are listed in the environmental information processing section, which is further subdivided into membrane transport, signal transduction, and signaling molecules and interaction. Beside human databases, information from other model organisms such as chimpanzee, mouse, rat, dogs, cows, and pigs is also available. KEGG pathway maps can be manipulated through the KEGG Markup Language (KGML) files, which provide graphical information to customize pathways.
The Cancer Cell Map
The Cancer Cell Map (http://cancer.cellmap.org) is the only database that focuses on signaling pathways implicated in cancer. This resource contains ten cancer-related pathways and each pathway has approximately 100 to 400 interactions. Interactions are manually curated and reviewed for biological validity. Extensive information is provided in each pathway, including the cellular locations of the proteins, the types of physical interactions including molecular interaction, biochemical reaction, catalysis and transport, and post-translational protein modifications. The original citations, experimental evidences, and links to other databases are also listed. Gene expression data can also be visualized in the context of Cancer Cell Map pathways using the Cytoscape network visualization and analysis software 60 .
Human Protein Reference Database
The HPRD (Human Protein Reference Database) contains ten cancer signaling pathways and ten immune signaling pathways which are graphically visualized in GenMAPP pathway maps54,61. The HPRD also offers the flexibility for investigators to refine their search of interested protein by multiple criteria, including molecular class from GO, domain name, motif, site of expression, length of protein sequence, molecular weight, and disease association (e.g. ovarian cancer and breast cancer). The protein domain architecture is graphically visualized with description of the domains and motifs within the queried protein. Post-translational modifications, protein interactions, and disease type are linked to PubMed, OMIM, Swiss Prot, Gene-Prot, Entrez Gene, or pathways within the HPRD. Individual genes within the pathway map are also linked to biologically relevant databases. Results from pathway analysis and HPRD entries can be readily exported. The use of XML (extensible markup language) for HPRD entries makes this database interoperable with other public databases. As with Cytoscape, the development of GenMAPP allows users to map microarray data onto pathway maps 61 .
Reactome
Reactome is a publicly available, peer-reviewed resource of human biological pathways. 62 Although the primary focus is on H. sapiens, it is now extending human pathways onto other organisms via putative orthologs to make them applicable to 21 model organisms, including mouse, rat, chicken, puffer fish, worm, fly, yeast, and E. coli. All the information in Reactome is cross-referenced with PubMed, GO, and the sequence databases at NCBI, Ensembl, and UniProt. Small molecules are linked to ChEBI (http://www.ebi.ac.uk/chebi), catalyst activities to the GO molecular function ontology, and sub-cellular locations to the GO cellular compartment ontology. The OMIM morbid map can be overlaid into the reaction map to see which genes have been implicated in the disease in the literature. Reactions from direct evidence in the literature and indirect evidence that are inferred via orthology in other species are indicated by color-coding. The Reactome Sky Painter tool facilitates the labeling of genes or proteins in the reaction maps. Thus, quantitative data from microarray experiments can be superimposed on Reactome maps to provide visualization and exploration in a pathway context.
Visualization Tools
Cross-talk between pathways can complicate the graphical representation of observed biological interactions. Therefore, visualization tools such as Cytoscape 60 and GenMAPP 61 have been developed to illustrate molecular interactions intuitively.
Cytoscape
Cytoscape is a software tool for the integration of pathways with expression profiles. It allows the querying of networks by using several filtering tools, and linking a given network to public databases for functional annotations 60 . An important feature of Cytoscape is its extensible software framework which allows users to implement new algorithms and network computations. In addition to its use by the Cancer Cell Map (described above), Cytoscape can also be used in conjunction with other protein interaction databases or genetic interaction databases63,64. Molecular species are represented as nodes and intermolecular interactions are linked as edges. Different visual properties such as node color, shape, and size can be chosen, and subsets of nodes and edges can be displayed based on the criteria that are selected by the user. Visualization properties and analysis parameters are customizable.
GenMAPP
GenMAPP (Gene Map Annotator and Pathway Profiler; previously called Gene MicroArray Pathway Profiler) is a computer program designed to display gene expression data in the context of biological pathways 61 . Based on the quantitative data that is loaded, the program will map genes onto relevant pathways and the user can set up criteria to color code the genes accordingly. GenMAPP visualize data in a file format called “MAPPs”, which allow users to organize the genes by their functional component. The user has the choice to download specific pathways or from the archive of MAPPs at www.netpath.org. The MAPPs database is manually curated, with interactions derived from textbooks, review articles, and public pathway databases. GenMAPP also has the feature to construct and modify the pathways by the user, a quality that is not possible if analyzing pre-existing pathway databases like EcoCyc, MetaCyc, and KEGG. Gene identification (ID) from GenBank, SWISS-PROT, Gene Ontology, or other known databases is used to link the gene object on the MAPP to public databases like SWISS-PROT or Entrez Gene by selecting the gene of interest. In addition, GenMAPP displays gene expression levels and provides statistical analysis based on the representation of altered genes in a given pathway MAPP.
Software Tools to Analyze HTP Data
GoMiner8,9, MAPPFinder 10 , and EASE 65 are software tools developed to correlate gene expression changes with GO terms to categorize the biological processes, cellular components, or molecular functions that are statistically affected. However, visualization of the pathway networks is challenging and complicated. Many software tools have been developed for microarray researchers to analyze large scale high-throughput data within the context of biological pathways, including the above mentioned Cytoscape and GenMAPP. Some of the most commonly used software tools are listed in Table 6. Here, we describe some of the freely available software tools that provide graphical representations of gene networks.
Software tools.
Pathway Processor
Pathway Processor is designed to visualize whole genome microarray data in the framework of metabolic networks and provides statistical significance of the reliability of each differentially expressed gene 66 . This program displays data based on the information from the KEGG pathway database. Pathway Processor is implemented as two programs: Pathway Analyzer and Expression Mapper. Pathway Analyzer is the portion responsible for the statistical analysis of pathway significance, while Expression Mapper facilitates the visualization of this data on KEGG pathway maps.
Whole Pathway Scope
Whole Pathway Scope (WPS) is a software tool to analyze high-throughput microarray experiments by referencing pathway or gene information from KEGG, BioCarta, and Gene Ontology 67 . The internal database also includes information from the Genetic Association Database and MedGene Database to allow users to rapidly identify disease-associated genes and highlight them inside their network diagram or select them for further network manipulation. One of the key features is the ability to view multiple experiments simultaneously and color-code the expression value with its p-value. In addition, this software allows users to customize their own metabolic pathway and gene groupings with the option of using statistical analysis.
Pathway Explorer
Pathway Explorer is a web-based service available at https://pathwayexplorer.genome.tugraz.at to map expression profiles of genes onto pathway maps extracted from KEGG, BioCarta, and GenMAPP 68 . This web-based service reduces the local requirement for computational resources. It offers customizable analysis of the data by analyzing in a single or multiple pathways, and a right-tailed Fisher's exact test and false discovery rate analysis were applied to determine the significance of the different pathways. Multiple experiments can also be displayed simultaneously on a single pathway with corresponding expression values. Data is linked to publicly available biological databases (e.g. the NCBI Entrez cross-database search, OMIM, KEGG pathways). The online accessibility of PathwayExplorer enables visualization of DNA or gene expression profiles within the context of biological pathways in a rapid manner.
Future Considerations
The development of various computational tools to interrogate biological databases is accelerating the process to understand high-throughput genomic studies. However, these new tools pose new challenges, and one must be cautious about the limitations and errors associated with various databases. For example, it has been reported that when a partial Enzyme Commission (EC) number, which is a combination of four digits to annotate enzymatic activities without the fourth digit, is assigned to a gene, several pathway databases have used partial EC number annotations and inaccurately assigned them to a set of reactions that are associated with the same partial EC number under each orthology group 69 . Pathway database users should be aware of the possible inherent problems associated with any databases due to the variable quality of the published data. Comprehensive examination of the literature, as well as additional experimental validation, should be used to confirm any findings. Cross-platform integrative analysis of genomics, epigenomics, and transcriptional profiling will offer a deeper understanding of the biological complexity underlying disease processes (Fig. 5) 70 . The current challenge is to incorporate these data together for direct comparison, visualization, and analysis in order to identify salient gene candidates 71 . Once this is accomplished, the next step will be to place these candidates in the context of their proper signaling pathways for a given cancer type. Ultimately, the software programs used to do this should be intuitive to use, provide accurate information, allow customizable analyses, and offer sophisticated statistical tools. All of these features will be essential for characterization of disrupted gene networks in cancer. This will set the stage for rational therapeutic selection based on the underlying genetic realties of a specific tumor38,41.

Genome-wide integrative analysis to identify pathways disrupted in cancer. Genome-wide analyses including copy number profiling, epigenetic profiling, and transcription profiling performed on the same cancer sample could narrow down the number of candidate genes, which would in turn help to pinpoint disrupted pathway involved in cancer.
Footnotes
Acknowledgments
The authors thank Bradley P. Coe and Ian M. Wilson for critical discussion. This work was supported by funds from National Institute of Dental and Craniofacial Research (NIDCR) grants R01 DE15965 and R01 DE13124, Genome Canada/Genome British Columbia, and the Canadian Institutes of Health Research. IFLT is supported by the Roman M. Babicki fellowship. RC and TPHB are supported by scholarships from the Canadian Institutes of Health Research and Michael Smith Foundation for Health Research.
