Global gene expression analysis reveals novel transcription factors associated with long-term low-level exposure of EA.hy926 human endothelial cells to bisphenol A

Abstract

Bisphenol A (BPA) is an endocrine disruptor that binds to estrogen receptors (ER); however, studies have shown that the ER pathway was not always the primary molecular mechanism of BPA's action in cells and that gene transcription could be altered by different exposure times and doses. Here, we sought to understand the correlation between the BPA-responsive genes that have associated biological functions and the transcription factors (TFs) involved in their regulation by repeatedly exposing human endothelial cells EA.hy926 to three nanomolar concentrations of BPA (10−9 M, 10−8 M, and 10−7 M) for 14 weeks, after which changes in global gene expression were determined by RNA sequencing. Cytoscape plug-in iRegulon was used to infer TFs involved in the control of BPA-deregulated genes. The results show a minimal overlap in deregulated genes between three concentrations of BPA, with 10−9 M BPA having the highest number of deregulated genes. TF analysis suggests that all three concentrations of BPA were active in the absence of an ER-mediated pathway. A unique set of TFs (NES≥4) has been identified for each BPA concentration, including the NFκB family and CEBPB for 10−9 M BPA, MEF family, AHR/ARNT, and ZBTB33 for 10−8 M BPA, and IRF1-7 and OVOL1/OVOL2 for 10−7 M BPA, whereas STAT1/STAT2 were common TFs for 10−9 M and 10−7 M BPA. Overall, our data suggest that long-term low-level exposure of EA.hy926 cells to BPA leads to concentration-specific changes in gene expression that are not controlled by the ER-mediated signaling but rather by other mechanisms.