- Open Access
KPP: KEGG Pathway Painter
BMC Systems Biologyvolume 9, Article number: S3 (2015)
High-throughput technologies became common tools to decipher genome-wide changes of gene expression (GE) patterns. Functional analysis of GE patterns is a daunting task as it requires often recourse to the public repositories of biological knowledge. On the other hand, in many cases researcher's inquiry can be served by a comprehensive glimpse. The KEGG PATHWAY database is a compilation of manually verified maps of biological interactions represented by the complete set of pathways related to signal transduction and other cellular processes. Rapid mapping of the differentially expressed genes to the KEGG pathways may provide an idea about the functional relevance of the gene lists corresponding to the high-throughput expression data.
Here we present a web based graphic tool KEGG Pathway Painter (KPP). KPP paints pathways from the KEGG database using large sets of the candidate genes accompanied by "overexpressed" or "underexpressed" marks, for example, those generated by microarrays or miRNA profilings.
KPP provides fast and comprehensive visualization of the global GE changes by consolidating a list of the color-coded candidate genes into the KEGG pathways. KPP is freely available and can be accessed at http://web.cos.gmu.edu/~gmanyam/kegg/
High-throughput technologies became common tools to decipher genome-wide changes of gene expression (GE) patterns or relative protein abundance. Typical output of these large-scale studies is represented by the list comprised of hundreds of gene candidates with attached quantitative labels. Functional analysis of these gene lists is a daunting task as it requires regular recourse to the public repositories of biological knowledge or use of expensive databases of manually curated biological annotation [1, 2]. On the other hand, in many cases researcher's inquiry can be successfully served by a comprehensive glimpse.
Functional analysis of markers identified from large-scale datasets can be performed using a wide variety of bioinformatics tools. As microarrays became a common tool to decipher global gene expression, centralized systems like Gene Expression Omnibus (GEO), ArrayExpress was developed to congregate the valuable profile data [3, 4]. An analysis of combined datasets generated in independent microarray experiments (so-called "microarray meta-analysis"), is often being employed , for example, to develop biomarker panels or to extract insights into the pathogenesis of various chronic diseases  including human malignancies . Meta-analysis lead to an increase of the complexity in microarray analysis; therefore, sophistication of subsequent functional analysis also increased. Gene Ontology (GO) and other pathway-centered types of analysis became indispensable .
KEGG (Kyoto Encyclopedia of Genes and Genomes) is a compendium of databases covering both annotated genomes and protein interaction networks for all sequenced organisms. Its integral part, KEGG PATHWAY, is a compilation of manually verified pathway maps displaying both the molecular interactions and the biochemical reactions . The recent version of this database includes a complete set of pathways related to signal transduction and other cellular processes . The extensive collection of the pathways at KEGG can be utilized for the rapid graphical evaluation of the functional relevance of the observed changes in GE patterns. This will save the precious time of the expert biologists and bioinformatics specialists.
Pathways assembled into the KEGG database are displayed as semi-static objects that can be manipulated using tools like KGML and KEGG application programmable interface (API) [11, 12]. KEGG API provides a routine that highlights specified genes within the particular metabolic pathway (http://www.genome.jp/kegg/tool/color_pathway.html). Similar task may be also executed using G-language Genome Analysis Environment . Both approaches work on the pathway by pathway basis. Another tool, Pathway Express, calculates the pathway-wise impact of differentially expressed genes based on normalized fold change and depicts the pathways with differentially expressed genes . However, the fold-change approach and its associated standard t-test statistics usually produce severely over-fitted models. A number of recently developed approaches generate gene rankings dissociated from the fold change estimates [15, 16]. An analysis of these gene lists may benefit from the binary graphical mapping of upregulated and downregulated elements within the complete collection of pathway maps. Resulting graphical pictures may be helpful both as tool for a quick assessment of the functional relevance of a gene list and as a set of the snapshots easily convertible into the illustrative material for presentations or manuscript figures.
With this notion, here we present a web-based tool, KEGG Pathway Painter (KPP). KPP performs a batch painting of relevant pathways according to the uploaded lists of up-regulated and down-regulated genes in KEGG. KPP returns a set of images that give a holistic perspective to the functional importance of the change in the GE patterns revealed by a given high-throughput experiment and facilitate the extraction of the biological insights.
KPP was implemented using PERL/CGI. Pathways assembled into the KEGG database are displayed as semi-static objects that can be manipulated using tools like KGML (KEGG Markup Language) and KEGG API (Application Programming Interface). The API allows access to the resources stored in KEGG system in an interactive and user-friendly way (http://www.genome.jp/kegg/rest/).
KEGG Pathway Painter (KPP) accepts the up-regulated and down-regulated gene lists as two different text files containing the gene identifiers of any sequenced organism. Permitted identifiers include GenBank id, NCBI GENE id, NCBI GI accession, Unigene ID and Uniprot ID. Conversion of the gene identifiers, extraction of the corresponding pathway and their painting is performed by specific API routines. The KPP processes data through direct interface to the KEGG database, and therefore, the KPP painted pathways are always up-to date with reference to KEGG knowledgebase. In KPP, genes of interest can be also highlighted with user-specified foreground and background colors allowing easy visual differentiating of up- and down-regulated genes.
Mapped genes are automatically consolidated within each pathway. The number of the KPP returned pathways could be filtered by either the total number of the painted genes in a given pathway or the ratio of painted genes to the total number of genes in a given pathway. The chosen pathways passing the criteria on filter are color coded according to users' preferences. Users can browse through these high-resolution pathway images along with gene information and an archive of the painted pathways can also be saved for future reference. After completion of the query, the URL to the index of resulting output images is e-mailed to the user along with the job summary.
Results and Discussion
The motivation for the development of KPP came up from the idea to build a user-friendly, platform-independent and simple tool to visualize the placement of genes in their associated pathways. The simplicity of KPP is due to the acceptance of gene identifiers without reference to respective microarray platform. This isolation enhances its utility for the studies of the data from RealTime-PCR or medium-throughput platforms or even for validation of the various hypotheses concerning an involvement of the groups of genes in one or another biological process.
This utility of KPP was demonstrated by highlighting of cell cycle related events using the publicly available prostate carcinoma dataset (GDS1439)  from the NCBI GEO database (see Figure 1), by aiding the selection of the mutations and epigenetic events to be tested as a companion diagnostics of treatment susceptibility and resistance in non-small lung carcinoma patients (not shown) and an analysis of the host-associated risk factors associated with lack of sustained virological response (SVR) in various cohorts of HCV patients [18, 19].
In one of these examples, KPP-aided visual parsing the pathways encompassing molecular components relevant to HCV pathogenesis allowed to pinpoint the Janus kinase-signal transducers and activators of transcription signaling cascade as the major pathogenetic component responsible for not achieving SVR , a conclusion that was later confirmed in in vitro experiments with blocking antibodies, a pharmacological inhibitor, and siRNAs .
In another example, KPP allowed to visualize a sustained pattern of treatment-induced gene expression in patients carrying interferon/ribavirin-responding IL28B genotype C/C, while in patients with therapy-resistant IL28B T* genotype, the background pre-activation of interferon-dependent genes precluded further therapeutic boost . Thus, KPP provided a critical insight into the lower rate of SVR observed in these patients. Furthermore, KPP analysis revealed LI28B genotype independent role of SOCS1 in therapeutic response . This KPP-aided hypothesis was later investigated both in vitro experiments showing that SOCS1 acts as a suppressor of type I IFN function against HCV  and in serum samples interferon/ribavirin-treated Hepatitis C patients who demonstrated that methylation-based silencing of SOCS-1 is associated with better therapeutic response . Thus, KPP was indispensable in acquiring mechanistic insights into the differential therapeutic response in Hepatitis C infected patients.
The major fetching point of the KPP tool lies in its tight connection with the KEGG database, as this will allow for the pathway visualization of every sequenced organism. However this flexibility comes at the cost of possible KEGG-attributed delay of the data transfer, the resultant tool is substantially more convenient for the user than the tools embed into existing pathway analysis environment, for example, Cytoskape (http://www.cytoscape.org/). Another commonly used pathway parsing tool, Reactome Skypainter (http://www.reactome.org/), is restricted to underlying knowledge base and, therefore, limits the potential set of insights to be extracted.
It is important to note that the painting of individual pathways can be performed through by the KEGG website itself (http://www.genome.jp/kegg/), however, the practicality of KPP is in its comprehensive visual representation of up- and downregulated genes in the KEGG dataset as a whole. In other words, KPP allows one to extract immediate and visual insights about cumulative change in each pathway under scrutiny. Users can browse through high-resolution pathway images and download an archive of the painted pathways that may be used as figures for upcoming manuscripts.
In summary, KPP provides fast and comprehensive visualization of the global GE changes by consolidating a list of the color-coded candidate genes into the KEGG pathways.
- KEGG Pathway Painter
- Gene Expression
- KEGG Markup Language
- Application Programming Interface
- Sustainer Virological Response
Ganter B, Giroux CN: Emerging applications of network and pathway analysis in drug discovery and development. Curr Opin Drug Discov Devel. 2008, 11 (1): 86-94.
Chen G, Cairelli MJ, Kilicoglu H, Shin D, Rindflesch TC: Augmenting microarray data with literature-based knowledge to enhance gene regulatory network inference. PLoS Comput Biol. 2014, 10 (6): e1003666-10.1371/journal.pcbi.1003666.
Kolesnikov N, Hastings E, Keays M, Melnichuk O, Tang YA, Williams E, Dylag M, Kurbatova N, Brandizi M, Burdett T, Megy K, Pilicheva E, Rustici G, Tikhonov A, Parkinson H, Petryszak R, Sarkans U, Brazma A: ArrayExpress update-simplifying data submissions. Nucleic Acids Res. 2015, 43 (Database): D1113-6.
Barrett T, Troup DB, Wilhite SE, Ledoux P, Rudnev D, Evangelista C, Kim IF, Soboleva A, Tomashevsky M, Marshall KA, Phillippy KH, Sherman PM, Muertter RN, Edgar R: NCBI GEO: archive for high-throughput functional genomic data. Nucleic Acids Res. 2009, 37 (Database): D885-90. 10.1093/nar/gkn764.
Tseng GC, Ghosh D, Feingold E: Comprehensive literature review and statistical considerations for microarray meta-analysis. Nucleic Acids Res. 2012, 40 (9): 3785-99. 10.1093/nar/gkr1265.
Mayburd A, Baranova A: Knowledge-based compact disease models identify new molecular players contributing to early-stage Alzheimer's disease. BMC Syst Biol. 2013, 7: 121-10.1186/1752-0509-7-121.
Yang Z, Chen Y, Fu Y, Yang Y, Zhang Y, Chen Y, Li D: Meta-analysis of differentially expressed genes in osteosarcoma based on gene expression data. BMC Med Genet. 2014, 15: 80-
Carnielli CM, Winck FV, Paes Leme AF: Functional annotation and biological interpretation of proteomics data. Biochim Biophys Acta. 2015, 1854 (1): 46-54. 10.1016/j.bbapap.2014.10.019.
Tanabe M, Kanehisa M: Using the KEGG database resource. Curr Protoc Bioinformatics. 2012, Chapter 1 (Unit1.12):
Kanehisa M, Goto S, Sato Y, Furumichi M, Tanabe M: KEGG for integration and interpretation of large-scale molecular data sets. Nucleic Acids Res. 2012, 40 (Database): D109-14.
Kawashima S, Katayama T, Sato Y, Kanehisa M: KEGG API: A web service using SOAP/WDSL to access the KEGG system. Genome Informatics. 2003, 14: 673-674.
Klukas C, Schreiber F: Dynamic exploration and editing of KEGG pathway diagrams. Bioinformatics. 2007, 23 (3): 344-50. 10.1093/bioinformatics/btl611.
Arakawa K, Mori K, Ikeda K, Matsuzaki T, Kobayashi Y, Tomita M: G-language Genome Analysis Environment: a workbench for nucleotide sequence data mining. Bioinformatics. 2003, 19 (2): 305-6. 10.1093/bioinformatics/19.2.305.
Khatri P, Sellamuthu S, Malhotra P, Amin K, Done A, Draghici S: Recent additions and improvements to the Onto-Tools. Nucleic Acids Res. 2005, 33 (Web Server): W762-5. 10.1093/nar/gki472.
Simon R: Microarray-based expression profiling and informatics. Curr Opin Biotechnol. 2008, 19 (1): 26-9. 10.1016/j.copbio.2007.10.008.
Emmert-Streib F, Glazko GV, Altay G, de Matos Simoes R: Statistical inference and reverse engineering of gene regulatory networks from observational expression data. Front Genet. 2012, 3: 8-
Varambally S, Yu J, Laxman B, Rhodes DR, Mehra R, Tomlins SA, Shah RB, Chandran U, Monzon FA, Becich MJ, Wei JT, Pienta KJ, Ghosh D, Rubin MA, Chinnaiyan AM: Integrative genomic and proteomic analysis of prostate cancer reveals signatures of metastatic progression. Cancer Cell. 2005, 8 (5): 393-406. 10.1016/j.ccr.2005.10.001.
Birerdinc A, Afendy A, Stepanova M, Younossi I, Manyam G, Baranova A, Younossi ZM: Functional pathway analysis of genes associated with response to treatment for chronic hepatitis C. J Viral Hepat. 2010, 17 (10): 730-6. 10.1111/j.1365-2893.2009.01235.x.
Younossi ZM, Birerdinc A, Estep M, Stepanova M, Afendy A, Baranova A: The impact of IL28B genotype on the gene expression profile of patients with chronic hepatitis C treated with pegylated interferon alpha and ribavirin. J Transl Med. 2012, 10: 25-10.1186/1479-5876-10-25.
Zhang L, Jilg N, Shao RX, Lin W, Fusco DN, Zhao H, Goto K, Peng LF, Chen WC, Chung RT: IL28B inhibits hepatitis C virus replication through the JAK-STAT pathway. J Hepatol. 2011, 55 (2): 289-98. 10.1016/j.jhep.2010.11.019.
Shao RX, Zhang L, Hong Z, Goto K, Cheng D, Chen WC, Jilg N, Kumthip K, Fusco DN, Peng LF, Chung RT: SOCS1 abrogates IFN's antiviral effect on hepatitis C virus replication. Antiviral Res. 2013, 97 (2): 101-7. 10.1016/j.antiviral.2012.12.001.
Tseng KC, Chou JL, Huang HB, Tseng CW, Wu SF, Chan MW: SOCS-1 promoter methylation and treatment response in chronic hepatitis C patients receiving pegylated-interferon/ribavirin. J Clin Immunol. 2013, 33 (6): 1110-6. 10.1007/s10875-013-9903-4.
The authors express gratitude to the general support provided by College of Science, George Mason University, a State Contract 14.607.21.0098 dated November 27th, 2014 (Ministry of Science and Education, Russia) and by the Human Proteome Scientific Program of the Federal Agency of Scientific Organizations, Russia.
Availability of support data: The data supporting the results of this article are included within the article and on the publicly available website http://web.cos.gmu.edu/~gmanyam/kegg/kpp.html
Open access fees were covered by funds of College of Science, George Mason University.
This article has been published as part of BMC Systems Biology Volume 9 Supplement 2, 2015: Selected articles from the IX International Conference on the Bioinformatics of Genome Regulation and Structure\Systems Biology (BGRS\SB-2014): Systems Biology. The full contents of the supplement are available online at http://www.biomedcentral.com/bmcsystbiol/supplements/9/S2.
The authors declare that they have no competing interests.
All three authors contributed to the study design, interpretation of results and producing the manuscript. All authors read and approved the final manuscript.