Transcriptomic analysis of bisphenol AF on early growth and development of zebrafish (Danio rerio) larvae

Bisphenol AF (BPAF), an alternative to bisphenol A, is widely detected in aquatic environments. Owing to health concerns, the toxic effects of BPAF on organisms are drawing attention. The present study aims to evaluate the toxicity of BPAF, combining the results of omics techniques and experiment. Employing transcriptome sequencing (RNA-seq), we obtained 391, 648, 512, and 545 differentially expressed genes (DEGs) in 0.1, 1, 10, and 100 μg/L BPAF-exposed zebrafish larvae, respectively. Gene ontology (GO) analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment revealed the early development, stimulus-response, and MAPK signaling pathway were significantly affected by BPAF. In addition, five hub genes (fgf3, fgf4, map2k1, myca, and casp3b) were highlighted as the key genes in MAPK signaling pathway using the protein-protein interaction network. Therefore, the RNA-seq results showed that early development and stimulus-response were the main processes affected by BPAF, which was consistent with our morphological and pathological results. The hatching rate of zebrafish embryos in 1 and 10 μg/L BPAF groups was significantly inhibited, and the oxidative stress indexes, including the level of total antioxidant capacity (T-AOC), superoxide dismutase (SOD), and lipid peroxidation (LPO), were significantly increased by the 100 μg/L BPAF treatment. Moreover, the activity of alkaline phosphatase (AKP) was significantly decreased in all BPAF exposure groups. In conclusion, exposure to BPAF at environmental relevant concentrations affected the early development and immune system of zebrafish larvae by modulating MAPK signaling pathway, and our results provide solid evidence for the future studies on the toxicity of bisphenols.