Investigation of autism-related transcription factors underlying sex differences in the effects of bisphenol A on transcriptome profiles and synaptogenesis in the offspring hippocampus

Thongkorn, Surangrat; Kanlayaprasit, Songphon; Kasitipradit, Kasidit; Lertpeerapan, Pattanachat; Panjabud, Pawinee; Hu, Valerie W.; Jindatip, Depicha; Sarachana, Tewarit

doi:10.1186/s13293-023-00496-w

Research
Open access
Published: 20 February 2023

Investigation of autism-related transcription factors underlying sex differences in the effects of bisphenol A on transcriptome profiles and synaptogenesis in the offspring hippocampus

Surangrat Thongkorn¹,
Songphon Kanlayaprasit²,
Kasidit Kasitipradit¹,
Pattanachat Lertpeerapan¹,
Pawinee Panjabud¹,
Valerie W. Hu³,
Depicha Jindatip^2,4 &
…
Tewarit Sarachana ORCID: orcid.org/0000-0003-4518-9399²

Biology of Sex Differences volume 14, Article number: 8 (2023) Cite this article

3124 Accesses
3 Citations
3 Altmetric
Metrics details

Abstract

Background

Bisphenol A (BPA) has been linked to susceptibility to autism spectrum disorder (ASD). Our recent studies have shown that prenatal BPA exposure disrupted ASD-related gene expression in the hippocampus, neurological functions, and behaviors associated with ASD in a sex-specific pattern. However, the molecular mechanisms underlying the effects of BPA are still unclear.

Methods

Transcriptome data mining and molecular docking analyses were performed to identify ASD-related transcription factors (TFs) and their target genes underlying the sex-specific effects of prenatal BPA exposure. Gene ontology analysis was conducted to predict biological functions associated with these genes. The expression levels of ASD-related TFs and targets in the hippocampus of rat pups prenatally exposed to BPA were measured using qRT-PCR analysis. The role of the androgen receptor (AR) in BPA-mediated regulation of ASD candidate genes was investigated using a human neuronal cell line stably transfected with AR-expression or control plasmid. Synaptogenesis, which is a function associated with genes transcriptionally regulated by ASD-related TFs, was assessed using primary hippocampal neurons isolated from male and female rat pups prenatally exposed to BPA.

Results

We found that there was a sex difference in ASD-related TFs underlying the effects of prenatal BPA exposure on the transcriptome profiles of the offspring hippocampus. In addition to the known BPA targets AR and ESR1, BPA could directly interact with novel targets (i.e., KDM5B, SMAD4, and TCF7L2). The targets of these TFs were also associated with ASD. Prenatal BPA exposure disrupted the expression of ASD-related TFs and targets in the offspring hippocampus in a sex-dependent manner. Moreover, AR was involved in the BPA-mediated dysregulation of AUTS2, KMT2C, and SMARCC2. Prenatal BPA exposure altered synaptogenesis by increasing synaptic protein levels in males but not in females, but the number of excitatory synapses was increased in female primary neurons only.

Conclusions

Our findings suggest that AR and other ASD-related TFs are involved in sex differences in the effects of prenatal BPA exposure on transcriptome profiles and synaptogenesis in the offspring hippocampus. These TFs may play an essential role in an increased ASD susceptibility associated with endocrine-disrupting chemicals, particularly BPA, and the male bias of ASD.

Highlights

Prenatal BPA exposure altered the transcriptome profiles of genes associated with ASD in the offspring hippocampus through several ASD-related transcription factors, including AR and ESR1, in a sex-dependent manner.
In addition to the known BPA targets (i.e., AR and ESR1), other ASD-related transcription factors, including KDM5B, SMAD4, and TCF7L2, are also involved in the BPA effects in the offspring hippocampus, possibly by direct interactions with BPA.
Prenatal BPA exposure also caused a significant reduction in the expression of ASD-related TFs in the hippocampus of female offspring but not in that of males.
ASD candidate genes (i.e., AUTS2, SMARCC2, and KMT2C), which are known AR targets, were differentially expressed in response to BPA in the human neuronal cell line with AR overexpression compared to the cell line with low AR protein levels, indicating that AR is involved in BPA-mediated effects on these ASD candidate genes.
BPA-responsive genes that are transcriptional targets of AR and other ASD-related transcription factors are associated with synaptogenesis. Prenatal BPA exposure caused an increase in synaptic protein levels in primary hippocampal neurons and in hippocampal tissues of male offspring but not in females.

Introduction

Bisphenol A (BPA) is an endocrine-disrupting chemical (EDC) that is widely used in polycarbonate plastics and epoxy resins. BPA can be found in several consumable products used in daily life, including the linings of food and beverage containers, dental sealants, thermal paper receipts, and recycled toilet paper [1, 2]. Moreover, BPA can be found in micro/nanoplastics, which have become a major environmental problem around the world [3]. In humans, the most common route of BPA exposure is oral ingestion, but BPA exposure can also occur through dermal contact or inhalation. Recent studies have shown that BPA can be detected in the serum and urine of pregnant women and children [4, 5]. BPA can cross the placenta [6] and the blood–brain barrier [7]. Exposure to BPA disrupts brain developmental processes, including neurogenesis [8], neuron migration [9], neurite outgrowth [10], and synaptogenesis [11, 12]. Our recent studies revealed that BPA is an environmental risk factor for neurological disorders, including autism spectrum disorder (ASD) [10, 13, 14] and Alzheimer's disease (AD) [15].

ASD is an early-onset neurodevelopmental disorder that is characterized by two main symptoms: (1) social interaction and communication deficits and (2) repetitive behaviors and restricted interests. ASD is a significant public health problem worldwide, with a prevalence of 1 in 44 children in the United States, as reported by the Centers for Disease Control and Prevention (CDC) [16]. ASD is a male-biased disorder that is at least four times more common in males than in females [17]. The exact cause of ASD is still unclear, but it is thought that both genetic and environmental factors influence ASD [10, 13, 14, 18,19,20,21,22]. Moreover, recent studies have reported that epigenetic regulatory mechanisms, including DNA methylation [19, 20, 22], histone modifications [23, 24], and noncoding RNA-associated gene silencing [25], are associated with ASD. Environmental factors that are thought to cause or increase the susceptibility of ASD include BPA and bisphenol-related compounds [10, 13, 14, 18, 26,27,28], polychlorinated biphenyls (PCBs), phthalates [29, 30], and cigarette smoke [31]. BPA can be detected in the serum and urine of individuals with ASD [26,27,28]. Kardas et al. investigated the serum BPA levels of children with ASD using high-performance liquid chromatography (HPLC) and reported that serum BPA was elevated in children with ASD compared to healthy control children [26]. In addition, Kondolot et al. measured BPA levels in ASD children and found a significant increase in plasma BPA levels in children with pervasive developmental disorder not otherwise specified (PDD-NOS), which is a subtype of ASD, compared to typically developing children [28]. Moreover, Hansen et al. measured BPA levels in maternal urine samples collected in gestational week 28 and assessed ASD symptoms of the children at the ages of 2 and 5 years [5]. Elevated odds ratios were observed among 5-year-old children within the 3rd tertile of BPA exposure with an ASD score above the 75th percentile, suggesting that prenatal BPA exposure may increase the risk of ASD symptoms [5]. These findings support that BPA exposure is a risk factor contributing to ASD.

Recent studies have reported that exposure to BPA disrupts the expression of genes in several brain regions, including the hippocampus [10, 13, 32], prefrontal cortex [14], hypothalamus [32], cerebellum [33], and amygdala [34]. The hippocampus is a brain region that plays an important role in learning and memory and is one of the brain regions known to be negatively impacted in ASD [35, 36]. Arambula et al. studied the effects of prenatal BPA exposure on the transcriptome profiles of the hypothalamus and hippocampus of neonatal rat offspring and found that prenatal BPA exposure induced sex-specific changes in Esr1 and Esr2 expression in the hypothalamus and oxytocin (Oxt) in the hippocampus [32]. Moreover, exposure to BPA during pregnancy caused apoptosis of neurons in the hippocampus of offspring [37]. Both the alterations in the expression of these genes and the apoptosis of hippocampal neurons induced by prenatal BPA exposure have been demonstrated to be involved in susceptibility to ASD [38,39,40].

Recently, we performed transcriptome profiling of hippocampal tissues isolated from neonatal rat pups prenatally exposed to BPA using RNA-seq analysis and found that prenatal BPA exposure caused dysregulation of genes in the hippocampus of offspring in a sex-dependent manner [10, 13]. Significantly differentially expressed genes (DEGs) in male and female rat pups showed a significant association with ASD, with a stronger association in the male group than in the female group [13]. In addition, we investigated the effects of prenatal BPA exposure on ASD-related genes (i.e., Arhgap32, Atp8a1, Cux1, Eif3h, Itga4, Kdm5c, and Mief2) known to regulate neuronal viability, neuritogenesis, and learning/memory and assessed these functions in the offspring. We found that the expression of these genes showed a male-specific relationship with abnormal neurological functions [10], suggesting that these genes may play an important role in BPA-induced adverse effects on hippocampus-related functions and on the learning/memory of male pups.

In addition to neuronal viability and neuritogenesis, exposure to BPA has been reported to affect synaptogenesis [11, 12, 41, 42]. Several studies have investigated the effects of BPA on adult synapse density [41, 42], but the effects of prenatal BPA exposure on synaptogenesis in the offspring brain remain unclear. Zhang et al. investigated the effects of prenatal BPA exposure on synaptogenesis in the offspring hippocampus and found that exposure to BPA decreased spine density in the hippocampus of the offspring [12]. Kawato et al. reported that perinatal exposure to BPA at 30 µg/kg·body weight of the mother reduced dendritic spine density in male offspring [11]. In female offspring, however, the total spine density was significantly increased at the estrus stage but decreased at the proestrus stage [11]. These findings indicate that there is a sex difference in the effects of prenatal BPA exposure on transcriptome profiles of genes associated with ASD and on neurological functions in the offspring hippocampus. However, the molecular mechanisms through which BPA exerts its sex-specific effects on the transcriptome profiles in the hippocampus of the offspring are still unclear.

In this study, we, therefore, sought to identify potential transcription factors that may be responsible for the sex difference in the effects of prenatal BPA exposure on the transcriptome profiles of ASD candidate genes in the hippocampus of offspring. First, we obtained transcriptome profiling data from RNA-seq analysis of hippocampal tissues isolated from neonatal rat pups prenatally exposed to BPA or vehicle control and identified differentially expressed genes (DEGs) in response to prenatal BPA exposure. Next, to identify transcription factors that may serve as upstream regulators of these DEGs, we obtained a list of transcription factors that have been associated with ASD and the targets of each ASD-related transcription factor. Then, hypergeometric distribution analysis was performed to identify transcription factors, whose transcriptional targets were significantly enriched in the list of BPA-responsive DEGs in the offspring hippocampus. Moreover, to determine whether BPA can directly interact with these ASD-related transcription factors, molecular docking analysis was performed using three-dimensional structures of BPA and these transcription factor molecules.

Since androgen receptor protein (AR) is one of the transcription factors known to interact with BPA and is thought to play an important role in the sex difference in the expression of ASD candidate genes [43,44,45,46,47,48,49], we further investigated whether AR is involved in the effects of BPA on the expression of ASD candidate genes. The human neuroblastoma cell line SH-SY5Y with low expression of AR protein was stably transfected with AR-overexpression plasmid or control plasmid and then treated with BPA or vehicle control. The expression of AR transcriptional targets which are also known ASD candidate genes was assessed using qRT-PCR analysis.

To further investigate biological functions, diseases, and pathways associated with the BPA-responsive DEGs that are the targets of ASD-related transcription factors, gene ontology analysis was performed using Ingenuity Pathway Analysis (IPA) software. As synaptogenesis was predicted to be one of the biological functions significantly associated with these DEGs, we further investigated the effects of prenatal BPA exposure on the synaptogenesis of hippocampal neurons isolated from male and female rat pups.

Methods

Collection of transcriptome profile data and the lists of DEGs in response to prenatal BPA exposure

RNA-seq data from our transcriptome profiling analysis of hippocampal tissues isolated from rat offspring prenatally exposed to BPA (n = 6 pups from independent litters; 3 males and 3 females) or vehicle control (n = 6 pups from independent litters; 3 males and 3 females) were obtained from the NCBI GEO Data sets database (accession number GSE140298) [13]. The details of animal husbandry, BPA treatment, and RNA-seq analysis were described in a previously published study [13]. Briefly, BPA was dissolved in molecular biology grade absolute ethanol and further diluted with corn oil to a final concentration of 5000 µg/kg·maternal body weight of BPA to treat each rat. The vehicle control treatment was prepared by mixing absolute ethanol with corn oil in amounts equivalent to those used for preparing BPA. After mating, each female rat was intragastrically administered either BPA or the vehicle control daily from GD1 until parturition. Total RNA was isolated from the hippocampus, and RNA-seq analysis was performed by BGI Genomics Co., Ltd using the Illumina HiSeq 4000 next-generation sequencing platform with 4 G reads (Illumina, Inc.) as previously described [13]. Poisson distribution analysis was performed to identify DEGs in response to BPA. Comparisons were performed with all male and female pups with the same treatment condition combined into one group and separately for each sex. DEGs with a p value < 0.05 and FDR < 0.05 were considered statistically significant.

Identification of TFs that are upstream regulators of DEGs

