Abstract:The embryonic period is a sensitive window for bisphenol A (BPA) exposure. However, embryonic development is a highly dynamic process with changing cell populations. The heterogeneity effects of BPA on fish embryo cells during development remain unclear. We applied single-cell RNA sequencing to analyze the impact of BPA exposure on transcriptome heterogeneity of 64 683 cells from zebrafish embryos at 8, 12, and 30 h postfertilization (hpf). Thirty-eight cell populations were identified and gene expression profiles of 16 cell populations were significantly altered by BPA. At 8 hpf, BPA mainly influenced the outer layer cell populations of embryos, such as neural plate border and enveloping layer cells. At 12 and 30 hpf, nervous system formation and heart morphogenesis were disturbed. The altered differential processes of the neural plate border, neural crest, and neuronal cells were found to lead to increased neurogenesis in zebrafish larvae. In the forebrain, midbrain, neurons, and optic cells, pathways related to cell division and DNA replication and repair were altered. Moreover, BPA also changed transforming growth factor (TGF) β signaling and heart tube morphogenesis in heart cells, leading to a decreased heartbeat in zebrafish larvae. Our study provides a comprehensive understanding of BPA toxicity on fish embryonic development at a single-cell level.The Supporting Information is available free of charge at <a class="ext-link" href="/doi/10.1021/acs.est.0c02428?goto=supporting-info">https://pubs.acs.org/doi/10.1021/acs.est.0c02428</a>.Additional experiment methods; experiment design of BPA; summary of sequencing results; expression of marker genes; cell lineages of zebrafish embryo; change of abundances of cell population; heat map of significant DEGs; cell-type-specific functional changes; time-course analysis; significantly changed time-course clusters of DEGs; average fluorescence intensity of apoptotic cells; images of apoptotic cells (Figures S1–S18); summary of cell types and their marker genes (Tables S1–S3) (<a class="ext-link" href="/doi/suppl/10.1021/acs.est.0c02428/suppl_file/es0c02428_si_001.pdf">PDF</a>)8 hpf significant pathways-KEGG; 8 hpf significant pathways-GO; and shared pathways (Supporting Dataset 1) (<a class="ext-link" href="/doi/suppl/10.1021/acs.est.0c02428/suppl_file/es0c02428_si_002.xlsx">XLSX</a>)8 hpf-kegg-network node; 12 hpf-kegg-network node; 30 hpf-kegg-network node; 8 hpf-go-network node; 12 hpf-go-network node; 30 hpf-go-network node; shared pathways in all groups (Supporting Dataset 2) (<a class="ext-link" href="/doi/suppl/10.1021/acs.est.0c02428/suppl_file/es0c02428_si_003.xlsx">XLSX</a>)12 hpf significant pathways-KEGG; 12 hpf significant pathways-GO; shared pathways (Supporting Dataset 3) (<a class="ext-link" href="/doi/suppl/10.1021/acs.est.0c02428/suppl_file/es0c02428_si_004.xlsx">XLSX</a>)30 hpf significant pathways-KEGG; 30 hpf significant pathways-GO; shared pathways (Supporting Dataset 4) (<a class="ext-link" href="/doi/suppl/10.1021/acs.est.0c02428/suppl_file/es0c02428_si_005.xlsx">XLSX</a>)This article has not yet been cited by other publications.

Single-Cell Sequencing Reveals Heterogeneity Effects of Bisphenol A on Zebrafish Embryonic Development