Walter Lichtensteiger1*
Catherine Bassetti-Gaille1
Hubert Rehrauer2
Jelena Kühn Georgijevic2Jesus A.F. Tresguerres3Margret Schlumpf1- 1GREEN Tox and Institute of Veterinary Pharmacology and Toxicology, University of Zurich, Zurich, Switzerland
- 2Functional Genomics Center, Swiss Federal Institute of Technology and University of Zurich, Zurich, Switzerland
- 3Department of Physiology, Universidad Complutense, Madrid, Spain
Endocrine disrupting chemicals (EDCs) can impair hippocampus-dependent behaviors in rat offspring and in children. In search for key processes underlying this effect, we compared the transcriptomes of rat hippocampus on postnatal day 6 after gestational and lactational exposure to three different EDCs at doses known to impair development of learning and memory. Aroclor 1254, a commercial PCB mixture (5 mg/kg or 0.5 mg/kg), or bisphenol A (5 mg/kg or 0.5 mg/kg) were administered in chow, chlorpyrifos (3 mg/kg or 1 mg/kg) was injected subcutaneously. Male hippocampus exhibited a common effect of all three chemicals on genes involved in cell-autonomous processes, Sox6, Sox11, Pou2f2/Oct2, and Pou3f2/Brn2, all upregulated at the high dose. Additional genes of the Sox and Pou families were affected by only one or two of the chemicals. Real time RT PCR showed a comparable expression change for bisphenol A also at the lower dose. Female hippocampus exhibited much fewer genes with expression changes (almost none with false discovery rate <0.05), and none of the genes of the Sox and Pou families was affected. Since gene network analyses in male hippocampus suggested a link between Sox6 and miR-24, known to be repressed by activation of ER-alpha and to repress Sox6 in other tissues, this microRNA was measured. miR-24 was downregulated by all chemicals at the high dose in males. Values of Sox6 mRNA and miR-24 were inversely correlated in individual male hippocampus samples, supporting the hypothesis that the change in Sox6 expression resulted from an action of miR-24. In contrast, miR-24 levels remained unchanged in hippocampus of females. A sexually dimorphic response of miR-24 may thus be at the basis of the sex difference in Sox6 expression changes following exposure to the three chemicals. ER-alpha expression was also sex-dependent, but the expression changes did not parallel those of potential downstream genes such as Sox6. Sox6 is known to suppress differentiation of Parvalbumin (Pvalb)-expressing interneurons. Individual Sox6 levels (FPKM) were inversely correlated with levels of Pvalb, but not with markers of Sox6-independent interneuron subpopulations, Nos1 and 5HT3aR. Effects on interneuron development are further suggested, in males, by expression changes of Nrg1 and its receptor Erbb4, controlling interneuron migration. Our study disclosed new types of EDC-responsive morphogenetic genes, and illustrated the potential relevance of microRNAs in sexually dimorphic EDC actions.
Introduction
The development of cognitive and emotional behaviors represents one of the best documented targets of endocrine disrupting chemicals (EDCs) in humans (Boucher et al., 2009; Rauh et al., 2011; Braun et al., 2017; Ejaredar et al., 2017; Bornehag et al., 2021). Effects on learning and memory and on anxiety have been observed as well in rodent models, for example after perinatal exposure to bisphenol A (Kubo et al., 2001; Ryan and Vandenbergh, 2010; Xu et al., 2010, 2012; Jašarevic et al., 2011; Matsuda et al., 2013; Kumar and Thakur, 2014; Johnson et al., 2016), polychlorinated biphenyls (PCBs) (Gilbert et al., 2000; Roegge et al., 2000; Widholm et al., 2001; Sugawara et al., 2006; Colciago et al., 2009), or phthalates (Zhao et al., 2020). Analogous effects have also been reported for perinatal exposure to the organophosphate pesticide chlorpyrifos (Levin et al., 2001, 2002; Aldridge et al., 2005; Braquenier et al., 2010; Haviland et al., 2010; Turgeman et al., 2011; Chen et al., 2012). Recent data indicate that chlorpyrifos can also interfere with endocrine mechanisms (see below).
From a regulatory perspective, the identification of molecular processes that are common to different types of EDCs with similar adverse effects on brain development, would be of great interest, since knowledge of critical molecular processes could help to design in vitro models and molecular markers for introduction into existing in vivo OECD test guidelines, and to establish links with the human situation. We compared gene expression patterns in developing rat hippocampus after pre- and postnatal exposure to three EDCs at dose levels known to affect learning and memory in rodents, the PCB mixture, Aroclor 1254 (Aro), bisphenol A (BPA), and the pesticide chlorpyrifos (CPF). Impaired memory function was chosen because it represents a well-documented apical adverse endpoint. Impairment of memory function in rats and mice has been observed after perinatal administration of Aroclor 1254 at oral doses between 1 mg/kg and 10 mg/kg, and after BPA at oral doses between 0.05 mg/kg and 50 mg/kg (References above). BPA doses between 0.0025 mg/kg and 0.02 mg/kg were ineffective (Fujimoto et al., 2006; Ferguson et al., 2012; Nakamura et al., 2012). Chlorpyrifos was effective at subcutaneous doses between 0.1 mg/kg (one report), 1 mg/kg and 5 mg/kg (see Supplemental Table S1). BPA is estrogenic but also interferes with thyroid hormone actions (Rubin, 2011). Aro interacts with thyroid hormones in developing brain (Gauger et al., 2004), and contains estrogenic congeners (Decastro et al., 2006). The acetylcholinesterase inhibitor CPF exerts non-endocrine effects on nerve cells, but has also been found to exhibit estrogenic activity, and to affect thyroid hormone levels in brain (Andersen et al., 2002; Slotkin and Seidler, 2009; Ventura et al., 2012; Slotkin et al., 2013; Hazarika et al., 2020).
Hippocampus is a hormone target during development and postnatal life in animals as well as humans (Ivanova and Beyer, 2000; González et al., 2007; Mogi et al., 2015; Gothié et al., 2017). A special feature of developing hippocampus is that there are two possible sources of sex steroids, the general circulation and local de novo synthesis from progesterone (Hojo et al., 2004; Prange-Kiel et al., 2009). The relative role of the two sources has not been clarified, but there is evidence for hippocampal estrogen synthesis during the neonatal period (Amateau et al., 2004; Kretz et al., 2004; Nuñez and McCarthy, 2009). The investigation focused on postnatal day (PND) 6, a stage allowing to investigate a number of different developmental processes, including cell proliferation, differentiation (in particular interneurons), and synaptogenesis (Danglot et al., 2006; Bowers et al., 2010). The early postnatal stage is covered by several OECD test guidelines addressing reproductive and neurodevelopmental toxicity, hence, molecular markers with predictive value could be included into these guidelines.
Transcriptomics analysis revealed analogous effects of the three chemicals on a small number of gene ontology processes related to early neurodevelopment, with sexually dimorphic expression changes of genes involved in cell-autonomous control of development, in particular Sox6, and of downstream genes of Sox6 linked to interneuron development. Sexually dimorphic effects of all three chemicals on an estrogen-regulated microRNA, miR-24, may be involved.
Methods
Chemicals
Aroclor 1254 (Lot 124-191-B) was obtained from CromLab S.L. (Barcelona, Spain), bisphenol A (2,2-bis-(4-hydroxyphenyl)-propane, 99%) from Sigma-Aldrich Chemie (Buchs, Switzerland), and chlorpyrifos (chlorpyrifos PESTANAL, analytical standard) from Sigma-Aldrich Quimica (Tres Cantos. Madrid, Spain).
Animals and Treatment
Animals
Long Evans rats were purchased from Janvier Labs (Le Genest Saint Isle Saint Berthevin, France). Experiments were conducted at the Dept. of Endocrinology, Complutense University, Madrid, and approved by the ethical committee of Complutense University. Females (9–10 weeks old) were housed at a 12 h light-dark cycle (lights on 07.00-19.00), 22°C, in groups of 4 during premating, two females with one male (12–14 weeks old) during the 1 week mating period, and in groups of two during pregnancy, being separated on gestational day (GD) 22, 1 day before parturition (GD 1 = day of sperm-positive vaginal smear, GD 23 = postnatal day (PND) 1 = day of birth).
Treatment
Aroclor 1254 (Aro) and bisphenol A (BPA) were administered in the chow at two dose levels, 5.0 mg/kg day (Aro5, BPA5) and 0.5 mg/kg day (Aro0.5, BPA0.5), to female F0 rats 2 weeks before mating, during mating, pregnancy and lactation, until PND 6, with water ad libitum. Both chemicals were admixed to phytoestrogen-low (soy-free) rat chow (IRTA) by the Dept. of Endocrinology (Madrid) and further diluted to the final concentrations, 50 mg/kg chow and 5 mg/kg chow, by IRTA (Institute of Agrifood Research and Technology, Torre Marimon, Barcelona), assuming a daily food intake of 10% of body weight (Historical data of female Long Evans rats in the laboratories of Zurich and Madrid: 8.2 and 10% of body weight, respectively). Feeding controls (CONF) received soy-free chow.
Chlorpyrifos (CPF) was dissolved in dimethyl sulfoxide (DMSO) and injected subcutaneously (1 ml DMSO/kg body weight), at 1 mg/kg (CPF1) and 3 mg/kg body weight per day (CPF3), to female F0 rats between GD 2 and GD 20 and again between PND2 and PND 6, in order to avoid possible interactions with parturition. Injection controls (CONSC) received 1 ml DMSO/kg body weight per day. The s. c. route was chosen because the majority of developmental studies with analysis of learning and memory, on which dose selection was based, had used s. c. injection (see Supplementary Table S1).
Dose selection
The higher dose level has consistently been found to impair learning/memory in adult rats and mice after pre- and/or early postnatal administration (Supplementary Table S1). In the case of BPA and CPF, a few studies on anxiety were also included in view of the possible role of hippocampus in anxiety. The lower dose of BPA and CPF was effective in some studies, whereas the lower dose of Aro had not been reported to affect learning or anxiety.
Tissue Preparation and Dissection of Hippocampus
Tissue preparation
Offspring were weighed and decapitated at postnatal day (PND) 6. Brain was excised and deep-frozen on dry ice by staff of GREEN Tox assisted by staff of Madrid. Sex was determined by gonadal inspection. Tissues were stored at −80°C before being transferred to Zurich on dry ice and again stored at −80°C. Time-window: Brains used for transcriptomics were dissected between 10.30 and 13.20 h. Real time RT PCR was also performed on brain tissue from three Aro5-exposed male pups dissected between 15.45 and 17.50 h. Levels of selected mRNAs analyzed in the project did not differ from those of pups dissected within the time-window.
Brain dissection
Dissection of left plus right dorsal hippocampus from deep-frozen brain was performed in Zurich. A frontal 1.0 mm thick slice was cut at a level immediately posterior to the optic chiasm in a stainless steel rat brain matrix (Zivic Instruments) cooled to below −10°C by dry ice. From this, a rectangular tissue piece comprising hippocampus plus laterally adjacent cortical tissue was dissected under microscopic control on a Pelletier element (−10°C) (Supplementary Figure S1), transferred into an Eppendorf tube containing 200 μl RNAlater (QIAGEN), and stored at −20°C. Brain dissections were performed in a cold room (4°C). This had been found to result in RIN (RNA integrity number) values for RNA of 8–10, as compared to RIN 5–8 with dissections on cooled plates at room temperature. All dissections were done by the same person (MS), without knowledge of the treatment group.
Transcriptomics
Tissue Homogenization and RNA Extraction
Tissues were transferred to Eppendorf tubes containing RLT buffer, 10 µl ß-mercaptoethanol/1 ml buffer, homogenized, and stored at -80°C. Total RNA was extracted with QIAGEN RNeasy-mini kit, treated with DNase-I (QIAGEN). RIN values ranged between 7.1 and 9.3 in males (Mean ± SD (N): 8.7 ± 0.5 (42)) and between 8.4 and 9.3 in females (9.0 ± 0.2 (24)). Samples with A260/A280 ratio <1.8 were further purified by NucleoSpin® RNA Clean-up XS kit (Macherey-Nagel). RNA samples were stored at −80°C (Males: CONF N = 5, CONSC N = 5, Aro0.5 N = 4, Aro5 N = 6, BPA0.5 N = 4, BPA5 N = 7, CPF1 N = 4, CPF3 N = 7. Females: CONF N = 6, CONSC N = 4, Aro5 N = 5, BPA5 N = 5, CPF3 N = 4).
RNA-Seq (Illumina RNA Sequencing)
RNA-seq was performed in the Functional Genomics Center (FGCZ) of the ETH Zurich and University of Zurich (J. K. G.) on extracts from dorsal hippocampus of rat pups euthanized between 10.30 and 13.20 h.
Library preparation
The quantity and quality of the isolated RNA was determined with a Qubit® (1.0) Fluorometer (Life Technologies, California, United States) and a Bioanalyzer 2100 (Agilent, Waldbronn, Germany). The TruSeq Stranded mRNA Sample Prep Kit (Illumina, Inc, California, United States) was used in the succeeding steps. Briefly, total RNA samples (100 ng) were poly-A selected and then reverse-transcribed into double-stranded cDNA with Actinomycin added during first-strand synthesis. The cDNA samples was fragmented, end-repaired and adenylated before ligation of TruSeq adapters. The adapters contain the index for multiplexing. Fragments containing TruSeq adapters on both ends were selectively enriched with PCR. The quality and quantity of the enriched libraries were validated using Qubit® (1.0) Fluorometer and the Bioanalyzer 2100 (Agilent, Waldbronn, Germany). The product is a smear with an average fragment size of approximately 360 bp. The libraries were normalized to 10 nM in Tris-Cl 10 mM, pH8.5 with 0.1% Tween 20.
Cluster Generation and Sequencing
The TruSeq SR Cluster Kit v4-cBot-HS (Illumina, Inc, California, United States) was used for cluster generation using 8 pM of pooled normalized libraries on the cBOT. Sequencing were performed on the Illumina HiSeq 2500 single end 125 bp using the TruSeq SBS Kit v4-HS (Illumina, Inc, California, United States). Adapter sequences shown in Supplementary Table S2.
Analysis of RNA-Seq Data
The first part of RNA-seq data analysis was performed by H. R. at the Functional Genomics Center (FGCZ). RNA-seq reads were aligned with the STAR-aligner (Dobin et al., 2013; Love et al., 2014). As reference we used the NCBI genome build for rat Rnor_6.0. Gene expression values were computed with the function featureCounts from the R package Rsubread (Liao et al., 2013). Differential expression was computed using the generalized linear model implemented in the Bioconductor package DESeq2.
Differential expression was computed using ANOVA, followed by Tukey HSD post-hoc test. The two control groups, CONF and CONSC, did not differ significantly from each other. This was further confirmed in a comparison of CONSC and CONF of six target genes by real time RT PCR (Supplementary Figure S2). Thus, a possible difference in the condition of offspring from injection and feeding control groups suggested by body weight, was not reflected in target gene expression in hippocampus. Hence, CONF and CONSC were combined.
Deseq2 data files were further analyzed by enrichment analysis for “significantly affected genes” (different from controls), “differentially affected genes” (different expression in different treatment groups) and “similarly affected genes” (analogous, statistically significant expression changes in the three treatment groups), using MetaCore™ (Thomson Reuters). A group of gene ontology (GO) processes that were similarly affected by the treatments, was further investigated. Extensive gene network analyses (>100) of genes from similarly affected GO processes revealed gene clusters with relevance for neurodevelopment, affected by all three or by two treatments, and a link to a microRNA, miR-24. Literature references documenting links between two genes (listed by MetaCore) were checked in detail in order to assess the reliability of the link. Possible relationships between two genes were further examined by correlation analyses of FPKM (fragments per kilobase per million reads) values of individual hippocampus samples.
Analysis of Individual mRNAs by Real Time RT PCR
mRNAs of genes of interest identified in network analyses of the transcriptomics data were quantitated by real time RT PCR on extracts from male hippocampus of both dose groups. The males had been euthanized between 10.15 and 13.45, with the exception of three male pups treated with Aroclor 1254, 5 mg/kg, dissected at 15.45, 17.07 and 17.50. In these three pups, the expression levels of the mRNAs studied did not differ from those of the remaining pups. For statistical reasons (see below), two groups of samples were analyzed, samples where transcriptomics had also been done (samples with transcriptomics), and samples where transcriptomics had not been done (samples without transcriptomics). Real time RT PCR was not performed in females, since transcriptomics analysis did not show an expression change in genes of interest.
Tissue Homogenization and RNA Extraction
Tissues were transferred to Eppendorf tubes containing Lysis Buffer RLT buffer, 10 µl ß-mercaptoethanol/1 ml buffer, homogenized, and homogenates stored at -80°C. Total RNA was extracted with QIAGEN RNeasy-mini kit (Qiagen), treated with DNase-I (QIAGEN). RNA Samples with A260/A280 ratio >1.8 were stored at -80°C.
Real Time RT PCR
Real Time RT PCR was performed as previously described (Lichtensteiger et al., 2015). mRNA of hippocampus samples from all eight experimental groups of one sex was analyzed in duplicate and the RNA standard curve in triplicate on 384 wells plates by TaqMan quantitative real-time polymerase chain reaction (CFX384 real-time system, Bio-Rad Laboratories Srl, Milan, Italy). Primers and TaqMan probes were designed with PrimerExpress Software, version: 3.0 (Applied Biosystems), using NCBI Gene Bank and Ensemble database (http://www.ncbi.nlm.nih.gov, http://www.ensembl.org), and synthesized by Microsynth (Balgach, Switzerland) (List see Supplementary Table S3). Since investigations in a foregoing study (Lichtensteiger et al., 2015) had shown that expression of 18 possible reference genes was affected by treatment with endocrine active chemicals, quantitation of target mRNAs was based on the standard curve method. Total RNA, containing 80% ribosomal RNA, was considered to represent a stable reference, because expression of two of the four ribosomal RNAs (Rn18, Rn28) had previously been found to be unchanged (Lichtensteiger et al., 2015).
Analysis of Real Time RT PCR Data
mRNA levels of individual samples were expressed as percentage of the mean of the corresponding control group. One-way analysis of variance (ANOVA) was performed for each mRNA species, followed by Dunnett’s multiple comparisons test (GraphPad Prism 7). Three ANOVA analyses were conducted for each mRNA species: Samples also used in transcriptomics analysis, samples not used in transcriptomics analysis, and both groups combined (Males total N: CONSC = 8, CONF = 10, Aro0.5 = 11, Aro5 = 10, BPA0.5 = 10, BPA5 = 12, CPF1 = 11, CPF3 = 11). The group of samples without transcriptomics analysis represents a statistical sample which is independent of transcriptomics analysis and hence, allows to control the transcriptomics data. This is particularly important for genes with a false discovery rate in transcriptomics of Fdr >0.05. In the present investigation, this was only the case for Pou2f1/Oct1 with Fdr = 2.14E-01 for bisphenol A, and Fdr = 1.12E-01 for Aroclor 1254. Fdr values of all other genes of interest ranged between 2E-02 and 9E-07. Figures show results with samples not used in transcriptomics analysis (1 per litter).
Determination of microRNA-24 (miR-24-3p)
Tissue Homogenization and microRNA Extraction
miR-24-3p analysis was limited to the high-dose groups because there was not enough tissue material left for analysis of the low-dose groups. Frozen tissues were homogenized in Lysis Buffer RLT buffer (Qiagen) and miRNAs were extracted using mirVana miRNA Isolation Kit (Ambion) following the manufacturer’s protocol. The majority of extracts was obtained from hippocampus samples also extracted for real time RT PCR.
Determination of miR-24-3p
miR-24-3p originates from two gene clusters, miR 23b, miR 27b, miR-24-1 (rat chromosome 17), and miR-23a, miR-27a, miR-24-2 (rat chromosome 19), but the mature sequence of miR-24-3p-1 and miR-24-3p-2 is the same, miR-24-3p (Chhabra et al., 2010; Marchetti et al., 2020). miR-24-3p is documented for a number of species. Rat miRBase ID is rno-miR-24-3p, miRBase Accession Number MIMAT0000794, the mature miRNA Sequence UGGCUCAGUUCAGCAGGAACAG.
Poly(A) tailing, ligation, reverse transcription and amplification reactions were performed with the TaqMan Advanced miRNA cDNA Synthesis Kit (Applied Biosystems/ThermoFisher Scientific), using 10 ng total RNA input. qRT-PCR experiments were performed in duplicate on 384-well plates with TaqMan Fast Advanced Master Mix (Applied Biosystems/ThermoFisher Scientific) and TaqMan primer/probe (TaqMan Advanced miRNA Assay, Assay ID rno481011_mir, Applied Biosystems/ThermoFisher Scientific). PCR amplifications were performed on CFX384 real-time system (Bio-Rad Laboratories Srl, Milan, Italy). Quantitation of target mRNAs was based on the standard curve method (see above).
Analysis of miR-24-3p Data
microRNA levels of individual samples were expressed as percentage of the mean of the corresponding control group. One-way analysis of variance (ANOVA) was performed for each microRNA, followed by Dunnett’s multiple comparisons test (GraphPad Prism 7).
Results
General Developmental Endpoints
Male and female offspring born to control dams receiving subcutaneous injections of vehicle showed lower body weights than offspring of the feeding control (Table 1). Yet, body weights as well as total number of pups exposed to the three chemicals were at the level of the corresponding, feeding or injection controls, indicating that the chemicals did not affect body weight or total number of pups. There also were no significant differences in sex ratio.
TABLE 1. Body weight and number of rat offspring at postnatal day 6 (litter mean).
Gene Expression in Hippocampus of Male and Female Rat Offspring on Postnatal Day 6
Analyses were performed on dorsal hippocampus (plus a small part of laterally adjacent neocortical tissue, Supplementary Figure S1) on postnatal day (PND) 6, when pups were still exposed to the chemicals. In the high-dose groups, the dose range consistently linked with alterations of hippocampus-dependent behaviors (Aroclor 1254, 5 mg/kg, bisphenol A, 5 mg/kg, chlorpyrifos, 3 mg/kg), transcriptomics analysis in hippocampus of male offspring revealed significant changes in gene expression (Table 2). ANOVA of hippocampus of male offspring exposed to a 10 times lower dose (Aroclor 1254, 0.5 mg/kg, bisphenol A, 0.5 mg/kg) or to a three times lower dose (chlorpyrifos, 1 mg/kg) yielded an uncertain result, with high false discovery rates (Fdr) in the great majority of genes (Table 2). Real time RT PCR, where a larger number of samples was studied, showed signficant changes in some target genes also in the low-dose groups (see below). Transcriptomics analysis of high-dose groups in females showed a much lower number of significantly affected genes, and the great majority of genes showed a high false discovery rate, suggesting high interindividual variability (Table 2, Table 4). Methodological aspects, RNA quality and RNA-seq conditions were found to be very good, at the level of the analysis of male samples. The RNA-seq assay of female hippocampus was repeated, with the same result. Since the number of samples in the various experimental groups was similar to that of males (see Methods), expression changes of a magnitude comparable to the high-dose groups of males should have been detectable.
TABLE 2. Number of genes with significant change of expression rate in transcriptomics analyses of rat hippocampus at postnatal day 6, according to p value and false discovery rate (Fdr).
The developmental switch of GABAergic transmission from excitation to inhibition, characterized by an increase in the expression ratio of the Cl-extruder KCC2, a neuron-specific potassium-chloride-co-transporter, relative to NKCC1, the sodium-potassium-chloride-co-transporter increasing intracellular Cl-concentration, takes place in hippocampus at the end of the first postnatal week (Watanabe and Fukuda, 2015). In line with previous findings (Nuñez and McCarthy, 2009), the KCC2/NKCC1 ratio (expressed as fragments per kilobase per million reads) was higher in female hippocampus (Figure 1, p < 0.0001), however, this general developmental trait was not significantly influenced by the chemicals.
FIGURE 1. Epression ratio of potassium-chloride-co-transporter KCC2 (Slc12a5), neuron-specific chloride ion extruder, to sodium-potassium-chloride-co-transporter NKCC1 (Slc12a2), increasing intracellular chloride concentration, in individual hippocampus samples of male and female rat offspring at postnatal day 6. RNAseq data, ratio of FPKM (fragments per kilobase per million reads), mean ± SD and individual samples (M, male, F, female). Two way ANOVA, sex difference (p < 0001), but no treatment effects and no sex × treatment interaction.
Gen Ontology Processes and Gene Networks Related to Neurodevelopment
Enrichment analysis in male hippocampus (significance level p < 0.01, effect threshold Thr (Log2 ratio) = 0.5 ≅ 141% of control level) focused on male offspring exposed to the higher dose levels, since ANOVA was essentially negative for the lower dose groups. It revealed a number of significantly affected gene ontology (GO) processes related to brain development and function (Supplementary Table S4). Since all three treatments impair development of learning and memory in rodents, we reasoned that it should be possible to identify GO processes in hippocampus that were similarly affected by all three chemicals, and similarities in gene expression patterns. An enrichment analysis for similarity of changes (p < 0.01, Thr = 0.5, Supplementary Figure S3) identified several GO processes related to development of hippocampus and cerebral cortex, with −log (p value) between 3 and 5 for all three treatments (Table 3, Supplementary Table S5 with all significantly affected genes). Interestingly, these GO processes control early developmental processes: Neuroepithelial cell differentiation, regulation of stem cell proliferation, regulation of neuroblast proliferation, positive regulation of neuroblast proliferation, positive regulation of neural precursor cell proliferation, cerebral cortex radially oriented cell migration, gliogenesis, oligodendrocyte development. In contrast, enrichment analysis for similarity of processes in female hippocampus (p < 0.01, Thr (Log2 ratio) = 0.1) did not show a clear-cut focus on GO processes related to neurodevelopment. Only two very general terms, “generation of neurons” (38/100) and “nervous system development” (40/100) reached p values below 1E-05 (Supplementary Figure S4). Moreover, there was little overlap of significantly affected genes.
TABLE 3. GO processes similarly affected by all three treatments in male hippocampus, with genes selected for network analysesa,b.
The group of similarly affected GO processes in male hippocampus shows several genes involved in cell-autonomous regulation of developmental processes, Gli6, Sox6, Sox11, and Pou3f2 (Brn2), the non-canonical Wnt receptor, Fdz3, and the transcription factor Tcf7L2 (Table 3). Possible links with endocrine actions were analyzed for these genes as well as for Nrg1 and its receptor Erbb4, previously found to be targeted by EDC mixtures (Lichtensteiger et al., 2015). Gene networks (n = 154) with information on transcriptional regulation, receptor binding, and phosphorylation, were constructed in MetaCore for each of these genes in combination with sex steroid receptors (Esr1, Esr2, AR), thyroid hormone receptors (TRa, TRb1), arylhydrocarbon receptor (AhR), and components of signaling pathways, with significance level p < 0.05 and Thr (Log ratio) = 0.5 (sometimes 0.3). Attention was given to possible links with miRs.
Sox and Pou Domain Genes Involved in Cell-Autonomous Regulation of Development
Network analyses based on the aforementioned genes revealed a possible involvement of a number of genes of the Sox and Pou families in male hippocampus (Figure 2). Sox6, Sox11, Pou2f2/Oct2, and Pou3f2/Brn2 stand out as being similarly affected by all three chemicals, besides Gli3 (Table 4). Expression of Sox5, which interacts with Sox6 during certain developmental phases (Azim et al., 2009; Batista-Brito et al., 2009), was unchanged. Some other Sox and Pou genes were significantly affected by only one or two chemicals. Expression of these genes in female hippocampus was at control level (Supplementary Table S6). In an additional set of hippocampus samples of male offspring from the same litters but not used for transcriptomics, levels of Sox6, Sox11, and Pou2f1/Oct1 mRNA were determined by real time RT PCR (Figure 3). Pou2f1/Oct1 was chosen because transcriptomics results were uncertain with high Fdr values for Aroclor 1254 and bisphenol A (5 mgkg). This analysis confirmed the effects of all three chemicals at the high dose level on Sox6, Sox11, and of bisphenol A and chlorpyrifos on Pou2f1/Oct1. Real time RT PCR revealed significant effects of bisphenol A on all three genes also at the lower dose level (Figure 3), possibly because this analysis was based on a larger sample size than transcriptomics. The magnitude of expression changes detected by RNA-seq and real time RT PCR was very similar (Supplementary Table S7).
FIGURE 2. Gene network of male rat hippocampus at postnatal day (PND) 6, with estrogen receptor alpha (Esr1). Up-regulation (red circles) or down-regulation (blue circle) of hippocampal genes exposed to Aroclor 1254 (5 mg/kg, A), bisphenol A (5 mg/kg, B), or chlorpyrifos (3 mg/kg, C) (see Table 4). Pou2f1/Oct1 and Sox11 exhibited a false discovery rate >0.05 after bisphenol A, but the effect was confirmed by real time RT PCR (Figure 3). Literature data indicate a possible inhibitory action of Esr1 on transcription of miR-24-3p, and of miR-24-3p on transcription of Sox6 (MetaCore, Shortest Paths network. TR: transcription, B: binding, M: miR binding, green: activation, red: inhibition).
TABLE 4. Male hippocampus at postnatal day 6: Expression of similarly affected genes (transcriptomics)a.
FIGURE 3. mRNA levels in male hippocampus at postnatal day 6 after pre- and postnatal administration of Aroclor 1254 (ARO0.5: 0.5 mg/kg, ARO5: 5 mg/kg), bisphenol A (BPA0.5: 0.5 mg/kg, BPA5: 5 mg/kg), or chlorpyrifos (CPF1: 1 mg/kg, CPF3: 3 mg/kg). Quantitative RT PCR of hippocampus samples not used for transcriptomics. Mean ± SEM, number of samples (one/litter). Difference from control: *p < 0.05, **p < 0.01, ***p < 0.001.
Genes Encoding for Growth Factors and Receptors of Transcellular Signaling Pathways
In male hippocampus, neuregulin 1 (Nrg1), which controls migration and other aspects of development of cortical interneurons (Mei and Nave, 2014), was found to be up-regulated in transcriptomics by all three treatments and in real time RT PCR by the higher dose levels of Aroclor 1254 and bisphenol A (Table 4; Figure 3). Its receptor Erbb4, expressed in hippocampal interneurons (Mei and Nave, 2014), was up-regulated by the higher doses of Aroclor 1254 and chlorpyrifos, as shown by both methods of analysis. In female hippocampus, a slight reduction of Erbb4 mRNA to 88.2% CON by 5 mg/kg of Aroclor 1254 is uncertain because of an increased Fdr. Nrg1 expression was unchanged (Supplementary Table S6). Interestingly, transcriptomics did not reveal expression changes of a number of growth factors and their receptors in male hippocampus, including BDNF, GDNF, EGF, IGF, NGF, and Vegfa. Fgf14 was upregulated by all three treatments and Fgf12 by Aroclor 1254 (data not shown). The GO analyses also showed significant effects on Fzd3, involved in non-canonical Wnt signaling, and Notch2 (Table 3), but since no additional changes in these pathways or related genes were observed in network analyses, they were not pursued further. Transcriptomics of female hippocampus did not reveal significant expression changes of the aforementioned genes.
miR-24 and Sexually Dimorphic Regulation of Sox6
It is not known what might be the basis of the sex difference in the responsiveness of the Sox and Pou genes to the three chemicals. In the case of Sox6, our gene network analyses suggested a possible link between ER-alpha (Esr1), miR-24, and Sox6 (Figure 2). ER-alpha-mediated repression of miR-24 transcription, and repression of Sox6 transcripton by miR-24, has been demonstrated in non-neural human cell lines and mouse pancreatic beta cells (Maillot et al., 2009; Melkman-Zehavi et al., 2011; Ferraro et al., 2012; Klinge, 2012). miR-24 has previously been detected in rat hippocampus (Liu et al., 2016). It is derived from two miR clusters which may be differentially regulated, but the mature miR-24-3p sequence derived from the two clusters is the same (Chhabra et al., 2010). miR-24-3p was significantly down-regulated in hippocampus of male offspring by all three chemicals administered at the higher dose level, whereas no significant expression change was observed in hippocampus of females (Figure 4). A significant inverse relationship existed between miR-24-3p levels and Sox6 mRNA levels (real time RT PCR data) in individual hippocampus samples of male offspring (Figure 4).
FIGURE 4. A. Level of miR-24-3p in hippocampus of male and female rat offspring at postnatal day 6 after pre- and postnatal administration of Aroclor 1254 (5 mg/kg, Aro5), bisphenol A (5 mg/kg, BPA5), or chlorpyrifos (3 mg/kg, CPF3). Quantitative RT PCR, individual values with mean ± S.D, *p < 0.05 different from control group. CON: controls, ARO5: Aroclor 1254, 5 mg/kg, BPA5: bisphenol A, 5 mg/kg, CPF3: chlorpyrifos 3 mg/kg. B. Relationship between miR-24-3p (% of control) and Sox6 mRNA (% of control, real time RT PCR) in individual hippocampus samples of male offspring. Inverse correlation, p = 0.0097, r = −0.4645, slope = −0.264. Green symbols: control, brown: Aroclor 1254 5 mg/kg, blue: bisphenol A 5 mg/kg, red: chlorpyrifos 3 mgkg.
The transcriptomics data set indicated increased ER-alpha expression in males but not females after exposure to the higher dose of all three chemicals, and no significant change in the lower dose groups. ER-beta (Esr2) levels remained unchanged (Supplementary Table S8). An analysis of a larger number of hippocampus samples by real time RT PCR confirmed the increased expression of ER-alpha mRNA in male offspring exposed to the chemicals at the higher doses, but revealed a decrease of expression after exposure to the lower doses (Figure 5A). ER-alpha mRNA levels in female hippocampus remained unchanged, except for a significant increase after 0.5 mg/kg BPA. Thus, the two sexes showed opposite responses to BPA. Control levels of ER-alpha mRNA were significantly higher in male hippocampus (131.8% of female level, Figure 5B), in line with previous reports of this developmental phase (Ivanova and Beyer, 2000). The ratio of ER-alpha to ER-beta expression (FPKM values of transcriptomics) differed between sexes (p = 0.0002) (Figure 5C). It was increased in hippocampus of males by all three treatments at the higher dose level (p = 0.0060), but unchanged in females, resulting in a significant sex × treatment interaction (p = 0.0113). ER-alpha/ER-beta ratios of controls were similar in both sexes. No effect was seen in hippocampus of males in the lower dose-range (Supplementary Figure S5).
FIGURE 5. (A) Estrogen receptor-alpha (ER-alpha, Esr1) in hippocampus of male and female offspring at postnatal day 6. Quantitative RT PCR, mean ± S.E.M, number of samples. ANOVA, ****p < 0.0001, ***p < 0.001, *p < 0.05. (B) Esr1 mRNA level is higher (131.8%) in male than female hippocampus. Quantitative RT PCR, mean ± S.E.M, number of samples. Unpaired t test, p < 0.05. (C) Ratio of ER-alpha (Esr1)/ER-beta (Esr2) in hippocampus of individual male and female offspring is increased in males by treatment with Aro, BPA and CPF at the higher dose level. RNAseq data, FPKM (fragments per kilobase per million reads), two-way ANOVA, treatment p = 0.0060, sex p = 0.0002, interaction treatment x sex p = 0.0113. CON: controls, ARO0.5, ARO5: Aroclor 1254 0.5 mg/kg and 5 mg/kg, BPA0.5, BPA5: bisphenol A 0.5 mg/kg and 5 mg/kg, CPF1, CPF3: chlorpyrifos 1 mg/kg and 3 mg/kg.
Target Genes of Sox6 and Sexually Dimorphic Responses in Interneurons
The developmental stage studied here corresponds to a period, when Sox6 is involved in differentiation of postmitotic interneuron subtypes, specifically, of medial ganglionic eminence (MGE)-derived parvalbumin (Pvalb)- and somatostatin (Sst)-expressing interneurons (Tricoire et al., 2011; Kelsom and Lu, 2013; Hu et al., 2017). FPKM values of individual hippocampus samples from male offspring of controls and high-dose treatment groups revealed a highly significant inverse correlation between Pvalb and Sox6, whereas no such correlation was seen between individual values of Sox6 and neuronal nitric oxide synthase (Nos1) or between Sox6 and the serotonergic ionotropic receptor Htr3a (Figure 6). The Pvalb/Nos1 expression ratio (FPKM values) was significantly reduced by Aroclor 1254 and chlorpyrifos, and showed a similar tendency with bisphenol A. Nos1 is expressed by some MGE-derived hippocampal interneurons, but is unaffected by conditional removal of Sox6 (Jaglin et al., 2012). Htr3a is a marker of CGE-derived, Sox6-independent interneurons (Vucurovic et al., 2010; Touzot et al., 2016). Transcriptomics analysis indicated a significant down-regulation of Pvalb by the higher doses of Aroclor 1254 and chlorpyrifos (Table 5), while real time RT PCR analysis of a different set of samples showed significantly reduced Pvalb mRNA levels after the higher doses of Aroclor 1254 and bisphenol A (Figure 3). According to real time RT PCR analysis, the lower dose of the three treatments may also have been active, but the effect did not reach statistical significance (Figure 3). Nos1 was slightly but significantly up-regulated by Aroclor 1254 (Table 5), possibly also by chlorpyrifos (Fdr >0.1). Sst expression was not significantly changed (Table 5).
FIGURE 6. Relationship between Sox6 and interneuron markers expressed in individual hippocampus samples of male offspring at postnatal day 6 (RNAseq data, FPKM values). (A) Inverse correlation between expression of Parvalbumin (Pvalb) expressed in medial ganglionic eminence (MGE)-derived interneurons and Sox6, p = 0.0039, r = −0.6155, slope = −6.065. No correlation between Sox6 and levels of neuronal nitric oxide synthase (Nos1) expressed in subgroup of MGE-derived interneurons, or serotonergic ionotropic receptor Htr3a expressed in caudal ganglionic eminence-derived interneurons, both independent of Sox6. (B) Correlation between Pvalb and Sox6 with identification of treatment groups. C. Expression ratio of Pvalb/Nos1, mean ± S.E.M, N = 5-6.
TABLE 5. Genes involved in development of different types of GABAergic hippocampal interneurons in male offspring (transcriptomics)a.
Sox6-independent NPY, expressed in MGE-, CGE-, and preoptic area-derived hippocampal interneurons (Wonders and Anderson, 2006; Gelman et al., 2009; Kelsom and Lu, 2013), as well as markers of CGE-dependent interneurons (Tricoire et al., 2011; Kelsom and Lu, 2013; Touzot et al., 2016), which are Sox6-independent, did not exhibit significant expression changes in male hippocampus (Table 5; p < 0.05, Fdr <0.05)). In hippocampus of female offspring, the aforementioned interneuron markers remained unchanged (Supplementary Table S6).
Discussion
The PCB mixture Aroclor 1254, bisphenol A, and chlorpyrifos, which exhibit endocrine disrupting properties and partly overlap but also differ in their mechanisms of action and effect patterns, all interfere with development of learning and memory, when administered to rodents during pre- and/or neonatal life (Roegge et al., 2000; Levin et al., 2001, 2002; Widholm et al., 2001; Aldridge et al., 2005; Sugawara et al., 2006; Miyagawa et al., 2007; Colciago et al., 2009; Johnson et al., 2009; Xu et al., 2010; Jašarevic et al., 2011; Turgeman et al., 2011; Chen et al., 2012; Kumar and Thakur, 2014; Johnson et al., 2016). We asked the question whether the three chemicals might share common gene expression patterns and signaling pathways that might be linked with development of hippocampus, which is essential for memory functions. Transcriptomics in hippocampus of postnatal day (PND) 6 revealed significant similar effects of the three chemicals on GO processes related to early neurodevelopment, but only in male offspring. Of particular interest are similar effects of all three chemicals on several Sox and Pou genes involved in cell-autonomous control of development, Sox6, Sox11, Pou2f2/Oct2, and Pou3f2/Brn2. This class of genes seems to represent a new focus with respect to neurodevelopmental effects of EDCs. We are only aware of a report by Gauger et al. (2004), who observed upregulation of Oct1 (Pou2f1) in cerebral cortex of rat fetuses exposed to Aroclor 1254. In PND 6 hippocampus, Oct1 was not affected by Aroclor 1254, but was upregulated by bisphenol A and chlorpyrifos. The fact that all three chemicals affected expression of the four genes of the Sox and Pou families, and in the same manner, suggests that the effect may be relevant to alterations in hippocampus development leading to impaired memory functions. There also were GO processes and genes in developing hippocampus exhibiting differential responses to the three chemicals (Supplementary Table S4).
The higher dose of the three chemicals, leading to upregulation of Sox6, Sox11, Pou2f2/Oct2, and Pou3f2/Brn2 mRNAs, corresponds to a level impairing learning and memory in rodent offspring (Roegge et al., 2000; Levin et al., 2001, 2002; Widholm et al., 2001; Aldridge et al., 2005; Sugawara et al., 2006; Miyagawa et al., 2007; Johnson et al., 2009; Xu et al., 2010; Jašarevic et al., 2011; Turgeman et al., 2011; Chen et al., 2012; Kumar and Thakur, 2014; Johnson et al., 2016). Transciptomics analysis of the low-dose groups did not yield significant results, possibly because of high inter-individual variability (false discovery rates above 0.05). Yet, real time RT PCR analysis of Sox6, Sox11, and Pou2f1/Oct1 showed an analogous upregulation of mRNA levels by bisphenol A (but not Aroclor 1254 or chlorpyrifos) also at the lower dose level, which is about 10 times higher than the lowest levels consistently found to affect memory development in rodents (Nakamura et al., 2012; Kumar and Thakur, 2014, see Supplementary Table S1).
Transcriptomics revealed a marked sex difference: None of the aforementioned genes showed an expression change in female hippocampus. In the case of Sox6, a sexually dimorphic effect of the three chemicals on miR-24 may play a role. The notion of a link between ER-alpha, miR-24, and Sox6, indicated by our network analyses, postulating an effect of ER-alpha activation on miR-24 on the one hand, and a suppressive effect of miR-24 on Sox6 expression on the other hand, was based on data from non-neural human cell lines and mouse pancreatic beta cells (Maillot et al., 2009; Melkman-Zehavi et al., 2011; Ferraro et al., 2012; Klinge, 2012). In line with this concept, all three chemicals caused a decrease in miR-24-3p and an increase in Sox6 mRNA level in male hippocampus. The inverse correlation between miR-24-3p and Sox6 mRNA level of individual samples was highly significant (p = 0.0097), suggesting an interaction of miR-24-3p with Sox6 in the same cell. In mammalian cells, mRNA levels reflect the action of microRNAs (Guo et al., 2010). In female hippocampus, both, miR-24-3p and its potential target Sox6 were unchanged, indicating a sex difference in responsiveness of miR-24-3p. miR-24 has recently been reported to exhibit sex-dependent responses to diethylhexylphthalate, with suppression in male and induction in female hippocampus (Luu et al., 2017). Hormone-sensitive microRNAs may represent an important novel target of EDCs in developing brain.
ER-alpha also exhibited sex differences in baseline expression levels and response to treatment. In keeping with earlier data of the same postnatal stage (Ivanova and Beyer, 2000), control ER-alpha mRNA levels were higher in males than females. The increased ER-alpha expression caused by the higher dose level of the three chemicals in males corresponds to the effect of estrogen observed in neonatal rat hippocampal slices of undefined sex (Rune et al., 2002). It resulted in a sex-specific increase of the ER-alpha/ER-beta ratio in males. Provided the expression change was reflected at protein level, this effect might have facilitated estrogenic actions of bisphenol A, PCBs or their metabolites, and chlorpyrifos (Hazarika et al., 2020). Yet, the inverse dose-response exhibited by ER-alpha in males and, in the case of 0.5 mg/kg bisphenol A, also in females, makes a direct link to the changes in downstream genes doubtful. Because of insufficient tissue material, analysis of miR-24-3p was not possible in the low-dose groups. However, neither Sox6, the presumed target of miR-24-3p, nor its downstream target Pvalb, showed an inverse dose-response. Rather, we consider the changes in ER-alpha as parallel events; estrogenic actions would also be possible at baseline ER-alpha levels.
The basis for the general low responsiveness of genes in female hippocampus at PND 6 remains unclear. EDCs are known to influence hippocampal structures in both sexes, but data on the first postnatal week are scarce. Ten out of 11 publications reviewed by Rebuli and Patisaul (2016) and a microarray analysis (Royland et al., 2008) reporting neonatal data, were limited to males or did not distinguish between sexes, and the gene expression study in neonatal females perinatally exposed to bisphenol A included only that sex (Facciolo et al., 2005). At PND 6, certain developmental processes in hippocampus are more advanced in females, as indicated by the higher expression ratio of KCC2/NKCC1, that governs the switch of excitatory to inhibitory action of GABA (Nuñez and McCarthy, 2007; Watanabe and Fukuda, 2015). The three chemicals did not affect the KCC2/NKCC1 ratio, thus, did not appear to change the development stage reached by the females, but the sensitivity of hippocampus to treatment effects may have differed between the developmental stages reached by females and males. In marked contrast to males, gene expression in PND 6 females showed high false discovery rates of the majority of genes, suggesting a combination of high inter-individual variability and small magnitude of gene expression changes. A possible role of developmental stage is indicated by the report of Levin and coworkers (2001, 2002) that memory function was impaired by chlorpyrifos after treatment between postnatal days 1–4 in male rats and after treatment between gestational days 17–20 in females.
The differences in susceptibility of neonatal hippocampus in male and female offspring could also result from specific sex differences in estrogenic mechanisms. Expression of both estrogen receptors is higher in males than females (Ivanova and Beyer, 2000), and analyses of cell proliferation suggest that the hippocampus of male rat neonates may be more sensitive to a weak activation of estrogenic mechanisms (Bowers et al., 2010). The fact that developmental hippocampus is able of de novo synthesis of sex steroids (Amateau et al., 2004; Kretz et al., 2004; Nuñez and McCarthy, 2009), and at the same time, may as well be influenced by circulating steroid hormones, complicates an interpretation of EDC effects. At peripubertal stage, there again are prominent effects of prenatal exposure to EDCs in male rat offspring, such as changes in expression and DNA methylation of BDNF (brain-derived neurotrophic factor) after prenatal bisphenol A exposure (Kundakovic et al., 2015), but gene expression changes have also been observed in female hippocampus at that stage, specifically, genes related to serotonin turnover, after the same treatment (Matsuda et al., 2013). Adult female hippocampus shows distinct, estradiol-dependent, gene expression changes during the estrous cycle (Iqbal et al., 2020), illustrating estrogen-sensitivity of hippocampal genes in the female at later life stages.
Investigations of memory functions in adult rodent offspring of both sexes perinatally exposed to one of the three chemicals show a heterogeneous picture with effects only in males or females, or differential or similar effects in males and females (Supplementary Table S1). In children, two studies reported a link between prenatal bisphenol A exposure and preferential impairment of cognitive functions in boys (Braun et al., 2017; Bornehag et al., 2021). A higher susceptibility of boys has also been observed with respect to other types of behaviors in relation to prenatal phthalate exposure (Engel et al., 2010). Many investigations do not report differences between boys and girls (e.g., Boucher et al., 2009, PCBs). One reason for the difficulty to associate hormone-regulated effects in developing hippocampus with behavioral outcome in adulthood could be that the behavioral parameters are not sufficiently detailed. Sex-dependent features of adult hippocampal function result from the activity of hormone-sensitive genes at multiple developmental stages, and, in addition, are influenced by circulating as well as locally synthetized sex steroids in a stage-dependent manner. This cannot be understood along a scheme of sexual brain differentiation leading to either male-or female-type behaviours, but rather, may lead to a mosaic of multiple functional processes in hippocampus, as well as in cortical and subcortical regions interacting with hippocampus, that are more or less influenced by sex steroids during development.
The change in Sox6 expression in male hippocampus was accompanied by an expression change of the Sox6-regulated Pvalb gene. Sox6 exerts different functions depending on developmental stage. It maintains stemness of neural precursor/progenitor cells (Ohta et al., 2013) and oligodendrocyte progenitors (Baroti et al., 2016), and, at later stages, plays a role in specification of cortical and hippocampal interneurons (Batista-Brito et al., 2009; Tricoire et al., 2011; Kelsom and Lu, 2013). Specifically, Sox6 regulates subtype determination MGE-derived postmitotic interneurons by suppressing differentiation of parvalbumin expressing interneurons and inducing differentiation of somatostatin expressing interneurons (Hu et al., 2017). The highly significant (p = 0.0039) inverse correlation between Sox6 and Pvalb indicates that the increased expression of Sox6 was functionally effective to downregulate Pvalb expression at the level of individual hippocampus samples. No relationship was seen between levels of Sox6 and of Nos1 (nNos), which is expressed by MGE-derived interneurons but is Sox6-independent in hippocampus (Jaglin et al., 2012), and between Sox6 and Htr3a, a marker of CGE-derived interneurons (Vucurovic et al., 2010). Developmental exposure to the three chemicals thus appeared to elicit a cascade of events involving a specific interneuron subtype in male but not female hippocampus, with downregulation of miR-24, upregulation of Sox6, and downregulation of Pvalb. There were no expression changes of markers of CGE-derived interneurons, or genes expressed by MGE-derived interneurons but not regulated by Sox6. An imbalance between interneurons is indicated by an increased Pvalb/Nos1 expression ratio. In the absence of morphological data, we do not know whether this represents a change in the number of Pvalb expressing cells, or in intracellular level of Pvalb mRNA. Since it had not been possible to perform perfusion-fixation, we tried to do immunocytochemistry on cryostat sections of brains frozen on dry-ice, but because of the ice crystals formed during freezing, tissue quality turned out to be insufficient for morphometric analysis. If the number of Pvalb interneurons was affected, this might disturb the balance between Pvalb interneurons (fast spiking basket cells and axo-axonic chandelier cells) (Donato et al., 2013), and other types of interneurons.
Our investigation revealed yet another target mechanism, by which hippocampal interneuron development may be disturbed: Expression of Nrg1 and its receptor Erbb4, which control interneuron migration (Faux et al., 2012), was altered, Nrg1 by Aroclor 1254 and bisphenol A, Erbb4 by the PCB mixture and chlorpyrifos. Thus, hippocampal interneuron development emerges as a potentially important target of EDCs.
Conclusion and Regulatory Perspectives
Our investigation disclosed a potential role of genes of the Sox and Pou families involved in cell-autonomous regulation of development, which so far seem to have received little attention in neurodevelopmental effects of EDCs. These genes were found to be indirect but common targets of three different EDCs known to impair hippocampal functions in adult offspring, suggesting that they occupy a central position in the pathogenetic process leading to adverse behavioral outcome. At the developmental stage studied, these genes exhibited a marked sexual dimorphism. In the case of Sox6, the sex difference was correlated with sex-dependent regulation of miR-24. This underscores the potential relevance of microRNAs in developmental effects of EDCs. Developing interneurons were found to represent an important EDC target, as demonstrated by the inverse correlation between mRNA expression levels of Sox6 and of parvalbumin, the marker of a Sox6-regulated subpopulation of hippocampal interneurons. A role of interneurons in developmental effects of EDCs is further suggested by expression changes of Nrg1 and its receptor Erbb4, regulating interneuron migration. It is evident that our descriptive findings need to be further corroborated by mechanistic studies analyzing the proposed EDC–miR–target gene interactions at the cellular level, and by morphological analyses of interneuron populations.
There is a great need to better link in vitro models designed for regulatory purposes, with adverse outcomes in vivo. Omics data could represent such links, but for that purpose, we need more information on the relationship between effects of EDCs on molecular processes during early stages of brain development, and delayed adverse outcome, because in vitro models focus on early developmental processes. Information on omics obtained at later postnatal and adult stages is of interest for the interpretation of functional alterations, but such analyses miss effects on early possible targets such as morphogenetic genes, which are visible only during restricted time-windows. Transcriptomics represents only one part of the approach. Since gene expression patterns change continuously throughout ontogeny, transcriptomics analyses provide information on signaling pathways targeted by a chemical at a given developmental stage, which can then be compared to in vitro systems, and they disclose potential target genes that could be studied by epigenetics. It is the information obtained with such a combined approach, complemented by metabolomics, which should increase the predictive value of in vitro and in vivo models.
Preliminary Reports
Preliminary reports were presented at Copenhagen Workshop on Endocrine Disrupters (COW 2017), Copenhagen, 2–5 May 2017, EUROTOX 2017, Bratislava, 10–13 September 2017, and 10th International Meeting on Steroids and Nervous System, Torino, Italy, 16–20 February 2019.
Data Availability Statement
The RNAseq data have been uploaded to GEO/ENA, ID: GSE180054.
Ethics Statement
The animal study was reviewed and approved by Ethical Committee of Complutense University, Madrid.
Author Contributions
CB-G: Tissue processing for RNA-seq, Real Time RT PCR, microRNA analysis JG: RNA-seq assay WL: Project design, animal experiments, transcriptome and data analysis, manuscript HR: RNA-seq assay, transcriptome analysis MS: Project design, brain preparation and dissection of hippocampus JT: Conduct of animal experiments.
Funding
The investigation was supported by a grant of the Swiss Federal Office for Public Health (Bundesamt für Gesundheit, Bern, Verfügung 13.005616).
Conflicts of Interest
Authors WL and MS are owners of GREEN Tox GmbH and worked without remuneration for the project. GREEN Tox GmbH did not make profit of the project. Funding by the Swiss Federal Office for Public Health covered exclusively the salary of CB-G, costs of animal experiments, and laboratory expenses.
The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
We thank Rocio Campon Montilla and Daniel Campon Montilla (Madrid) for their excellent technical assistance.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2021.718796/full#supplementary-material
References
Aldridge, J. E., Levin, E. D., Seidler, F. J., and Slotkin, T. A. (2005). Developmental Exposure of Rats to Chlorpyrifos Leads to Behavioral Alterations in Adulthood, Involving Serotonergic Mechanisms and Resembling Animal Models of Depression. Environ. Health Perspect. 113, 527–531. doi:10.1289/ehp.7867
Amateau, S. K., Alt, J. J., Stamps, C. L., and McCarthy, M. M. (2004). Brain Estradiol Content in Newborn Rats: Sex Differences, Regional Heterogeneity, and Possible De Novo Synthesis by the Female Telencephalon. Endocrinology 145, 2906–2917. doi:10.1210/en.2003-1363
Andersen, H. R., Vinggaard, A. M., Høj Rasmussen, T., Gjermandsen, I. M., and Cecilie Bonefeld-Jørgensen, E. (2002). Effects of Currently Used Pesticides in Assays for Estrogenicity, Androgenicity, and Aromatase Activity In Vitro. Toxicol. Appl. Pharmacol. 179, 1–12. doi:10.1006/taap.2001.9347
Azim, E., Jabaudon, D., Fame, R. M., and Macklis, J. D. (2009). SOX6 Controls Dorsal Progenitor Identity and Interneuron Diversity during Neocortical Development. Nat. Neurosci. 12, 1238–1247. doi:10.1038/nn.2387
Baroti, T., Zimmermann, Y., Schillinger, A., Liu, L., Lommes, P., Wegner, M., et al. (2016). Transcription Factors Sox5 and Sox6 Exert Direct and Indirect Influences on Oligodendroglial Migration in Spinal Cord and Forebrain. Glia 64, 122–138. doi:10.1002/glia.22919
Batista-Brito, R., Rossignol, E., Hjerling-Leffler, J., Denaxa, M., Wegner, M., Lefebvre, V., et al. (2009). The Cell-Intrinsic Requirement of Sox6 for Cortical Interneuron Development. Neuron 63, 466–481. doi:10.1016/j.neuron.2009.08.005
Bornehag, C.-G., Engdahl, E., Unenge Hallerbäck, M., Wikström, S., Lindh, C., Rüegg, J., et al. (2021). Prenatal Exposure to Bisphenols and Cognitive Function in Children at 7 Years of Age in the Swedish SELMA Study. Environ. Int. 150, 106433. doi:10.1016/j.envint.2021.106433
Boucher, O., Muckle, G., and Bastien, C. H. (2009). Prenatal Exposure to Polychlorinated Biphenyls: A Neuropsychologic Analysis. Environ. Health Perspect. 117, 7–16. doi:10.1289/ehp.11294
Bowers, J. M., Waddell, J., and McCarthy, M. M. (2010). A Developmental Sex Difference in Hippocampal Neurogenesis Is Mediated by Endogenous Oestradiol. Biol. Sex Dif 1, 8–20. doi:10.1186/2042-6410-1-8
Braquenier, J.-B., Quertemont, E., Tirelli, E., and Plumier, J.-C. (2010). Anxiety in Adult Female Mice Following Perinatal Exposure to Chlorpyrifos. Neurotoxicol. Teratol. 32, 234–239.
Braun, J. M., Muckle, G., Arbuckle, T., Bouchard, M. F., Fraser, W. D., Ouellet, E., et al. (2017). Associations of Prenatal Urinary Bisphenol A Concentrations with Child Behaviors and Cognitive Abilities. Environ. Health Perspect. 125, 067008. doi:10.1289/EHP984
Chen, X.-P., Chen, W.-Z., Wang, F.-S., and Liu, J.-X. (2012). Selective Cognitive Impairments Are Related to Selective hippocampus and Prefrontal Cortex Deficits after Prenatal Chlorpyrifos Exposure. Brain Res. 1474, 19–28. doi:10.1016/j.brainres.2012.07.036
Chhabra, R., Dubey, R., and Saini, N. (2010). Cooperative and Individualistic Functions of the microRNAs in the miR-23a∼27a∼24-2 Cluster and its Implication in Human Diseases. Mol. Cancer 9, 232–247. doi:10.1186/1476-4598-9-232
Colciago, A., Casati, L., Mornati, O., Vergoni, A. V., Santagostino, A., Celotti, F., et al. (2009). Chronic Treatment with Polychlorinated Biphenyls (PCB) during Pregnancy and Lactation in the Rat. Toxicol. Appl. Pharmacol. 239, 46–54. doi:10.1016/j.taap.2009.04.023
Danglot, L., Triller, A., and Marty, S. (2006). The Development of Hippocampal Interneurons in Rodents. Hippocampus 16, 1032–1060. doi:10.1002/hipo.20225
Decastro, B. R., Korrick, S. A., Spengler, J. D., and Soto, A. M. (2006). Estrogenic Activity of Polychlorinated Biphenyls Present in Human Tissue and the Environment. Environ. Sci. Technol. 40, 2819–2825. doi:10.1021/es051667u
Dobin, A., Davis, C. A., Schlesinger, F., Drenkow, J., Zaleski, C., Jha, S., et al. (2013). STAR: Ultrafast Universal RNA-Seq Aligner. Bioinformatics 29, 15–21. doi:10.1093/bioinformatics/bts635
Donato, F., Rompani, S. B., and Caroni, P. (2013). Parvalbumin-expressing Basket-Cell Network Plasticity Induced by Experience Regulates Adult Learning. Nature 504, 272–276. doi:10.1038/nature12866
Ejaredar, M., Lee, Y., Roberts, D. J., Sauve, R., and Dewey, D. (2017). Bisphenol A Exposure and Children's Behavior: A Systematic Review. J. Expo. Sci. Environ. Epidemiol. 27, 175–183. doi:10.1038/jes.2016.8
Engel, S. M., Miodovnik, A., Canfield, R. L., Zhu, C., Silva, M. J., Calafat, A. M., et al. (2010). Prenatal Phthalate Exposure Is Associated with Childhood Behavior and Executive Functioning. Environ. Health Perspect. 118, 565–571. doi:10.1289/ehp.0901470
Facciolo, R. M., Madeo, M., Alò, R., Canonaco, M., and Dessì-Fulgheri, F. (2005). Neurobiological Effects of Bisphenol A May Be Mediated by Somatostatin Subtype 3 Receptors in Some Regions of the Developing Rat Brain. Toxicol. Sci. 88, 477–484. doi:10.1093/toxsci/kfi322
Faux, C., Rakic, S., Andrews, W., and Britto, J. M. (2012). Neurons on the Move: Migration and Lamination of Cortical Interneurons. Neurosignals 20, 168–189. doi:10.1159/000334489
Ferguson, S. A., Law, C. D., and Abshire, J. S. (2012). Developmental Treatment With Bisphenol A Causes Few Alterations on Measures of Postweaning Activity and Learning. Neurotoxicol. Teratol. 34, 598–606.
Ferraro, L., Ravo, M., Nassa, G., Tarallo, R., De Filippo, M. R., Giurato, G., et al. (2012). Effects of Oestrogen on microRNA Expression in Hormone-Responsive Breast Cancer Cells. Horm. Canc. 3, 65–78. doi:10.1007/s12672-012-0102-1
Fujimoto, T., Kubob, K., and Aoua, S. (2006). Prenatal Exposure to Bisphenol A Impairs Sexual Differentiation of Exploratory Behavior and Increases Depression-Like Behavior in Rats. Brain Res. 1068, 49–55.
Gauger, K. J., Kato, Y., Haraguchi, K., Lehmler, H.-J., Robertson, L. W., Bansal, R., et al. (2004). Polychlorinated Biphenyls (PCBs) Exert Thyroid Hormone-like Effects in the Fetal Rat Brain but Do Not Bind to Thyroid Hormone Receptors. Environ. Health Perspect. 112, 516–523. doi:10.1289/ehp.6672
Gelman, D. M., Martini, F. J., Nóbrega-Pereira, S., Pierani, A., Kessaris, N., and Marín, O. (2009). The Embryonic Preoptic Area Is a Novel Source of Cortical GABAergic Interneurons. J. Neurosci. 29, 9380–9389. doi:10.1523/jneurosci.0604-09.2009
Gilbert, M. E., Mundy, W. R., and Crofton, K. M. (2000). Spatial Learning and Long-Term Potentiation in the Dentate Gyrus of the Hippocampus in Animals Developmentally Exposed to Aroclor 1254. Toxicol. Sci. 57, 102–111.
González, M., Cabrera-Socorro, A., Pérez-García, C. G., Fraser, J. D., López, F. J., Alonso, R., et al. (2007). Distribution Patterns of Estrogen Receptor α and β in the Human Cortex and hippocampus during Development and Adulthood. J. Comp. Neurol. 503, 790–802. doi:10.1002/cne.21419
Gothié, J.-D., Demeneix, B., and Remaud, S. (2017). Comparative Approaches to Understanding Thyroid Hormone Regulation of Neurogenesis. Mol. Cell Endocrinol. 459, 104–115. doi:10.1016/j.mce.2017.05.020
Guo, H., Ingolia, N. T., Weissman, J. S., and Bartel, D. P. (2010). Mammalian microRNAs Predominantly Act to Decrease Target mRNA Levels. Nature 466, 835–840. doi:10.1038/nature09267
Haviland, J. A., Butz, D. E., and Porter, W. P. (2010). Long-Term Sex Selective Hormonal and Behavior Alterations in Mice Exposed to Low Doses of Chlorpyrifos in Utero. Reprod. Toxicol. 29, 74–79.
Hazarika, J., Ganguly, M., Borgohain, G., Baruah, I., Sarma, S., Bhuyan, P., et al. (2020). Disruption of Androgen Receptor Signaling by Chlorpyrifos (CPF) and its Environmental Degradation Products: a Structural Insight. J. Biomol. Struct. Dyn. 39, 1–12. doi:10.1080/07391102.2021.1875885
Hojo, Y., Hattori, T.-a., Enami, T., Furukawa, A., Suzuki, K., Ishii, H.-T., et al. (2004). Adult Male Rat hippocampus Synthesizes Estradiol from Pregnenolone by Cytochromes P45017 and P450 Aromatase Localized in Neurons. Proc. Natl. Acad. Sci. 101, 865–870. doi:10.1073/pnas.2630225100
Hu, J. S., Vogt, D., Lindtner, S., Sandberg, M., Silberberg, S. N., and Rubenstein, J. L. R. (2017). Coup-TF1 and Coup-TF2 Control Subtype and Laminar Identity of MGE-Derived Neocortical Interneurons. Development 144, 2837–2851. doi:10.1242/dev.150664
Iqbal, J., Tan, Z.-N., Li, M.-X., Chen, H.-B., Ma, B., Zhou, X., et al. (2020). Estradiol Alters Hippocampal Gene Expression during the Estrous Cycle. Endocr. Res. 45, 84–101. doi:10.1080/07435800.2019.1674868
Ivanova, T., and Beyer, C. (2000). Ontogenetic Expression and Sex Differences of Aromatase and Estrogen Receptor-Α/β mRNA in the Mouse hippocampus. Cel Tissue Res. 300, 231–237. doi:10.1007/s004410000199
Jaglin, X. H., Hjerling-Leffler, J., Fishell, G., and Batista-Brito, R. (2012). The Origin of Neocortical Nitric Oxide Synthase-Expressing Inhibitory Neurons. Front. Neural Circuits 6, 44–16. doi:10.3389/fncir.2012.00044
Jašarevic, E., Sieli, P. T., Twellman, E. E., Welsh, T. H., Schachtman, T. R., Roberts, R. M., et al. (2011). Disruption of Adult Expression of Sexually Selected Traits by Developmental Exposure to Bisphenol A. Proc. Natl. Acad. Sci. 108, 11715–11720. doi:10.1073/pnas.1107958108
Johnson, F. O., Chambers, J. E., Nail, C. A., Givaruangsawat, S., and Carr, R. L. (2009). Developmental Chlorpyrifos and Methyl Parathion Exposure Alters Radial-Arm Maze Performance in Juvenile and Adult Rats. Toxicol. Sci. 109, 132–142. doi:10.1093/toxsci/kfp053
Johnson, S. A., Javurek, A. B., Painter, M. S., Ellersieck, M. R., Welsh, T. H., Camacho, L., et al. (2016). Effects of Developmental Exposure to Bisphenol A on Spatial Navigational Learning and Memory in Rats: A CLARITY-BPA Study. Horm. Behav. 80, 139–148. doi:10.1016/j.yhbeh.2015.09.005
Kelsom, C., and Lu, W. (2013). Development and Specification of GABAergic Cortical Interneurons. Cell Biosci 3, 19. doi:10.1186/2045-3701-3-19
Klinge, C. M. (2012). miRNAs and Estrogen Action. Trends Endocrinol. Metab. 23, 223–233. doi:10.1016/j.tem.2012.03.002
Kretz, O., Fester, L., Wehrenberg, U., Zhou, L., Brauckmann, S., Zhao, S., et al. (2004). Hippocampal Synapses Depend on Hippocampal Estrogen Synthesis. J. Neurosci. 24, 5913–5921. doi:10.1523/jneurosci.5186-03.2004
Kubo, K., Arai, O., Ogata, R., Omura, M., Hori, T., and Aou, S. (2001). Exposure to Bisphenol A During the Fetal and Suckling Periods Disrupts Sexual Differentiation of the Locus Coeruleus and of Behavior in the Rat. Neurosci. Lett. 304, 73–76.
Kumar, D., and Thakur, M. K. (2014). Perinatal Exposure to Bisphenol-A Impairs Spatial Memory through Upregulation of Neurexin1 and Neuroligin3 Expression in Male Mouse Brain. PLoS ONE 9 (10), e110482. doi:10.1371/journal.pone.0110482
Kundakovic, M., Gudsnuk, K., Herbstman, J. B., Tang, D., Perera, F. P., and Champagne, F. A. (2015). DNA Methylation of BDNF as a Biomarker of Early-Life Adversity. Proc. Natl. Acad. Sci. USA 112, 6807–6813. doi:10.1073/pnas.1408355111
Levin, E. D., Addy, N., Baruah, A., Elias, A., Christopher, N. C., Seidler, F. J., et al. (2002). Prenatal Chlorpyrifos Exposure in Rats Causes Persistent Behavioral Alterations. Neurotoxicology and Teratology 24, 733–741. doi:10.1016/s0892-0362(02)00272-6
Levin, E. D., Addy, N., Nakajima, A., Christopher, N. C., Seidler, F. J., and Slotkin, T. A. (2001). Persistent Behavioral Consequences of Neonatal Chlorpyrifos Exposure in Rats. Develop. Brain Res. 130, 83–89. doi:10.1016/s0165-3806(01)00215-2
Liao, Y., Smyth, G. K., and Shi, W. (2013). The Subread Aligner: Fast, Accurate and Scalable Read Mapping by Seed-And-Vote. Nucleic Acids Res. 41 (10), e108. doi:10.1093/nar/gkt214
Lichtensteiger, W., Bassetti-Gaille, C., Faass, O., Axelstad, M., Boberg, J., Christiansen, S., et al. (2015). Differential Gene Expression Patterns in Developing Sexually Dimorphic Rat Brain Regions Exposed to Antiandrogenic, Estrogenic, or Complex Endocrine Disruptor Mixtures: Glutamatergic Synapses as Target. Endocrinology 156, 1477–1493. doi:10.1210/en.2014-15
Liu, W., Chen, X., and Zhang, Y. (2016). Effects of microRNA-21 and microRNA-24 Inhibitors on Neuronal Apoptosis in Ischemic Stroke. Am. J. Transl. Res. 8, 3179–3187.
Love, M. I., Huber, W., and Anders, S. (2014). Moderated Estimation of Fold Change and Dispersion for RNA-Seq Data with DESeq2. Genome Biol. 15, 550–570. doi:10.1186/s13059-014-0550-8
Luu, B. E., Green, S. R., Childers, C. L., Holahan, M. R., and Storey, K. B. (2017). The roles of hippocampal microRNAs in response to acute postnatal exposure to di(2-ethylhexyl) phthalate in female and male rats. NeuroToxicology 59, 98–104. doi:10.1016/j.neuro.2017.02.004
Maillot, G., Lacroix-Triki, M., Pierredon, S., Gratadou, L., Schmidt, S., Bénès, V., et al. (2009). Widespread Estrogen-dependent Repression of microRNAs Involved in Breast Tumor Cell Growth. Cancer Res. 69, 8332–8340. doi:10.1158/0008-5472.can-09-2206
Marchetti, M., Meloni, M., Anwar, M., Al-Haj-Zen, A., Sala-Newby, G., Slater, S., et al. (2020). MicroRNA-24-3p Targets Notch and Other Vascular Morphogens to Regulate post-ischemic Microvascular Responses in Limb Muscles. Ijms 21, 1733. doi:10.3390/ijms21051733
Matsuda, S., Matsuzawa, D., Ishii, D., Tomizawa, H., Sajiki, J., and Shimizu, E. (2013). Perinatal Exposure to Bisphenol A Enhances Contextual Fear Memory and Affects the Serotoninergic System in Juvenile Female Mice. Horm. Behav. 63, 709–716. doi:10.1016/j.yhbeh.2013.03.016
Mei, L., and Nave, K.-A. (2014). Neuregulin-ERBB Signaling in the Nervous System and Neuropsychiatric Diseases. Neuron 83, 27–49. doi:10.1016/j.neuron.2014.06.007
Melkman-Zehavi, T., Oren, R., Kredo-Russo, S., Shapira, T., Mandelbaum, A. D., Rivkin, N., et al. (2011). miRNAs Control Insulin Content in Pancreatic β-cells via Downregulation of Transcriptional Repressors. EMBO J. 30, 835–845. doi:10.1038/emboj.2010.361
Miyagawa, K., Narita, M., Narita, M., Akama, H., and Suzuki, T. (2007). Memory Impairment Associated with a Dysfunction of the Hippocampal Cholinergic System Induced by Prenatal and Neonatal Exposures to Bisphenol-A. Neurosci. Lett. 418, 236–241. doi:10.1016/j.neulet.2007.01.088
Mogi, K., Takanashi, H., Nagasawa, M., and Kikusui, T. (2015). Sex Differences in Spatiotemporal Expression of AR, ERα, and ERβ mRNA in the Perinatal Mouse Brain. Neurosci. Lett. 584, 88–92. doi:10.1016/j.neulet.2014.10.028
Nakamura, K., Itoh, K., Dai, H., Han, L., Wang, X., Kato, S., et al. (2012). Prenatal and Lactational Exposure to Low-Doses of Bisphenol A Alters Adult Mice Behavior. Brain Develop. 34, 57–63. doi:10.1016/j.braindev.2010.12.011
Nuñez, J. L., and McCarthy, M. M. (2009). Resting Intracellular Calcium Concentration, Depolarizing Gamma-Aminobutyric Acid and Possible Role of Local Estradiol Synthesis in the Developing Male and Female hippocampus. Neuroscience 158, 623–634. doi:10.1016/j.neuroscience.2008.09.061
Ohta, S., Misawa, A., Lefebvre, V., Okano, H., Kawakami, Y., and Toda, M. (2013). Sox6 Up-Regulation by Macrophage Migration Inhibitory Factor Promotes Survival and Maintenance of Mouse Neural Stem/progenitor Cells. PLoS ONE 8 (9), e74315. doi:10.1371/journal.pone.0074315
Prange-Kiel, J., Fester, L., Zhou, L., Jarry, H., and Rune, G. M. (2009). Estrus Cyclicity of Spinogenesis: Underlying Mechanisms. J. Neural Transm. 116, 1417–1425. doi:10.1007/s00702-009-0294-x
Rauh, V., Arunajadai, S., Horton, M., Perera, F., Hoepner, L., Barr, D. B., et al. (2011). Seven-year Neurodevelopmental Scores and Prenatal Exposure to Chlorpyrifos, a Common Agricultural Pesticide. Environ. Health Perspect. 119, 1196–1201. doi:10.1289/ehp.1003160
Rebuli, M. E., and Patisaul, H. B. (2016). Assessment of Sex Specific Endocrine Disrupting Effects in the Prenatal and Pre-pubertal Rodent Brain. J. Steroid Biochem. Mol. Biol. 160, 148–159. doi:10.1016/j.jsbmb.2015.08.021
Roegge, C. S., Byung-Woun, S., Crofton, K. M., and Schantz, S. L. (2000). Gestational-lactational Exposure to Aroclor 1254 Impairs Radial-Arm Maze Performance in Male Rats. Toxicol. Sci. 57, 121–130. doi:10.1093/toxsci/57.1.121
Royland, J. E., Wu, J., Zawia, N. H., and Kodavanti, P. R. S. (2008). Gene Expression Profiles in the Cerebellum and hippocampus Following Exposure to a Neurotoxicant, Aroclor 1254: Developmental Effects☆. Toxicol. Appl. Pharmacol. 231, 165–178. doi:10.1016/j.taap.2008.04.022
Rubin, B. S. (2011). Bisphenol A: An Endocrine Disruptor with Widespread Exposure and Multiple Effects. J. Steroid Biochem. Mol. Biol. 127, 27–34. doi:10.1016/j.jsbmb.2011.05.002
Rune, G. M., Wehrenberg, U., Prange-Kiel, J., Zhou, L., Adelmann, G., and Frotscher, M. (2002). Estrogen Up-Regulates Estrogen Receptor α and Synaptophysin in Slice Cultures of Rat hippocampus. Neuroscience 113, 167–175. doi:10.1016/s0306-4522(02)00152-5
Ryan, B. C., and Vandenbergh, J. G. (2006). Developmental Exposure to Environmental Estrogens Alters Anxiety and Spatial Memory in Female Mice.. Horm. Behav. 50, 85–93.
Slotkin, T. A., Cooper, E. M., Stapleton, H. M., and Seidler, F. J. (2013). Does Thyroid Disruption Contribute to the Developmental Neurotoxicity of Chlorpyrifos? Environ. Toxicol. Pharmacol. 36, 284–287. doi:10.1016/j.etap.2013.04.003
Slotkin, T. A., and Seidler, F. J. (2009). Oxidative and Excitatory Mechanisms of Developmental Neurotoxicity: Transcriptional Profiles for Chlorpyrifos, Diazinon, Dieldrin, and Divalent Nickel in PC12 Cells. Environ. Health Perspect. 117, 587–596. doi:10.1289/ehp.0800251
Sugawara, N., Nakai, K., Nakamura, T., Ohba, T., Suzuki, K., Kameo, S., et al. (2006). Developmental and Neurobehavioral Effects of Perinatal Exposure to Polychlorinated Biphenyls in Mice. Arch. Toxicol. 80, 286–292. doi:10.1007/s00204-005-0042-4
Touzot, A., Ruiz-Reig, N., Vitalis, T., and Studer, M. (2016). Molecular Control of Two Novel Migratory Paths for CGE-Derived Interneurons in the Developing Mouse Brain. Development 143, 1753–1765. doi:10.1242/dev.131102
Tricoire, L., Pelkey, K. A., Erkkila, B. E., Jeffries, B. W., Yuan, X., and McBain, C. J. (2011). A Blueprint for the Spatiotemporal Origins of Mouse Hippocampal Interneuron Diversity. J. Neurosci. 31, 10948–10970. doi:10.1523/jneurosci.0323-11.2011
Turgeman, G., Pinkas, A., Slotkin, T. A., Tfilin, M., Langford, R., and Yanai, J. (2011). Reversal of Chlorpyrifos Neurobehavioral Teratogenicity in Mice by Allographic Transplantation of Adult Subventricular Zone-Derived Neural Stem Cells. J. Neurosci. Res. 89, 1185–1193. doi:10.1002/jnr.22631
Ventura, C., Núñez, M., Miret, N., Martinel Lamas, D., Randi, A., Venturino, A., et al. (2012). Differential Mechanisms of Action Are Involved in Chlorpyrifos Effects in Estrogen-dependent or -independent Breast Cancer Cells Exposed to Low or High Concentrations of the Pesticide. Toxicol. Lett. 213, 184–193. doi:10.1016/j.toxlet.2012.06.017
Vucurovic, K., Gallopin, T., Ferezou, I., Rancillac, A., Chameau, P., van Hooft, J. A., et al. (2010). Serotonin 3A Receptor Subtype as an Early and Protracted Marker of Cortical Interneuron Subpopulations. Cereb. Cortex 20, 2333–2347. doi:10.1093/cercor/bhp310
Watanabe, M., and Fukuda, A. (2015). Development and Regulation of Chloride Homeostasis in the central Nervous System. Front. Cel. Neurosci. 9, 371. doi:10.3389/fncel.2015.00371
Widholm, J. J., Clarkson, G. B., Strupp, B. J., Crofton, K. M., Seegal, R. F., and Schantz, S. L. (2001). Spatial Reversal Learning in Aroclor 1254-exposed Rats: Sex-specific Deficits in Associative Ability and Inhibitory Control. Toxicol. Appl. Pharmacol. 174, 188–198. doi:10.1006/taap.2001.9199
Wonders, C. P., and Anderson, S. A. (2006). The Origin and Specification of Cortical Interneurons. Nat. Rev. Neurosci. 7, 687–696. doi:10.1038/nrn1954
Xu, X.-H., Zhang, J., Wang, Y.-M., Ye, Y.-P., and Luo, Q.-Q. (2010). Perinatal Exposure to Bisphenol-A Impairs Learning-Memory by Concomitant Down-Regulation of N-Methyl-D-Aspartate Receptors of hippocampus in Male Offspring Mice. Horm. Behav. 58, 326–333. doi:10.1016/j.yhbeh.2010.02.012
Xu, X., Hong, X., Xie, L., Li, T., Yang, Y., Zhang, Q., et al. (2012). Gestational and Lactational Exposure to Bisphenol-A Affects Anxiety- and Depression-like Behaviors in Mice. Horm. Behav. 62, 480–490. doi:10.1016/j.yhbeh.2012.08.005
Keywords: endocrine disrupter, development, hippocampus, sox6, pou gene, microRNA, sex difference, transcriptomics
Citation: Lichtensteiger W, Bassetti-Gaille C, Rehrauer H, Georgijevic JK, Tresguerres JA and Schlumpf M (2021) Converging Effects of Three Different Endocrine Disrupters on Sox and Pou Gene Expression in Developing Rat Hippocampus: Possible Role of microRNA in Sex Differences. Front. Genet. 12:718796. doi: 10.3389/fgene.2021.718796
Received: 01 June 2021; Accepted: 20 September 2021;
Published: 11 November 2021.
Edited by:
Steffen H. Keiter, Orebro University, SwedenReviewed by:
Laura N Vandenberg, University of Massachusetts Amherst, United StatesAjay Pradhan, Örebro University, Sweden
Copyright © 2021 Lichtensteiger, Bassetti-Gaille, Rehrauer, Georgijevic, Tresguerres and Schlumpf. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Walter Lichtensteiger, walter.lichtensteiger@access.uzh.ch