Part 4: Visualization of transcriptomics expression datasets in the enriched AOP network part 1

The AOP project ► Key objective 2

Author: Shakira Agata

This Jupyter Notebook describes the steps needed for the mapping of transcriptomics datasets: GSE109565, E-MEXP-2599 and E-MEXP-3583 in the enriched AOP network. This notebook is subdivided into the following seven sections:

• Section 1: System preparation
• Section 2: Retrieval of molecular inflammation-process related AOP network
• Section 3: Adaptation of gene node color of molecular inflammation-process related AOP network
• Section 4: Mapping of dataset:GSE109565
  • Section 4.1 PCB126 concentration 1
  • Section 4.2 PCB126 concentration 2
  • Section 4.3 PCB126 concentration 3
  • Section 4.4 Roundup
• Section 5: Mapping of dataset:E-MEXP-2599
  • Section 5.1 CdCl2 exposure time 1
  • Section 5.2 CdCl2 exposure time 2
  • Section 5.3 CsA exposure time 1
  • Section 5.4 CsA exposure time 2
  • Section 5.5 Diquat dibromide exposure time 1
  • Section 5.6 Diquat dibromide exposure time 2
• Section 6: Mapping of dataset:E-MEXP-3583
  • Section 6.1 Ag+ exposure time 1
  • Section 6.2 Ag+ exposure time 2
  • Section 6.3 AgNP exposure time 1
  • Section 6.4 AgNP exposure time 2
• Section 7: Metadata

Section 1: System preparation

In this section, you will import the required packages and tools you need for this Jupyternotebook.

step 1: You imported Pandas, Py4cytoscape and style mapping functions of Py4cytoscape.

import pandas as pd

import py4cytoscape as p4c
p4c.cytoscape_ping()
p4c.cytoscape_version_info()

You are connected to Cytoscape!





{'apiVersion': 'v1',
 'cytoscapeVersion': '3.10.1',
 'automationAPIVersion': '1.9.0',
 'py4cytoscapeVersion': '1.9.0'}

from py4cytoscape import get_node_color
from py4cytoscape import set_node_color_mapping
from py4cytoscape import gen_node_color_map
from py4cytoscape import set_edge_color_default
from py4cytoscape import set_node_color_default
from py4cytoscape import set_edge_source_arrow_shape_default
from py4cytoscape import set_edge_target_arrow_shape_default
from py4cytoscape import get_arrow_shapes
from py4cytoscape import get_edge_target_arrow_shape
from py4cytoscape import set_edge_target_arrow_shape_mapping
from py4cytoscape import gen_edge_arrow_map
from py4cytoscape import select_nodes
from py4cytoscape import get_table_value
from py4cytoscape import get_network_suid
from py4cytoscape import clear_selection
from py4cytoscape import set_node_color_bypass
from py4cytoscape import set_edge_color_bypass
from py4cytoscape import set_edge_target_arrow_color_default
from py4cytoscape import set_node_size_bypass
from py4cytoscape import create_subnetwork

Section 2: Retrieval of molecular inflammation-process related AOP network

In this section, you will open your Cytoscape molecular inflammatory-process related AOP network using the py4cytoscape session function. This allows you to open the complete network with previously established style settings and bypass settings.

step 2: You open the session you saved in the previous Jupyternotebook.

p4c.open_session('Agata,S.-Part4-Complete Molecular inflammation-process related AOP network.cys')

Opening C:\Users\shaki\Downloads\Agata,S.-Part4-Complete Molecular inflammation-process related AOP network.cys...





{}

Section 3: Adaptation of gene node color of molecular inflammation-process related AOP network

In this section, you will change the node color of genes and adapt the style for easier intepretation of the upcoming results. This is needed in preparation for the mapping of transcriptomics datasets. These datasets may contain genes that are not present in the build AOP network and so therefore should receive a distinct color to correctly inform user.

step 3: You can change the style with the following command.

style_name = "default"
defaults = {'NODE_SHAPE': "ELLIPSE", 'NODE_SIZE': 20, 'EDGE_TRANSPARENCY': 140, 'NODE_LABEL_POSITION': "C,C,c,0.00,0.00"}
nodeLabels = p4c.map_visual_property('node label', 'name', 'p') 
edgeWidth = p4c.map_visual_property('edge width', 'weight', 'p') 
arrowShapes = p4c.map_visual_property('Edge Target Arrow Shape','interaction', 'd')
p4c.create_visual_style(style_name, defaults, [nodeLabels, edgeWidth])
p4c.set_visual_style(style_name)

