In this section we will focus on pathways associated to calcium channels. This interaction module has been resolved by modularity analysis in the SARS CoV pathways analysis page in this web site.
First, we can have a look at the pathway/gene network graph. Although some calcium channels subunits in this module are associated to two generic pathways related to membrane potential, they also point to NCAM and LGI-ADAM interaction pathways, both related to neurite outgrowth and nerve development.
Gene-pathways association map for this module. Node size by predicted activity against coronavirus aggregated at the node level.
Next, observe how pathways in this module score in this aggregation scatter plot.
Pathways scores for coronavirus obtained from aggregation of the target/pathways association map values. Log10SumPositiveInteractions.- Logged addition of all active interactions of predicted SARS active molecules with all the proteins belonging each pathway. AggregatedSarsScore.- Average value of predicted activity against SARS viruses aggregated at the pathway level from the values assigned to each protein/gene. [[activeTargetsInPathway]] / [proteinsInPathway]]].- Number of proteins considered relevant for SARS coronavirus infection in a particular pathway vs the total number of protein of that pathway referred in the ChEMBL DB. CombinedScore.- Linear combination of the number of evidences supporting a score (log10SumPositiveInteractions) and the predicted score from aggregation for a pathway (aggregatedSarsScore).
And how pathways in this module score in the volcano enrichment plot.
Volcano plot displaying the relationship between enrichmentScore and negative value of the flagged false discovery rate (-log10FDR). Size and colors by the aggregation scores described for the aggregation chart. Dotted line marks the 0.05 pValue threshold. Pathways in this module flagged.
Aggregation and enrichment scores are ranked in the bar chart below.
Best scored pathways bar chart. Pathways ranked by its overall combination of scores.
We can also visualize GSEA Enrichment plots for some representative pathways.
For full reference, see “Gene set enrichment analysis: A knowledge-based approach for interpreting genome-wide expression profiles.” and “ReactomePA: an R/Bioconductor package for reactome pathway analysis and visualization.“.
GSEA enrichment plots for four representative pathways in this module.
By using Reactome pathways hierarchies, we can identify additional biological roles to be played with associations of pathways in this module with same and higher level pathways.
high-low level pathways associations graph set by pathways in this module. Reactome low level pathways. Node size by predicted activity against coronavirus aggregated at the pathway level.
This hierarchical association with different level pathways can also be visualized in a tree map, where parent pathways in headings embrace their child partners, sized and colored by their respective scores as described in the legend.
