In this hands-on session you will go through some basic applications in network biology.
You will learn
- Where to find biological interaction data?
- How to use Cytoscape?
- How to visualize experimental data on a network?
- How to open a pathway as a network?
- How to extend a network with regulatory information?
The questions in part 2, 3 and 4 are meant to guide you through the practical session. You will receive all the answers at the end of the course. Feel free to ask the instructors if you need help with any steps in the practical.
Part 1: Learn how to use Cytoscape
We will go through this part together with you to show you how Cytoscape works.
Cytoscape (www.cytoscape.org) is a popular and commonly used open source tool for network biology. The program has lots of different functionalities making it possible to perform a variety of network analysis, from the basics to advanced analysis. Since it is as open community project, many different researchers and developers contribute to Cytoscape. This is demonstrated by the wide variety of Cytoscape apps.
For the following two parts you will need different Cytoscape functionality that is explained in one of their tutorials online. You don’t need to go through the tutorials when you followed the introduction session.
Part 2: Using networks to analyze transcriptomics datasets
In this session we will continue using the dataset that you analyzed in the pathway analysis hands-on.
You can find the “dataset.xlsx” file in the network-analysis folder (in the zip file).
We will use the statistical results of the dataset from the paper “Gene Expression Patterns Induced by HPV-16 L1 Virus-Like Particles in Leukocytes from Vaccine Recipients” by Garcia-Pineres et al. published in The Journal of Immunology in 2008. We compare gene expression before and after vaccination.
In the next step you will further investigate how the top 10 up-regulated genes are connected with each other. You will use GeneMania (www.genemania.org) to find the interactions between the genes.
- Copy the gene symbols of the genes (external_gene_id) in answer 2A and use them as the query on GeneMania.
- Click on ‘Show advanced options’.
- Make sure that only ‘Genetic interactions, Pathway and Physical Interactions’ are selected.
- Make sure that you increase the number of related genes to 50.
Go to the “Functions” tab on the right side. GeneMania performs an over-repesentation analysis and provides a list of GO terms.
In the next step, we want to import this network in Cytoscape to visualize our experimental data on the gene nodes. Therefore, we need to download the network as text.
Then, we go to Cytoscape and import a network from the downloaded file (genemania_network.txt) – File → Import → Network → File.
Now, we want to import and visualize the gene expression dataset (network-analysis/dataset.xlsx) in the network. Follow the same instructions as in the introduction tutorial for importing data.
Make sure that you link the correct ID column to the correct key column:
Part 3: Biological pathways as networks
Open a new session in Cytoscape: File → New → Session.
In this part, we are going to open a pathway from WikiPathways as a network and extend it with regulatory interactions. We will use two different Cytoscape apps: WikiPathways and CyTargetLinker. The installation instructions (sent a week ago) explain how to install apps in Cytoscape.
First we are going to open the Type II Interferon Signaling Pathway in Cytoscape. You can find the GPML file in the network-analysis folder.
- File → Import → Network → File…
- Select network-analysis/WP619_71168.gpml
- Make sure you select “Network” in the first properties dialog:
Use the NetworkAnalyzer to calculate the different network properties that were discussed during the lecture. Go to Tools → NetworkAnalyzer → Analyze Network. Treat the network as undirected.
Next, we are going to use intuitive visualization of the node degree to easily find hub proteins in the network.
Part 4: Regulatory interactions
In this part, we will use the pathway network from Part 3. Now, we will extend the network with regulatory information. We will investigate which microRNAs and drugs are regulating which proteins in the pathway. For this step, we will use the CyTargetLinker app for Cytoscape.
The regulatory information is provided as networks, called Regulatory Interaction Networks (RegINs). For this practical the miRTarBase and the DrugBank RegIN are stored in the downloaded zip file in the directory “network-analysis/RegINs“.
We are going to extend the Type II Interferon Signaling Pathway with regulatory information.
- Click on Apps → CyTargetLinker → Extend Network
- Select User Network → choose the network name of the pathway network
- Select Network Attribute → choose the name of the column containing a biological identifier (“GeneID”)
- Select Datasource Networks → Browse for the directory which contains the RegINs (network-analysis/RegINs)
- Select Direction → “Add regulators”
For the next questions, it is easier if you hide the microRNA-target interactions from miRTarBase. Go to the CyTargetLinker tab on the left and use the drop-down box next to the miRTarBase RegIN to hide the interactions.
For the next questions, it is easier if you hide the drug-target interactions from DrugBank. Go to the CyTargetLinker tab on the left and use the drop-down box next to the DrugBank RegIN to hide the interactions. Make sure that the miRTarBase interactions are shown!