{'message': 'Visual Style applied.'}

Section 4: Mapping of dataset: GSE109565

In this section, you will map the transcriptomics expression data of datasets:GSE109565. GSE109565 focused on the alteration in gene expression profiles of HepaRG cells following exposure to three concentrations of polychlorinated biphenyl (PCB) 126 and one concentration of Roundup.

4.1 PCB concentration 1

step 4: You first import the expression table into a new dataframe.

PCB1_DEG= pd.read_csv('PCB concentration 1-GSE109565.top.table.tsv',sep='\t')
PCB1_DEG

step 5: You will now integrate the expression table into the nodetable of the AOP network. This will be done using the function: p4c.load_table_data_from_file where you select ‘True’ for the second variable as the first row of the PCB1_DEG table has the needed column names. You also select ’node’ for table and ‘CTL.GeneID’ for table_key_column as you want the expression data to be matched to the Gene ID’s that were added by CyTargetLinker extension in the previous notebook.

p4c.load_table_data_from_file('PCB concentration 1-GSE109565.top.table.tsv', first_row_as_column_names=True,table='node', table_key_column='CTL.GeneID')

{'mappedTables': [171586, 171624]}

step 6: You define the style for the mapping so that the expression values (log2FoldChange) are mapped to the gene nodes. This is done by first retrieving the log2FoldChange column.

Log2Foldchange_column = p4c.get_table_columns(table='node', columns='log2FoldChange')

Log2Foldchange_column

step 7: This is followed by definition of the color scheme so that low and high expression values receive distinct color.For the color scheme, you will use the following colors:

• Low expression value (minimum) = blue node color
• No expression value = white node color
• High expression value (maximum) = red node color

This color scheme was also described in the official py4cytoscape documentation.

Blue_expression_color= Log2Foldchange_column.min().values[0]
Red_expression_color= Log2Foldchange_column.max().values[0]
White_expression_color= Blue_expression_color + (Red_expression_color -Blue_expression_color)/2

step 8: You apply this color scheme to the network.

p4c.set_node_color_mapping('log2FoldChange', [Blue_expression_color,White_expression_color,Red_expression_color], ['#0000FF', '#FFFFFF', '#FF0000'],mapping_type='c', style_name='Sample1')

''

4.2 PCB concentration 2

step 9: You first import the expression table into a new dataframe.

PCB2_DEG= pd.read_csv('PCB concentration 2-GSE109565.top.table.tsv',sep='\t')
PCB2_DEG

step 10: You will now integrate the expression table into the nodetable of the AOP network. This will be done using the function: p4c.load_table_data_from_file where you select ‘True’ for the second variable as the first row of the PCB2_DEG table has the needed column names. You also select ’node’ for table and ‘CTL.GeneID’ for table_key_column as you want the expression data to be matched to the Gene ID’s that were added by CyTargetLinker extension in the previous notebook.

p4c.load_table_data_from_file('PCB concentration 2-GSE109565.top.table.tsv', first_row_as_column_names=True,table='node', table_key_column='CTL.GeneID')

{'mappedTables': [171586, 171624]}

step 11: You define the style for the mapping so that the expression values (log2FoldChange) are mapped to the gene nodes. This is done by first retrieving the log2FoldChange column.

Log2Foldchange_column = p4c.get_table_columns(table='node', columns='log2FoldChange')

Log2Foldchange_column

step 12: This is followed by definition of the color scheme so that low and high expression values receive distinct color.For the color scheme, you will use the following colors:

• Low expression value (minimum) = blue node color
• No expression value = white node color
• High expression value (maximum) = red node color

This color scheme was also described in the official py4cytoscape documentation (1).

Blue_expression_color= Log2Foldchange_column.min().values[0]
Red_expression_color= Log2Foldchange_column.max().values[0]
White_expression_color= Blue_expression_color + (Red_expression_color -Blue_expression_color)/2

step 13: You apply this color scheme to the network.

