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

Part 5: Visualization of transcriptomics expression datasets in the enriched AOP network part 2

The AOP project ► Key objective 2

Author: Shakira Agata

This Jupyter notebook describes the steps needed for the mapping of transcriptomics datasets in the constructed enriched AOP network. For this notebook, open license transcriptomics datasets were derived from ArrayExpress and Gene Expression Omnibus (GEO). These datasets were preprocessed followed by execution of statistical analysis to identify differential expression genes (DEG). The tabulation of differential gene expression data was subsequently mapped/integrated into the network. This notebook is subdivided into the following six sections:

• Section 1: System preparation
• Section 2: Retrieval of the enriched AOP network
• Section 3: Adaptation of gene node color within the enriched AOP network
• Section 4: Mapping of dataset: E-GEOD-69851
  • Section 4.1 Bisphenol A
  • Section 4.2 Farnesol
  • Section 4.3 Tetrachlorodibenzo p-dioxin
  • Section 4.4 Troglitazone
  • Section 4.5 Valproic acid
• Section 5: Mapping of dataset: E-GEOD-69851
  • Section 5.1 ACR exposure time 1
  • Section 5.2 ACR exposure time 2
  • Section 5.3 MA exposure time 1
  • Section 5.4 MA exposure time 2
  • Section 5.5 CP exposure time 1
  • Section 5.6 CP exposure time 2
• Section 6: 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 glob

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 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 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 within the enriched 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 commands.

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.'}

set_node_color_default('#a7a5a5',style_name='default')
set_edge_color_default('#01e735', style_name='default')

''

p4c.clone_network()
p4c.rename_network('clone-GSE69844')

{'network': 235771, 'title': 'clone-GSE69844'}

Section 4: Mapping of dataset:GSE69844

In this section, you will map the transcriptomics expression data of dataset:GSE69844.This will be done in similar fashion as in previous section, but will be streamlined due to high number of datafiles. In preparation for this section, you must first download the datafiles of the chemicals into separate folders.

4.1 Bisphenol A

4.1.1 Concentration 1

step 4: You first import the expression table into a new dataframe followed by data manipulation. This is done so that the final table can be used as input for the py4cytoscape function.

BisphenolA_concentration1= pd.read_csv('GSE69844.BisphenolA-1uM.tsv',sep='\t')
Adjusted_BisphenolA_concentration1= BisphenolA_concentration1[BisphenolA_concentration1['padj'] < 0.05]
BisphenolA_concentration_1=Adjusted_BisphenolA_concentration1.drop('ID', axis=1)
BisphenolA_Concentration1 = BisphenolA_concentration_1[['Entrez.Gene','Gene.Symbol','padj', 't', 'B','logFC','GB_LIST','SPOT_ID']]
BisphenolA_Concentration1.to_excel('GSE69844.BisphenolA-concentration1.xlsx',index=False)
BisphenolA_Concentration1

**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('GSE69844.BisphenolA-concentration1.xlsx', first_row_as_column_names=True,table='node', table_key_column='CTL.GeneID',network='Agata,S.-Part4-Molecular inflammation-process related AOP network')

{'mappedTables': [388893, 388931]}

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='logFC')

**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 (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 8: You apply this color scheme to the network.

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

''

3.1.2 Concentration 2

step 9: You first import the expression table into a new dataframe followed by data manipulation. This is done so that the final table can be used as input for the py4cytoscape function.

BisphenolA_concentration2= pd.read_csv('GSE69844.BisphenolAConcentration2.tsv',sep='\t')
Adjusted_BisphenolA_concentration2= BisphenolA_concentration2[BisphenolA_concentration2['adj.P.Val'] < 0.05]
BisphenolA_concentration_2=Adjusted_BisphenolA_concentration2.drop('ID', axis=1)
BisphenolA_Concentration2 = BisphenolA_concentration_2[['Gene.Symbol', 'adj.P.Val', 't', 'B','logFC','GB_LIST','SPOT_ID']]
BisphenolA_Concentration2.to_excel('GSE69844.BisphenolA-concentration2.xlsx',index=False)
BisphenolA_Concentration2

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 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('GSE69844.BisphenolA-concentration2.xlsx', first_row_as_column_names=True,table='node', table_key_column='CTL.GeneID')

{'mappedTables': [388893, 388931]}

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='logFC')

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('logFC', [Blue_expression_color,White_expression_color,Red_expression_color], ['#0000FF', '#FFFFFF', '#FF0000'],mapping_type='c', style_name='default')

