Network architecture of transcriptomic stress responses in zebrafish embryos

Protein-protein interaction (PPI) network topology can contribute to explaining fundamental properties of genes, from expression levels to evolutionary constraints. Genes central to a network are more likely to be both conserved and highly expressed, whereas genes that are able to evolve in response to selective pressures but expressed at lower levels are located on the periphery of the network. The stress response is likewise thought to be conserved and its associated genes highly expressed, however, experimental evidence for these patterns is limited. Therefore, we examined here whether the transcriptomic response to two environmental stressors (heat, UV, and their combination) is related to PPI architecture in zebrafish (Danio rerio) embryos. We show that stress response genes are situated more centrally in the PPI network. The transcriptomic response to heat was located in both central and peripheral positions, whereas UV response transcripts occupied central to intermediate positions. Network position was further linked to the magnitude of fold changes of genes and number of their associated phenotype GO terms. Across treatments, differentially expressed genes in different parts of the network affected identical phenotypes. Our results indicate that the zebrafish stress response is considered conserved but also have stressor-specific aspects. These properties can aid in better understanding the organismal response to diverse and co-occurring stressors. Given the speed of contemporary changes in aquatic ecosystems, our approach can aid in identifying novel key regulators of the systemic response to specific stressors.