p4c.set_node_color_mapping('log2FoldChange', [Blue_expression_color,White_expression_color,Red_expression_color], ['#0000FF', '#FFFFFF', '#FF0000'],mapping_type='c', style_name='Sample1')

''

4.3 PCB concentration 3

step 14: You first import the expression table into a new dataframe.

PCB3_DEG= pd.read_csv('PCB concentration 3-GSE109565.top.table.tsv',sep='\t')
PCB3_DEG

step 15: You will now integrate the expression table into the nodetable of the AOP network. This will be done using the function: p4c.load_table_data_from_file where you select ‘True’ for the second variable as the first row of the PCB1_DEG table has the needed column names. You also select ’node’ for table and ‘CTL.GeneID’ for table_key_column as you want the expression data to be matched to the Gene ID’s that were added by CyTargetLinker extension in the previous notebook.

p4c.load_table_data_from_file('PCB concentration 3-GSE109565.top.table.tsv', first_row_as_column_names=True,table='node', table_key_column='CTL.GeneID')

{'mappedTables': [171586, 171624]}

step 16: You define the style for the mapping so that the expression values (log2FoldChange) are mapped to the gene nodes. This is done by first retrieving the log2FoldChange column.

Log2Foldchange_column = p4c.get_table_columns(table='node', columns='log2FoldChange')

Log2Foldchange_column

step 17: This is followed by definition of the color scheme so that low and high expression values receive distinct color.For the color scheme, you will use the following colors:

• Low expression value (minimum) = blue node color
• No expression value = white node color
• High expression value (maximum) = red node color

This color scheme was also described in the official py4cytoscape documentation (1).

Blue_expression_color= Log2Foldchange_column.min().values[0]
Red_expression_color= Log2Foldchange_column.max().values[0]
White_expression_color= Blue_expression_color + (Red_expression_color -Blue_expression_color)/2

step 18: You apply this color scheme to the network.

p4c.set_node_color_mapping('log2FoldChange', [Blue_expression_color,White_expression_color,Red_expression_color], ['#0000FF', '#FFFFFF', '#FF0000'],mapping_type='c', style_name='default')

''

p4c.save_session('PCB3network.cys')

This file has been overwritten.





{}

4.4 Roundup

step 19: You first import the expression table into a new dataframe.

Roundup_DEG= pd.read_csv('GSE109565.top.table-Roundup.tsv',sep='\t')
Roundup_DEG= Roundup_DEG[Roundup_DEG['padj'] < 0.05]
Roundup_DEG.to_excel('GSE109565.top.table-Roundup.xlsx')

step 20: You will now integrate the expression table into the nodetable of the AOP network. This will be done using the function: p4c.load_table_data_from_file where you select ‘True’ for the second variable as the first row of the PCB1_DEG table has the needed column names. You also select ’node’ for table and ‘CTL.GeneID’ for table_key_column as you want the expression data to be matched to the Gene ID’s that were added by CyTargetLinker extension in the previous notebook.

p4c.load_table_data_from_file('GSE109565.top.table-Roundup.xlsx', first_row_as_column_names=True,table='node', table_key_column='CTL.GeneID')

{'mappedTables': [171586, 171624]}

step 21: You define the style for the mapping so that the expression values (log2FoldChange) are mapped to the gene nodes. This is done by first retrieving the log2FoldChange column.

Log2Foldchange_column = p4c.get_table_columns(table='node', columns='log2FoldChange')

Log2Foldchange_column

step 22: This is followed by definition of the color scheme so that low and high expression values receive distinct color.For the color scheme, you will use the following colors:

• Low expression value (minimum) = blue node color
• No expression value = white node color
• High expression value (maximum) = red node color

This color scheme was also described in the official py4cytoscape documentation (1).

Blue_expression_color= Log2Foldchange_column.min().values[0]
Red_expression_color= Log2Foldchange_column.max().values[0]
White_expression_color= Blue_expression_color + (Red_expression_color -Blue_expression_color)/2

step 23: You apply this color scheme to the network.

p4c.set_node_color_mapping('log2FoldChange', [Blue_expression_color,White_expression_color,Red_expression_color], ['#0000FF', '#FFFFFF', '#FF0000'],mapping_type='c', style_name='Sample1')