''

3.1.3 Concentration 3

step 14: You first import the expression table into a new dataframe followed by data manipulation. This is done so that the final table can be used as input for the py4cytoscape function.

BisphenolA_concentration3= pd.read_csv('GSE69844.BisphenolAConcentration3.tsv',sep='\t')
Adjusted_BisphenolA_concentration3= BisphenolA_concentration3[BisphenolA_concentration3['adj.P.Val'] < 0.05]
BisphenolA_concentration_3=Adjusted_BisphenolA_concentration3.drop('ID', axis=1)
BisphenolA_Concentration3 = BisphenolA_concentration_3[['Gene.Symbol', 'adj.P.Val', 't', 'B','logFC','GB_LIST','SPOT_ID']]
BisphenolA_Concentration3.to_excel('GSE69844.BisphenolA-concentration3.xlsx',index=False)
BisphenolA_Concentration3

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('GSE69844.BisphenolA-concentration3.xlsx', first_row_as_column_names=True,table='node', table_key_column='CTL.geneID')

{'mappedTables': [388893, 388931]}

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='logFC')

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('logFC', [Blue_expression_color,White_expression_color,Red_expression_color], ['#0000FF', '#FFFFFF', '#FF0000'],mapping_type='c', style_name='default')

''

4.2 Farnesol

3.2.1 Concentration 1

step 19: You first import the expression table into a new dataframe and filtered for the significant genes. It turns out that concentration 1 of Farnesol does not yield significant genes and thus can’t be mapped in the AOP network.

Farnesol_concentration1= pd.read_csv('GSE69844.FarnesolConcentration1.tsv',sep='\t')
Adjusted_Farnesol_concentration1= Farnesol_concentration1[Farnesol_concentration1['adj.P.Val'] < 0.05]
Farnesol_concentration_1=Adjusted_Farnesol_concentration1.drop('ID', axis=1)
Farnesol_Concentration1 = Farnesol_concentration_1[['Gene.Symbol', 'adj.P.Val', 't', 'B','logFC','GB_LIST','SPOT_ID']]
Farnesol_Concentration1.to_excel('GSE69844.Farnesol-concentration1.xlsx',index=False)
Farnesol_Concentration1

3.2.2 Concentration 2

step 20: You first import the expression table into a new dataframe and filtered for the significant genes. It turns out that concentration 2 of Farnesol does not yield significant genes and thus can’t be mapped in the AOP network.

Farnesol_concentration2= pd.read_csv('GSE69844.FarnesolConcentration2.tsv',sep='\t')
Adjusted_Farnesol_concentration2= Farnesol_concentration2[Farnesol_concentration2['adj.P.Val'] < 0.05]
Farnesol_concentration_2=Adjusted_Farnesol_concentration2.drop('ID', axis=1)
Farnesol_Concentration2 = Farnesol_concentration_2[['Gene.Symbol', 'adj.P.Val', 't', 'B','logFC','GB_LIST','SPOT_ID']]
Farnesol_Concentration2.to_excel('GSE69844.Farnesol-concentration2.xlsx',index=False)
Farnesol_Concentration2

3.2.3 Concentration 3

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

Farnesol_concentration3= pd.read_csv('GSE69844.FarnesolConcentration3.tsv',sep='\t')
Adjusted_Farnesol_concentration3= Farnesol_concentration3[Farnesol_concentration3['adj.P.Val'] < 0.05]
Farnesol_concentration_3=Adjusted_Farnesol_concentration3.drop('ID', axis=1)
Farnesol_Concentration3 = Farnesol_concentration_3[['Gene.Symbol', 'adj.P.Val', 't', 'B','logFC','GB_LIST','SPOT_ID']]
Farnesol_Concentration3.to_excel('GSE69844.Farnesol-concentration3.xlsx',index=False)
Farnesol_Concentration3

step 22: 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('GSE69844.Farnesol-concentration3.xlsx', first_row_as_column_names=True,table='node', table_key_column='CTL.geneID')

