Visualized gene network reveals the novel target transcripts sox2 and pax6 of neuronal development in trans-placental exposure to bisphenol A

Background: Bisphenol A (BPA) is a ubiquitous endocrine disrupting chemical in our daily life, and its health effect in response to prenatal exposure is still controversial. Early-life BPA exposure may impact brain development and contribute to childhood neurological disorders. The aim of the present study was to investigate molecular target genes of neuronal development in trans-placental exposure to BPA. Methodology: A meta-analysis of three public microarray datasets was performed to screen for differentially expressed genes (DEGs) in exposure to BPA. The candidate genes of neuronal development were identified from gene ontology analysis in a reconstructed neuronal sub-network, and their gene expressions were determined using real-time PCR in 20 umbilical cord blood samples dichotomized into high and low BPA level groups upon the median 16.8 nM. Principal Findings: Among 36 neuronal transcripts sorted from DAVID ontology clusters of 457 DEGs using the analysis of Bioconductor limma package, we found two neuronal genes, sex determining region Y-box 2 (Sox2) and paired box 6 (Pax6), had preferentially down-regulated expression (Bonferroni correction p-value <10^-4 and log2-transformed fold change ≤-1.2) in response to BPA exposure. Fetal cord blood samples had the obviously attenuated gene expression of Sox2 and Pax6 in high BPA group referred to low BPA group. Visualized gene network of Cytoscape analysis showed that Sox2 and Pax6 which were contributed to neural precursor cell proliferation and neuronal differentiation might be downregulated through sonic hedgehog (Shh), vascular endothelial growth factor A (VEGFA) and Notch signaling. Conclusions: These results indicated that trans-placental BPA exposure down-regulated gene expression of Sox2 and Pax6 potentially underlying the adverse effect on childhood neuronal development.