''

Section 5: Mapping of dataset: E-MEXP-2599

In this section, you will map the transcriptomics expression data of datasets:E-MEXP-2599. This dataset aimed to identify transcriptional alterations in renal proximal tubular cell cultures due to nephrotoxin exposure.

5.1 CdCl2 exposure time 1

step 24: You first import the expression table into a new dataframe.

Cdcl2_X12= pd.read_csv('topTable_X12_CdCl2_.5 - X12_vehicle...control_.0.TSV.tsv',sep='\t')
Cdcl2_X12

step 25: You will now integrate the expression table into the nodetable of the AOP network. This will be done using the function: p4c.load_table_data_from_file where you select ‘True’ for the second variable as the first row of the PCB1_DEG table has the needed column names. You also select ’node’ for table and ‘CTL.GeneID’ for table_key_column as you want the expression data to be matched to the Gene ID’s that were added by CyTargetLinker extension in the previous notebook.

p4c.load_table_data_from_file('topTable_X12_CdCl2_.5 - X12_vehicle...control_.0.TSV.tsv', first_row_as_column_names=True,table='node', table_key_column='CTL.GeneID')

{'mappedTables': [171586, 171624]}

step 26: You define the style for the mapping so that the expression values (log2FoldChange) are mapped to the gene nodes. This is done by first retrieving the log2FoldChange column.

Log2Foldchange_column = p4c.get_table_columns(table='node', columns='log2FC')

Log2Foldchange_column

step 27: This is followed by definition of the color scheme so that low and high expression values receive distinct color.For the color scheme, you will use the following colors:

• Low expression value (minimum) = blue node color
• No expression value = white node color
• High expression value (maximum) = red node color

This color scheme was also described in the official py4cytoscape documentation (1).

Blue_expression_color= Log2Foldchange_column.min().values[0]
Red_expression_color= Log2Foldchange_column.max().values[0]
White_expression_color= Blue_expression_color + (Red_expression_color -Blue_expression_color)/2

step 28: You apply this color scheme to the network.

p4c.set_node_color_mapping('log2FC', [Blue_expression_color,White_expression_color,Red_expression_color], ['#0000FF', '#FFFFFF', '#FF0000'],mapping_type='c', style_name='Sample1')

''

5.2 CdCl2 exposure time 2

step 29: You first import the expression table into a new dataframe.

Cdcl2_X48= pd.read_csv('topTable_X48_CdCl2_.5 - X48_vehicle...control_.0.TSV.tsv',sep='\t')
Cdcl2_X48

step 30: You will now integrate the expression table into the nodetable of the AOP network. This will be done using the function: p4c.load_table_data_from_file where you select ‘True’ for the second variable as the first row of the PCB1_DEG table has the needed column names. You also select ’node’ for table and ‘CTL.GeneID’ for table_key_column as you want the expression data to be matched to the Gene ID’s that were added by CyTargetLinker extension in the previous notebook.

p4c.load_table_data_from_file('topTable_X48_CdCl2_.5 - X48_vehicle...control_.0.TSV.tsv', first_row_as_column_names=True,table='node', table_key_column='CTL.GeneID')

{'mappedTables': [171586, 171624]}

step 31: You define the style for the mapping so that the expression values (log2FoldChange) are mapped to the gene nodes. This is done by first retrieving the log2FoldChange column.

Log2Foldchange_column = p4c.get_table_columns(table='node', columns='log2FC')

Log2Foldchange_column

step 32: This is followed by definition of the color scheme so that low and high expression values receive distinct color. For the color scheme, you will use the following colors:

• Low expression value (minimum) = blue node color
• No expression value = white node color
• High expression value (maximum) = red node color

This color scheme was also described in the official py4cytoscape documentation (1).

Blue_expression_color= Log2Foldchange_column.min().values[0]
Red_expression_color= Log2Foldchange_column.max().values[0]
White_expression_color= Blue_expression_color + (Red_expression_color -Blue_expression_color)/2

step 33: You apply this color scheme to the network.

p4c.set_node_color_mapping('log2FC', [Blue_expression_color,White_expression_color,Red_expression_color], ['#0000FF', '#FFFFFF', '#FF0000'],mapping_type='c', style_name='Sample1')