{'mappedTables': [388893, 388931]}

step 23: 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='logFC')

Log2Foldchange_column

step 24: 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 25: You apply this color scheme to the network.

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

''

4.4 Tetrachlorodibenzo p-dioxin

3.4.1 Concentration 1

step 26: You first import the expression table into a new dataframe followed by data manipulation. This is done so that the final table can be used as input for the py4cytoscape function.

Tpdioxin_concentration1= pd.read_csv('GSE69844.TpdioxinConcentration1.tsv',sep='\t')
Adjusted_Tpdioxin_concentration1= Tpdioxin_concentration1[Tpdioxin_concentration1['adj.P.Val'] < 0.05]
Tpdioxin_concentration_1=Adjusted_Tpdioxin_concentration1.drop('ID', axis=1)
Tpdioxin_Concentration1 = Tpdioxin_concentration_1[['Gene.Symbol', 'adj.P.Val', 't', 'B','logFC','GB_LIST','SPOT_ID']]
Tpdioxin_Concentration1.to_excel('GSE69844.Tpdioxin-concentration1.xlsx',index=False)
Tpdioxin_Concentration1

step 27: 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.GeneMame’ for table_key_column as you want the expression data to be matched to the Gene names that were added by CyTargetLinker extension in the previous notebook.

p4c.load_table_data_from_file('GSE69844.Tpdioxin-concentration1.xlsx', first_row_as_column_names=True,table='node', table_key_column='CTL.geneID')

{'mappedTables': [388893, 388931]}

step 28: 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='logFC')

step 29: 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 30: You apply this color scheme to the network.

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

''

3.4.2 Concentration 2

step 31: You first import the expression table into a new dataframe followed by data manipulation. This is done so that the final table can be used as input for the py4cytoscape function.

Tpdioxin_concentration2= pd.read_csv('GSE69844.TpdioxinConcentration2.tsv',sep='\t')
Adjusted_Tpdioxin_concentration2= Tpdioxin_concentration2[Tpdioxin_concentration2['adj.P.Val'] < 0.05]
Tpdioxin_concentration_2=Adjusted_Tpdioxin_concentration2.drop('ID', axis=1)
Tpdioxin_Concentration2 = Tpdioxin_concentration_2[['Gene.Symbol', 'adj.P.Val', 't', 'B','logFC','GB_LIST','SPOT_ID']]
Tpdioxin_Concentration2.to_excel('GSE69844.Tpdioxin-concentration2.xlsx',index=False)
Tpdioxin_Concentration2

step 32: 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.GeneMame’ for table_key_column as you want the expression data to be matched to the Gene names that were added by CyTargetLinker extension in the previous notebook.

p4c.load_table_data_from_file('GSE69844.Tpdioxin-concentration2.xlsx', first_row_as_column_names=True,table='node', table_key_column='CTL.geneID')

{'mappedTables': [388893, 388931]}

step 33: 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='logFC')

step 34: 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 35: You apply this color scheme to the network.

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

''

3.4.3 Concentration 3

step 36: You first import the expression table into a new dataframe followed by data manipulation. This is done so that the final table can be used as input for the py4cytoscape function.

Tpdioxin_concentration3= pd.read_csv('GSE69844.TpdioxinConcentration3.tsv',sep='\t')
Adjusted_Tpdioxin_concentration3= Tpdioxin_concentration3[Tpdioxin_concentration3['adj.P.Val'] < 0.05]
Tpdioxin_concentration_3=Adjusted_Tpdioxin_concentration3.drop('ID', axis=1)
Tpdioxin_Concentration3 = Tpdioxin_concentration_3[['Gene.Symbol', 'adj.P.Val', 't', 'B','logFC','GB_LIST','SPOT_ID']]
Tpdioxin_Concentration3.to_excel('GSE69844.Tpdioxin-concentration3.xlsx',index=False)
Tpdioxin_Concentration3

step 37: 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.GeneMame’ for table_key_column as you want the expression data to be matched to the Gene names that were added by CyTargetLinker extension in the previous notebook.

p4c.load_table_data_from_file('GSE69844.Tpdioxin-concentration3.xlsx', first_row_as_column_names=True,table='node', table_key_column='CTL.geneID')

