Prediction of zebrafish embryonic developmental toxicity by integrating omics with adverse outcome pathway

New approach methodologies (NAMs), especially omics-based high-throughput bioassays have been developed rapidly, providing rich mechanistic information such as molecular initiation events (MIEs) and (sub)cellular key events (KEs) in adverse outcome pathways (AOPs). However, how to apply the knowledge of MIEs/KEs to predict adverse outcomes (AOs) induced by chemicals represents a new challenge for computational toxicology. Here, an integrated method named ScoreAOP was developed and evaluated to predict chemicals’ developmental toxicity for zebrafish embryos by integrating four related AOPs and dose-dependent reduced zebrafish transcriptome (RZT). The rules of ScoreAOP included 1) sensitivity of responsive KEs demonstrated by point of departure of KEs (PODKE), 2) evidence reliability and 3) distance between KEs and AOs. Moreover, eleven chemicals with different modes of action (MoAs) were tested to evaluate ScoreAOP. Results showed that eight of the eleven chemicals caused developmental toxicity at tested concentration in apical tests. All the tested chemicals’ developmental defects were predicted using ScoreAOP, whereas eight out of the eleven chemicals predicted by ScoreMIE which was developed to score MIEs disturbed by chemicals based on in vitro bioassays data. Finally, in terms of mechanism explanation, ScoreAOP clustered chemicals with different MoAs while ScoreMIE failed, and ScoreAOP revealed the activation of aryl hydrocarbon receptor (AhR) plays a significant role in dysfunction of cardiovascular system, resulting in zebrafish developmental defects and mortality. In conclusion, ScoreAOP represents a promising approach to apply mechanism information obtained from omics to predict AOs induced by chemicals.