Fig. 3

Signed co-expression network analysis for P. orbicularis. a Correlation matrix for P. orbicularis. Values in the cells represent significant (P < 0.01) Pearson’s correlation of module eigenvalue to physiological parameters (top panel). Names (left panel) are arbitrary color-coded names for each module; values in parenthesis represent the number of genes per module. Empty cells indicate non-significant correlations (P ≥ 0.01). Individual cortisol [log (pg.mg^− 1)] and weight (g) are continuous values. Time (24 hpi and 96 hpi) and Status (coded 1 for control_24h, control_96h and resistant_96h groups and 2 for infected_24h) are discrete numeric values. b Heatmaps of the top 30 genes in module_blue-host and module_{turquoise-host.} Scales represent Log2 (prior 2) of the individual expression levels. Individuals were clustered using hierarchical clustering procedures implemented in the pheatmap R package [53]