Concentration-dependent transcriptome of zebrafish larvae for environmental bisphenol S assessment

The widespread use of bisphenol S (BPS) as an alternative to bisphenol A has captured attention due to its potential toxicity to aquatic organisms. In the present study, the zebrafish was used as a model to evaluate the toxicity of BPS and determine the underlying mechanisms. The environmental concentration-dependent (0, 0.1, 1, 10, 100, and 1000 μg/L BPS) transcriptome approach was employed in combination with toxicity assays to address the problem. Based on a weighted correlation network analysis, we speculated that excess reactive oxygen species (ROS) may initiate cellular events in BPS-exposed zebrafish, leading to multiple toxic effects. Furthermore, we used pathway enrichment analysis to identify key pathways (MAPK signalling pathway and metabolic pathways) that link the molecular mechanisms with different toxic effects. In addition, we performed protein-protein network and shortest path analyses to identify six hub genes (erbb2, rrm2, rps27a, his2h3c, cdk1, and mcm5) and their interactions. Moreover, we suggest that BPS may interact with erbb2 by molecular docking. Thus, the BPS-erbb2 interaction may activate the MAPK signalling and metabolic pathways, resulting in ROS production and then caused multiple toxic effects in zebrafish. This study provides information for characterising the mechanisms of BPS exposure in aquatic environments.