''

5.3 CsA exposure time 1

step 34: You first import the expression table into a new dataframe.

CsA_X12= pd.read_csv('topTable_X12_CsA_.5 - X12_vehicle...control_.0.TSV.tsv',sep='\t')
CsA_X12

step 35: You will now integrate the expression table into the nodetable of the AOP network. This will be done using the function: p4c.load_table_data_from_file where you select ‘True’ for the second variable as the first row of the PCB1_DEG table has the needed column names. You also select ’node’ for table and ‘CTL.GeneID’ for table_key_column as you want the expression data to be matched to the Gene ID’s that were added by CyTargetLinker extension in the previous notebook.

p4c.load_table_data_from_file('topTable_X12_CsA_.5 - X12_vehicle...control_.0.TSV.tsv', first_row_as_column_names=True,table='node', table_key_column='CTL.GeneID')

{'mappedTables': [171586, 171624]}

step 36: You define the style for the mapping so that the expression values (log2FoldChange) are mapped to the gene nodes. This is done by first retrieving the log2FoldChange column.

Log2Foldchange_column = p4c.get_table_columns(table='node', columns='log2FC')

Log2Foldchange_column

step 37: This is followed by definition of the color scheme so that low and high expression values receive distinct color.For the color scheme, you will use the following colors:

• Low expression value (minimum) = blue node color
• No expression value = white node color
• High expression value (maximum) = red node color

This color scheme was also described in the official py4cytoscape documentation (1).

Blue_expression_color= Log2Foldchange_column.min().values[0]
Red_expression_color= Log2Foldchange_column.max().values[0]
White_expression_color= Blue_expression_color + (Red_expression_color -Blue_expression_color)/2

step 38: You apply this color scheme to the network.

p4c.set_node_color_mapping('log2FC', [Blue_expression_color,White_expression_color,Red_expression_color], ['#0000FF', '#FFFFFF', '#FF0000'],mapping_type='c', style_name='Sample1')

''

5.4 CsA exposure time 2

step 39: You first import the expression table into a new dataframe.

CsA_X48= pd.read_csv('topTable_X48_CsA_.5 - X48_vehicle...control_.0.TSV.tsv',sep='\t')
CsA_X48

step 40: You will now integrate the expression table into the nodetable of the AOP network. This will be done using the function: p4c.load_table_data_from_file where you select ‘True’ for the second variable as the first row of the PCB1_DEG table has the needed column names. You also select ’node’ for table and ‘CTL.GeneID’ for table_key_column as you want the expression data to be matched to the Gene ID’s that were added by CyTargetLinker extension in the previous notebook.

p4c.load_table_data_from_file('topTable_X48_CsA_.5 - X48_vehicle...control_.0.TSV.tsv', first_row_as_column_names=True,table='node', table_key_column='CTL.GeneID')

{'mappedTables': [171586, 171624]}

step 41: You define the style for the mapping so that the expression values (log2FoldChange) are mapped to the gene nodes. This is done by first retrieving the log2FoldChange column.

Log2Foldchange_column = p4c.get_table_columns(table='node', columns='log2FC')

Log2Foldchange_column

step 42: This is followed by definition of the color scheme so that low and high expression values receive distinct color.For the color scheme, you will use the following colors:

• Low expression value (minimum) = blue node color
• No expression value = white node color
• High expression value (maximum) = red node color

This color scheme was also described in the official py4cytoscape documentation (1).

Blue_expression_color= Log2Foldchange_column.min().values[0]
Red_expression_color= Log2Foldchange_column.max().values[0]
White_expression_color= Blue_expression_color + (Red_expression_color -Blue_expression_color)/2

step 43: You apply this color scheme to the network.

p4c.set_node_color_mapping('log2FC', [Blue_expression_color,White_expression_color,Red_expression_color], ['#0000FF', '#FFFFFF', '#FF0000'],mapping_type='c', style_name='Sample1')

''

5.5 Diquat dibromide exposure time 1

step 44: You first import the expression table into a new dataframe.

Diquat_X12= pd.read_csv('topTable_X12_Diquat_30 - X12_vehicle...control_.0.TSV.tsv',sep='\t')
Diquat_X12

