I would like to visualize this interaction list. Is there an online/web based way to do it? Is there a way to analyze the data?
Assuming I have exported the list into lines such as:
2 A2M 348 APOE
Bioinformatics Stack Exchange is a question and answer site for researchers, developers, students, teachers, and end users interested in bioinformatics. It only takes a minute to sign up.Sign up to join this community
since I am not yet allowed to comment, I will have to pots this as an answer.
I guess Cytoscape (http://cytoscape.org/) would be a nice way to analyse such interaction lists/networks (still depends a bit on what you plan to do).
However, for this you will have to preprocess the downloadable excel sheet a bit.
python3 script.py <path-to-tsv-file> <delimeter for cytoscape>. The new format is written to std out on the console, therefore maybe you need to write the output to a new file. A sample call could look like
python3 testbench.py human_interactome_2017-06-12.tsv.csv ";". Make sure to add the delimeter in quotes.
import sys HIIIidx = 4 HII5idx = 5 venkatesanidx = 6 yuidx = 7 Litidx = 8 print(sys.argv) if len(sys.argv) <= 1: print("Usage: script.py file delimeter") delim = ';' inFile = sys.argv if len(sys.argv) > 2: delim = sys.argv def printInteraction(syma, symb, inter): print(syma + delim + inter + delim + symb) with open( inFile , 'r') as file: file.readline() for line in file: aline = line.strip('\n').split('\t') if aline[HIIIidx] == '1': printInteraction(aline, aline, 'HIII') if aline[HII5idx] == '1': printInteraction(aline, aline, 'HII5') if aline[venkatesanidx] == '1': printInteraction(aline, aline, 'VEN') if aline[yuidx] == '1': printInteraction(aline, aline, 'YU') if aline[Litidx] == '1': printInteraction(aline, aline, 'LIT')
Your list has two identifiers for the same node per line. In order to use it, you will need to change that. If you want to use the gene name (2nd and 4th fields, in your example), just run:
awk 'print $2,$4' netw.txt > netw.gr
If you want to use the Entrez geneIDs instead, run:
awk 'print $1,$3' netw.txt > netw.gr
Then, as others already mentioned, install
Cytoscape, launch it and import your gene list:
In the welcome screen, select "From Network File...":
Then, select your network file (
netw.gr in the example above) and choose "Advanced Options":
Set the delimiter to SPACE and uncheck "Use first line as column names":
Finally, click on the header of each column and set the first to "source node" and the second to "target node":
Now, click "OK" and you will have successfully imported your network into Cytoscape.
Usually interactions are represented in a graph, with edges representing the interactions (colors, thickness or values on the edges represent different properties of those edges).