{'mappedTables': [388893, 388931]}

step 38: 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='logFC')

step 39: 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 40: You apply this color scheme to the network.

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

''

4.5 Valproic acid

3.5.1 Concentration 1

step 41: You first import the expression table into a new dataframe and filtered for the significant genes. It turns out that concentration 1 of Valproic acid does not yield significant genes and thus can’t be mapped in the AOP network.

Valproicacid_concentration1= pd.read_csv('GSE69844.ValproicacidConcentration1.tsv',sep='\t')
Adjusted_Valproicacid_concentration1= Valproicacid_concentration1[Valproicacid_concentration1['adj.P.Val'] < 0.05]
Valproicacid_concentration_1=Adjusted_Valproicacid_concentration1.drop('ID', axis=1)
Valproicacid_Concentration1 = Valproicacid_concentration_1[['Gene.Symbol', 'adj.P.Val', 't', 'B','logFC','GB_LIST','SPOT_ID']]
Valproicacid_Concentration1.to_excel('GSE69844.Valproicacid-concentration1.xlsx',index=False)
Valproicacid_Concentration1

step 42: 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('GSE69844.Valproicacid-concentration1.xlsx', first_row_as_column_names=True,table='node', table_key_column='CTL.geneID')

{'mappedTables': [388893, 388931]}

step 43: 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='logFC')

step 44: 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 45: You apply this color scheme to the network.

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

''

3.5.2 Concentration 2

step 46: You first import the expression table into a new dataframe and filtered for the significant genes. It turns out that concentration 2 of Valproic acid does not yield significant genes and thus can’t be mapped in the AOP network.

Valproicacid_concentration2= pd.read_csv('GSE69844.ValproicacidConcentration2.tsv',sep='\t')
Adjusted_Valproicacid_concentration2= Valproicacid_concentration2[Valproicacid_concentration2['adj.P.Val'] < 0.05]
Valproicacid_concentration_2=Adjusted_Valproicacid_concentration2.drop('ID', axis=1)
Valproicacid_Concentration2 = Valproicacid_concentration_2[['Gene.Symbol', 'adj.P.Val', 't', 'B','logFC','GB_LIST','SPOT_ID']]
Valproicacid_Concentration2.to_excel('GSE69844.Valproicacid-concentration2.xlsx',index=False)
Valproicacid_Concentration2

3.5.3 Concentration 3

step 47: You first import the expression table into a new dataframe and filtered for the significant genes. It turns out that concentration 3 of Valproic acid does not yield significant genes and thus can’t be mapped in the AOP network.

Valproicacid_concentration3= pd.read_csv('GSE69844.ValproicacidConcentration3.tsv',sep='\t')
Adjusted_Valproicacid_concentration3= Valproicacid_concentration3[Valproicacid_concentration3['adj.P.Val'] < 0.05]
Valproicacid_concentration_3=Adjusted_Valproicacid_concentration3.drop('ID', axis=1)
Valproicacid_Concentration3 = Valproicacid_concentration_3[['Gene.Symbol', 'adj.P.Val', 't', 'B','logFC','GB_LIST','SPOT_ID']]
Valproicacid_Concentration3.to_excel('GSE69844.Valproicacid-concentration3.xlsx',index=False)
Valproicacid_Concentration3

4.7 Troglitazone

4.7.1 Concentration 1

step 48: You first import the expression table into a new dataframe followed by data manipulation. This is done so that the final table can be used as input for the py4cytoscape function.

Troglitazone_concentration1= pd.read_csv('GSE69844.TroglitazoneConcentration1.tsv',sep='\t')
Adjusted_Troglitazone_concentration1= Troglitazone_concentration1[Troglitazone_concentration1['adj.P.Val'] < 0.05]
Troglitazone_concentration_1=Adjusted_Troglitazone_concentration1.drop('ID', axis=1)
Troglitazone_Concentration1 = Troglitazone_concentration_1[['Gene.Symbol', 'adj.P.Val', 't', 'B','logFC','GB_LIST','SPOT_ID']]
Troglitazone_Concentration1.to_excel('GSE69844.Troglitazone-concentration1.xlsx',index=False)
Troglitazone_Concentration1