step 45: You will now integrate the expression table into the nodetable of the AOP network. This will be done using the function: p4c.load_table_data_from_file where you select ‘True’ for the second variable as the first row of the PCB1_DEG table has the needed column names. You also select ’node’ for table and ‘CTL.GeneID’ for table_key_column as you want the expression data to be matched to the Gene ID’s that were added by CyTargetLinker extension in the previous notebook.

p4c.load_table_data_from_file('topTable_X12_CsA_.5 - X12_vehicle...control_.0.TSV.tsv', first_row_as_column_names=True,table='node', table_key_column='CTL.GeneID')

{'mappedTables': [171586, 171624]}

step 46: You define the style for the mapping so that the expression values (log2FoldChange) are mapped to the gene nodes. This is done by first retrieving the log2FoldChange column.

Log2Foldchange_column = p4c.get_table_columns(table='node', columns='log2FC')

Log2Foldchange_column 

step 47: This is followed by definition of the color scheme so that low and high expression values receive distinct color.For the color scheme, you will use the following colors:

• Low expression value (minimum) = blue node color
• No expression value = white node color
• High expression value (maximum) = red node color

This color scheme was also described in the official py4cytoscape documentation (1).

Blue_expression_color= Log2Foldchange_column.min().values[0]
Red_expression_color= Log2Foldchange_column.max().values[0]
White_expression_color= Blue_expression_color + (Red_expression_color -Blue_expression_color)/2

step 48: You apply this color scheme to the network.

p4c.set_node_color_mapping('log2FC', [Blue_expression_color,White_expression_color,Red_expression_color], ['#0000FF', '#FFFFFF', '#FF0000'],mapping_type='c', style_name='Sample1')

''

5.6 Diquat dibromide exposure time 2

step 49: You first import the expression table into a new dataframe.

Diquat_X48= pd.read_csv('topTable_X48_Diquat_30 - X48_vehicle...control_.0.TSV.tsv',sep='\t')
Diquat_X48

step 50: You will now integrate the expression table into the nodetable of the AOP network. This will be done using the function: p4c.load_table_data_from_file where you select ‘True’ for the second variable as the first row of the PCB1_DEG table has the needed column names. You also select ’node’ for table and ‘CTL.GeneID’ for table_key_column as you want the expression data to be matched to the Gene ID’s that were added by CyTargetLinker extension in the previous notebook.

p4c.load_table_data_from_file('topTable_X48_Diquat_30 - X48_vehicle...control_.0.TSV.tsv', first_row_as_column_names=True,table='node', table_key_column='CTL.GeneID')

{'mappedTables': [171586, 171624]}

step 51: You define the style for the mapping so that the expression values (log2FoldChange) are mapped to the gene nodes. This is done by first retrieving the log2FoldChange column.

Log2Foldchange_column = p4c.get_table_columns(table='node', columns='log2FC')

Log2Foldchange_column

step 52: This is followed by definition of the color scheme so that low and high expression values receive distinct color.For the color scheme, you will use the following colors:

• Low expression value (minimum) = blue node color
• No expression value = white node color
• High expression value (maximum) = red node color

This color scheme was also described in the official py4cytoscape documentation (1).

Blue_expression_color= Log2Foldchange_column.min().values[0]
Red_expression_color= Log2Foldchange_column.max().values[0]
White_expression_color= Blue_expression_color + (Red_expression_color -Blue_expression_color)/2

step 53: You apply this color scheme to the network.

p4c.set_node_color_mapping('log2FC', [Blue_expression_color,White_expression_color,Red_expression_color], ['#0000FF', '#FFFFFF', '#FF0000'],mapping_type='c', style_name='Sample1')

''

Section 6: Mapping of dataset: E-MEXP-3583

In this section, you will map the transcriptomics expression data of datasets:E-MEXP-3583. In this study, the effect of AgNPs and Ag+ on the transcriptome of human lung epithelial cell line A549 was determined by comparing two exposure times (24 and 48 hours).

6.1 Ag+ exposure time 1

step 54: You first import the expression table into a new dataframe.

Ag_X24= pd.read_csv('topTable_Ag._.1.3_24 - H2O.control_.0.0_24.tsv',sep='\t')
Ag_X24