To identify transcription factors that have been associated with ASD and may serve as the upstream regulators of these DEGs, we obtained the list of human TFs from the Human Transcription Factors database (http://humantfs.ccbr.utoronto.ca/allTFs.php) [50] and overlapped them with the list of genes associated with ASD from the SFARI database (https://gene-archive.sfari.org/database/human-gene/) [51] using Venny 2.1 software (https://bioinfogp.cnb.csic.es/tools/venny/) [52]. Transcriptional targets of these ASD-related transcription factors were predicted using the TRANSFAC Curated Transcription Factor Targets database [53, 54], the TRANSFAC Predicted Transcription Factor Targets database [53, 54], the JASPAR Predicted Transcription Factor Targets database [55, 56], the CHEA Transcription Factor Targets database [57], the ENCODE Transcription Factor Targets database [58, 59], and the MotifMap Predicted Transcription Factor Targets database [60] through the Harmonizome database [61]. To identify TFs whose transcriptional targets were significantly enriched in the lists of DEGs, hypergeometric distribution analysis between the list of DEGs in response to prenatal BPA exposure and the list of transcriptional targets of each TF was then performed using the Keisan Online Calculator package (http://keisan.casio.com/exec/system/1180573201) [62]. A q value < 0.05 was considered a significant association.

Prediction of biological functions, disorders, and pathways associated with DEGs that are regulated by ASD-related TFs

Biological functions, disorders, and pathways associated with DEGs that were transcriptional targets of ASD-related TFs were predicted using Ingenuity Pathway Analysis (IPA) software (QIAGEN Inc., Hilden, Germany, https://www.qiagenbioinformatics.com/products/ingenuity-pathway-analysis/) [63]. The list of DEGs overlapped with the list of genes experimentally validated to be associated with each function/disorder/canonical pathway in the Ingenuity Knowledge Base database. Fisher’s exact test was then performed to calculate p values, and a p value < 0.05 was considered statistically significant.

Expression of ASD-related TFs in the developing hippocampus

To assess the expression of ASD-related TFs in the developing hippocampus, the expression profile data of each TF in the hippocampus of mice at embryonic days 11.5, 13.5, 15.5, 18.5, and postnatal day (PND) 4 were obtained from the Allen Brain Atlas database (https://mouse.brain-map.org/) [64]. For rat offspring (PND1), the expression of ASD-related TFs in the hippocampus was also investigated using qRT-PCR analyses.

In silico virtual screening of ASD-related TFs that are potential BPA targets

To predict whether BPA can directly interact with ASD-related TFs whose target genes are significantly differentially expressed in the hippocampus of rat offspring prenatally exposed to BPA, molecular docking analysis of BPA and those TFs was performed. For TFs whose ligands have been identified, molecular docking analysis of TFs and their known ligand(s) was also conducted for comparison with BPA. The three-dimensional structures of BPA and known ligands of each TF were obtained from the PubChem database (https://pubchem.ncbi.nlm.nih.gov/) [65] and processed to generate pdb files with clean geometry using BIOVIA Discovery Studio Visualizer 2020 [66, 67]. The processed BPA and ligand structures were then converted to pdbqt files for molecular docking analysis. For TFs, the three-dimensional structures of TFs were retrieved from the RCSB PDB Protein Data Bank database (https://www.rcsb.org/) [68]. TF structures obtained by the X-ray diffraction method with a resolution of less than 3 Å were selected for molecular docking analysis. Structures of water or other molecules surrounding TFs were removed, and all missing hydrogens and Kollman charges were added to the TF structure using AutoDockTools 4.2 software [69]. The processed TF structures were then converted to pdbqt files. Molecular docking analysis was performed as previously described [14]. Briefly, rigid docking of the TFs with BPA and their known ligands was performed using AutoDockTools 4.2 software with the Lamarckian GA algorithm and default parameters recommended by the software. For each run, a total of ten conformations of the BPA/ligand–TF complex were generated, and the conformation with the lowest binding free energy was selected. The mean and standard deviation values of the lowest binding-free energies from three independent runs were calculated and reported in this study. The interactions (i.e., hydrogen bonds and others) between each TF and BPA or known ligands were visualized using BIOVIA Discovery Studio Visualizer 2020 [66, 67].

Animal husbandry and treatment

Eight-week-old male and female Wistar rats were obtained from the National Laboratory Animal Center (NLAC) in Thailand. All experimental procedures were approved by the Chulalongkorn University Animal Care and Use Committee (Animal Use Protocol No. 1673007, 1773011, and 2073011). All procedures were performed in accordance with relevant guidelines and regulations. All rats were housed at the Chulalongkorn University Laboratory Animal Center (CULAC) under standard temperature (21 ± 1 °C) and humidity (30–70%) conditions on a 12-h light/dark cycle with food and RO-UV water available ad libitum. Rats were treated as previously described in the transcriptome profiling study [10, 13,14,15]. Briefly, we synchronized the estrous cycle and stimulated the mating behaviors of female rats by adding the bedding from a male rat cage to the cage of female rats for 5 days prior to mating. After that, one male and one female were placed together in the same cage overnight. The presence of a vaginal plug was observed on the following day and considered gestational day 1 (GD1).

After mating, female rats were divided into two groups: the BPA treatment group and the vehicle control group. Each female rat was weighed daily to monitor pregnancy and calculate the dose of BPA or vehicle control and was intragastrically administered either BPA or vehicle control daily from GD1 until parturition. For BPA treatment, BPA (Sigma-Aldrich, USA) was dissolved in molecular biology grade absolute ethanol (Merck Millipore, USA) to a final concentration of 250 mg/mL and then further diluted with corn oil to a final concentration of 5000 µg/kg·maternal body weight of BPA to treat each rat. The vehicle control treatment was prepared by mixing molecular biology grade absolute ethanol (Merck Millipore, USA) with corn oil in amounts equivalent to those used for preparing BPA. The dose of BPA used in this study is equal to the No-Observed-Adverse-Effect Level (NOAEL) of BPA determined by the U.S. Food and Drug Administration [70]. To prevent cross-contamination of the treatment conditions, rats in the BPA and control groups were raised separately in individually ventilated cages in a biohazard containment housing system. Separate sets of stainless-steel needles and all consumable products were used for oral gavage. All reusable materials were cleaned with ethanol and rinsed with copious amounts of Milli-Q deionized water before use.

The sex of every neonatal rat pup was determined on PND1 by visual inspection of the external genitalia. Male and female pups were distinguished by the distance between the external genitalia and the anus of the pup by two trained researchers [71]. Male pups exhibited a greater distance between the external genitalia and the anus than females (Additional file 1). For hippocampus isolation, neonatal rat pups were anesthetized by intraperitoneal injection of 100 mg/kg·BW sodium pentobarbital and euthanized by decapitation on ice. The brain was immediately placed in ice-cold, freshly prepared dissection medium (Additional file 2), and the hippocampus was quickly dissected on ice under a Nikon SMZ18 Stereo Microscope (Nikon, Japan). Meninges were removed completely. The hippocampus was then stored in RNAlater (Ambion, USA) at − 80 °C for subsequent qRT-PCR analysis or placed in a prechilled 35-mm dish containing ice-cold dissection medium for isolation of primary hippocampal neurons.

RNA isolation

Total RNA was extracted using the mirVana miRNA Isolation Kit (Ambion, USA) according to the manufacturer’s protocol. Briefly, hippocampal tissues from male and female rat pups (n = 3 pups/sex/treatment group from three different dams) were obtained and lysed in denaturing lysis buffer, which stabilized RNA and inactivated RNases. Total RNA was then extracted from the hippocampal tissue lysates using acid–phenol:chloroform and further purified over a glass–fiber filter. RNA was eluted from the purification column using the elution solution. The purity of total RNA was assessed using a NanoDrop spectrophotometer (Thermo Fisher Scientific, USA) and quantified using an Invitrogen Qubit 2.0 Fluorometer (Thermo Fisher Scientific, USA).

Quantitative RT–PCR analysis

Purified total RNA samples isolated from the hippocampus of rat offspring prenatally exposed to BPA or vehicle control were reverse transcribed to cDNA using a Maxime RT PreMix Kit (iNtRON Biotechnology, South Korea) according to the manufacturer’s protocol. Briefly, total RNA (1 µg) was added to a Maxime RT PreMix tube, and distilled water was then added to the tube to a total volume of 20 µL. The cDNA synthesis reaction was performed by incubating the tube at 45 °C for 60 min, followed by 95 °C for 5 min. After incubation, cDNA was diluted with nuclease-free water to obtain a final concentration of 10 ng/µL and used as the template for subsequent qPCR analysis. Quantitative PCR analysis was conducted using RealMOD Green W 2X qPCR Master Mix (iNtRON Biotechnology, South Korea) according to the manufacturer’s protocol. Briefly, 1 µL of cDNA was mixed with RealMOD Green W 2 × qPCR mix, forward primer, reverse primer, and nuclease-free water. Each cDNA sample was run in triplicate reactions. qPCR was performed in a Bio-Rad CFX Connect Real-Time System (Bio-Rad, USA) using the following conditions: an initial activation step at 95 °C for 10 min, followed by 40 cycles of denaturing at 95 °C for 20 s and annealing/extension at 55 °C for 30 s. A no-template control was included for every gene to check for contamination in the reaction. A melting curve analysis was set from 65 to 95 °C to assess the specificity of primers and the formation of the PCR product. The expression levels were calculated by the 2^−ΔΔCt method using rat 18S ribosomal RNA (Rn18s) or human 18S ribosomal RNA (RN18S) as an endogenous control gene. Forward and reverse primers were designed for rat TFs (i.e., Ar, Esr1, Egr2, Kdm5b, Sox5, Smad4, Tcf7l2, and Yy1), rat DEGs (i.e., Adnp, Chd2, Cic, Elavl2, Esrrb, Kat6a, Kmt2c, Nfib, and Snap25), and human genes (i.e., AUTS2, DICER1, ESR1, ESR2, KMT2C, and SMARCC2). The sequences of all primers are shown in Additional file 3.

Overexpression of AR in a human neuroblastoma cell line and BPA treatment

To investigate the role of AR in the effects of BPA on the expression of ASD candidate genes, the human neuroblastoma cell line SH-SY5Y (ATCC, USA) stably transfected with AR-expression plasmid or negative control plasmid was used. These cell lines were kindly provided by Associate Professor Ray-Chang Wu and Professor Valerie Hu, George Washington University, USA. For the AR-expression plasmid, the DNA sequence encoding the DNA binding and transcription activation domains of human AR was inserted into pSG5.HA expression vector containing a neomycin resistance gene. An empty pSG5.HA expression vector containing a neomycin resistance gene was used as a negative control vector. SH-SY5Y cells were stably transfected with the AR-expression plasmid or negative control plasmid using Lipofectamine LTX reagent with PLUS reagents (Invitrogen, USA) according to the manufacturer’s protocol. These cells were cultured in 1:1 MEM/EBSS with l-glutamine (GE Healthcare, USA) and Ham’s F12 media (GE Healthcare, USA) supplemented with 15% fetal bovine serum (Sigma-Aldrich, USA) and 1000 µg/mL geneticin (G418; Invitrogen, USA) in a 37 °C incubator with 5% CO₂. Geneticin was added to the media to select only the cells containing the neomycin-resistant plasmid. Western blot analysis was conducted to assess the expression of AR protein in these stably transfected cells (Additional file 4).

For BPA treatment, SH-SY5Y cells stably transfected with AR-expression or negative control plasmid were trypsinized and washed with PBS (GE Healthcare, USA). Then, cells were resuspended in 1:1 phenol red-free DMEM/F12 media (Invitrogen, USA) supplemented with 15% fetal bovine serum stripped with charcoal–dextran (Invitrogen, USA) and 1000 µg/mL geneticin (G418, Invitrogen, USA), and seeded at 1.5 × 10⁶ cells per well in a 6-well plate precoated with 10 µg/cm² collagen type I (Sigma-Aldrich, USA). Phenol red-free media and charcoal–dextran stripped serum were used in this experiment to minimize confounding effects from the phenol-related compound in the media and endogenous hormones in the serum. At 24 h after seeding, cells were induced to differentiate by replacing the medium with 1:1 phenol red-free DMEM/F12 medium (Invitrogen, USA) supplemented with 1% fetal bovine serum stripped with charcoal–dextran (Invitrogen, USA), 1,000 µg/ml geneticin (G418; Invitrogen, USA), and 10 µM retinoic acid (RA; Sigma-Aldrich, USA). BPA (Sigma-Aldrich, USA) was dissolved in molecular biology grade absolute ethanol (Merck Millipore, USA) to a final concentration of 250 mg/mL, and then added to the differentiation media to a final concentration of 1 ng/mL. For the vehicle control treatment, molecular biology grade absolute ethanol (Merck Millipore, USA) in amounts equivalent to those used for preparing BPA was added to the differentiation media. The cells were differentiated for 7 days and then lysed in GENEzol reagent (Geneaid Biotech, Taiwan) for RNA extraction. Total RNA was then isolated from the cell lysates using the phenol:chloroform method according to the manufacturer’s protocol. The purity and concentration of total RNA were assessed using a NanoDrop spectrophotometer (Thermo Fisher Scientific, USA).

Synaptogenesis assay

Primary hippocampal neurons were isolated from neonatal rat pups that were prenatally exposed to BPA (5000 µg/kg·maternal BW per day) or vehicle control. Primary hippocampal neurons were cultured according to Beaudoin et al. (2012) with slight modifications [72]. Briefly, the hippocampus of neonatal rat pups on PND1 was dissected and kept in ice-cold, freshly prepared dissection medium (Additional file 2). Then, the medium was removed, and the tissues were digested using 5 ml 2.5% w/v trypsin solution (Thermo Fisher Scientific, USA) and incubated at 37 °C for 20 min. Then, 25 µl 1% DNase I solution (Sigma-Aldrich, USA) was added to the tube and incubated at room temperature for 5 min. The hippocampal tissue chunks were gently washed twice in ice-cold, freshly prepared dissection medium. Hippocampal tissues were resuspended in plating medium (Additional file 2) and gently triturated using a fire-polished glass Pasteur pipette to dissociate the cells. The cells were counted and seeded in a 35-mm cell culture dish containing a 22-mm coverslip precoated with poly-l-lysine (Sigma-Aldrich, USA) at 2 × 10⁵ cells/coverslip. The cells were allowed to attach to the coverslip for 4 h in a CO₂ incubator. Then, the plating medium was gently removed and replaced with fresh maintenance medium (Additional file 2). After two days of incubation at 37 °C and 5% CO₂, half of maintenance medium was replaced with the maintenance medium supplemented with cytosine arabinoside (ara-C, 1-b-d-arabinofuranosylcytosine; Sigma-Aldrich, USA) at a final concentration of 2 µM to inhibit the proliferation of nonneuronal cells. Cells were further incubated for a total of 14 days after seeding. Half of the maintenance medium was replaced with fresh maintenance medium without ara-C once every 3 days.

To investigate the effects of prenatal BPA treatment on synaptogenesis, mature primary hippocampal neurons (DIV14) were gently washed with ice-cold PBS and then fixed with freshly prepared 4% paraformaldehyde in PBS for 15 min. The cells were washed with PBS for 5 min three times and blocked using 3% BSA (Capricorn Scientific, Germany) at room temperature for 30 min in a dark chamber. The cells were then immunostained for synaptic markers, including Syn1 (presynaptic marker), Psd95 (postsynaptic marker), and Map2 (mature neuron marker). Rabbit anti-Syn1 (ab8; Abcam, UK), mouse anti-Psd95 (ab2723; Abcam, UK), and chicken anti-Map2 (ab5392; Abcam, UK) antibodies were added to the cells and incubated at 4 °C overnight. On the following day, donkey anti-rabbit Alexa 647 (ab150063; Abcam, UK), donkey anti-mouse Alexa 488 (ab150109; Abcam, UK), and donkey anti-chicken Alexa 405 (ab175675; Abcam, UK) secondary antibodies were added to the cells and incubated for 1 h at room temperature. The coverslip with cells was mounted to a glass slide using ProLong Diamond Antifade reagent (Invitrogen, USA). The experiment was performed using three litters of pups/sex/treatment group. For each litter of pups, a total of 15 differentiated primary hippocampal neurons were captured using a LSM800 confocal laser scanning microscope (Carl Zeiss, Germany). The fluorescence intensities of Syn1, Psd95, and Map2 staining were measured using ZEN Blue software (Carl Zeiss, Germany). The number of Syn1 and Psd95 colocalization sites on mature neurons was counted as synaptic puncta using Imaris software (Bitplane, UK).

Western blot analysis

Western blot analysis was performed to investigate the expression of synaptic proteins in the hippocampus of rat offspring prenatally exposed to BPA and vehicle control. Hippocampal tissues from the pups in the BPA group (n = 6 pups; 3 males and 3 females from 3 different dams) and the vehicle control group (n = 6 pups; 3 males and 3 females from 3 different dams) were used. One male and one female pup were randomly selected from each dam. The hippocampal tissues were dissected from the rat pups (PND1), and protein was isolated using Genezol reagent (Geneaid Biotech, Taiwan) according to the manufacturer’s protocol. Protein pellets were resuspended in lysis buffer consisting of 7 M urea, 2 M thiourea, 4% w/v CHAPS, and 100 mM dithiothreitol (DTT). The concentrations of the protein samples were measured using Bradford protein assays using Bradford reagent (Bio-Rad, USA), with bovine serum albumin (BSA; Capricorn Scientific, Germany) as a standard. A total of 15 µg of each protein sample was separated via 10% SDS-PAGE and then transferred to 0.2 µm Immu-Blot PVDF membranes (Bio-Rad, USA) using a Mini-PROTEAN Tetra system (Bio-Rad, USA). The membrane was blocked with 5% nonfat dry milk (Bio-Rad, USA) in TBST for 1 h at room temperature. After that, the membrane was incubated overnight at 4 °C with rabbit anti-Syn1 antibody (1:10,000, 106011; Synaptic Systems, Germany), rabbit anti-Syp antibody (1:1,000, ab52636; Abcam, UK), mouse anti-Psd95 antibody (1:1,000, ab2723; Abcam, UK), and rabbit anti-Map2 antibody (1:1,000, #4542; Cell Signaling Technology, USA). The membrane was then washed and incubated with donkey anti-rabbit (ab7083, Abcam, UK) conjugated with HRP or donkey anti-mouse conjugated with HRP antibody (ab205724; Abcam, UK) for 1 h at room temperature. The membrane was visualized by Amersham ECL Western Blotting Detection Reagent (Cytiva, USA), and band images were taken using an Amersham ImageQuant 800 Western Blot Imaging System (Cytiva, USA). Band intensity was analyzed using ImageJ software [73].

Prediction of TFs involved in the effects of prenatal BPA exposure on the transcriptome profiles of the offspring hippocampus using data from other independent studies

To determine whether the ASD-related TFs predicted to be BPA targets and regulate DEGs in the offspring hippocampus were reproducible when other transcriptome profiling data were used, transcriptome profiling data of hippocampal tissues from rats or mice exposed to BPA in utero from previously published studies of other groups were obtained from the NCBI PubMed database [74] in a search performed on September 12, 2022, using the keywords “bisphenol A” and “hippocampus.” Studies that met the following criteria were included in this study: (i) the animal model was rodents and exposed to BPA in utero; (ii) transcriptome profiling analysis of the hippocampus of offspring was performed using RNA-seq or microarray techniques; and (iii) raw RNA-seq or microarray data or the complete list of all DEGs in response to prenatal/perinatal BPA exposure in the hippocampus were provided. When raw transcriptome profiling data were available, the data were obtained and reanalyzed using the Galaxy platform (https://usegalaxy.org/) [75]. DEGs were then identified by DESeq2 analysis, and p-value < 0.05 and FDR < 0.05 were considered statistically significant. If raw transcriptome data were not available, complete lists of all DEGs reported in the study were obtained. Hypergeometric distribution analysis was then performed to determine the association between the list of DEGs from these studies and the list of transcriptional targets of ASD-related TFs.

Statistical analyses

The statistical analyses were performed using the SPSS software package for Windows [76]. Two-tailed Student’s t test analysis was performed to determine the statistical significance of the differences in the mean values. A p value < 0.05 was considered statistically significant.

Results

Dysregulated genes in the hippocampus of offspring prenatally exposed to BPA were transcriptional targets of ASD-related TFs

Our recent studies have shown that prenatal BPA exposure altered the transcriptome profiles in the hippocampus of offspring, and changes in the expression of ASD candidate genes in response to BPA showed a male-specific correlation with hippocampal functions [10, 13]. Moreover, we found that BPA-responsive genes in the prefrontal cortex of rat offspring were potentially regulated by several ASD-related TFs, including AR, ESR1, and RORA [14]. To determine whether genes that were differentially expressed in the hippocampus of offspring prenatally exposed to BPA were transcriptional targets of ASD-related TFs, we obtained lists of DEGs from our previous RNA-seq analysis of hippocampal tissues isolated from rat offspring prenatally exposed to BPA (n = 6 pups; 3 males and 3 females) or vehicle control (n = 6 pups; 3 males and 3 females) for subsequent analysis. When both male and female pups in the same treatment group were combined into one group, a total of 4525 genes were reported to be significantly dysregulated in response to prenatal BPA exposure [13]. When each sex of the pups was analyzed separately, a total of 2078 and 3055 genes were found to be differentially expressed in the male and female groups, respectively [13]. These lists of DEGs were used for subsequent analysis in this study.

To identify transcription factors that have been associated with ASD and may serve as upstream regulators of these DEGs, we obtained a list of 1639 proteins identified as transcription factor proteins from the Human Transcription Factors database. Among these TFs, a total of 96 TFs have been determined to be ASD candidate genes/proteins by the SFARI database. The list of these ASD-related TFs is shown in Additional file 5. Since TFs are known to regulate the transcription of multiple genes, we obtained the list of transcriptional targets of each ASD-related TF from multiple databases through the Harmonizome database. The transcriptional targets of 34 ASD-related TFs have been identified and provided in the Harmonizome database.

To predict whether the DEGs in the hippocampus of offspring prenatally exposed to BPA were regulated by ASD-related TFs, the lists of DEGs and the list of transcriptional targets of each ASD-related TF were compared. Hypergeometric distribution analysis was then performed to assess the significance of the association between DEGs and transcriptional targets of each TF (Additional file 6). The hypergeometric distribution analysis revealed eight ASD-related TFs whose transcriptional targets were significantly enriched in the lists of DEGs (q value < 0.05; Table 1). These TFs were AR, ESR1, EGR2, KDM5B, SOX5, SMAD4, TCF7L2, and YY1. Interestingly, the targets of seven TFs (i.e., AR, ESR1, EGR2, KDM5B, SOX5, SMAD4, and YY1) were significantly associated with BPA-responsive genes in a sex-dependent pattern. The targets of AR, ESR1, EGR2, KDM5B, and SOX5 were manually curated by the TRANSFAC Curated Transcription Factor Targets database. This finding suggests that these ASD-related TFs may serve as potential upstream regulators of BPA-responsive genes, and BPA may disrupt the transcriptome profiles in the hippocampus of male and female offspring through different ASD-related TFs.

Table 1 Hypergeometric distribution analysis of transcriptional targets of ASD-related TFs and BPA-responsive genes in the hippocampus of male and female offspring

Full size table

BPA-responsive genes that are transcriptional targets of ASD-related TFs are also associated with ASD

Transcriptional targets of ASD-related TFs that were significantly dysregulated in the hippocampus of pups exposed to BPA in utero were analyzed by Ingenuity Pathway Analysis (IPA) software to predict biological functions, disorders, and pathways associated with these TF targets (Tables 2, 3, Additional files 7, 8, 9, 10, 11, 12, 13, 14). As AR and ESR1 are known BPA targets and are thought to be involved in the male bias of ASD [46, 48, 49], we, therefore, focused on these two TFs in this study. Interestingly, AR targets that were differentially expressed in male pups prenatally exposed to BPA were significantly associated with “autism spectrum disorder or intellectual disability” (p value = 1.04E−07) and “pervasive developmental disorder” (p value = 1.42E−04) (Table 2). Moreover, neurological functions known to be disrupted in ASD, including “proliferation of neuronal cells” (p value = 1.87E−05), “developmental process of synapse” (p value = 1.97E−04), and “excitatory postsynaptic potential” (p value = 8.25E−04) were also highlighted. ESR1 targets that were significantly altered in the hippocampus of male pups exposed to BPA were associated with “abnormal morphology of neurons” (p value = 3.22E−05), “neuritogenesis” (p value = 4.35E-05), and “proliferation of neuronal cells” (p value = 5.14E−05), all of which have been implicated in ASD (Table 3).

Table 2 Neurological diseases and functions associated with AR target genes that were differentially expressed in the hippocampus of male offspring prenatally exposed to BPA

Full size table

Table 3 Neurological diseases and functions associated with ESR1 target genes that were differentially expressed in the hippocampus of male offspring prenatally exposed to BPA

Full size table

Transcriptional target genes of EGR2, KDM5B, SOX5, SMAD4, TCF7L2, and YY1 that were differentially expressed in the BPA group were also analyzed by IPA (Additional files 7, 8, 9, 10, 11, 12, 13, 14). In addition to AR targets, KDM5B targets dysregulated in male and female pups, and SMAD4 targets dysregulated in male pups were also significantly associated with “autism spectrum disorder or intellectual disability” (Additional files 8, 9, 10). These results suggest that transcriptional targets of these ASD-related TFs that were differentially expressed in the offspring hippocampus are associated with ASD, and BPA may increase the susceptibility of ASD in male and female offspring by disrupting the expression of genes related to ASD through different TFs. AR, KDM5B, and SMAD4 may be involved in the effects of prenatal BPA exposure on the expression of ASD candidate genes in male offspring, whereas only KDM5B may be involved in females.

Molecular docking analysis of BPA and ASD-related TFs whose targets were differentially expressed in the hippocampus of offspring prenatally exposed to BPA

To examine whether BPA can directly interact with ASD-related TFs, the binding affinity between BPA and each ASD-related TF molecule was assessed using molecular docking analysis with Discovery Studio 2020 and AutoDock 4.2 software. Known ligands of each TF, when available, were also included for comparisons. Molecular docking analysis of BPA and the 8 ASD-related TFs whose targets were differentially expressed in response to prenatal BPA exposure showed that BPA directly interacted with AR, ESR1, EGR2, KDM5B, SOX5, SMAD4, TCF7L2, and YY1 (Table 4). Notably, SMAD4, KDM5B, and TCF7L2 have not been reported to interact with BPA, suggesting that these ASD-related TFs may serve as novel targets of BPA. When compared to known ligands, BPA interacted with AR at the same binding sites with 5α-dihydrotestosterone and testosterone and interacted with ESR1 at the same binding sites with 17β-estradiol (Fig. 1a, b). Unlike AR and ESR1, BPA interacted with KDM5B and TCF7L2 at different locations from their ligands, GSK-J1 and 6-[2-(2H-tetrazol-5-yl)ethyl]-1,2-benzoxazol-3-one, respectively (Fig. 2a, b). It is noteworthy that three TFs (i.e., AR, ESR1, and SMAD4) that interacted with BPA with low binding free energies also showed a male-specific enrichment of their transcriptional targets in the list of BPA-responsive genes (Table 1), suggesting that BPA may exert male-specific effects on the expression of ASD candidate genes in the offspring hippocampus and related functions through direct interactions with these TFs.

Table 4 Molecular docking analysis of BPA and ASD-related TFs

Full size table

Effects of prenatal BPA exposure on the expression of ASD-related TFs

In addition to direct interactions with TFs, BPA may disrupt the expression of ASD candidate genes in the offspring hippocampus through other pathways that affect the expression of TFs [32, 34, 77]. To examine the effects of prenatal BPA exposure on the expression of ASD-related TFs in the hippocampus of offspring, we first assessed whether these ASD-related TFs were expressed in the developing brain of rodents by obtaining the expression data of these TFs in embryonic and early postnatal stages of rodents from the Allen Brain Atlas database (Additional file 15). Seven of these TFs (i.e., Ar, Esr1, Egr2, Kdm5b, Sox5, Smad4, and Tcf7l2) were expressed in the rostral secondary prosencephalon and the telencephalic vesicle of the developing brain (Additional file 15). The expression data of Yy1 were not available. Moreover, qRT-PCR analysis of hippocampal tissues isolated from rat pups (PND1) prenatally exposed to BPA (n = 6 pups; 3 males and 3 females from 3 different dams) or vehicle control (n = 6 pups; 3 males and 3 females from 3 different dams) was conducted (Fig. 3). When male and female pups in the same treatment group were combined into one group, the expression levels of ASD-related TFs were not significantly changed. However, when each sex of pups was analyzed separately, the qRT-PCR analysis revealed a sex difference in the effects of prenatal BPA exposure on the expression of these TFs. In male pups, the expression of these TFs in the hippocampus was not impacted by prenatal BPA exposure, but in females, the expression levels of several TFs (i.e., Ar, Esr1, Egr2, Sox5, Smad4, and Tcf7l2) were significantly reduced (Fig. 3). This result suggests that, unlike in males, BPA may also alter the transcriptome profiles in the hippocampus of female offspring by suppressing the expression of these ASD-related TFs in addition to direct interactions with TFs.

Quantitative RT-PCR analysis of transcriptional targets of ASD-related TFs in the hippocampus of male and female offspring

To confirm the sex difference in the expression of genes that were transcriptional targets of ASD-related TFs and were dysregulated in the hippocampus of offspring, qRT-PCR analysis was performed using hippocampal tissues from rat pups (PND1) prenatally exposed to BPA (n = 6 pups; 3 males and 3 females from 3 different dams) or vehicle control (n = 6 pups; 3 males and 3 females from 3 different dams). A total of nine genes (i.e., Adnp, Chd2, Cic, Elavl2, Esrrb, Kat6a, Kmt2c, Nfib, and Snap25) were selected for qRT-PCR analysis. These genes have been identified by the TRANSFAC Curated database to be transcriptional targets of AR, ESR1, EGR2, SOX5, and YY1, which showed a binding free energy of less than − 5 kcal/mol when interacting with BPA (Table 4). Moreover, the expression of Ar, Esr1, Egr2, and Sox5 was downregulated in the hippocampus of female offspring in response to prenatal BPA exposure but not in males (Fig. 3). Information about each ASD-related TF and its transcriptional targets that were selected for qRT-PCR analysis in this study is provided in Additional file 16. When both sexes of pups in the same treatment group were combined into one group, we found that there was no significant difference in the expression of these TF target genes in the BPA group (Fig. 4a). However, when each sex of pups was analyzed separately, we found a sex difference in the effects of prenatal BPA exposure on the expression of these genes. The expression of Kmt2c, Esrrb, and Kat6a was significantly upregulated in males but not in females (Fig. 4b, c), whereas the expression of Snap25 and Chd2 was significantly decreased in females only (Fig. 4c). Notably, Kmt2c, Esrrb, Snap25, and Chd2 are transcriptional targets of AR. Kat6a is a transcriptional target of SOX5. In addition to these nine genes, we conducted a qRT-PCR analysis of the hippocampus of rat offspring prenatally exposed to BPA and vehicle control in our previous study and found that Auts2, Dicer1, Foxp2, and Smarcc2 were differentially expressed in response to BPA [13]. These genes are also transcriptional targets of the ASD-related TFs and, therefore, included in the subsequent analysis (Additional file 16). Notably, Auts2 and Foxp2 are transcriptional targets of AR.

To examine the relationships between the expression of ASD-related TFs and transcription targets that were differentially expressed in the offspring hippocampus, a correlation analysis was performed using the expression levels of ASD-related TFs (i.e., Ar, Esr1, Egr2, Sox5, and Yy1) and their transcriptional targets (i.e., Auts2, Chd2, Dicer1, Esrrb, Foxp2, Kmt2c, Smarcc2, and Snap25) from the qRT-PCR analysis (Fig. 5). Interestingly, the correlation analysis revealed sex-dependent relationships between these TFs and the transcriptional targets. The expression levels of AR and its targets (i.e., Auts2, Chd2, Esrrb, Foxp2, Kmt2c, and Snap25) were strongly correlated in a sex-specific pattern in response to prenatal BPA exposure, suggesting that AR is involved in the sex-specific effects of BPA on the expression levels of these ASD candidate genes.

AR is involved in the effects of BPA on the expression of ASD-candidate genes in the human neuroblastoma cell line

To further determine whether AR is also involved in the effects of BPA on the expression of ASD-candidate genes in human neuronal cells, the human neuroblastoma cell line SH-SY5Y stably transfected with AR-expression plasmid (pAR) or negative control plasmid (pNeg) was used. The expression levels of AR in the transfected cells were assessed using Western blot analysis (Additional file 4). The SH-SY5Y cells stably transfected with the negative control plasmid exhibited very low expression of AR proteins, while the cells transfected with the AR-expression plasmid showed a significant increase in AR protein expression. We also examined the expression of estrogen receptors in the transfected cells and found that the levels of ESR1 and ESR2 were not significantly changed (Additional file 17). These cells were differentiated into mature neurons and were treated with 1 ng/mL BPA (n = 5) or the vehicle control (n = 5) during differentiation. The concentration of BPA used for treatment in this experiment was calculated from an estimated BPA level in the fetus's brain when 5000 µg/kg·maternal body weight of BPA was given to a pregnant rat [78]. After 7 days of differentiation, RNA was extracted from the cells, and qRT-PCR analysis of selected ASD candidate genes (i.e., AUTS2, ESRRB, KMT2C, SMARCC2, and DICER1) was performed (Fig. 6). ESRRB expression in the SH-SY5Y cells transfected with pAR or pNeg plasmid was undetectable by qRT-PCR analysis and thus excluded from data analysis.

The expression patterns of AUTS2, KMT2C, and SMARCC2 in the cells with AR overexpression (pAR) and the cells with low AR levels (pNeg) were changed in response to BPA treatment (Fig. 6). Without BPA, the expression levels of AUTS2 and SMARCC2 tended to decrease but not significantly when AR was overexpressed (Fig. 6a, c). However, in the presence of BPA, the expression levels of these genes in the cells with AR overexpression were significantly reduced when compared to the cells with low AR expression levels. In contrast, KMT2C expression in cells stably transfected with the pAR plasmid was lower than that in cells transfected with the negative control plasmid (Fig. 6b). When treated with BPA, the expression of KMT2C in the cells with AR overexpression was significantly increased compared to cells treated with the vehicle control and exhibited no significant difference when compared to cells with low AR expression treated with BPA. These results suggest that AR, together with other mechanisms, regulates the expression of these ASD candidate genes, and BPA suppresses AUTS2 and SMARCC2 expression and enhances KMT2C expression through AR.

Sex differences in the effects of prenatal BPA exposure on synaptogenesis

Several AR targets (i.e., Adnp, Auts2, Snap25) that were selected for confirmation by qRT-PCR analysis and differentially expressed in the BPA group are known to play an important role in synaptogenesis and synaptic functions associated with ASD [79,80,81]. For instance, the reduction of the AUTS2 level leads to increased excitatory synapses in the brain [80]. To investigate the effects of prenatal BPA exposure on synaptogenesis, we isolated primary hippocampal neurons at PND1 from neonatal rat pups prenatally exposed to BPA or vehicle control. Primary hippocampal neurons were cultured for 14 days and immunostained for the neuronal marker protein Map2, the presynaptic protein Syn1, and the postsynaptic protein Psd95 (Fig. 7a). The fluorescence intensity of each marker and the colocalization of the presynaptic and postsynaptic markers were examined. We found that the expression levels of the presynaptic marker Syn1 and the postsynaptic marker Psd95 were significantly increased in primary hippocampal neurons from males prenatally exposed to BPA but not in those from females (Fig. 7b). Conversely, the number of colocalized Syn1 and Psd95 foci was significantly increased only in female primary hippocampal neurons (Fig. 7d). To further examine the expression of synaptic proteins in hippocampal tissues, Western blot analysis of hippocampal tissues from rat pups prenatally exposed to BPA (n = 6 pups; 3 males and 3 females from 3 different dams) or vehicle control (n = 6 pups; 3 males and 3 females from 3 different dams) was performed (Fig. 7e–n). When both sexes of pups in the same treatment group were combined, no changes in the expression of Map2 and these synaptic proteins were observed (data not shown). However, when each sex of pups was analyzed separately, the Western blot analysis showed significant upregulation of Syp, Syn1, and Psd95 proteins only in male hippocampal tissues (Fig. 7e–h) but not in female hippocampal tissues (Fig. 7j–m). No significant change in the expression of Map2 was observed in either sex (Fig. 7i, n). These results indicate that prenatal exposure to BPA causes sex-specific effects on the expression of synaptic proteins and synaptogenesis.

Reanalysis of transcriptome data from other BPA studies revealed dysregulation of ASD-related TF targets in the hippocampus of offspring

To determine whether ASD-related TFs might be involved in the effects of prenatal BPA exposure on transcriptome profiles in the offspring hippocampus in independent studies, we obtained transcriptome data from other BPA studies conducted by independent groups of researchers from the NCBI PubMed database and reanalyzed them as described in the Methods. The details of each study, including the title, sample size, and treatment, are described in Additional file 18. The lists of DEGs in the offspring hippocampus in response to prenatal/perinatal BPA exposure from other BPA studies are shown in Additional file 19. Hypergeometric distribution analysis was performed to assess the enrichment of transcriptional targets of ASD-related TFs in the lists of BPA-responsive genes in the offspring hippocampus from each transcriptomic study (Additional file 20). Transcriptional targets of AR, ESR1, EGR2, KDM5B, SOX5, SMAD4, TCF7L2, and YY1 were also found to be differentially expressed in the hippocampus of rat/mouse pups prenatally/perinatally exposed to BPA in independent studies, supporting that these TFs may be involved in the effects of BPA on gene expression. It is noteworthy that transcriptional targets of AR, ESR1, SMAD4, and YY1 were also found to be significantly enriched in the list of BPA-responsive genes in the hippocampus of male pups reported by Zhang et al. [82].

Discussion

The prevalence of ASD in males is approximately four times higher than that in females [17]. The exact cause and the mechanisms underlying the male bias of ASD are still unclear, but there is evidence that sex chromosomes and hormones are involved in ASD susceptibility. In the developing brain, sex hormones play an essential role in sex differences in brain structures, behaviors, cognition, and neurological functions associated with ASD [46,47,48,49]. Elevated concentrations of sex hormones, including testosterone, androstenedione, progesterone, and 17α-hydroxy-progesterone, were found in amniotic fluid samples of males who were later diagnosed with ASD [83]. The expression levels of estrogen receptors, aromatase (an enzyme that converts androgens to estrogen), and several estrogen receptor coactivators were also reduced in the frontal cortex of ASD individuals [40]. Moreover, our previous studies have shown that the expression of retinoic acid-related orphan receptor-alpha (RORA) was reduced in blood-derived cell lines and postmortem brain tissues from ASD individuals [84, 85]. This protein is differentially regulated by male and female sex hormones through AR and ESR1 and serves as an upstream regulator of several ASD-related genes, including CYP19A1, a gene encoding the aromatase enzyme that converts androgens to estrogen [25, 46, 48, 49]. Androgen negatively regulates the expression of RORA, whereas estrogen enhances the levels of RORA [46]. The findings from these studies serve as supporting evidence that sex hormones and their signaling are associated with ASD.

Both genetic and environmental factors can lead to abnormalities in sex hormones and their molecular signaling cascades during pregnancy, which may cause or increase the susceptibility of ASD in children [10, 13, 14, 18]. There is accumulating evidence that exposure to endocrine-disrupting chemicals, including BPA, phthalates, polychlorinated biphenyls (PCBs), and polybrominated diphenyl ethers (PBDEs), is an environmental risk factor for ASD [18, 86,87,88]. Our recent studies in rats have shown that maternal exposure to BPA during the gestation period altered the transcriptome profiles of genes associated with ASD and related neurological functions in the hippocampus and the frontal cortex of offspring in a sex-dependent manner [10, 13, 14]. Several ASD candidate genes, including Autism Susceptibility Gene 2 (Auts2) and Forkhead Box Protein P2 (Foxp2), were found to be significantly reduced in the hippocampus of male offspring but not in females in response to prenatal BPA exposure [13]. Moreover, we found that in utero BPA exposure increased neuritogenesis and reduced the size of the hippocampal cell body in both sexes of the offspring but decreased neuronal viability, the density of neuronal cells in the hippocampus, and learning/memory only in the male offspring [10]. Similarly, prenatal BPA exposure disrupted the transcriptome-interactome profiles of ASD candidate genes, including Auts2, Ankrd11, and Ntng1, in the offspring frontal cortex in a sex-specific pattern, possibly through AR, ESR1, and RORA [14]. These findings strongly suggest that BPA alters the expression of ASD-related genes and neurological functions in the offspring brain through sex-dependent mechanisms. However, the molecular mechanisms underlying the sex-specific effects of BPA have not been investigated. This is the first study to identify ASD-related TFs potentially involved in the sex differences in the effects of prenatal BPA exposure on the transcriptome profiles and neurological functions associated with ASD in the hippocampus of offspring. The major findings of our study are summarized in Fig. 8.

In this study, the rat dams were treated with BPA at 5000 µg/kg·maternal body weight or vehicle control from gestational day 1 until parturition by oral gavage. This dose of BPA is equal to the no-observed-adverse-effect level (NOAEL) for BPA exposure determined by the United States Food and Drug Administration (FDA) and the European Food Safety Authority (EFSA) [70]. Moreover, this dose was also used in our previous studies on the effects of prenatal BPA exposure on the transcriptome profiles of the hippocampus and hippocampal functions [10, 13], allowing for integrated interpretations of the results. The transcriptome profiling data from an RNA-seq analysis of hippocampal tissues isolated from male and female rat pups prenatally exposed to BPA or vehicle control were obtained from our previous study [13]. To determine whether BPA-responsive genes in the offspring hippocampus can be transcriptionally regulated by ASD-related TFs, we identified ASD-related TFs by overlapping the list of human TFs from the Human Transcription Factors database with the list of ASD candidate genes from the SFARI database (Additional file 5). In this study, we focused on ASD-related TFs, because there is evidence that several of them, including AR and ESR1, play an important role in the male bias of ASD [46,47,48,49] and thus serve as promising candidates for an investigation of the molecular mechanisms underpinning sex-specific effects of prenatal BPA exposure in the context of ASD. However, it is possible that, in addition to ASD-related TFs, BPA can also exert its effects through other TFs that have not been associated with ASD, which could be further investigated in the future.

Hypergeometric distribution analysis revealed that the transcriptional targets of AR, ESR1, KDM5B, SOX5, SMAD4, and YY1 were significantly enriched in the list of BPA-responsive genes in the hippocampus of male offspring. In contrast, the targets of EGR2, KDM5B, and SOX5 were significantly enriched in the list of BPA-responsive genes in females (Table 1). This result suggests that BPA may exert its effects on the transcriptome profiles in the hippocampus of male and female offspring through different TFs depending on the sex of the offspring. Although these TFs were obtained from a human transcription factor database, they are known to be highly conserved in mammals [89]. According to the Allen Brain Atlas database, all of these TFs, except Yy1, were also expressed in the developing brain of rodents during the prenatal stage (Additional file 15). Although there were no expression data for Yy1 in the database, our qRT-PCR analysis and results from other studies have confirmed that it is expressed in the developing brain (Fig. 3) [90, 91]. Moreover, molecular docking analysis showed that BPA can directly interact with these ASD-related TFs (Table 4). BPA is known to be an AR antagonist and an ESR1 agonist [92,93,94] and has been predicted to bind EGR2, SOX5, and YY1 directly [14]. However, direct interactions between BPA and KDM5B, SMAD4, or TCF7L2 have not been reported elsewhere, suggesting that these TFs are novel targets of BPA.

Gene ontology analysis by IPA software revealed that transcriptional targets of AR, KDM5B, and SMAD4 that were differentially expressed in response to BPA treatment were significantly associated with “autism spectrum disorder or intellectual disability” (Table 2, Additional files 8, 9, 10). Interestingly, the transcriptional targets of AR and SMAD4 were significantly enriched in the lists of BPA-responsive genes only in male offspring but not females, suggesting that BPA may exert its male-specific effects on neurological functions associated with ASD and increase the risk of ASD in male offspring through AR and SMAD4, which deserves further investigation. Moreover, the transcriptional targets of AR that were differentially expressed in the BPA group are significantly associated with “proliferation of neuronal cells,” “developmental process of synapse,” and “excitatory postsynaptic potential,” all of which have been linked to ASD [95,96,97]. Examples of AR targets associated with these neurological functions include Adnp, Chd2, Cic, Kmt2c, and Snap25. Interestingly, these genes are also ASD candidate genes in the SFARI database [98,99,100,101,102]. In addition, genetic abnormalities of these ASD-related TFs have been found in ASD individuals [44, 103,104,105,106,107]. Moreover, YY1 interacts with the methyl-CpG binding protein 2 (MECP2) protein and regulates the expression of ANT1 encoding a mitochondrial adenine nucleotide translocase enzyme associated with Rett syndrome which is an ASD-related neurodevelopmental disorder with loss of motor skills and speech [108]. Taken together, these findings suggest that altered functions of these TFs, either caused by genetic mutations or BPA exposure, may lead to dysregulations of ASD-related genes and altered neurological functions, which, in turn, cause or increase ASD susceptibility. The effects of BPA binding on the activities and functions of these TFs should be further investigated.

In addition to direct interactions between BPA and these ASD-related TFs, we also investigated the effects of prenatal BPA exposure on the expression of genes encoding these ASD-related TFs in the hippocampus of offspring (Fig. 3). The qRT-PCR analysis revealed a sex difference in the effects of prenatal BPA exposure on the expression levels of Ar, Esr1, Egr2, Sox5, Smad4, and Tcf7l2. These genes were significantly suppressed in the hippocampus of female offspring prenatally exposed to BPA but not in males. It is noteworthy that BPA exposure has been reported to reduce the expression of Esr1 in the hippocampus of offspring [109]. Chang et al. found that Sprague–Dawley rat pups exposed to BPA throughout gestation and lactation exhibited decreased Esr1 expression in the hippocampus, along with elevated DNA methylation of the Esr1 gene promoter [109]. In addition, Arambula et al. found that gestational exposure to BPA at 25,000 µg/kg maternal body weight tended to reduce the expression of Esr1 in the hippocampus of female offspring but not in males, albeit not statistically significant [32]. These findings suggest that, in addition to direct interaction with these TFs, BPA also alters the expression of these TFs through other transcriptional regulatory mechanisms (e.g., DNA methylation), which are sex dependent. However, such effects may also depend on the dose and time of BPA exposure.

We further investigated the expression of transcriptional targets of these TFs in the offspring hippocampus by qRT-PCR analysis and found that prenatal BPA exposure altered the expression levels of AR target genes in a sex-dependent manner (Fig. 4). The expression levels of Kmt2c and Esrrb were increased in males but not in females, whereas the expression levels of Snap25 and Chd2 were reduced in females only. Moreover, our previous study found that Auts2 and Foxp2, which are also transcriptional targets of AR, were also reduced in the hippocampus of male offspring prenatally exposed to BPA but not in females [13]. The correlation analysis revealed sex-specific relationships between the expression levels of AR and its targets (i.e., Auts2, Chd2, Esrrb, Foxp2, Kmt2c, and Snap25) (Fig. 5), suggesting that AR is involved in the sex-specific effects of BPA on the expression levels of these ASD candidate genes. Interestingly, the patterns of the correlations between the expression levels of TFs and their targets in the male offspring exposed to BPA tended to be more similar to those in the female control group than to the male control group. This finding suggests that prenatal BPA exposure may demasculinize the hippocampus of male pups through these TFs.

We then further investigated the role of AR in BPA-mediated effects on the expression of the AUTS2, KMT2C, and SMARCC2 genes in the human neuronal cell line SH-SY5Y. This cell line was used as a model, because it is a human cell line and has been widely used as a model for neuroscience research. Moreover, our western blot analysis showed low expression of AR protein in SH-SY5Y cells stably transfected with the negative control plasmid (Additional file 4), making it suitable for studying the role of AR in the expression of ASD candidate genes in the presence or absence of BPA. We found that BPA significantly suppressed the expression levels of AUTS2 and SMARCC2 but enhanced KMT2C in cells with high AR expression when compared to cells with low AR expression. These results suggest that BPA suppresses AUTS2 and SMARCC2 expression and enhances KMT2C expression in human neuronal cells through AR. In addition to sex hormone receptors, sex chromosomes are thought to play an important role in sex differences in the brain [110]. As SH-SY5Y cells are derived from a female patient with XX chromosomes, future studies should include human neuronal cell lines with XY chromosomes for comparisons.

Notably, AUTS2, SMARCC2, and KMT2C are known to be ASD-candidate genes [101, 111, 112]. AUTS2 (Autism Susceptibility Gene 2 or Activator of Transcription and Developmental Regulator) has been associated with ASD [111]. Mutations in the Auts2 gene cause developmental delay, intellectual disability, microcephaly, and short stature [113]. Li et al. revealed that mice with low expression of Auts2 showed hypoplasia in the dentate gyrus of the hippocampus [114]. Studies in zebrafish showed that Auts2 gene disruptions caused an increase in cell death during the neuronal differentiation stage [115]. SMARCC2 (SWI/SNF Related, Matrix Associated, Actin Dependent Regulator of Chromatin Subfamily C Member 2) encodes a member of the SWI/SNF family of proteins involved in transcriptional activation and repression by the chromatin remodeling process [116]. De novo variants in the SMARCC2 gene cause intellectual disability and developmental delay and have been associated with ASD [117]. KMT2C (lysine methyltransferase 2C) is an autism candidate gene and downstream target gene of AR. Kmt2c is involved in the histone methyltransferase process that catalyzes the methyl group to the amino acid group of Lys-4' of histone H3 (H3K4) [118]. Mutations in the Kmt2c gene have been reported in neurodevelopmental disorders, including ASD, Kleefstra syndrome, and intellectual disability [119, 120]. Frega et al. demonstrated that the deficiency in KMT2C, along with other genes (i.e., EHMT1, SMARCB1, and MDB5), led to hyperactivity in neuronal network communication and altered the excitatory and inhibitory balance [121].

Since transcriptional targets of ASD-related TFs, including AR, were associated with synaptic functions, we further performed synaptogenesis assays and found a sex difference in the effects of prenatal BPA exposure on synaptic protein expression and synaptogenesis (Fig. 7). Prenatal BPA exposure caused an increase in the expression levels of synaptic proteins in the hippocampal tissues of male offspring but not in those of females. Notably, this finding is consistent with previous studies, which reported that the loss of Auts2 led to increased excitatory synapses, disrupting the balance between excitatory and inhibitory synapses and causing cognitive and social communication deficits [80]. Our previous study also showed that prenatal BPA exposure reduced neuronal viability and neuronal density in the hippocampus and impaired learning/memory only in male offspring, while females were unaffected [10]. Some studies found that BPA exposure through gestational and postnatal stages reduced the total spine density in the hippocampus of male offspring [11, 12]. Moreover, BPA causes neuroinflammation by microglial activation [122], suggesting that BPA may negatively impact the microglia-mediated synaptic pruning process after birth. The effects of BPA on changes in synaptic formation and synaptic pruning processes in the offspring brain from embryonic to adult stages should be investigated further.

Differences in the expression of Ar and Esr1 in the hippocampus have also been investigated. Simerly et al. examined the distribution of Ar and Esr1-expressing neurons in the hippocampus and found that the expression of Ar was high, while the expression of Esr1 was low in the CA1 and CA2/3 of the hippocampus. In the dentate gyrus, Esr1 showed high expression, whereas the expression of Ar was low [123]. According to our previous study, we found that neuronal viability and neuronal density were significantly reduced in CA2/3 neurons of male pups prenatally exposed to BPA [10], suggesting that BPA might deregulate the expression of genes through AR, resulting in a reduced number of neurons in the hippocampus of male offspring. Moreover, previous studies in rodents have shown that young adult males performed better in the Morris water maze and radial arm maze than females (see reviews in [124]). When male offspring were exposed to the androgen receptor antagonist flutamide prenatally, they exhibited female-like spatial learning, suggesting that direct effects of androgen signaling are required [125]. Inasmuch as BPA is also known to be an AR antagonist similar to flutamide, it is possible that BPA reduced the learning/memory of male offspring through AR.

Perspectives and significance

In summary, we propose that maternal BPA exposure during gestation alters the transcriptome profiles in the hippocampus of offspring and neurological functions associated with ASD through sex-dependent mechanisms (Fig. 8). In male offspring, BPA disrupts the transcriptome profiles of genes associated with ASD in the hippocampus through direct interactions with AR and other ASD-related TFs (i.e., Esr1, Egr2, Kdm5b, Sox5, Smad4, Tcf7l2, and Yy1). In females, BPA not only directly interacts with ASD-related TFs but also dysregulates the gene expression of these TFs (i.e., Ar, Esr1, Egr2, Sox5, Smad4, and Tcf7l2) in the hippocampus. Altered TF-mediated signaling, together with other mechanisms (e.g., indirect signaling, genetic/sex chromosome-linked factors, sex hormones, and epigenetic mechanisms), causes sex-specific dysregulation of ASD-linked genes in males (i.e., Auts2, Foxp2, Kmt2c, Esrrb, and Kat6a) and in females (i.e., Chd2, Snap25, and Smarcc2) which, in turn, leads to altered neurological functions associated with ASD, particularly synaptogenesis. Male-specific downregulation of certain AR target genes, including Auts2 and Foxp2, can disrupt excitatory synapses in the hippocampus, causing deficits in learning/memory and other ASD-related symptoms observed in male offspring. This study leads to a better understanding of sex differences in BPA-mediated mechanisms in the offspring hippocampus, which provides a plausible explanation of the role of environmental factors, particularly endocrine-disrupting chemicals, on the male bias of ASD and may lead to the development of therapeutic targets for BPA-related ASD in the future.

Availability of data and materials

The transcriptome profiling data used in this study have been deposited in the NCBI GEO data set database (GSE140298).

References

Biedermann S, Tschudin P, Grob K. Transfer of bisphenol A from thermal printer paper to the skin. Anal Bioanal Chem. 2010;398(1):571–6.
Article CAS PubMed Google Scholar
Noonan GO, Ackerman LK, Begley TH. Concentration of bisphenol A in highly consumed canned foods on the US market. J Agric Food Chem. 2011;59(13):7178–85.
Article CAS PubMed Google Scholar
Prüst M, Meijer J, Westerink RH. The plastic brain: neurotoxicity of micro-and nanoplastics. Part Fibre Toxicol. 2020;17(1):1–16.
Article Google Scholar
Calafat AM, Ye X, Wong L-Y, Reidy JA, Needham LL. Exposure of the US population to bisphenol A and 4-tertiary-octylphenol: 2003–2004. Environ Health Perspect. 2008;116(1):39–44.
Article CAS PubMed Google Scholar
Hansen JB, Bilenberg N, Timmermann CAG, Jensen RC, Frederiksen H, Andersson A-M, et al. Prenatal exposure to bisphenol A and autistic-and ADHD-related symptoms in children aged 2 and5 years from the Odense Child Cohort. Environ Health. 2021;20(1):1–12.
Article Google Scholar
Mao J, Jain A, Denslow ND, Nouri M-Z, Chen S, Wang T, et al. Bisphenol A and bisphenol S disruptions of the mouse placenta and potential effects on the placenta–brain axis. Proc Natl Acad Sci. 2020;117(9):4642–52.
Article CAS PubMed PubMed Central Google Scholar
Engdahl E, van Schijndel MD, Voulgaris D, Di Criscio M, Ramsbottom KA, Rigden DJ, et al. Bisphenol A inhibits the transporter function of the blood–brain barrier by directly interacting with the ABC transporter breast cancer resistance protein (BCRP). Int J Mol Sci. 2021;22(11):5534.
Article CAS PubMed PubMed Central Google Scholar
Kim K, Son TG, Park HR, Kim SJ, Kim HS, Kim HS, et al. Potencies of bisphenol A on the neuronal differentiation and hippocampal neurogenesis. J Toxicol Environ Health A. 2009;72(21–22):1343–51.
Article CAS PubMed Google Scholar
Ling W, Endo T, Kubo K-I, Nakajima K, Kakeyama M, Tohyama C. In utero bisphenol A exposure induces abnormal neuronal migration in the cerebral cortex of mice. Front Endocrinol. 2016;7:7.
Article Google Scholar
Thongkorn S, Kanlayaprasit S, Panjabud P, Saeliw T, Jantheang T, Kasitipradit K, et al. Sex differences in the effects of prenatal bisphenol A exposure on autism-related genes and their relationships with the hippocampus functions. Sci Rep. 2021;11(1):1241.
Article CAS PubMed PubMed Central Google Scholar
Kawato S, Ogiue-Ikeda M, Soma M, Yoshino H, Kominami T, Saito M, et al. Perinatal exposure of bisphenol A differently affects dendritic spines of male and female grown-up adult hippocampal neurons. Front Neurosci. 2021;15:1040.
Article Google Scholar
Zhang H, Wang Z, Meng L, Kuang H, Liu J, Lv X, et al. Maternal exposure to environmental bisphenol A impairs the neurons in hippocampus across generations. Toxicology. 2020;432: 152393.
Article CAS PubMed Google Scholar
Thongkorn S, Kanlayaprasit S, Jindatip D, Tencomnao T, Hu VW, Sarachana T. Sex differences in the effects of prenatal bisphenol A exposure on genes associated with autism spectrum disorder in the hippocampus. Sci Rep. 2019;9(1):1–14.
Article CAS Google Scholar
Kanlayaprasit S, Thongkorn S, Panjabud P, Jindatip D, Hu VW, Kikkawa T, et al. Autism-related transcription factors underlying the sex-specific effects of prenatal bisphenol A exposure on transcriptome-interactome profiles in the offspring prefrontal cortex. Int J Mol Sci. 2021;22(24):13201.
Article CAS PubMed PubMed Central Google Scholar
Sukjamnong S, Thongkorn S, Kanlayaprasit S, Saeliw T, Hussem K, Warayanon W, et al. Prenatal exposure to bisphenol A alters the transcriptome-interactome profiles of genes associated with Alzheimer’s disease in the offspring hippocampus. Sci Rep. 2020;10(1):1–15.
Article Google Scholar
Maenner MJ, Shaw KA, Bakian AV, Bilder DA, Durkin MS, Esler A, et al. Prevalence and characteristics of autism spectrum disorder among children aged 8 years—autism and developmental disabilities monitoring network, 11 sites, United States, 2018. MMWR Surveill Summ. 2021;70(11):1.
Article PubMed PubMed Central Google Scholar
Maenner MJ, Shaw KA, Baio J. Prevalence of autism spectrum disorder among children aged 8 years—autism and developmental disabilities monitoring network, 11 sites, United States, 2016. MMWR Surveill Summ. 2020;69(4):1.
Article PubMed PubMed Central Google Scholar
Moosa A, Shu H, Sarachana T, Hu VW. Are endocrine disrupting compounds environmental risk factors for autism spectrum disorder? Horm Behav. 2018;101:13–21.
Article CAS PubMed Google Scholar
Saeliw T, Tangsuwansri C, Thongkorn S, Chonchaiya W, Suphapeetiporn K, Mutirangura A, et al. Integrated genome-wide Alu methylation and transcriptome profiling analyses reveal novel epigenetic regulatory networks associated with autism spectrum disorder. Mol Autism. 2018;9(1):1–19.
Article Google Scholar
Tangsuwansri C, Saeliw T, Thongkorn S, Chonchaiya W, Suphapeetiporn K, Mutirangura A, et al. Investigation of epigenetic regulatory networks associated with autism spectrum disorder (ASD) by integrated global LINE-1 methylation and gene expression profiling analyses. PLoS ONE. 2018;13(7): e0201071.
Article PubMed PubMed Central Google Scholar
Pichitpunpong C, Thongkorn S, Kanlayaprasit S, Yuwattana W, Plaingam W, Sangsuthum S, et al. Phenotypic subgrouping and multi-omics analyses reveal reduced diazepam-binding inhibitor (DBI) protein levels in autism spectrum disorder with severe language impairment. PLoS ONE. 2019;14(3): e0214198.
Article CAS PubMed PubMed Central Google Scholar
Saeliw T, Permpoon T, Iadsee N, Tencomnao T, Hu VW, Sarachana T, et al. LINE-1 and Alu methylation signatures in autism spectrum disorder and their associations with the expression of autism-related genes. Sci Rep. 2022;12(1):1–16.
Article Google Scholar
Lasalle JM. Autism genes keep turning up chromatin. OA autism. 2013;1(2):14.
Article PubMed PubMed Central Google Scholar
Sun W, Poschmann J, DelRosario RC-H, Parikshak NN, Hajan HS, Kumar V, et al. Histone acetylome-wide association study of autism spectrum disorder. Cell. 2016;167(5):1385-97e11.
Article CAS PubMed Google Scholar
Sarachana T, Zhou R, Chen G, Manji HK, Hu VW. Investigation of post-transcriptional gene regulatory networks associated with autism spectrum disorders by microRNA expression profiling of lymphoblastoid cell lines. Genome Med. 2010;2(4):1–18.
Article Google Scholar
Kardas F, Bayram AK, Demirci E, Akin L, Ozmen S, Kendirci M, et al. Increased serum phthalates (MEHP, DEHP) and bisphenol A concentrations in children with autism spectrum disorder: the role of endocrine disruptors in autism etiopathogenesis. J Child Neurol. 2016;31(5):629–35.
Article PubMed Google Scholar
Stein TP, Schluter MD, Steer RA, Guo L, Ming X. Bisphenol A exposure in children with autism spectrum disorders. Autism Res. 2015;8(3):272–83.
Article PubMed PubMed Central Google Scholar
Kondolot M, Ozmert EN, Ascı A, Erkekoglu P, Oztop DB, Gumus H, et al. Plasma phthalate and bisphenol a levels and oxidant-antioxidant status in autistic children. Environ Toxicol Pharmacol. 2016;43:149–58.
Article CAS PubMed Google Scholar
Oulhote Y, Lanphear B, Braun JM, Webster GM, Arbuckle TE, Etzel T, et al. Gestational exposures to phthalates and folic acid, and autistic traits in Canadian children. Environ Health Perspect. 2020;128(2): 027004.
Article CAS PubMed PubMed Central Google Scholar
Carter CJ, Blizard R. Autism genes are selectively targeted by environmental pollutants including pesticides, heavy metals, bisphenol A, phthalates and many others in food, cosmetics or household products. Neurochem Int. 2016;101:83–109.
Article CAS Google Scholar
Kalkbrenner AE, Braun JM, Durkin MS, Maenner MJ, Cunniff C, Lee L-C, et al. Maternal smoking during pregnancy and the prevalence of autism spectrum disorders, using data from the autism and developmental disabilities monitoring network. Environ Health Perspect. 2012;120(7):1042–8.
Article PubMed PubMed Central Google Scholar
Arambula SE, Belcher SM, Planchart A, Turner SD, Patisaul HB. Impact of low dose oral exposure to bisphenol A (BPA) on the neonatal rat hypothalamic and hippocampal transcriptome: a CLARITY-BPA consortium study. Endocrinology. 2016;157(10):3856–72.
Article CAS PubMed PubMed Central Google Scholar
Mathisen GH, Yazdani M, Rakkestad KE, Aden PK, Bodin J, Samuelsen M, et al. Prenatal exposure to bisphenol A interferes with the development of cerebellar granule neurons in mice and chicken. Int J Dev Neurosci. 2013;31(8):762–9.
Article CAS PubMed Google Scholar
Arambula SE, Jima D, Patisaul HB. Prenatal bisphenol A (BPA) exposure alters the transcriptome of the neonate rat amygdala in a sex-specific manner: a CLARITY-BPA consortium study. Neurotoxicology. 2018;65:207–20.
Article CAS PubMed Google Scholar
Schumann CM, Hamstra J, Goodlin-Jones BL, Lotspeich LJ, Kwon H, Buonocore MH, et al. The amygdala is enlarged in children but not adolescents with autism; the hippocampus is enlarged at all ages. J Neurosci. 2004;24(28):6392–401.
Article CAS PubMed PubMed Central Google Scholar
Banker SM, Gu X, Schiller D, Foss-Feig JH. Hippocampal contributions to social and cognitive deficits in autism spectrum disorder. Trends Neurosci. 2021;44(10):793–807.
Article CAS PubMed PubMed Central Google Scholar
Wang Y, Du X, Wang D, Wang J, Du J. Effects of bisphenol A exposure during pregnancy and lactation on hippocampal function in newborn rats. Int J Med Sci. 2020;17(12):1751.
Article CAS PubMed PubMed Central Google Scholar
Dong D, Zielke HR, Yeh D, Yang P. Cellular stress and apoptosis contribute to the pathogenesis of autism spectrum disorder. Autism Res. 2018;11(7):1076–90.
Article PubMed PubMed Central Google Scholar
Winslow JT, Insel TR. The social deficits of the oxytocin knockout mouse. Neuropeptides. 2002;36(2–3):221–9.
Article CAS PubMed Google Scholar
Crider A, Thakkar R, Ahmed AO, Pillai A. Dysregulation of estrogen receptor beta (ERβ), aromatase (CYP19A1), and ER co-activators in the middle frontal gyrus of autism spectrum disorder subjects. Mol Autism. 2014;5(1):1–10.
Article Google Scholar
Hyun S-A, Ko MY, Jang S, Lee B-S, Rho J, Kim KK, et al. Bisphenol-A impairs synaptic formation and function by RGS4-mediated regulation of BDNF signaling in the cerebral cortex. Disease Models Mech. 2022;15(7):49177.
Article Google Scholar
Liu Z-H, Ding J-J, Yang Q-Q, Song H-Z, Chen X-T, Xu Y, et al. Early developmental bisphenol-A exposure sex-independently impairs spatial memory by remodeling hippocampal dendritic architecture and synaptic transmission in rats. Sci Rep. 2016;6(1):1–11.
Google Scholar
Lombardo MV, Auyeung B, Pramparo T, Quartier A, Courraud J, Holt RJ, et al. Sex-specific impact of prenatal androgens on social brain default mode subsystems. Mol Psychiatry. 2020;25(9):2175–88.
Article CAS PubMed Google Scholar
Henningsson S, Jonsson L, Ljunggren E, Westberg L, Gillberg C, Råstam M, et al. Possible association between the androgen receptor gene and autism spectrum disorder. Psychoneuroendocrinology. 2009;34(5):752–61.
Article CAS PubMed Google Scholar
Hu VW, Nguyen A, Kim KS, Steinberg ME, Sarachana T, Scully MA, et al. Gene expression profiling of lymphoblasts from autistic and nonaffected sib pairs: altered pathways in neuronal development and steroid biosynthesis. PLoS ONE. 2009;4(6): e5775.
Article PubMed PubMed Central Google Scholar
Sarachana T, Xu M, Wu R-C, Hu VW. Sex hormones in autism: androgens and estrogens differentially and reciprocally regulate RORA, a novel candidate gene for autism. PLoS ONE. 2011;6(2): e17116.
Article CAS PubMed PubMed Central Google Scholar
Sarachana T, Hu VW. Genome-wide identification of transcriptional targets of RORA reveals direct regulation of multiple genes associated with autism spectrum disorder. Mol Autism. 2013;4(1):1–19.
Article Google Scholar
Sarachana T, Hu VW. Differential recruitment of coregulators to the RORA promoter adds another layer of complexity to gene (dys) regulation by sex hormones in autism. Mol Autism. 2013;4(1):1–17.
Article Google Scholar
Hu VW, Sarachana T, Sherrard RM, Kocher KM. Investigation of sex differences in the expression of RORA and its transcriptional targets in the brain as a potential contributor to the sex bias in autism. Mol Autism. 2015;6(1):1–19.
Article Google Scholar
Lambert SA, Jolma A, Campitelli LF, Das PK, Yin Y, Albu M, et al. The human transcription factors. Cell. 2018;172(4):650–65.
Article CAS PubMed Google Scholar
Abrahams BS, Arking DE, Campbell DB, Mefford HC, Morrow EM, Weiss LA, et al. SFARI Gene 2.0: a community-driven knowledgebase for the autism spectrum disorders (ASDs). Mol Autism. 2013;4(1):1–3.
Article Google Scholar
Oliveros JC. Venny. An interactive tool for comparing lists with Venn's diagrams. 2007–2015. Available from: https://bioinfogp.cnb.csic.es/tools/venny/index.html.
Matys V, Fricke E, Geffers R, Gößling E, Haubrock M, Hehl R, et al. TRANSFAC^®: transcriptional regulation, from patterns to profiles. Nucleic Acids Res. 2003;31(1):374–8.
Article CAS PubMed PubMed Central Google Scholar
Matys V, Kel-Margoulis OV, Fricke E, Liebich I, Land S, Barre-Dirrie A, et al. TRANSFAC^® and its module TRANSCompel^®: transcriptional gene regulation in eukaryotes. Nucleic Acids Res. 2006;34(Suppl 1):D108–10.
Article CAS PubMed Google Scholar
Mathelier A, Zhao X, Zhang AW, Parcy F, Worsley-Hunt R, Arenillas DJ, et al. JASPAR 2014: an extensively expanded and updated open-access database of transcription factor binding profiles. Nucleic Acids Res. 2014;42(D1):D142–7.
Article CAS PubMed Google Scholar
Sandelin A, Alkema W, Engström P, Wasserman WW, Lenhard B. JASPAR: an open-access database for eukaryotic transcription factor binding profiles. Nucleic Acids Res. 2004;32(Suppl 1):D91–4.
Article CAS PubMed PubMed Central Google Scholar
Lachmann A, Xu H, Krishnan J, Berger SI, Mazloom AR, Ma’ayan A. ChEA: transcription factor regulation inferred from integrating genome-wide ChIP-X experiments. Bioinformatics. 2010;26(19):2438–44.
Article CAS PubMed PubMed Central Google Scholar
Feingold E, Pachter L. The ENCODE (ENCyclopedia of DNA elements) project. Science. 2004;306(5696):636–40.
Article Google Scholar
Consortium EP. A user’s guide to the encyclopedia of DNA elements (ENCODE). PLoS Biol. 2011;9(4): e1001046.
Article Google Scholar
Xie X, Rigor P, Baldi P. MotifMap: a human genome-wide map of candidate regulatory motif sites. Bioinformatics. 2009;25(2):167–74.
Article CAS PubMed Google Scholar
Rouillard AD, Gundersen GW, Fernandez NF, Wang Z, Monteiro CD, McDermott MG, et al. The harmonizome: a collection of processed datasets gathered to serve and mine knowledge about genes and proteins. Database. 2016;2016:baw100.
Article PubMed PubMed Central Google Scholar
Calculator KO. Hypergeometric distribution Calculator 2022. Available from: http://keisan.casio.com/exec/system/1180573201.
Krämer A, Green J, Pollard J Jr, Tugendreich S. Causal analysis approaches in ingenuity pathway analysis. Bioinformatics. 2014;30(4):523–30.
Article PubMed Google Scholar
Lein ES, Hawrylycz MJ, Ao N, Ayres M, Bensinger A, Bernard A, et al. Genome-wide atlas of gene expression in the adult mouse brain. Nature. 2007;445(7124):168–76.
Article CAS PubMed Google Scholar
Kim S, Chen J, Cheng T, Gindulyte A, He J, He S, et al. PubChem in 2021: new data content and improved web interfaces. Nucleic Acids Res. 2021;49(D1):D1388–95.
Article CAS PubMed Google Scholar
Dassault Systèmes. BIOVIA Discovery Studio Visualizer 2020. San Diego: Dassault Systèmes; 2022.
BIOVIA Discovery Studio Visualizer 2020 [Internet]. 2022.
Berman HM, Westbrook J, Feng Z, Gilliland G, Bhat TN, Weissig H, et al. The protein data bank. Nucleic Acids Res. 2000;28(1):235–42.
Article CAS PubMed PubMed Central Google Scholar
Morris GM, Huey R, Lindstrom W, Sanner MF, Belew RK, Goodsell DS, et al. AutoDock4 and AutoDockTools4: automated docking with selective receptor flexibility. J Comput Chem. 2009;30(16):2785–91.
Article CAS PubMed PubMed Central Google Scholar
Organization WH. Food and Agriculture Organization of United Nations: Bisphenol A (BPA) Current state of knowledge and future actions by WHO and FAO. International Food Safety Authorities Network (INFOSAN). 2009.
McCarthy MM. Incorporating sex as a variable in preclinical neuropsychiatric research. Schizophr Bull. 2015;41(5):1016–20.
Article PubMed PubMed Central Google Scholar
Beaudoin GM, Lee S-H, Singh D, Yuan Y, Ng Y-G, Reichardt LF, et al. Culturing pyramidal neurons from the early postnatal mouse hippocampus and cortex. Nat Protoc. 2012;7(9):1741–54.
Article CAS PubMed Google Scholar
Schneider CA, Rasband WS, Eliceiri KW. NIH Image to ImageJ: 25 years of image analysis. Nat Methods. 2012;9(7):671–5.
Article CAS PubMed PubMed Central Google Scholar
Wheeler DL, Barrett T, Benson DA, Bryant SH, Canese K, Chetvernin V, et al. Database resources of the national center for biotechnology information. Nucleic Acids Res. 2007;36(Suppl 1):D13–21.
Article PubMed PubMed Central Google Scholar
Jalili V, Afgan E, Gu Q, Clements D, Blankenberg D, Goecks J, et al. The Galaxy platform for accessible, reproducible and collaborative biomedical analyses: 2020 update. Nucleic Acids Res. 2020;48(W1):W395–402.
Article CAS PubMed PubMed Central Google Scholar
SPSS Inc. IBM SPSS Statistics for Windows. Armonk, NY: IBM Corp2022.
Xu XB, He Y, Song C, Ke X, Fan SJ, Peng WJ, et al. Bisphenol A regulates the estrogen receptor alpha signaling in developing hippocampus of male rats through estrogen receptor. Hippocampus. 2014;24(12):1570–80.
Article CAS PubMed Google Scholar
Domoradzki J, Pottenger L, Thornton C, Hansen S, Card T, Markham D, et al. Metabolism and pharmacokinetics of bisphenol A (BPA) and the embryo-fetal distribution of BPA and BPA-monoglucuronide in CD Sprague-Dawley rats at three gestational stages. Toxicol Sci. 2003;76(1):21–34.
Article CAS PubMed Google Scholar
Sragovich S, Malishkevich A, Piontkewitz Y, Giladi E, Touloumi O, Lagoudaki R, et al. The autism/neuroprotection-linked ADNP/NAP regulate the excitatory glutamatergic synapse. Transl Psychiatry. 2019;9(1):1–14.
Article CAS Google Scholar
Hori K, Yamashiro K, Nagai T, Shan W, Egusa SF, Shimaoka K, et al. AUTS2 regulation of synapses for proper synaptic inputs and social communication. Iscience. 2020;23(6): 101183.
Article CAS PubMed PubMed Central Google Scholar
Fossati G, Morini R, Corradini I, Antonucci F, Trepte P, Edry E, et al. Reduced SNAP-25 increases PSD-95 mobility and impairs spine morphogenesis. Cell Death Differ. 2015;22(9):1425–36.
Article CAS PubMed PubMed Central Google Scholar
Zhang Z, Wang H, Lei X, Ommati MM, Tang Z, Yuan J. Bisphenol a exposure decreases learning ability through the suppression of mitochondrial oxidative phosphorylation in the hippocampus of male mice. Food Chem Toxicol. 2022;165:113167.
Article CAS PubMed Google Scholar
Baron-Cohen S, Auyeung B, Nørgaard-Pedersen B, Hougaard DM, Abdallah MW, Melgaard L, et al. Elevated fetal steroidogenic activity in autism. Mol Psychiatry. 2015;20(3):369–76.
Article CAS PubMed Google Scholar
Hu VW, Sarachana T, Kim KS, Nguyen A, Kulkarni S, Steinberg ME, et al. Gene expression profiling differentiates autism case–controls and phenotypic variants of autism spectrum disorders: evidence for circadian rhythm dysfunction in severe autism. Autism Res. 2009;2(2):78–97.
Article PubMed PubMed Central Google Scholar
Nguyen A, Rauch TA, Pfeifer GP, Hu VW. Global methylation profiling of lymphoblastoid cell lines reveals epigenetic contributions to autism spectrum disorders and a novel autism candidate gene, RORA, whose protein product is reduced in autistic brain. FASEB J. 2010;24(8):3036–51.
Article CAS PubMed PubMed Central Google Scholar
Lyall K, Croen LA, Sjödin A, Yoshida CK, Zerbo O, Kharrazi M, et al. Polychlorinated biphenyl and organochlorine pesticide concentrations in maternal mid-pregnancy serum samples: association with autism spectrum disorder and intellectual disability. Environ Health Perspect. 2017;125(3):474–80.
Article CAS PubMed Google Scholar
Braun JM, Kalkbrenner AE, Just AC, Yolton K, Calafat AM, Sjödin A, et al. Gestational exposure to endocrine-disrupting chemicals and reciprocal social, repetitive, and stereotypic behaviors in 4-and 5-year-old children: the HOME study. Environ Health Perspect. 2014;122(5):513–20.
Article PubMed PubMed Central Google Scholar
Sagiv SK, Kogut K, Gaspar FW, Gunier RB, Harley KG, Parra K, et al. Prenatal and childhood polybrominated diphenyl ether (PBDE) exposure and attention and executive function at 9–12 years of age. Neurotoxicol Teratol. 2015;52:151–61.
Article CAS PubMed PubMed Central Google Scholar
Nitta KR, Jolma A, Yin Y, Morgunova E, Kivioja T, Akhtar J, et al. Conservation of transcription factor binding specificities across 600 million years of bilateria evolution. Elife. 2015;4:e04837.
Article PubMed PubMed Central Google Scholar
Dong X, Kwan KM. Yin Yang 1 is critical for mid-hindbrain neuroepithelium development and involved in cerebellar agenesis. Mol Brain. 2020;13(1):1–18.
Article Google Scholar
Zurkirchen L, Varum S, Giger S, Klug A, Häusel J, Bossart R, et al. Yin Yang 1 sustains biosynthetic demands during brain development in a stage-specific manner. Nat Commun. 2019;10(1):1–16.
Article CAS Google Scholar
Lee HJ, Chattopadhyay S, Gong E-Y, Ahn RS, Lee K. Antiandrogenic effects of bisphenol A and nonylphenol on the function of androgen receptor. Toxicol Sci. 2003;75(1):40–6.
Article CAS PubMed Google Scholar
Wang H, Ding Z, Shi Q-M, Ge X, Wang H-X, Li M-X, et al. Anti-androgenic mechanisms of Bisphenol A involve androgen receptor signaling pathway. Toxicology. 2017;387:10–6.
Article CAS PubMed Google Scholar
Kurosawa T, Hiroi H, Tsutsumi O, Ishikawa T, Osuga Y, Fujiwara T, et al. The activity of bisphenol A depends on both the estrogen receptor subtype and the cell type. Endocr J. 2002;49(4):465–71.
Article CAS PubMed Google Scholar
Antoine MW, Langberg T, Schnepel P, Feldman DE. Increased excitation-inhibition ratio stabilizes synapse and circuit excitability in four autism mouse models. Neuron. 2019;101(4):648-61e4.
Article CAS PubMed PubMed Central Google Scholar
Sohal VS, Rubenstein JL. Excitation-inhibition balance as a framework for investigating mechanisms in neuropsychiatric disorders. Mol Psychiatry. 2019;24(9):1248–57.
Article PubMed PubMed Central Google Scholar
Marchetto MC, Belinson H, Tian Y, Freitas BC, Fu C, Vadodaria K, et al. Altered proliferation and networks in neural cells derived from idiopathic autistic individuals. Mol Psychiatry. 2017;22(6):820–35.
Article CAS PubMed Google Scholar
Yuen CRK, Merico D, Bookman M, Howe LJ, Thiruvahindrapuram B, Patel RV, et al. Whole genome sequencing resource identifies 18 new candidate genes for autism spectrum disorder. Nature Neurosci. 2017;20(4):602–11.
Article CAS Google Scholar
Guerini FR, Bolognesi E, Chiappedi M, Manca S, Ghezzo A, Agliardi C, et al. SNAP-25 single nucleotide polymorphisms are associated with hyperactivity in autism spectrum disorders. Pharmacol Res. 2011;64(3):283–8.
Article CAS PubMed Google Scholar
Iossifov I, Levy D, Allen J, Ye K, Ronemus M, Lee Y-H, et al. Low load for disruptive mutations in autism genes and their biased transmission. Proc Nat Acad Sci. 2015;112(41):E5600–7.
Article CAS PubMed PubMed Central Google Scholar
O’Roak BJ, Vives L, Girirajan S, Karakoc E, Krumm N, Coe BP, et al. Sporadic autism exomes reveal a highly interconnected protein network of de novo mutations. Nature. 2012;485(7397):246–50.
Article PubMed PubMed Central Google Scholar
Vandeweyer G, Helsmoortel C, Van Dijck A, Vulto‐van Silfhout AT, Coe BP, Bernier R, et al., editors. The transcriptional regulator ADNP links the BAF (SWI/SNF) complexes with autism. American Journal of Medical Genetics Part C: Seminars in Medical Genetics; 2014: Wiley Online Library.
Gabriele M, Vulto-van Silfhout AT, Germain P-L, Vitriolo A, Kumar R, Douglas E, et al. YY1 haploinsufficiency causes an intellectual disability syndrome featuring transcriptional and chromatin dysfunction. Am J Hum Genet. 2017;100(6):907–25.
Article CAS PubMed PubMed Central Google Scholar
Iossifov I, O’roak BJ, Sanders SJ, Ronemus M, Krumm N, Levy D, et al. The contribution of de novo coding mutations to autism spectrum disorder. Nature. 2014;515(7526):216–21.
Article CAS PubMed PubMed Central Google Scholar
Lamb AN, Rosenfeld JA, Neill NJ, Talkowski ME, Blumenthal I, Girirajan S, et al. Haploinsufficiency of SOX5 at 12p12.1 is associated with developmental delays with prominent language delay, behavior problems, and mild dysmorphic features. Hum Mutation. 2012;33(4):728–40.
Article CAS Google Scholar
Swanberg SE, Nagarajan RP, Peddada S, Yasui DH, LaSalle JM. Reciprocal co-regulation of EGR2 and MECP2 is disrupted in Rett syndrome and autism. Hum Mol Genet. 2009;18(3):525–34.
Article CAS PubMed Google Scholar
Zettergren A, Jonsson L, Johansson D, Melke J, Lundström S, Anckarsäter H, et al. Associations between polymorphisms in sex steroid related genes and autistic-like traits. Psychoneuroendocrinology. 2013;38(11):2575–84.
Article CAS PubMed Google Scholar
Forlani G, Giarda E, Ala U, Di Cunto F, Salani M, Tupler R, et al. The MeCP2/YY1 interaction regulates ANT1 expression at 4q35: novel hints for Rett syndrome pathogenesis. Hum Mol Genet. 2010;19(16):3114–23.
Article CAS PubMed PubMed Central Google Scholar
Chang H, Wang M, Xia W, Chen T, Huo W, Mao Z, et al. Perinatal exposure to low-dose bisphenol A disrupts learning/memory and DNA methylation of estrogen receptor alpha in the hippocampus. Toxicology Research. 2016;5(3):828–35.
Article CAS PubMed PubMed Central Google Scholar
Arnold AP. Sex chromosomes and brain gender. Nat Rev Neurosci. 2004;5(9):701–8.
Article CAS PubMed Google Scholar
Sultana R, Yu C-E, Yu J, Munson J, Chen D, Hua W, et al. Identification of a novel gene on chromosome 7q11.2 interrupted by a translocation breakpoint in a pair of autistic twins. Genomics. 2002;80(2):129–34.
Article CAS PubMed Google Scholar
Neale BM, Kou Y, Liu L, Ma’Ayan A, Samocha KE, Sabo A, et al. Patterns and rates of exonic de novo mutations in autism spectrum disorders. Nature. 2012;485(7397):242–5.
Article CAS PubMed PubMed Central Google Scholar
Beunders G, Van De Kamp J, Vasudevan P, Morton J, Smets K, Kleefstra T, et al. A detailed clinical analysis of 13 patients with AUTS2 syndrome further delineates the phenotypic spectrum and underscores the behavioural phenotype. J Med Genet. 2016;53(8):523–32.
Article CAS PubMed Google Scholar
Li J, Sun X, You Y, Li Q, Wei C, Zhao L, et al. Auts2 deletion involves in DG hypoplasia and social recognition deficit: The developmental and neural circuit mechanisms. Sci Adv. 2022;8(9):eabk1238.
Article CAS PubMed PubMed Central Google Scholar
Monderer-Rothkoff G, Tal N, Risman M, Shani O, Nissim-Rafinia M, Malki-Feldman L, et al. AUTS2 isoforms control neuronal differentiation. Mol Psychiatry. 2021;26(2):666–81.
Article PubMed Google Scholar
Kadam S, McAlpine GS, Phelan ML, Kingston RE, Jones KA, Emerson BM. Functional selectivity of recombinant mammalian SWI/SNF subunits. Genes Dev. 2000;14(19):2441–51.
Article CAS PubMed PubMed Central Google Scholar
Machol K, Rousseau J, Ehresmann S, Garcia T, Nguyen TTM, Spillmann RC, et al. Expanding the spectrum of BAF-related disorders: de novo variants in SMARCC2 cause a syndrome with intellectual disability and developmental delay. Am J Hum Genet. 2019;104(1):164–78.
Article CAS PubMed Google Scholar
Lavery WJ, Barski A, Wiley S, Schorry EK, Lindsley AW. KMT2C/D COMPASS complex-associated diseases [KCDCOM-ADs]: an emerging class of congenital regulopathies. Clin Epigenetics. 2020;12(1):1–20.
Article Google Scholar
Chen C-H, Huang A, Huang Y-S, Fang T-H. Identification of a rare novel KMT2C mutation that presents with schizophrenia in a multiplex family. J Personalized Med. 2021;11(12):1254.
Article CAS Google Scholar
Koemans TS, Kleefstra T, Chubak MC, Stone MH, Reijnders MR, de Munnik S, et al. Functional convergence of histone methyltransferases EHMT1 and KMT2C involved in intellectual disability and autism spectrum disorder. PLoS Genet. 2017;13(10): e1006864.
Article PubMed PubMed Central Google Scholar
Frega M, Selten M, Mossink B, Keller JM, Linda K, Moerschen R, et al. Distinct pathogenic genes causing intellectual disability and autism exhibit a common neuronal network hyperactivity phenotype. Cell Rep. 2020;30(1):173-86.e6.
Article CAS PubMed Google Scholar
Sillapachaiyaporn C, Chuchawankul S, Nilkhet S, Moungkote N, Sarachana T, Ung AT, et al. Ergosterol isolated from cloud ear mushroom (Auricularia polytricha) attenuates bisphenol A-induced BV2 microglial cell inflammation. Food Res Int. 2022;157:111433.
Article CAS PubMed Google Scholar
Simerly R, Swanson L, Chang C, Muramatsu M. Distribution of androgen and estrogen receptor mRNA-containing cells in the rat brain: an in situ hybridization study. J Comp Neurol. 1990;294(1):76–95.
Article CAS PubMed Google Scholar
Kight KE, McCarthy MM. Androgens and the developing hippocampus. Biol Sex Differ. 2020;11(1):1–14.
Article Google Scholar
Isgor C, Sengelaub DR. Prenatal gonadal steroids affect adult spatial behavior, CA1 and CA3 pyramidal cell morphology in rats. Horm Behav. 1998;34(2):183–98.
Article CAS PubMed Google Scholar

Download references

Acknowledgements

We wish to acknowledge Asst.Prof.Dr.Suwanakiet Sawangkoon, Assoc.Prof.Dr.Anusak Kijtawornrat, Dr.Nitira Anakkul, Dr.Choopet Nitsakulthong, Ms.Sornsawan Chumjai, Ms.Saifon Sreechomphoe, and Mr.Phanupong Dungkhokkruat of the Chulalongkorn University Laboratory Animal Center for their assistance with the ethical approval process and training ST, SK, KK, PL, PP and TS in the proper care and use of laboratory animals. We are grateful to the Research Affairs, Faculty of Medicine, Chulalongkorn University, for assistance regarding Imaris software. We would also like to acknowledge Asst.Prof.Dr.Tewin Tencomnao, Department of Clinical Chemistry, Faculty of Allied Health Sciences, Chulalongkorn University, for the antibodies used in this study. We wish to acknowledge Assoc.Prof.Dr.Ray-Chang Wu, The George Washington University, USA, for plasmid construction. We thank Dr.Bumpenporn Sanannam, Mr.Thanit Saeliw, and Ms.Kwanjira Songsritaya for their help in animal treatment and tissue collection. We thank Dr.Chanin Sillapachaiyaporn for his help training the molecular docking program.

Funding

This study was supported by a grant from the National Research University Project, Office of Higher Education Commission (NRU59-031-HR), Thailand Science Research and Innovation Fund Chulalongkorn University (HEA663700091), and Ratchadapisek Somphot Fund for Supporting Research Unit, Chulalongkorn University (GRU 6300437001-1, GRU_64_033_37_004, and GRU 6506537004–1) to TS. This research and innovation activity is also funded by the National Research Council of Thailand (NRCT) to ST. ST received a Royal Golden Jubilee Ph.D. Programme Scholarship (grant no. PHD/0029/2561) from the Thailand Research Fund and National Research Council of Thailand and “The 90th Anniversary Chulalongkorn University Fund (Ratchadaphiseksomphot Endowment Fund: GCUGR1125632108D-108)”, Graduate School, Chulalongkorn University. SK is a postdoctoral researcher and is financially supported by the Second Century Fund (C2F), Chulalongkorn University, Bangkok, Thailand. KK, PL, and PP were supported by “The 100th Anniversary Chulalongkorn University Fund for Doctoral Scholarship” and “The 90th Anniversary Chulalongkorn University Fund (Ratchadaphiseksomphot Endowment Fund: GCUGR1125632109D-109, GCUGR1125651062D-062, and GCUGR1125651060D-060)” from the Graduate School, Chulalongkorn University. Animal husbandry and housing were financially supported by a Chulalongkorn University Laboratory Animal Center (CULAC) Grant (Animal Use Protocol No. 2073011) to KK and TS.

Author information

Authors and Affiliations

Department of Clinical Chemistry, Faculty of Allied Health Sciences, Chulalongkorn University, Bangkok, Thailand
Surangrat Thongkorn, Kasidit Kasitipradit, Pattanachat Lertpeerapan & Pawinee Panjabud
SYstems Neuroscience of Autism and PSychiatric Disorders (SYNAPS) Research Unit, Department of Clinical Chemistry, Faculty of Allied Health Sciences, Chulalongkorn University, 154 Soi Chula 12, Rama 1 Road, Wangmai, Pathumwan, Bangkok, 10330, Thailand
Songphon Kanlayaprasit, Depicha Jindatip & Tewarit Sarachana
Department of Biochemistry and Molecular Medicine, The George Washington University School of Medicine and Health Sciences, The George Washington University, Washington, DC, USA
Valerie W. Hu
Department of Anatomy, Faculty of Medicine, Chulalongkorn University, Bangkok, Thailand
Depicha Jindatip

Authors

Surangrat Thongkorn
View author publications
You can also search for this author in PubMed Google Scholar
Songphon Kanlayaprasit
View author publications
You can also search for this author in PubMed Google Scholar
Kasidit Kasitipradit
View author publications
You can also search for this author in PubMed Google Scholar
Pattanachat Lertpeerapan
View author publications
You can also search for this author in PubMed Google Scholar
Pawinee Panjabud
View author publications
You can also search for this author in PubMed Google Scholar
Valerie W. Hu
View author publications
You can also search for this author in PubMed Google Scholar
Depicha Jindatip
View author publications
You can also search for this author in PubMed Google Scholar
Tewarit Sarachana
View author publications
You can also search for this author in PubMed Google Scholar

Contributions

ST performed all the experiments, analyzed the data, and drafted the manuscript under the supervision of TS and DJ. SK, KK, PL, and PP assisted ST with animal treatment and tissue collection. SK performed Western blot analysis. DJ provided the IMARIS software. VWH provided the Ingenuity Pathway Analysis program. TS conceived the study, designed the experiments, analyzed the data, interpreted the results, determined the conclusions, and participated in writing and editing this manuscript. All the authors read and approved the final manuscript.

Corresponding author

Correspondence to Tewarit Sarachana.

Ethics declarations

Ethics approval and consent to participate

All experimental animal procedures were approved by the Chulalongkorn University Animal Care and Use Committee (Animal Use Protocol No. 1673007, No. 1773011, and No. 2073011), Chulalongkorn University.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1

. Sex determination in neonatal rat pups by observing the distance between the external genitalia and anus.

Additional file 2

. Compositions of the media used for primary hippocampal cell culture in this study.

Additional file 3

. List of primers for qRT-PCR analyses.

Additional file 4.

Western blot analysis of AR-overexpressing human neuroblastoma SH-SY5Y cells. The expression of AR protein in stably transfected cells was significantly upregulated compared to that in negative control plasmid transfected cells. The differences between the two groups were analyzed using a two-tailed Student’s t-test. A p-value < 0.05 was considered significant.

Additional file 5

. A list of ASD-related transcription factors.

Additional file 6

. Hypergeometric distribution analysis results between DEGs and transcriptional targets of each TF. A p-value < 0.05 was considered significant. NA = Not applicable

Additional file 7

. Biological functions, disorders, and pathways associated with the transcriptional targets of EGR2 that were dysregulated in the female hippocampus predicted by IPA software. Statistical significance was determined using Fisher’s exact test. A p-value < 0.05 was considered significant.

Additional file 8

. Biological functions, disorders, and pathways associated with the transcriptional targets of KDM5B that were dysregulated in the female hippocampus predicted by IPA software. Statistical significance was determined using Fisher’s exact test. A p-value < 0.05 was considered significant.

Additional file 9

. Biological functions, disorders, and pathways associated with the transcriptional targets of KDM5B that were dysregulated in the male hippocampus predicted by IPA software. Statistical significance was determined using Fisher’s exact test. A p-value < 0.05 was considered significant.

Additional file 10

. Biological functions, disorders, and pathways associated with the transcriptional targets of SMAD4 that were dysregulated in the male hippocampus predicted by IPA software. Statistical significance was determined using Fisher’s exact test. A p-value < 0.05 was considered significant.

Additional file 11

. Biological functions, disorders, and pathways associated with the transcriptional targets of SOX5 that were dysregulated in the female hippocampus predicted by IPA software. Statistical significance was determined using Fisher’s exact test. A p-value < 0.05 was considered significant.

Additional file 12

. Biological functions, disorders, and pathways associated with the transcriptional targets of SOX5 that were dysregulated in the male hippocampus predicted by IPA software. Statistical significance was determined using Fisher’s exact test. A p-value < 0.05 was considered significant.

Additional file 13

. Biological functions, disorders, and pathways associated with the transcriptional targets of TCF7L2 that were dysregulated in both sexes predicted by IPA software. Statistical significance was determined using Fisher’s exact test. A p-value < 0.05 was considered significant

Additional file 14

. Biological functions, disorders, and pathways associated with the transcriptional targets of YY1 that were dysregulated in the male hippocampus predicted by IPA software. Statistical significance was determined using Fisher’s exact test. A p-value < 0.05 was considered significant.

Additional file 15

. ASD-related transcription factor expression data from The Allen Brain Atlas by In Situ Hybridization at several stages of brain development. Expression data of ASD-related transcription factors in developing hippocampus-related areas (rostral secondary prosencephalon; RSP and telencephalic vesicle; Tel). All data processing was performed by R software, and raw expression data were calculated from the sum of expressing pixel intensities divided by the total number of pixels that intersect. NA = Not applicable.

Additional file 16

. ASD-related TFs and their transcriptional targets were selected for qRT-PCR analysis.

Additional file 17

. Expression of estrogen receptors in pAR-overexpressing neuroblastoma SH-SY5Y cell line.

Additional file 18

. Previously published BPA transcriptome studies obtained from NCBI GEO DataSets were used to reanalyze BPA-responsive genes.

Additional file 19

. List of differentially expressed genes in the hippocampus of rodent offspring prenatally/perinatally exposed to BPA from other BPA studies.

Additional file 20

. Hypergeometric distribution analysis results between DEGs in other independent BPA studies and the transcriptional target of ASD-related transcription factors. A p-value < 0.05 was considered significant. NA = Not applicable.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Cite this article

Thongkorn, S., Kanlayaprasit, S., Kasitipradit, K. et al. Investigation of autism-related transcription factors underlying sex differences in the effects of bisphenol A on transcriptome profiles and synaptogenesis in the offspring hippocampus. Biol Sex Differ 14, 8 (2023). https://doi.org/10.1186/s13293-023-00496-w

Download citation

Received: 08 June 2022
Accepted: 07 February 2023
Published: 20 February 2023
DOI: https://doi.org/10.1186/s13293-023-00496-w

Investigation of autism-related transcription factors underlying sex differences in the effects of bisphenol A on transcriptome profiles and synaptogenesis in the offspring hippocampus

Abstract

Background

Methods

Results

Conclusions

Highlights

Introduction

Methods

Collection of transcriptome profile data and the lists of DEGs in response to prenatal BPA exposure

Identification of TFs that are upstream regulators of DEGs

Prediction of biological functions, disorders, and pathways associated with DEGs that are regulated by ASD-related TFs

Expression of ASD-related TFs in the developing hippocampus

In silico virtual screening of ASD-related TFs that are potential BPA targets

Animal husbandry and treatment

RNA isolation

Quantitative RT–PCR analysis

Overexpression of AR in a human neuroblastoma cell line and BPA treatment

Synaptogenesis assay

Western blot analysis

Prediction of TFs involved in the effects of prenatal BPA exposure on the transcriptome profiles of the offspring hippocampus using data from other independent studies

Statistical analyses

Results

Dysregulated genes in the hippocampus of offspring prenatally exposed to BPA were transcriptional targets of ASD-related TFs

BPA-responsive genes that are transcriptional targets of ASD-related TFs are also associated with ASD

Molecular docking analysis of BPA and ASD-related TFs whose targets were differentially expressed in the hippocampus of offspring prenatally exposed to BPA

Effects of prenatal BPA exposure on the expression of ASD-related TFs

Quantitative RT-PCR analysis of transcriptional targets of ASD-related TFs in the hippocampus of male and female offspring

AR is involved in the effects of BPA on the expression of ASD-candidate genes in the human neuroblastoma cell line

Sex differences in the effects of prenatal BPA exposure on synaptogenesis

Reanalysis of transcriptome data from other BPA studies revealed dysregulation of ASD-related TF targets in the hippocampus of offspring

Discussion

Perspectives and significance

Availability of data and materials

References

Acknowledgements

Funding

Author information

Authors and Affiliations

Contributions

Corresponding author

Ethics declarations

Ethics approval and consent to participate

Consent for publication

Competing interests

Additional information

Publisher's Note

Supplementary Information

Rights and permissions

About this article

Cite this article

Share this article

Keywords

Biology of Sex Differences

Contact us