step 49: 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('GSE69844.Troglitazone-concentration1.xlsx', first_row_as_column_names=True,table='node', table_key_column='CTL.geneID')

{'mappedTables': [388893, 388931]}

step 50: 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='logFC')

step 51: 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 52: You apply this color scheme to the network.

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

''

3.7.2 Concentration 2

step 53: You first import the expression table into a new dataframe followed by data manipulation. This is done so that the final table can be used as input for the py4cytoscape function.

Troglitazone_concentration2= pd.read_csv('GSE69844.TroglitazoneConcentration2.tsv',sep='\t')
Adjusted_Troglitazone_concentration2= Troglitazone_concentration2[Troglitazone_concentration2['adj.P.Val'] < 0.05]
Troglitazone_concentration_2=Adjusted_Troglitazone_concentration2.drop('ID', axis=1)
Troglitazone_Concentration2 = Troglitazone_concentration_2[['Gene.Symbol', 'adj.P.Val', 't', 'B','logFC','GB_LIST','SPOT_ID']]
Troglitazone_Concentration2.to_excel('GSE69844.Troglitazone-concentration2.xlsx',index=False)
Troglitazone_Concentration2

step 54: 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 names that were added by CyTargetLinker extension in the previous notebook.

p4c.load_table_data_from_file('GSE69844.Troglitazone-concentration2.xlsx', first_row_as_column_names=True,table='node', table_key_column='CTL.geneID')

{'mappedTables': [388893, 388931]}

step 55: 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='logFC')

step 56: 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 57: You apply this color scheme to the network.

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

''

3.7.3 Concentration 3

step 58: You first import the expression table into a new dataframe followed by data manipulation. This is done so that the final table can be used as input for the py4cytoscape function.

Troglitazone_concentration3= pd.read_csv('GSE69844.TroglitazoneConcentration3.tsv',sep='\t')
Adjusted_Troglitazone_concentration3= Troglitazone_concentration3[Troglitazone_concentration3['adj.P.Val'] < 0.05]
Troglitazone_concentration_3=Adjusted_Troglitazone_concentration3.drop('ID', axis=1)
Troglitazone_Concentration3 = Troglitazone_concentration_3[['Gene.Symbol', 'adj.P.Val', 't', 'B','logFC','GB_LIST','SPOT_ID']]
Troglitazone_Concentration3.to_excel('GSE69844.Troglitazone-concentration3.xlsx',index=False)
Troglitazone_Concentration3

step 59: 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 names that were added by CyTargetLinker extension in the previous notebook.

p4c.load_table_data_from_file('GSE69844.Troglitazone-concentration3.xlsx', first_row_as_column_names=True,table='node', table_key_column='CTL.geneID')

{'mappedTables': [388893, 388931]}

step 60: 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='logFC')

step 61: 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 62: You apply this color scheme to the network.

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

''

Section 5: Mapping of dataset:GSE44729

In this section, you will map the transcriptomics expression data of datasets: GSE44729. This dataset aimed to transcriptonally profile BEAS-2B cells for the comparison between controls and skin sensitizers, controls and respiratory sensitizers and controls and non-sensitizing irritants.

5.1 ACR exposure time 1

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

ACR_10h= pd.read_csv('adaptedACR10h.tsv',sep='\t')
ACR_10h

step 64: Unfortunately, the PTAFR gene is not included in the AOP network and thus can’t be mapped.

p4c.load_table_data_from_file('adaptedACR10h.tsv', first_row_as_column_names=True,table='node', table_key_column='CTL.geneID')

{'mappedTables': [388893, 388931]}

5.2 ACR exposure time 2

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

ACR_24h= pd.read_csv('adaptedACR24h.tsv',sep='\t')
ACR_24h

step 66: 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.

ACR24h= pd.read_csv('adaptedACR24h.tsv',sep='\t')
ACR24h_version1= ACR24h[ACR24h['padj'] < 0.05]
ACR24h_version2=ACR24h_version1.drop('ID', axis=1)
ACR24h_version3= ACR24h_version2[['GENE_SYMBOL', 'padj', 't', 'B','logFC','SPOT_ID']]
ACR24h_version3.to_excel('ACR24h-adjusted.xlsx',index=False)
ACR24h_version3