step 55: You will now integrate the expression table into the nodetable of the AOP network. This will be done using the function: p4c.load_table_data_from_file where you select ‘True’ for the second variable as the first row of the PCB1_DEG table has the needed column names. You also select ’node’ for table and ‘CTL.GeneID’ for table_key_column as you want the expression data to be matched to the Gene ID’s that were added by CyTargetLinker extension in the previous notebook.

p4c.load_table_data_from_file('topTable_Ag._.1.3_24 - H2O.control_.0.0_24.tsv', first_row_as_column_names=True,table='node', table_key_column='CTL.GeneID')

{'mappedTables': [171586, 171624]}

step 56: You define the style for the mapping so that the expression values (log2FoldChange) are mapped to the gene nodes. This is done by first retrieving the log2FoldChange column.

Log2Foldchange_column = p4c.get_table_columns(table='node', columns='log2FC')

Log2Foldchange_column

step 57: This is followed by definition of the color scheme so that low and high expression values receive distinct color.For the color scheme, you will use the following colors:

• Low expression value (minimum) = blue node color
• No expression value = white node color
• High expression value (maximum) = red node color

This color scheme was also described in the official py4cytoscape documentation (1).

Blue_expression_color= Log2Foldchange_column.min().values[0]
Red_expression_color= Log2Foldchange_column.max().values[0]
White_expression_color= Blue_expression_color + (Red_expression_color -Blue_expression_color)/2

step 58: You apply this color scheme to the network.

p4c.set_node_color_mapping('log2FC', [Blue_expression_color,White_expression_color,Red_expression_color], ['#0000FF', '#FFFFFF', '#FF0000'],mapping_type='c', style_name='Sample1')

''

6.2 Ag+ exposure time 2

step 59: You first import the expression table into a new dataframe.

Ag_X48= pd.read_csv('topTable_Ag._.1.3_48 - H2O.control_.0.0_48.tsv',sep='\t')
Ag_X48

step 60: You will now integrate the expression table into the nodetable of the AOP network. This will be done using the function: p4c.load_table_data_from_file where you select ‘True’ for the second variable as the first row of the PCB1_DEG table has the needed column names. You also select ’node’ for table and ‘CTL.GeneID’ for table_key_column as you want the expression data to be matched to the Gene ID’s that were added by CyTargetLinker extension in the previous notebook.

p4c.load_table_data_from_file('topTable_Ag._.1.3_48 - H2O.control_.0.0_48.tsv', first_row_as_column_names=True,table='node', table_key_column='CTL.GeneID')

{'mappedTables': [171586, 171624]}

step 61: You define the style for the mapping so that the expression values (log2FoldChange) are mapped to the gene nodes. This is done by first retrieving the log2FoldChange column.

Log2Foldchange_column = p4c.get_table_columns(table='node', columns='log2FC')

Log2Foldchange_column 

step 62: This is followed by definition of the color scheme so that low and high expression values receive distinct color.For the color scheme, you will use the following colors:

• Low expression value (minimum) = blue node color
• No expression value = white node color
• High expression value (maximum) = red node color

This color scheme was also described in the official py4cytoscape documentation (1).

Blue_expression_color= Log2Foldchange_column.min().values[0]
Red_expression_color= Log2Foldchange_column.max().values[0]
White_expression_color= Blue_expression_color + (Red_expression_color -Blue_expression_color)/2

step 63: You apply this color scheme to the network.

p4c.set_node_color_mapping('log2FC', [Blue_expression_color,White_expression_color,Red_expression_color], ['#0000FF', '#FFFFFF', '#FF0000'],mapping_type='c', style_name='Sample1')

''

6.3 AgNP exposure time 1

step 64: You first import the expression table into a new dataframe.

Agnp_X24= pd.read_csv('topTable_AgNP_12.1_24 - H2O.control_.0.0_24.tsv',sep='\t')
Agnp_X24

step 65: You will now integrate the expression table into the nodetable of the AOP network. This will be done using the function: p4c.load_table_data_from_file where you select ‘True’ for the second variable as the first row of the PCB1_DEG table has the needed column names. You also select ’node’ for table and ‘CTL.GeneID’ for table_key_column as you want the expression data to be matched to the Gene ID’s that were added by CyTargetLinker extension in the previous notebook.

