R version: 4.2.1
Bioconductor version: 3.16
Cy3D (Cytoscape app): 1.1.3
KEGGscape (Cytoscape app): 0.9.1
Visualization of Trans-omic networks helps biological interpretation by illustrating pathways where the signals are transmitted (Gehlenborg et al., 2010).
To characterize signals that go across multiple omic layers, Yugi and colleagues have proposed a method for network visualization (Yugi et al., 2014) by stacking multiple 2D pathways in a 3D space.
The 3D network visualization was realized by VANTED (Rohn et al., 2012). However, the visualization relies on time-consuming manual operation. Here we propose transomics2cytoscape, an R package that automatically creates 3D network visualization in combination with Cytoscape (Shannon, 2003), Cy3D App, and Cytoscape Automation (Otasek et al., 2019).
This package requires Cytoscape to be installed and you need to run Cytoscape before running the following R code.
transomics2cytoscape has 2 main functions
Below is a diagram of the transomics2cytoscape workflow.
create3Dnetwork has 3 arguments.
The 1st one is a directory path where you put the network files to be layered in 3D space. The 2nd one is a file path of TSV for the Z-axis layout of the network files (called “Layer definition file”). The last one is a file path of XML used to style Cytoscape.
suid <- create3Dnetwork(networkDataDir, networkLayers, stylexml)
createTransomicEdges has 2 arguments.
The 1st one is the SUID of the network created by
The 2nd one is a file path of TSV for the transomic interactions
(called “Transomic interaction file”).
suid <- createTransomicEdges(suid, layer1to2)
transomics2cytoscape can layer all the networks that Cytoscape can import.
You need to put these files in the directory of the 1st argument of
You don’t need to put files for the KEGG pathway.
For KEGG pathway, you can import the network just by writing the KEGG
pathway ID in the “Layer definition file” described later.
“Layer definition file” is a TSV file for the Z-axis layout of the network files.
A file that defines network layer index and the Z-height of the network in 3D space. The format is as follows.
kinase-layer rno04910 1500 false enzyme-layer rno01100 900 true metabolite-layer rno01100 1 false
The 1st column is the network layer index.
This information is added to the node table column
The 2nd column is the KEGG pathway ID or the network file name in the directory
of the 1st argument of
You don’t need to prepare a network file for the KEGG pathway.
You can import the KEGG pathway simply by writing the KEGG pathway ID.
The 3rd column is the Z-height of the network.
The fourth column specifies whether you want to connect the interaction to the “edge” of the network in that row. If this is true, transomics2cytoscape will create a node at the midpoint of all edges of the network. This is useful when you want to represent an interaction that activates and inactivates a reaction in the network. For example, KEGG global metabolic pathway network does not have an enzyme node, unlike a normal metabolic pathway. This column should be set to “true” in such cases.
A Cytoscape style file. For more information about Cytoscape style file, see the Cytoscape user manual. Note that you can only use style properties that are supported by Cy3D.
“Trans-omic interaction file” is a TSV file that defines the edges that connect the different network layers. The format is as follows.