step 67: You now set the working directory to the clone of the network

p4c.set_current_network('clone-GSE69844')

{}

step 68: 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 table has the needed column names. You also select ’node’ for table and ‘CTL.GeneMame’ for table_key_column as you want the expression data to be matched to the Gene names that were added by CyTargetLinker extension in the previous notebook.

p4c.load_table_data_from_file('ACR24h-adjusted.xlsx', first_row_as_column_names=True,table='node', table_key_column='CTL.geneID',network='clone-GSE69844')

{'mappedTables': [235742, 235780]}

step 69: 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='logFC',network='clone-GSE69844')

step 70: 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 71: You apply this color scheme to the network.

p4c.set_node_color_mapping('logFC', [Blue_expression_color,White_expression_color,Red_expression_color], ['#0000FF', '#FFFFFF', '#FF0000'],mapping_type='c', style_name='default',network='clone-GSE69844')

''

5.3 MA exposure time 1

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

MA_10h= pd.read_csv('adaptedMA10h.tsv',sep='\t')
MA_10h

**step 73:**Unfortunately, these four genes are not included in the AOP network and thus can’t be mapped.

5.4 MA exposure time 2

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

MA_24h= pd.read_csv('adaptedMA24h.tsv',sep='\t')
MA_24h

step 75: 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.

MA24h= pd.read_csv('adaptedMA24h.tsv',sep='\t')
MA24h_version1= MA24h[MA24h['padj'] < 0.05]
MA24h_version2=MA24h_version1.drop('ID', axis=1)
MA24h_version3= MA24h_version2[['GENE_SYMBOL', 'padj', 't', 'B','logFC','SPOT_ID']]
MA24h_version3.to_excel('MA24h-adjusted.xlsx',index=False)
MA24h_version3

step 76: You now set the working directory to the clone of the network

p4c.set_current_network('clone-GSE69844')

{}

step 77: 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 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 names that were added by CyTargetLinker extension in the previous notebook.

p4c.load_table_data_from_file('MA24h-adjusted.xlsx', first_row_as_column_names=True,table='node', table_key_column='CTL.GeneID',network='clone-GSE69844')

{'mappedTables': [235742, 235780]}

step 78: 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='logFC',network='clone-GSE69844')

step 79: 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 80: You apply this color scheme to the network.

p4c.set_node_color_mapping('logFC', [Blue_expression_color,White_expression_color,Red_expression_color], ['#0000FF', '#FFFFFF', '#FF0000'],mapping_type='c', style_name='default',network='clone-GSE69844')

''

5.5 CP exposure time 1

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

CP_10h= pd.read_csv('adaptedCP10h.tsv',sep='\t')
CP_10h

step 82: 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.

CP10h= pd.read_csv('adaptedCP10h.tsv',sep='\t')
CP10h_version1= CP10h[CP10h['padj'] < 0.05]
CP10h_version2=CP10h_version1.drop('ID', axis=1)
CP10h_version3= CP10h_version2[['GENE_SYMBOL', 'padj', 't', 'B','logFC','SPOT_ID']]
CP10h_version3.to_excel('CP10h-adjusted.xlsx',index=False)
CP10h_version3

step 83: You now set the working directory to the clone of the network

p4c.set_current_network('clone-GSE69844')

{}

step 84: 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 table has the needed column names. You also select ’node’ for table and ‘CTL.GeneMame’ for table_key_column as you want the expression data to be matched to the Gene names that were added by CyTargetLinker extension in the previous notebook.

p4c.load_table_data_from_file('CP10h-adjusted.xlsx', first_row_as_column_names=True,table='node', table_key_column='CTL.geneID',network='clone-GSE69844')

{'mappedTables': [235742, 235780]}

step 85: 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='logFC',network='clone-GSE69844')

step 86: 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 87: You apply this color scheme to the network.

p4c.set_node_color_mapping('logFC', [Blue_expression_color,White_expression_color,Red_expression_color], ['#0000FF', '#FFFFFF', '#FF0000'],mapping_type='c', style_name='default',network='clone-GSE69844')

''

5.6 CP exposure time 2

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

CP_24h= pd.read_csv('adaptedCP24h.tsv',sep='\t')
CP_24h