- Open Access
NetworkPainter: dynamic intracellular pathway animation in Cytobank
BMC Bioinformaticsvolume 16, Article number: 172 (2015)
High-throughput technologies such as flow and mass cytometry have the potential to illuminate cellular networks. However, analyzing the data produced by these technologies is challenging. Visualization is needed to help researchers explore this data.
We developed a web-based software program, NetworkPainter, to enable researchers to analyze dynamic cytometry data in the context of pathway diagrams. NetworkPainter provides researchers a graphical interface to draw and “paint” pathway diagrams with experimental data, producing animated diagrams which display the activity of each network node at each time point.
NetworkPainter enables researchers to more fully explore multi-parameter, dynamical cytometry data.
Cellular signaling is enormously complex, arising from dynamic interactions among thousands of molecules. Understanding signaling at the molecular level therefore requires comprehensive, dynamic measurements of individual cells . Flow and mass cytometry are two of the most informative technologies, capable of measuring multiple properties of thousands of individual cells per second . Fluorescence-based flow cytometry enables up to 12 simultaneous measurements . Mass cytometry promises as many as 100 . Flow cytometry has provided crucial information about DNA copy number , protein expression , and phenotypes . Mass cytometry has been used to study leukemia , hematopoiesis , and immune signaling .
Recently, Kotecha and colleagues developed Cytobank, a web-based platform for storing, exploring, and sharing cytometry data [http://www.cytobank.org; . Cytobank stores primary flow and mass cytometry data and provides annotation, gating, and visualization tools. Cytobank also enables researchers to easily share data and analysis with collaborators. Cytobank’s vertical integration enables researchers to drill down from publication quality figures to the underlying analysis, processing, filtering, and ultimately to the raw data, making data analysis fully transparent and reproducible. Additional tools are needed to help researchers analyze dynamic data.
Visualization is an effective technique for investigating complex data [12-14]. Flow cytometry researchers commonly use several graphical analysis software programs and packages including FlowJo [http://www.flowjo.com], flowViz , flowCore , iFlow , and SPICE  to browse and analyze cytometry data. Recently, researchers have developed several new techniques for visually analyzing mass cytometry data. Qui et al. developed spanning-tree progression analysis of density-normalized events (SPADE) to cluster and visualize measurements across multiple cell types . Amir et al. developed viSNE to analyze mass cytometry data using visuals similar to scatter plots .
In addition, several tools are available for visualizing temporal data in the context of the underlying molecular network. Secrier and Schneider have extensively reviewed these dynamical data visualization tools . These tools use three main strategies to visualize temporal data on top of networks. First, several tools including TVNViewer  and VisANT  can depict temporal dynamics using multiple vertically stacked or spatially tiled layers corresponding to different time points. Each layer displays the same network structure and colors and/or boldens the observed nodes and edges according to their observed activity. Second, several tools including the Cytoscape  plugins MODAM , MultiColoredNodes , SpotXplore , and VistaClara , and VANTED  place small sub-plots next to each observed node and edge to indicate its observed temporal dynamics. Third, several tools including BioTapestry , DynNetwork [http://code.google.com/p/dynnetwork], and the Pathway Tools Cellular Overview  animate network diagrams by coloring or boldening each observed node and edge at each time point. To optimally display large networks and high-dimensional data, these existing network visualization tools depict networks using streamlined node-link or circuit diagrams.
However, cell signaling researchers often prefer classical textbook-style pathway diagrams for medium size networks and multi-parameter data. Unfortunately, none of the textbook-style pathway drawing tools such as ChemBioDraw [http://www.cambridgesoft.com/software/overview.aspx] and Cell Illustrator  are capable of displaying temporal dynamics. These software programs are also proprietary and cannot be integrated or embedded within other software or databases. CellDesigner, which provides a graphical style intermediate between the streamlined node-link and classical styles, also cannot display dynamics directly on top of a network .
We developed a web-based software program, NetworkPainter, to enable cell signaling researchers to visualize and communicate multi-parameter data in the context of classical textbook-style pathway diagrams. NetworkPainter enables researchers to draw and “paint” classical pathway diagrams with data, creating animated diagrams. Optionally, NetworkPainter also displays heatmaps next to each observed node to indicate its underlying activity distribution across different experimental conditions or patients. Furthermore, we have integrated NetworkPainter with Cytobank to encourage experimentalists to try this visualization without any complicated importing or exporting. In addition, we developed a standalone version of NetworkPainter available at http://covert.stanford.edu/networkpainter which is capable of visualizing any dynamic data. Both versions also enable researchers to share diagrams with collaborators.
NetworkPainter was motivated by the unique needs of flow- and mass-cytometry researchers to visualize multi-parameter data. NetworkPainter fills a dynamical data visualization need unmet both by existing network visualization software programs which are optimized for larger networks and by existing classical pathway drawing tools which are not capable of displaying temporal dynamics. NetworkPainter combines the capabilities of classical pathway drawing tools with the capabilities of network visualization tools to visualize dynamics through animation and multiple sub-plots. We believe that NetworkPainter can help scientists interpret and communicate multi-parameter dynamic biological data.
Here, we describe the features and implementation of NetworkPainter. We present a mass cytometry timecourse of the human peripheral blood mononuclear cell (PBMC) immune signaling network  as an example use case. We conclude by discussing our future plans for NetworkPainter.
NetworkPainter provides researchers a web-based graphical interface to visually analyze cytometry data. NetworkPainter only requires a web browser with Adobe Flash Player. Below we briefly describe how to use NetworkPainter. Additional file 1 provides further information about how to use the software.
First, researchers use the graphical interface to draw a pathway diagram, formalizing their prior biological knowledge about their pathway, including its molecular components and their interactions. Users simply drag and drop molecules to add them to the diagram and assign them to subcellular compartments. Users can draw arrows by selecting two molecules and selecting “Draw arrow” from the right-click context menu. Users have full control over the graphical appearance of each molecule and compartment including its color, shape, and size. Twenty-six shapes are available including polygons, as well as several commonly used graphical representations of DNA, RNA, and protein. Figure 1 depicts a PBMC immune signaling diagram created with NetworkPainter.
To help users create visually pleasing diagrams, NetworkPainter provides an automatic layout tool to improve the arrangement and spacing of molecules in diagrams. In addition, users can import diagrams from the KEGG PATHWAY database  or from diagrams published by other users. Furthermore, NetworkPainter provides basic editing features including cut, copy, paste, undo, redo, auto save, find, zoom, and pan.
Second, users upload and annotate experimental data. Cytobank users use Cytobank’s existing graphical interface to create an experiment, upload flow or mass cytometry data, and annotate the individual, cell type, condition, dosage, channel, and timepoint of each measurement. Standalone users can upload data using the JSON format described in the online help.
Third, users link diagram nodes with experimental measurements and perturbations. Researchers use drop-down lists to select the experimental channel corresponding to each observed signaling molecule. Optionally, researchers can also use drop-down lists to enter the enhansive or repressive effect of each experimental perturbation.
Fourth, researchers can specify how to display the experimental individual, cell type, condition, and dosage dimensions of their data. Researchers can select two dimensions to display in small heatmaps below each experimentally linked signaling molecule. These heatmaps display the measurement of the corresponding experimental channel for each value of the two selected heatmap dimensions at each timepoint, averaged over the unselected heatmap dimensions. Optionally, NetworkPainter can cluster and order the selected heatmap dimensions using hierarchical agglomerative clustering  and optimal leaf ordering . A legend and tooltips indicate the data plotted in each heatmap cell. These heatmaps are designed to help researchers visualize more of their data simultaneously, as well as compare measurements across individuals, cell types, conditions, and dosages.
Fifth, researchers can specify the averaging algorithm NetworkPainter uses to collapse the experimental data displayed in the signal molecules and heatmaps. Users can select among all of the statistics (e.g. mean, median, min, max), equations (e.g. raw, log, fold, arcsinh, difference), and controls (e.g. first row, first column, table min, table max) available in Cytobank’s illustrations. Researchers can also select the color mapping between the experimentally observed values and painted colors. The online help lists all of the available statistics, equations, controls, and colormaps.
Sixth, researchers can use the playback controls at the bottom-left of the diagram to dynamically paint each experimentally perturbed or observed signaling molecule over the measured timecourse. After clicking the play button, NetworkPainter colors each experimentally linked signaling molecule according to the measured value of the corresponding experimental channel at the current timepoint and the selected colormap, averaged over all other experimental dimensions. NetworkPainter interpolates timepoints to produce smooth animations. Unobserved nodes are colored grey and optionally, the compartments and membranes are colored grey to emphasize the dynamics of the observed nodes. Optionally, small heatmaps below the signaling molecules enumerate the values of the experimental channel across the two selected heatmap dimensions. Users can also control the playback speed and looping. Figure 2 depicts eight static snapshots of a human PBMC immune signaling diagram painted by NetworkPainter with the first reported mass cytometry timecourse containing fourteen measured signaling nodes . An interactive, animated version of Figure 2 is available at http://covert.stanford.edu/networkpainter/KarrEtAl2014Fig2. The case study below describes how NetworkPainter can be used to analyze this data.
Seventh, researchers can save diagrams to either the Cytobank or standalone server for later use and/or sharing with collaborators. Cytobank users can give collaborators permission to view or edit diagrams by granting permission to their associated experiments. Standalone users can use a simple graphical interface to grant read, write, or administer privileges to other users. Once collaborators are granted permission, they can edit a diagram. However, each diagram can only be edited by a single user at a time. Both Cytobank and standalone users can also publish diagrams to share them with all users.
Lastly, researchers can export diagrams to several image (GIF, JPEG, PNG, SVG) and animation (GIF, SWF) formats for use in papers, presentations, and websites. NetworkPainter also exports diagrams to the Biological Pathway Exchange (BioPax), CellML, and Systems Biology Markup Language (SBML) standards for use with other network analysis software programs, as well as to a JSON format which can be subsequently imported back into NetworkPainter. In addition, NetworkPainter can generate MATLAB scripts for Boolean simulations.
NetworkPainter was implemented as a web-based software program to enable platform independence, no installation, collaboration, and cloud-based computation. NetworkPainter is composed of a web-based graphical user interface, and a back-end server. The graphical interface provides users a diagram editor, paints signaling molecules with experimental data, and exports diagrams. The back-end server stores diagrams for later use, coordinates diagram permissions, helps export animated diagrams, and helps import diagrams from pathway databases.
The user interface was implemented using Adobe Flex [http://www.adobe.com/products/flex.html] using the Degrafa declarative graphics framework [http://code.google.com/p/degrafa]. Animation was implemented using Tweener [http://code.google.com/p/tweener]; automatic graph layout was implemented using GraphViz ; graphical export was implemented using Inkscape [http://www.inkscape.org]; and animation export was implemented using Adobe Flex.
The Cytobank server was implemented using JRuby [http://jruby.org], MySQL [http://www.mysql.com], Apache [http://httpd.apache.org], and Tomcat [http://tomcat.apache.org]. The standalone server was implemented using PHP [http://www.php.net], MySQL, and Apache. The servers and user interfaces communicate using JSON.
The NetworkPainter source code is freely available at http://github.com/CovertLab/NetworkPainter, including all of the code for the user interface and all of the code for the standalone server.
Case Study: mass cytometry cell signaling timecourse
We illustrate the functionality of NetworkPainter using the first reported mass cytometry timecourse of human peripheral blood mononuclear cell (PBMC) signaling . PBMCs are a diverse population of cells with distinct functions including B cells, T cells, NK cells, and macrophages. They are critical to both the innate and adaptive immune systems, helping the body identify and fight pathogens including viruses and bacteria.
Bodenmiller et al. used mass cytometry to investigate how PBMC subpopulations differentially respond to stimuli and drugs. They measured the effects of twelve stimuli on fourteen signaling nodes and ten surface markers in fourteen PBMC cell types at eight timepoints. In total, Bodenmiller et al. collected data from over 2,000 conditions.
Using SPADE, Bodenmiller et al. found complex dynamics, as well as significant cell type variability in the phosphorylation responses of the fourteen observed signaling nodes. In particular, they found that LPS rapidly induces p38, ERK, and NF-κB within 15–30 min in monocytes, followed by slower S6 induction at 2 h. In contrast, they also found that LPS slowly induces STAT3, STAT5, and ITK after 2 h in T cells and NK cells, and STAT1 in B cells after 4 h. Bodenmiller et al. attributed these differences in physiology to monocyte IL-6 secretion. Bodenmiller concluded that mass cytometry can be used to rapidly screen drug candidates against multiple cell types.
We used NetworkPainter to reexamine the Bodenmiller et al. dataset. First, we used NetworkPainter to draw a signaling diagram containing the fourteen observed signaling molecules (Figure 1). Next, we linked the diagram nodes to the observed channels, and used NetworkPainter to color the diagram with mass cytometry measurements (Figure 2). This enabled us to quickly contextualize, browse, and analyze the Bodenmiller et al. dataset. The animated diagram illustrated the same temporal signaling profile and cell type differences among monocytes, B cells, T cells, and NK cells reported by Bodenmiller et al. Furthermore, NeworkPainter's simple graphical interface allowed us to explore alternate network architectures and explanations for the observed dynamics and cell type variability. This example highlights the utility of NetworkPainter to help researchers quickly discern patterns from complex cytometry data.
Results and Discussion
NetworkPainter is a web-based program for visualizing dynamic cytometry data in the context of animated pathway diagrams. NetworkPainter provides a simple graphical interface for drawing and painting pathway diagrams with experimental cytometry data either stored in the Cytobank platform or uploaded to the standalone version. NetworkPainter stores diagrams on its server to enable users to easily share diagrams with collaborators. NetworkPainter exports diagrams to a variety of image and animation formats, enabling users to create high-quality graphics for papers, presentations, and websites.
Overall, high-throughput technologies such as mass cytometry have the potential to provide valuable insights into the molecular circuits that govern biological behavior and disease. Data visualization software is needed to help researchers explore and analyze the large data sets created by these technologies. We believe that NetworkPainter is a valuable tool for analyzing multi-parameter dynamical data. Furthermore, we believe that NetworkPainter is a valuable communication tool. NetworkPainter can help researchers communicate complex data to the broader scientific community in a simple, intuitive graphical format.
Availability and requirements
Project name: NetworkPainter
Operating system(s): Platform independent
Programming language: Actionscript, Ruby, PHP
Other requirements: Web browser, Flash player
License: Attribution Assurance License
Any restrictions to use by non-academics: None
NetworkPainter is available through the Cytobank cytometry platform at http://www.cytobank.org/networkpainter.html to analyze data stored in Cytobank. NetworkPainter is also freely available at http://covert.stanford.edu/networkpainter for use with any dynamic data. Source code is available at http://github.com/CovertLab/NetworkPainter under the Attribution Assurance License.
peripheral blood mononuclear cell
PHP: hypertext preprocessor
Przytycka TM, Singh M, Slonim DK. Toward the dynamic interactome: it's about time. Brief Bioinform. 2010;11(1):15–29.
Chattopadhyay PK, Roederer M. Cytometry: today’s technology and tomorrow’s horizons. Methods. 2012;57(3):251–8.
Perfetto SP, Chattopadhyay PK, Roederer M. Seventeen-colour flow cytometry: unravelling the immune system. Nat Rev Immunol. 2004;4(8):648–55.
Ornatskya O, Bandura D, Baranova V, Nitza M, Winnika MA, Tanner S. Highly multiparametric analysis by mass cytometry. J Immun Meth. 2010;361(1–2):1–20.
Darzynkiewicz Z. Critical aspects in analysis of cellular DNA content. Curr Protoc Cytom. 2011;56:7.2.1–8.
Sigal A, Danon T, Cohen A, Milo R, Geva-Zatorsky N, Lustig G, et al. Generation of a fluorescently labeled endogenous protein library in living human cells. Nat Protoc. 2007;2(6):1515–27.
Biancotto A, Fuchs JC, Williams A, Dagur PK, McCoy Jr JP. High dimensional flow cytometry for comprehensive leukocyte immunophenotyping (CLIP) in translational research. J Immunol Methods. 2011;363(2):245–61.
Bandura DR, Baranov VI, Ornatsky OI, Antonov A, Kinach R, Lou X, et al. Mass cytometry: technique for real time single cell multitarget immunoassay based on inductively coupled plasma time-of-flight mass spectrometry. Anal Chem. 2009;81(16):6813–22.
Bendall SC, Simonds EF, Qiu P, Amir e-AD, Krutzik PO, Finck R, et al. Single-cell mass cytometry of differential immune and drug responses across a human hematopoietic continuum. Science. 2011;332(6030):687–96.
Bodenmiller B, Zunder ER, Finck R, Chen TJ, Savig ES, Bruggner RV, et al. Multiplexed mass cytometry profiling of cellular states perturbed by small-molecule regulators. Nat Biotechnol. 2012;30(9):858–67.
Chen TJ, Kotecha N. Cytobank: providing an analytics platform for community cytometry data analysis and collaboration. Curr Top Microbiol Immunol. 2014;377:127–57.
O'Donoghue SI, Gavin AC, Gehlenborg N, Goodsell DS, Hériché JK, Nielsen CB, et al. Visualizing biological data–now and in the future. Nat Meth. 2010;7 Suppl 3:S2–4.
Gehlenborg N, O'Donoghue SI, Baliga NS, Goesmann A, Hibbs MA, Kitano H, et al. Visualization of omics data for systems biology. Nat Meth. 2010;7 Suppl 3:S56–68.
Secrier M, Schneider R. Visualizing time-related data in biology, a review. Brief Bioinform. 2013;15(5):771–82.
Sarkar D, Le Meur N, Gentleman R. Using flowViz to visualize flow cytometry data. Bioinformatics. 2008;24(6):878–879
Hahne F, LeMeur N, Brinkman RR, Ellis B, Haaland P, Sarkar D, et al. FlowCore: a bioconductor package for high throughput flow cytometry. BMC Bioinformatics. 2009;10:106.
Lee K, Hahne F, Sarkar D, Gentleman R. iFlow: A graphical user interface for flow cytometry tools in bioconductor. Adv Bioinformatics. 2009;2009(10):38–39.
Roederer M, Nozzi JL, Nason MC. SPICE: exploration and analysis of post-cytometric complex multivariate datasets. Cytometry A. 2011;79(2):167–74.
Qiu P, Simonds EF, Bendall SC, Gibbs Jr KD, Bruggner RV, Linderman MD, et al. Extracting a cellular hierarchy from high-dimensional cytometry data with SPADE. Nat Biotechnol. 2011;29(10):886–91.
Amir e-AD, Davis KL, Tadmor MD, Simonds EF, Levine JH, Bendall SC, et al. viSNE enables visualization of high dimensional single-cell data and reveals phenotypic heterogeneity of leukemia. Nat Biotechnol. 2013;31(6):545–52.
Curtis RE, Yuen A, Song L, Goyal A, Xing EP. TVNViewer: an interactive visualization tool for exploring networks that change over time or space. Bioinformatics. 2011;27(13):1880–1.
Hu Z, Chang YC, Wang Y, Huang CL, Liu Y, Tian F, et al. VisANT 4.0: Integrative network platform to connect genes, drugs, diseases and therapies. Nucleic Acids Res. 2013;41(Web server issue):W225–31.
Smoot ME, Ono K, Ruscheinski J, Wang PL, Ideker T. Cytoscape 2.8: new features for data integration and network visualization. Bioinformatics. 2011;27(3):431–2.
Enjalbert B, Jourdan F, Portais JC. Intuitive visualization and analysis of multi-omics data and application to Escherichia coli carbon metabolism. PLoS One. 2011;6(6), e21318.
Warsow G, Greber B, Falk SS, Harder C, Siatkowski M, Schordan S, et al. ExprEssence–revealing the essence of differential experimental data in the context of an interaction/regulation net-work. BMC Syst Biol. 2010;4:164.
Westenberg MA, Roerdink JB, Kuipers OP, van Hijum SA. SpotXplore: a cytoscape plugin for visual exploration of hotspot expression in gene regulatory networks. Bioinformatics. 2010;26(22):2922–3.
Kincaid R, Kuchinsky A, Creech M. VistaClara: an expression browser plug-in for cytoscape. Bioinformatics. 2008;24(18):2112–4.
Klukas C, Schreiber F. Integration of -omics data and networks for biomedical research with VANTED. J Integer Bioinform. 2010;7(2):112.
Longabaugh WJ. BioTapestry: a tool to visualize the dynamic properties of gene regulatory networks. Methods Mol Biol. 2012;786:359–94.
Latendresse M, Karp PD. Web-based metabolic network visualization with a zooming user interface. BMC Bioinformatics. 2011;12:176.
Nagasaki M, Saito A, Jeong E, Li C, Kojima K, Ikeda E, et al. Cell illustrator 4.0: a computational platform for systems biology. Stud Health Technol Inform. 2011;162:160–81.
Matsuoka Y, Funahashi A, Ghosh S, Kitano H. Modeling and simulation using cell designer. Methods Mol Biol. 2014;1164:121–45.
Kanehisa M, Goto S. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28(1):27–30.
Eisen MB, Spellman PT, Brown PO, Botstein D. Cluster analysis and display of genome-wide expression patterns. Proc Natl Acad Sci U S A. 1998;95(25):14863–8.
Bar-Joseph Z, Gifford DK, Jaakkola TS. Fast optimal leaf ordering for hierarchical clustering. Bioinformatics. 2001;17 Suppl 1:S22–9.
Ellson J, Gansner E, Koutsofios L, North SC, Woodhull G. Graphviz– open source graph drawing tools. Lect Note Comput Sci. 2002;2265:483–4.
We thank Tiffany Chen, Jennifer Davis, Christina Hall, Angela Landrigan, and Chad Rosenberg for helping integrate NetworkPainter with Cytobank; and Bernd Bodenmiller and Rachel Finck for providing the case study data.
This work was supported by National Institute of Health grants (1P50GM107615, 5DP1LM01150-05, CA125994-01A1) and an Allen Foundation Distinguished Investigator Award to MWC; a National Defense Science and Engineering Graduate Fellowship, a National Science Foundation Graduate Fellowship (DGE-1147470), a Stanford Graduate Fellowship and a James S. McDonnell Foundation Postdoctoral Fellowship Award in Studying Complex Systems to JRK; a Stanford Electrical Engineering Graduate Fellowship to HG; and a National Institutes of Health K99/R00 grant (CA143231) to JMI. The Cytobank project has been funded in part by National Institutes of Health grants and contracts (HHSN268201300037C, GM096579, AI094929).
NK and JMI are founders of Cytobank Inc. SB is an employee of Cytobank Inc. All other authors declare that they have no competing interests.
JRK, JMI, and MWC conceived the project. JRK, HG, and EYC developed the NetworkPainter user interface and standalone server. JRK and SLB integrated NetworkPainter with Cytobank. NK supervised the integration of NetworkPainter with Cytobank. JRK developed the case study. MC supervised all aspects of the project. JRK wrote the manuscript. All authors approved the final manuscript.
NetworkPainter user instructions.