p4c.load_table_data_from_file('topTable_AgNP_12.1_24 - H2O.control_.0.0_24.tsv', first_row_as_column_names=True,table='node', table_key_column='CTL.GeneID')

{'mappedTables': [171586, 171624]}

step 66: You define the style for the mapping so that the expression values (log2FoldChange) are mapped to the gene nodes. This is done by first retrieving the log2FoldChange column.

Log2Foldchange_column = p4c.get_table_columns(table='node', columns='log2FC')

Log2Foldchange_column 

step 67: This is followed by definition of the color scheme so that low and high expression values receive distinct color.For the color scheme, you will use the following colors:

• Low expression value (minimum) = blue node color
• No expression value = white node color
• High expression value (maximum) = red node color

This color scheme was also described in the official py4cytoscape documentation (1).

Blue_expression_color= Log2Foldchange_column.min().values[0]
Red_expression_color= Log2Foldchange_column.max().values[0]
White_expression_color= Blue_expression_color + (Red_expression_color -Blue_expression_color)/2

step 68: You apply this color scheme to the network.

p4c.set_node_color_mapping('log2FC', [Blue_expression_color,White_expression_color,Red_expression_color], ['#0000FF', '#FFFFFF', '#FF0000'],mapping_type='c', style_name='Sample1')

''

6.4 AgNP exposure time 2

step 69: You first import the expression table into a new dataframe.

Agnp_X48= pd.read_csv('topTable_AgNP_12.1_48 - H2O.control_.0.0_48.tsv',sep='\t')
Agnp_X48

step 70: You will now integrate the expression table into the nodetable of the AOP network. This will be done using the function: p4c.load_table_data_from_file where you select ‘True’ for the second variable as the first row of the PCB1_DEG table has the needed column names. You also select ’node’ for table and ‘CTL.GeneID’ for table_key_column as you want the expression data to be matched to the Gene ID’s that were added by CyTargetLinker extension in the previous notebook.

p4c.load_table_data_from_file('topTable_AgNP_12.1_48 - H2O.control_.0.0_48.tsv', first_row_as_column_names=True,table='node', table_key_column='CTL.GeneID')

{'mappedTables': [171586, 171624]}

step 71: You define the style for the mapping so that the expression values (log2FoldChange) are mapped to the gene nodes. This is done by first retrieving the log2FoldChange column.

Log2Foldchange_column = p4c.get_table_columns(table='node', columns='log2FC')

Log2Foldchange_column

step 72: This is followed by definition of the color scheme so that low and high expression values receive distinct color.For the color scheme, you will use the following colors:

• Low expression value (minimum) = blue node color
• No expression value = white node color
• High expression value (maximum) = red node color

This color scheme was also described in the official py4cytoscape documentation (1).

Blue_expression_color= Log2Foldchange_column.min().values[0]
Red_expression_color= Log2Foldchange_column.max().values[0]
White_expression_color= Blue_expression_color + (Red_expression_color -Blue_expression_color)/2

step 73: You apply this color scheme to the network.

p4c.set_node_color_mapping('log2FC', [Blue_expression_color,White_expression_color,Red_expression_color], ['#0000FF', '#FFFFFF', '#FF0000'],mapping_type='c', style_name='Sample1')

''

p4c.notebook_export_show_image()

Section 7: Metadata

step 74. At last, the metadata belonging to this jupyternotebook is displayed which contains the version numbers of packages and system-set-up for interested users. This requires the usage of packages:Watermark and print_versions.

%load_ext watermark
!pip install print-versions

Requirement already satisfied: print-versions in c:\users\shaki\anaconda3\lib\site-packages (0.1.0)

%watermark

Last updated: 2025-06-03T20:22:47.796839+02:00

Python implementation: CPython
Python version       : 3.12.3
IPython version      : 8.25.0

Compiler    : MSC v.1938 64 bit (AMD64)
OS          : Windows
Release     : 11
Machine     : AMD64
Processor   : Intel64 Family 6 Model 140 Stepping 1, GenuineIntel
CPU cores   : 8
Architecture: 64bit

from print_versions import print_versions
print_versions(globals())

json==2.0.9
ipykernel==6.28.0
pandas==2.2.2
py4cytoscape==1.9.0