Persistent autism-relevant behavioral phenotype and social neuropeptide alterations in female mice offspring induced by maternal transfer of PBDE congeners in the commercial mixture DE-71

Polybrominated diphenyl ethers (PBDEs) are ubiquitous persistent organic pollutants (POPs) that are known neuroendocrine disrupting chemicals with adverse neurodevelopmental effects. PBDEs may act as risk factors for autism spectrum disorders (ASD), characterized by abnormal psychosocial functioning, although direct evidence is currently lacking. Using a translational exposure model, we tested the hypothesis that maternal transfer of a commercial mixture of PBDEs, DE-71, produces ASD-relevant behavioral and neurochemical deficits in female offspring. C57Bl6/N mouse dams (F0) were exposed to DE-71 via oral administration of 0 (VEH/CON), 0.1 (L-DE-71) or 0.4 (H-DE-71) mg/kg bw/d from 3 wk prior to gestation through end of lactation. Mass spectrometry analysis indicated in utero and lactational transfer of PBDEs (in ppb) to F1 female offspring brain tissue at postnatal day (PND) 15 which was reduced by PND 110. Neurobehavioral testing of social novelty preference (SNP) and social recognition memory (SRM) revealed that adult L-DE-71 F1 offspring display deficient short- and long-term SRM, in the absence of reduced sociability, and increased repetitive behavior. These effects were concomitant with reduced olfactory discrimination of social odors. Additionally, L-DE-71 exposure also altered short-term novel object recognition memory but not anxiety or depressive-like behavior. Moreover, F1 L-DE-71 displayed downregulated mRNA transcripts for oxytocin (Oxt) in the bed nucleus of the stria terminalis (BNST) and supraoptic nucleus, and vasopressin (Avp) in the BNST and upregulated Avp1ar in BNST, and Oxtr in the paraventricular nucleus. Our work demonstrates that developmental PBDE exposure produces ASD-relevant neurochemical, olfactory processing and behavioral phenotypes that may result from early neurodevelopmental reprogramming within central social and memory networks. Supplementary Information The online version contains supplementary material available at 10.1007/s00204-021-03163-4.


Introduction
Autism spectrum disorder (ASD) is a group of neurodevelopmental conditions defined clinically by deficits in social reciprocity and communication, and restricted interest and repetitive behaviors (American Psychiatric Association 2013). Hallmarks of ASD, as classified by the NIH Research Domain Criteria (RDoC) (Social Processes: Workshop Proceedings, 2012) include disturbances in the social cognition (SC) domain such as facial recognition ability, empathy and evaluation of emotions of others (Weigelt et al. 2012); (Ewbank et al. 2017). The prevalence of ASD has increased dramatically over the past 3 decades. In the United States, the Centers for Disease Control (CDC) estimates that ASD affects 1 in 54 neurotypical children (Maenner et al. 2020), while the worldwide prevalence is estimated to be 1-2% (Kim et al. 2011). While genetic heritability is an important factor in ASD etiology, the incremental incidence of autism over the last several decades, raises the possibility that environmental factors, such as xenobiotic chemicals, may contribute alongside genetic predisposition and influence ASD risk (Grandjean and Landrigan 2014); (Pelch et al. 2019). Although the incidence of autism is four times greater in boys, girls and women with autism are often undiagnosed, misdiagnosed or receive a diagnosis of autism at later age (Rynkiewicz et al. 2019) suggesting underestimation in females. According to the female protective model, females may also benefit from a higher threshold of genetic liability to manifest ASD phenotype (Werling and Geschwind 2013); (Zhang et al. 2020). Nevertheless, female ASD cases may display susceptibility to xenobiotic chemicals (Terasaki et al. 2016) that can potentially enhance the risk of neurodevelopmental disorders (NDDs). Indeed, we have found that female mice offspring exposed to PBDEs during prenatal and postnatal development exhibit endocrine and metabolic disruption, indicating that females may provide a susceptible substrate for studying xenobiotic effects on neurodevelopment .
Polybrominated diphenyl ethers (PBDEs) are a class of brominated flame retardants added to a wide range of products including consumer building material, electronics, textiles, plastics and foams including infant products (Ionas et al. 2016) since the 1970s (Stapleton et al. 2005). Three commercial formulations of PBDEs were prevalent in commerce, including penta-BDE, octa-BDE and deca-BDE. Two commercial PBDE mixtures, penta-and octa-BDEs, were banned in Europe in 2003 and all PBDEs were voluntarily phased out in the US by 2013, leading to a slow, but measurable, decrease in environmental levels as well as in human sera and breastmilk concentrations of some PBDE congeners (Drage et al. 2019); (Guo et al. 2016). Notwithstanding a commitment to a voluntary phase out of deca-BDE by 2013, PBDE contamination is predicted to remain an ongoing problem through the next several decades due to their long half-lives, persistence in e-waste (Ohajinwa et al. 2019), recycling into consumer products and inadvertent reappearance into environment (Abbasi et al. 2019). In an unprecedented action, the U.S. EPA formally banned the production, import and distribution of deca-BDE in February 2021. Nevertheless, PBDEs are still being detected in various tissue samples worldwide, including human breastmilk (Terry et al. 2017); (Hurley et al. 2017); (Lyche et al. 2015); (Chen et al. 2014); (Darrow et al. 2017).
Compared to adults, infants and toddlers are at greater risk of the adverse health effects resulting from PBDE exposure since they disproportionately accumulate 3-to 9-fold greater body burdens (Costa et al. 2014). Circulating levels of PBDEs in US children are 10-to 1000-fold higher than similar age populations in Mexico and Europe (Rose et al. 2010). Elevated exposures in infants are due to the maternal transfer of PBDEs via cord blood and breastmilk (Toms et al. 2008). After weaning in early childhood, an additional route of exposure is dust ingestion and inhalation associated with children's mouthing and crawling behaviors (Stapleton et al. 2008); (Johnson-Restrepo and Kannan 2009). Therefore, high PBDE exposure poses significant health risks during critical periods of development.
Major health effects associated with PBDE exposures are endocrine disruption, reproductive and developmental toxicity and neurotoxicity (Costa and Giordano 2007); (Darnerud 2008); (Kodavanti and Curras-Collazo 2010); ; (Dingemans et al. 2011). However, epidemiological studies examining the association between PBDE exposure and ASD show inconsistent findings. PBDE exposure (e.g., PBDE congeners BDE-153 and -47) during both pre-and post-natal development has been linked to adverse neurological outcomes such as impairments in executive function, poor attention and behavioral regulation, reduced social scores, and lower IQ. Early-life exposure to PBDEs  has been associated with externalizing behaviors such as hyperactivity and impulsivity (Roze et al. 2009); (Ding et al. 2015); (Herbstman et al. 2010); (Hoffman et al. 2012); (Vuong et al. 2018). With regard to the association of PBDEs with social behavior deficits and ASD, preschool-aged children with greater ΣPBDE exposures were rated as less assertive by their teachers (Lipscomb et al. 2017) or showed greater anxious behavior (Adgent et al. 2014). In the HOME prospective cohort study, serum levels of PBDEs in mothers were associated with greater (BDE-28) or fewer  autistic behaviors in their children (Braun et al. 2014). Similarly, significantly higher risk of poor social competence symptoms was shown as a consequence of postnatal BDE-47 exposure (Gascon et al. 2011). Although the possibility that environmental toxicants serve as risk factors for social neurodevelopmental disorders (NDDs) has not been established (Messer 2010), PBDEs may have deleterious effects on children's social development relevant to ASD (Ding et al. 2015); (Messer 2010;Gascon et al. 2011;Braun et al. 2014); (Gibson et al. 2018). Studies in experimental animals demonstrate that certain PBDE congeners produce adverse effects on behavior, learning, and memory in exposed offspring (Costa and Giordano 2007); (Kodavanti and Curras-Collazo 2010); (Pinson et al. 2016), but information about the negative impact of PBDEs on psycho-social behavior is limited ; (Kim et al. 2015). We hypothesized that developmental PBDE exposure produces ASDrelevant social behavioral and neurochemical phenotypes in a mouse toxicant model. Social recognition, or the ability to distinguish between familiar and novel conspecifics, is a fundamental process across species required for forming long-term attachments, hierarchies, and other complex social strategies that enhances survival (Brennan and Kendrick 2006). Disturbances in this capacity are present in individuals with ASD who have difficulties identifying faces of novel conspecifics from those previously encountered (Weigelt et al. 2012);(Ewbank et al. 2017). Rodents, because of their highly social nature, are used as proxies for studying autism-relevant social competence (Young et al. 2002). Mouse social behavior paradigms rely on the natural propensity of mice for investigation of social novelty compared to previously encountered conspecifics when given the choice (Moy et al. 2004). This preference for social novelty has been shown to be absent in monogenetic, idiopathic and environmental models of ASD (Silverman et al. 2010); (Sgritta et al. 2019); (Buffington et al. 2016)). In the current study, we used a toxicant exposure mouse model to characterize social recognition ability, repetitive behaviors and concomitant autism comorbidities such as anxiety, memory impairment and altered olfactory processing.
While the behavioral deficits in typical ASD rodent models are well established, the underlying neural mechanisms are not well understood. The neuropeptides oxytocin and vasopressin are considered major neurotransmitters implicated in social information processing and social cognition that are disrupted in ASD patients (Landgraf and Neumann 2004). Rodent studies have shown that these neuropeptidergic systems are involved in several social cognition domains such as social memory, social/emotional recognition and social reward (Ferguson et al. 2001); (Raam et al. 2017); (Bielsky et al. 2004); (Ferretti et al. 2019). Work by us and our collaborators has provided evidence that PBDEs (and the structural analogues, polychlorinated biphenyls (PCBs)) disrupt the magnocellular neuroendocrine system responsible for vasopressin production involved in osmoregulation, cardiovascular function and social behavior (Kodavanti and Curras-Collazo 2010); (Coburn et al. 2005); (Coburn et al. 2007); (Currás-Collazo 2011); (Coburn et al. 2015); (Mucio-Ramírez et al. 2017);(Alvarez-Gonzalez et al. 2020). We have shown that exposure to DE-71 during in utero and lactation via maternal transfer can nearly abolish vasopressin immunoreactivity in the activated supraoptic (SON) and paraventricular nuclei (PVN) of the hypothalamus (Coburn et al. 2007); (Mucio-Ramírez et al. 2017). Therefore, we also tested the hypothesis that PBDEs disrupt gene expression of prosocial neuropeptides such as vasopressin, oxytocin, PACAP and their receptors in regions of the social brain network, which may underlie deficient social behavior (Bicks et al. 2020); (Tanimizu et al. 2017); (Ferguson et al. 2000).
To lend insight to whether early-life exposure to PBDEs can produce ASD-relevant phenotypes, we exposed mouse dams to a commercial mixture of PBDEs, DE-71, at low doses to mimic chronic, low-level exposure to BDE congeners and doses encountered by infants and toddlers. We demonstrate that perinatal exposure to DE-71 produces dosedependent deficits in social recognition memory and general memory, altered olfactory function and altered neuromolecular phenotypes in brain regions that coordinate complex 1 3 social behaviors. To the best of our knowledge, this study is the first to show a comprehensive profile of autistic-relevant behavior, comorbidities and brain transcriptional changes in female offspring impacted by maternal transfer of PBDEs. Concomitant characterization of ASD-relevant behavioral and neurochemical phenotypes exhibited by offspring developmentally exposed to and reprogrammed by DE-71, provides an integrative framework for exploring environmental risk factors that may contribute to the increasing incidence of ASD. A portion of our findings has been published in preliminary form (Kozlova et al., 2019).

Animal housing and care
C57Bl/6 N mice were generated using breeders obtained from Charles River Labs (West Sacramento, CA). Mice were housed 2-4 per cage in standard polycarbonate plastic cages with corn-cob bedding in a specific pathogen-free vivarium and kept on a 12:12h light:dark cycle in a controlled temperature (21.1-22.8 °C) and humidity (20-70%) environment. Mice were provided rodent chow (Laboratory Rodent Diet 5001; LabDiet, USA) and water ad libitum. Care and treatment of animals was performed in compliance with NIH guidelines and approved by the University of California, Riverside Institutional Animal Care and Use Committee (AUP# 20170026 and 20200018). DE-71 (technical pentabromodiphenyl oxide;Lot no. 1550OI18A), was obtained from Great Lakes Chemical Corporation (West Lafayette, IN). Ninety-seven percent of this mixture contains the following congener composition (in %): BDE-17 (0.11), -28 (0.24), -47 (33.3), -85 (2.54), -99 (45.3), -100 (8.24), -138 (0.45), -139 (0.92), -153 (3.57), -154 (3.19) as described . DE-71 dosing solutions were prepared in corn oil vehicle (VEH/CON) to yield two doses: 0.1 mg/kg/d (L-DE-71) and 0.4 mg/kg/d (H-DE-71) using 2 mL of stock solution/kg body weight. The DE-71 doses were selected to contain the same molar concentrations of BDE-47 used in other mouse studies ; (Wang et al. 2018). BDE-47 has been the focus of many PBDE studies, in part, because it is the primary congener found in human breast milk (Guo et al. 2016); (Darnerud et al. 2015). Offspring were exposed to DE-71 via maternal transfer using a 10-week dosing regimen (Fig. 1a) as described previously . Mice were randomly assigned to one of the three exposure groups: corn oil vehicle control (VEH/CON), L-DE-71 or H-DE-71. This exposure paradigm was chosen to model chronic, low-level exposure to the mother and transfer of PBDEs to infant during gestation (1st, 2nd and 3rd trimester) and lactation as shown in humans (Toms et al. 2008); (Schecter et al. 2007); (Chao et al. 2011); (Zhao et al. 2013). After 3 weeks of pre-dosing, virgin females were paired with an untreated male using harem-style breeding. The presence of a vaginal plug was designated as gestational day (GD) 0. Females that failed to conceive within 10 d were removed from the study. The litters were not culled to avoid potential artificial equalization of variance that masks the developmental and reproductive effects of toxicants as justified previously (Suvorov and Vandenberg 2016). F1 offspring were weaned at PND21 and housed in same-sex cages (2-4/cage). Dams (F0) and their adult female offspring (F1) were subjected to behavioral testing and later sacrificed by exsanguination via cardiac puncture under terminal isoflurane anesthesia (5%) followed by cervical dislocation.

DE-71 exposure and experimental design
To reduce cross-over effects, behavioral tests were distributed across three different cohorts. Mice were run through a battery of behavioral tests in the following order for Cohort 1 (mean age): Suok (PND 46); social novelty preference test (SNP; PND 71); 3 chamber social novelty (PND 87); elevated plus maze (EPM; PND 72). The brains of Cohort 1 were collected at sacrifice on PND 108 and used in RT-qPCR. The following tests were performed on Cohort 2: Marble burying (MB; PND 81); olfactory habituation/dishabituation (OHT; PND 79); olfactory preference (OPT; PND 102); forced swim test (FST; PND 74). Cohort 3 was subjected to social recognition memory (SRMT; PND 30); juvenile open field (OFT; PND 31); juvenile MB (PND 35); novel object recognition (NORT; PND 111) tests. F1 and F0 were tested similarly, except that F0 did not get tested on the SRMT. Analytical characterization by mass spectrometry was performed on brains from Cohort 1 (PND 110) and a subset of Cohort 3 (PND 15). Enzyme-linked immunosorbent assays (ELISA) were performed on plasma from Cohorts 1-3. Whenever possible, the dam or litter was used as the statistical unit of analysis for F1 i.e. values of offspring in each litter were averaged to represent one sample. This allowed us to control for large differences between litters if they existed (Jiménez and Zylka 2021). Based on our previous findings we estimate 5.8, 5.8, 7.2 pups/litter for VEH/CON, L-DE-71 and H-DE-71, respectively, with no differences in the secondary sex ratio . Therefore, there were on average three females per litter per group. A total of 45 litters, distributed equally across experimental groups, were needed to cover the range of experiments conducted in these studies. The number of litters per Cohort is as follows: Cohort 1 (17), Cohort 2 (18), Cohort 3 (10). In addition, results were replicated in a minimum of three independent experiments.

Nest scoring
To test for possible effects of DE-71 on maternal parameters, nests of single-housed dams built from pressed cotton squares (5 × 5 cm; Nestlets) were evaluated at PND 0-1 using a modified scoring system (Hess et al. 2008). Manual scores were assigned by several experimenters according to the height and closure of the walls surrounding the nest cavity. Scores were assigned according to whether the nest contained a center (1) plus a 50% border, (2) 75% border (3) or 100% border (4). A score of 5 was given if the nest resembled a dome ( Supplementary Fig. 1). Nest scores were boosted by 0.5 if the nest was elevated. For interrater reliability, the Bland-Altman method was used to calculate bias as the mean of the differences (0 representing two judges were not producing different results) and precision as 95% limits of agreement (standard deviation of mean bias ± 1.96) (Supplementary Fig. 1).

Congener analysis in adult offspring brain
PBDE concentrations were measured in PND110 whole brain homogenate extracts by the Stapleton laboratory at Duke University using gas chromatography coupled with electron capture negative ion mass spectrometry (GC/ECNI-MS; Agilent 5975 N MS) as described previously . Briefly, approximately 0.2-0.5 g of tissue was first ground with clean sodium sulfate, spiked with two isotopically labeled standards (F-BDE-69 and 13C BDE-209) and then extracted using 50:50 DCM:hexane. Extracts were concentrated, measured for lipid content using a gravimetric analysis, and then purified using acidified silica before analysis for 26 different PBDE congeners ranging from BDE-30 to BDE 209. Laboratory processing blanks (clean sodium sulfate only) were analyzed alongside samples to monitor background contamination. Recoveries of F-BDE-69, and 13C BDE 209, averaged 91 (± 6.9%) and 106 (± 19.9%), respectively, in all samples. All samples were blank-corrected on a congener-specific basis using the average concentrations measured in the laboratory processing blanks. Method detection limits (MDLs) were estimated using either a signal to noise ratio of 10, or, if analytes were detected in laboratory blanks, by calculating three times the standard deviation of the laboratory blanks. MDLs differed by congener and ranged from 0.8 (BDE-47) to 6.6 ng/g (BDE-206).

PBDE congener analysis in postnatal offspring brain
Due to force majeure, i.e. SARS-CoV-2 pandemic, we were unable to carry out planned analytical characterization of PND 15 tissues in collaboration with the Stapleton lab, therefore, the mass spectrometry (MS) system in the Schramm lab was used. The performance of both methods was comparable, especially with regard to the limit of quantification. Using high-resolution gas chromatography-high-resolution mass spectrometry (HRGC/HRMS), PBDE concentrations were measured in PND15 whole brain homogenates (0.1-0.2 g) as described . PBDE analytes included 37 PBDE congeners 10,15,17,28,30,47,49,66,71,77,85,99,100,119,126,138,139,140,153,154,156,176,180,183,184,191,196,197,201,203,204,205,206,207,208,209). Samples were ground and homogenized to a fine powder under liquid nitrogen. Each sample (100-200 mg) was mixed with CHEM TUBE-Hydromatrix (Agilent Technologies) and spiked with 13 C-labelled PBDE standard mix (BFR-LCS, Wellington Laboratories). For pressurized liquid extraction (Dionex ASE 200) n-hexane/ acetone (3:1, v/v) was used at 120 °C and 12 MPa. The volume of the extract was reduced to ~ 5 mL using a vacuum rotary evaporator. Samples were purified using an automated system (DEXTech, LCTech, Germany), where the sample was passed and fractionated over an acidic silica, alumina and carbon column. Concentrated extracts were spiked with the recovery standard (BFR-SCS, Wellington Laboratories) and analyzed by HRGC/HRMS (Agilent 6890/Thermo MAT95 XL) using electron impact ionization (EI), in the selected ion monitoring mode. The instrumental parameters are listed in Supplementary Table 1. Average recovery for 13 C-labelled PBDE standards ranged between 40 and 120%. All samples were blank-corrected on a congener-specific basis using the average of three procedural blank samples. Analytes with concentrations after blank correction that were lower than three times the standard deviation of the blank values or that were not detected before blank correction were considered as not detectable (n.d.). The limit of quantification (LOQ) of the instrumental methodology was considered as a signal/noise ratio of 9:1 (Supplementary Table 2). Congener concentrations that were below the detection limit were assigned a randomly generated value of LOQ/2. The accuracy of our method was confirmed by successful participation in interlaboratory comparison studies.

Comparison of MS methods
The GC/ECNI-MS method used the ECNI ionization mode to improve sensitivity and provides equal sensitivity to HRGC/HRMS that uses electron impact.

Neurobehavioral testing paradigms
At least 30 min prior to testing, mice were moved to a designated behavior room. Ethanol (70%) was used to remove debris and odors from apparati between individual mouse trials. Unless stated otherwise, mouse behavior was scored using automated video-tracking software (Ethovision XT 15, Noldus) or manual scoring software (BORIS (Friard and Gamba 2016) or JWatcher), performed blind to treatment by trained observers. Mice were tested between 10 am and 4 pm during the light phase under bright light conditions, unless otherwise stated.

Social novelty preference
Social novelty preference (SNP) was conducted and analyzed according to methods adopted from published protocols (Moy et al. 2004). Briefly, mice were habituated for 30 min to a polycarbonate cage identical to their home cage (26 cm length × 16 width cm × 12 cm height), followed by 30 min to two wire interaction corrals (11 cm height × 10 cm diameter) placed on each side of the cage. During a 5-min training trial, a stimulus mouse was placed into one corral while the empty corral was removed. After a 30 min retention period, social recognition was assessed in the following 5 min test, during which the test mouse explored the same stimulus mouse (now familiar mouse) versus a novel stimulus mouse. Prior to testing days, sex-and age-matched conspecific stimulus mice were trained to stay in corrals for 15 min for 3 times per day for 7-14 d. Stimulus mice were single-housed to preserve their unique scent. Investigation by test mouse was measured as time spent sniffing (snout within 2 cm of stimulus). Test robustness was measured using an Investigation Index calculated as the ratio of time spent investigating the novel mouse to total investigation time during the training trial ( Supplementary Fig. 2). Social recognition is represented as time spent investigating novel stimulus as percent of total investigation time during test trial. To evaluate between group differences, a Recognition Index was calculated as the time spent investigating Novel minus Familiar/total investigation time in test period.

Three-chamber sociability task
Sociability was assessed as described (Yang et al. 2011). In brief, during the first habituation phase, test mice were habituated for 10 min to the center chamber of a Plexiglass three-chambered apparatus (22 × 40 × 23 cm). Next, the retractable doors partitioning the chambers were opened to permit exploration of all three chambers (second habituation phase). Sociability was tested in the following 10 min session, during which the test mouse was permitted to explore an empty 9 × 27 cm corral (novel object) versus a different corral occupied by a stimulus mouse (novel social object). Inherent side preference during the second habituation phase was evaluated as right chamber time minus left chamber time/right chamber time + left chamber time × 100. Mean values for test mice meeting the inclusion criterion (0 ± 15%) are shown in Supplementary Fig. 2. Sociability was analyzed during the subsequent testing phase both as time spent in chamber and time spent sniffing within 2 cm of stimulus.

Marble burying and nestlet shredding tests
The marble burying (MB) and nestlet shredding tests were utilized for analysis of elicited repetitive behavior in rodents that is considered analogous to that observed in autistic individuals (Silverman et al. 2010). During the marble burying test, the test mouse was placed in the corner of a polycarbonate cage (19 × 29 × 13 cm) containing 5 cm of bedding (Angoa-Pérez et al. 2013) and allowed to interact for 30 min with an array of equidistant marbles (8 × 4 for adults or 5 × 4 for juvenile). A minimum 2/3 of the marble was defined as being buried in the 32-marble array and 1/2 buried in the 20-marble array. Images of the cage were scored by 2-3 experimenters who were blind to treatment and a mean score obtained. For interrater reliability on marble burying the Bland-Altman method was used to calculate bias as the mean of the differences (0 representing that two judges were not producing different results) and precision as 95% limits of agreement (standard deviation of mean bias ± 1.96) (Supplementary Fig. 3). After a 5 min rest period, the test mouse was placed into another cage of the same size with 0.5 cm of bedding containing a pre-weighted square of cotton fiber (Nestlet). After 30 min, the remaining Nestlet was weighed and percent shredding calculated.

Social recognition memory test
A two-trial social recognition memory test (SRMT) was performed as previously described (Tanimizu et al. 2017) to assess long-term social recognition memory. Test mice (PND 28-40) were exposed to a juvenile sex-matched conspecific stimulus mouse (PND 15-32) during two 3 min trials following an intertrial delay of 24 h. For each experiment, test mice were individually placed into polycarbonate cages (26 × 16 × 12 cm) and allowed to habituate for 1 h under dim conditions. A juvenile sex-matched conspecific was then placed into the cage, and the mice were allowed to interact for 3 min during trial 1 (Day 1). In trial 2, performed 24 h later on Day 2, the same test mouse was exposed to the same stimulus mouse (familiar from Day 1). Each stimulus was not used more than 4 times per day. The tests were digitally recorded and scored for social investigation behavior. To evaluate the differences in the ability of mice to form a long-term social memory a Recognition Index (RI) was calculated as the ratio of the duration of investigation of stimulus mouse on Day 2 over Day 1. We used a one-sample t test to determine if the mean RI of each group was statistically different from the previously reported mean RI of 0.65 (Kogan et al. 2000); (Tanimizu et al. 2017). In a second set of experiments, test mice were exposed to a novel stimulus mouse on Day 2. In this context, we used a one-sample t test to determine if the sample mean RI was statistically different from 1, representing similar preference for mouse presented on Day 1 and new novel stimulus mouse presented on Day 2.

Novel object recognition test
The novel object recognition test (NORT) was used to assess non-social recognition memory. We adapted an established protocol to a two-day protocol that tested the same mouse using short-and long-term retention times. On Day 1, the test mouse was habituated to an empty square Plexiglas open field arena (39 × 39 × 38 cm) for 15 min as described (Murai et al. 2007), followed by a 20 min rest in its home cage. During the acquisition phase, the test mouse was placed in the open field containing two identical objects (F vs F`) and allowed to freely explore the environment and objects. During the short-term memory (30 min retention) testing session, the test mouse was again placed in the apparatus and allowed to explore a familiar and novel object (F vs N). After a 24 h retention time (Day 2), long-term memory was assessed by placing mice into the open field containing both the familiar and a new novel object (F vs N`). All test/train sessions lasted 5 min. Preference for the novel object was expressed as the ratio of time exploring the novel stimulus relative to the total exploration time. To evaluate the differences in ability to form NOR memory, a Discrimination Index was calculated as the difference in exploration time between novel and familiar objects relative to total exploration time, where 0 indicates equal preference. Test objects were first validated for intrinsic preference in untreated mice. After analysis of the data using Ethovision, we applied the following exclusion criteria: a) total distance travelled one or more standard deviation(s) lower than the group mean for any trial, or b) less than 6 visits to the familiar or novel target zones.

Innate olfactory preference test
To test the ability of mice to detect attractive or aversive odorants, the innate Olfactory Preference Test (OPT) was performed and analyzed as described (Kobayakawa et al. 2007). Mice were habituated to the experimental conditions. First, they were placed individually into an empty test cage (19 × 29 × 13 cm) and sequentially transferred to three other cages every 15 min. After the final habituation, mice were transferred into the test cage containing a filter paper (2 × 2 cm) infiltrated with either 500 uL of a fresh solution of test odorants: 10% peanut butter, 1% vanilla, 1% butyric acid, or deionized water. The four test odorants were presented to the test mouse in a randomized order. Time spent sniffing the filter paper during the 3-min odorant trials was video-recorded and later measured.

Olfactory habituation test
The ability of mice to detect and differentiate social and nonsocial odorants was examined using the olfactory habituation/dishabituation test (OHT) (Silverman et al. 2010). OHT involves presenting a test animal with various non-social and social odorants. Mice were acclimated for 45 min to an empty mouse cage (26 × 16 × 12 cm) and lid with a cotton-tipped applicator inserted through the water bottle hole to reduce the novelty of the applicator during test sessions. Non-social odors were prepared from extracts immediately before testing. They included: (1) deionized water; (2) almond (1:100; McCormick); (3) banana (1:100; McCormick). Two social odors were obtained the morning of test day by swiping applicator across the bottom of two different stimulus mouse cages containing soiled bedding from sex-matched conspecifics. Cages housed 3-4 mice and bedding was at least 3 d old. Stimuli were presented in 2-min triplicate trials in the following order: water, almond, banana, social odor 1, social odor 2. Time spent sniffing the applicator was recorded with a stopwatch. Parameters measured were habituation, defined as a decrement in olfactory investigation of the same odor after repeated presentations and dishabituation, defined as a reinstatement of olfactory investigation upon presentation of a new odorant.

Suok
Suok is an elevated platform behavioral paradigm used to analyze anxiety, anxiety-induced motor impairments and motor-vestibular anomalies in mice. The apparatus consists of a smooth aluminum beam (2 m long, 3 cm diameter) elevated to 20 cm and fixed to two clear acrylic walls as described (Kalueff et al. 2008). Bilateral to a central segment (38 cm) of the aluminum rod are 10 cm segments labeled by line markings. After acclimation to the dimly lit testing room, mice were placed in the center of the beam and several behaviors were scored over a 5 min trial: (1) horizontal and locomotor (normalized) activity, assessed by number of segments traveled, (2) sensorimotor coordination, measured by the number of hind leg slips and falls from the rod, (3) exploratory behavior such as side looks and head dips, (4) anxiogenic behaviors such as increased latency to leave the central zone and unprotected stretch attend postures (SAP), in which the mouse stretches forward and retracts without moving its feet (considered a non-social form of ambivalence), (5) vegetative responses (combined number of urinations and defecation boli), and (6) autogrooming behaviors.
Hyperactivity, loss of sensorimotor coordination, increased anxiety and displacement behavior are represented by elevated values for #1, 2, 4 and #5, and 6, respectively. Locomotor activity was calculated as total test time minus time spent immobile in center. Measures were recorded manually by stopwatch.

Open field test
The open field test allows rapid assessment of rodent locomotion, anxiety and habituation without a training requirement (Hall 1934). The open field apparatus, a Plexiglas square arena of 39 × 39 × 38 cm was designed as a large, brightly lit, open and aversive environment. Locomotor and other activity over a 1 h period was digitally recorded and scored using Ethovision for distance traveled, velocity and total time in periphery (10 cm adjacent to wall) and center.

RNA extraction from brain micropunches
At sacrifice, using isoflurane anesthesia and cervical dislocation, whole brains were rapidly dissected and snap frozen in 2-methylbutane over dry ice. Brains were cryosectioned (0.3 mm thick) coronally and sections mounted on sterile glass slides and stored at -80 0 C. Five regions of interest were punched out bilaterally from tissue sections under a stereomicroscope using a microdissection blunted needle (16-gauge) adapted from the Palkovits micropunch technique (Palkovits 1973). The anatomical precision was determined by comparing cryosections with cresyl violet stained sections of reference mouse brains and corresponding anatomical plates in the atlas of Paxinos and Franklin. Tissue punches were immediately homogenized in TRIzol Reagent (Thermo Fisher Scientific, USA) using a hand-held homogenizer. Total RNA was prepared via a modified partial phenol-methanol extraction protocol (RNeasy Micro Kit, Qiagen, USA). Purity and quantity of RNA were assessed by determining the optical density (OD) photometrically using 260/280 nm and 260/230 nm ratios (NanoDrop ND-2000, Thermo-Fisher Scientific Inc., Waltham, MA, USA). RNA integrity was assessed using an Agilent 2100 Bioanalyzer (Agilent Technologies Inc. Santa Clara, CA, USA) (Supplementary data 2).

Quantitative polymerase chain reaction (RT-qPCR) analysis of brain micropunches
RT-qPCR was used to quantitate mRNA transcripts for pro-social peptides, AVP, OXT, PACAP and their receptors. Custom-or predesigned DNA oligonucleotide PCR primers were obtained from Integrated DNA Technologies.
Primers were designed to meet several criteria using NCBI Primer Blast and then optimized by testing against complementary DNA generated using RT-PCR and gel electrophoresis. Only primers that gave single-band amplicons in the presence of RT and that matched the base length of the predicted target were selected. In addition, primers selected yielded 89% to 107% efficiency on RT-qPCR (Table 1). Oxtr and the reference gene, ActB, were multiplexed using hydrolysis probes with double-quenchers. For all other primers, intercalating dye chemistry was used. RT-qPCR was performed on RNA (1-4 ng) samples, run in triplicate, on a CFX Connect (Bio-Rad, USA) thermocycler with the Luna Universal or Probe one-step qPCR Master Mixes (New England Biolabs, Ispwich, MA). In each experiment, no-template controls (NTCs) without mRNA were run to rule out extraneous nucleic acid contamination and primer dimer formation. Negative RT controls, which contained the complete RNA synthesis reaction components without the addition of the enzyme reverse transcriptase (RT) were used to rule out presence of genomic DNA (gDNA). Fold-change gene expression was measured relative to the reference gene, ActB, and differential gene expression was determined compared to null group (VEH/CON) using the Pfaffl method (Pfaffl 2001). Molecular work was carried out in adherence to MIQE guidelines (Bustin et al. 2009) (Supplementary data 2).

Enzyme immunoassays of oxytocin and arginine-8-vasopressin
Blood was collected by cardiac puncture and the plasma separated at 2000 × g centrifugation for 20 min at 4℃. Plasma levels of the neuropeptides OXT and arginine-8-vasopressin were quantified using commercially available ELISA kits from Arbor Assays (Ann Arbor, MI USA; OXT, K048-H1, Arg-8-Vasopressin, K049-C1) and Enzo Life Sciences (Farmingdale, NY, USA; OXT, ADI901153A0001; Arg8-Vasopressin, ADI-900-017) following the manufacturer's instructions. For the Arbor Assay kits, samples were first treated using the acetonebased extraction solution followed by vacuum lyophilization of the resulting supernatant at 37℃ to reduce the non-specific binding. For oxytocin, the colorimetric reaction product was read as optical density at 450 nm on a plate reader (SpectraMax 190, Molecular Devices). The kit has a sensitivity of 1.7 pg/sample in a dynamic range of 16.38-10,000 pg/mL. Arg8-Vasopressin was detected using a luminescence plate reader (Victor3, Perkin Elmer). The Arg8-vasopressin kit has a sensitivity of 0.9 pg/mL in a dynamic range of 1.638-1,000 pg/mL. For the Enzo Life Sciences kits, samples underwent solid phase extraction using 200 mg C18 Sep-pak columns as previously described (Deol et al. 2020). Plasma oxytocin and arginine vasopressin were quantified by interpolating absorbance or luminosity values, respectively, using a 4-parameterlogarithmic standard curve (MyAssays).

Statistical analyses
Statistical analysis was performed using GraphPad Prism (version 8.4.3 San Diego, CA, USA). Within group comparisons were performed using paired Student's t test. Between groups comparisons were accomplished using one-way, twoway or mixed model ANOVA with or without a repeated measures design. Non-parametric statistical tests (i.e., Kruskal-Wallis H test) were used when normality and/or equal variances assumptions were not met as measured using the Shapiro-Wilk and F-tests. If an equal variance assumption was not met, a Brown-Forsythe ANOVA or Welch's correction was used. Post hoc comparisons were performed using appropriate tests. Biological outliers were excluded when animals were unable to perform behavioral tests. Type 1 error rate (α) was set at 0.05; F and P values are presented in the figure legends or Supplemental statistical information.
The data are expressed as the mean ± s.e.m, in bar or scatter plots or as median and inter quartile range representing minimum and maximum values in whisker plots.

DE-71 dosing paradigm and maternal parameters
C57Bl/6 mice dams were exposed to DE-71 and later investigated along with the F1 female offspring as shown in diagram (Fig. 1). Using this dosing paradigm, we have previously reported no differences in litter size at birth, secondary sex ratio, nor gestational maternal parameters . Moreover, dams exposed to DE-71 did not build inferior nests and F1 litters at PND 46 had normal body mass relative to VEH/CON ( Supplementary  Fig. 1). In combination, these data indicate that perinatal DE-71 exposure does not interfere with pup health, maternal nest quality nor related behaviors shown to be affected by exposure to PCBs, a structural/functional analogue class of PBDEs (Abu-Arafeh et al. 2016).

PBDE congener analysis in offspring brain
PBDE congener content was determined using HRGC/ HRMS or GC/ECNI-MS in F1 female brain from offspring during the lactational period (PND 15) or as adults ( 14 congeners detected showed significantly elevated concentrations in DE-71 exposed offspring at PND 15 relative to VEH/CON (P < 0.05 and P < 0.01). Of note, BDE-153 was ~ tenfold enriched and BDE-47 was slightly depleted (~ twofold) relative to the DE-71 mixture as reported previously . By PND 110, the BDE composition in F1 brain was limited to BDE-153 (Fig. 1f), which was significantly elevated in L-DE-71 and H-DE-71 relative to VEH/CON (P < 0.01 and P < 0.05). In comparison, BDE-153 at ppb (and an additional 6 congeners) has been reported in postmortem brain samples from 4-71 year-old neurotypical controls and autistic humans born in 1940-2000 (Mitchell et al. 2012).

Early-life exposure to DE-71 induces deficits relevant to core symptoms of autism in F1 female progeny
Social novelty preference. Testing mice on a social novelty preference (SNP) test has been suggested to be ethologically relevant to social domains affected in autistics (Moy et al. 2004). On this test, all F1 exposure groups except the L-DE-71 F1 group (P < 0.05) showed a preference for the novel over familiar stimulus (Fig. 2a), and this was also represented in the recognition index vs VEH/CON (Fig. 2b, P < 0.05). In contrast, there was no effect of exposure in F0; all groups showed a preference for novel stimulus (Fig. 2c, P < 0.0001) and no group differences were observed in the recognition index (Fig. 2d). The investigation index for F1 and F0 groups approached 1 (Supplementary Fig. 2) indicating that the reduced exploration of novel over familiar shown by L-DE-71 F1 was not due to a decrease in total investigation time indicating no lack of participation.
Sociability. To determine social interest, an independent social cognition domain, we examined mouse behavior on a 3-chamber sociability test. All F1 groups (VEH/CON, L-DE-71, H-DE-71) showed preference for a novel social stimulus relative to a non-social novel stimulus as measured by sniffing time (Fig. 2e, P < 0.05, 0.01, 0.05, respectively), indicating normal sociability. Using chamber time VEH/ CON and L-DE-71, but not H-DE-71 F1, showed a preference for social stimulus (Fig. 2g, P < 0.05, P < 0.01, ns, respectively). Sniffing time has been suggested to have superior validity over chamber time scores since active behaviors that are most directly related to social investigation are captured (Fairless et al. 2011) and the physical proximity allows for transmission of volatile and nonvolatile oderants (Brennan and Kendrick 2006); (Luo et al. 2003); (Fairless et al. 2011). For F0, chamber (Fig. 2f, P < 0.01) and sniffing time scores inidcated a preference for social stimulus and no effect of exposure was found (Fig. 2h P < 0.05 for VEH/ CON and L-DE-71, and P < 0.01 for H-DE-71). As a measure of test robustness there was no indication of side preference during training for F1 and F0 ( Supplementary Fig. 2).
Repetitive behavior. On the marble burying test, which measures repetitive and perseverative behavior in rodents (Thomas et al. 2009), L-DE-71 (but not H-DE-71) adult F1 Fig. 2 Early-life exposure to DE-71 induces deficits relevant to core symptoms of ASD in F1 female progeny. a, c Social Novelty Preference scores for dams and female offspring: unlike the F1 VEH/CON and F1 H-DE-71, F1 L-DE-71 females failed to spend more time with a novel relative to a familiar conspecific stimulus. F0 dams exposed to DE-71 did not show abnormal social recognition relative to VEH/ CON. b, d Recognition Index scores show decreased preference for novel stimulus in L-DE-71 F1 relative to VEH/CON but not in F0. e, f Time spent sniffing in sociability test. All exposure groups spent significantly more time sniffing social stimulus indicating normal sociability. g, h Chamber time scores in sociability. All groups show significantly greater time spent in the social chamber relative to non-social except for F1 H-DE-71. i Marble burying scores showed offspring L-DE-71 buried a greater amount of marbles as compared to VEH/CON and H-DE-71, but not in dams. j Nestlet shredding was not affected in exposed F1 but was reduced in L-DE-71 F0 relative to corresponding VEH/CON. *P < .05, **P < .01; ****P < .0001 compared to VEH/CON (b,d,i,j), familiar (a,c) or non-social chamber (e,f,g,h). n = 6-11 litters/group (a-b), 19-26 subjects/group (c-d), 6-9 litters/group (e), 8 litters/group (f), 16-33 subjects/group (g), 13-24 subjects/group (h), 19-37 subjects/group for F1 and 11-16 subjects/group for F0 (i), 18-36 subjects/group for F1 and 16-19 subjects/group for F0 (j). F familiar, N novel, N non-social, S social, N non-social, E empty, S social buried a significantly greater number of marbles relative to VEH/CON (Fig. 2i, P < 0.05). A subgroup of F1 was tested at PND 30, but no group differences were detected, possibly indicating age-related physical hypoactivity, reduced habituation to the test arena or a latently emerging impact of PBDEs ( Supplementary Fig. 3). Again, no group differences were seen in F0 (Fig. 2i). Mean values for nestlet shredding were not affected by DE-71 exposure in F1. However, in F0, the L-DE-71 group showed a mean reduction in nestlet shredding relative to VEH/CON (Fig. 2j, P < 0.05). However, less nestlet shredding did not translate into poorer maternal nest scores (Supplementary Fig. 1).

Exposure to L-DE-71 but not H-DE-71 reduces long-term social recognition memory in F1 female progeny
We determined that SNP scores requiring a 30 min memory retention were abnormal in exposed F1 but not F0. To test the hypothesis that DE-71 compromises consolidation of long-term social recognition memory, we subjected F1 to a social recognition memory test (SRMT) (Tanimizu et al. 2017). On this test, mice with intact memory exhibit less time investigating a familiar juvenile conspecific 24 h after a first exposure. Figure 3a shows that VEH/CON and H-DE-71 mice were able to form a social recognition memory of the stimulus by Day 2 since they spent significantly less time with a familiar stimulus mouse (P < 0.05 and P < 0.0001, respectively). Mean RI values for F1 were no different for VEH/CON (0.71) or less for H-DE-71 (0.56, P < 0.05), suggesting no deficits in recognition memory (Fig. 3b). In contrast, L-DE-71 F1 showed an apparently greater RI (mean RI, 0.85, P = 0.07), indicating potentially deficient long-term social recognition memory (Fig. 3b). Next, we examined investigation time with a second novel stimulus mouse on Day 2 to determine whether the reduction of investigation time on Day 2 is specific to social memory formation and not due to disengagement. In Fig. 3c, no significant reduction of investigation time was noted as expected for VEH/CON and H-DE-71, suggesting that the reduction in investigation of familiar mouse above was specific to recognition memory formation. In contrast, L-DE-71 F1 exhibited a significant reduction in investigation time of new novel mouse on Day 2 (P < 0.05), indicating deficient social memory. In Fig. 3d L-DE-71 showed a significantly lower RI than 1 (0.89, P < 0.05). During test optimization using untreated controls, we confirmed that our test was valid, i.e., that the 3 but not 1 min of social exposure on Day 1 was sufficient to form a memory on Day 2 (P < 0.01) as reported (Tanimizu et al. 2017) (Supplementary Fig. 4). In summary, these results indicate that developmental exposure to DE-71 at 0.1 mg/ kg/d but not 0.4 mg/kg significantly reduces long-term social recognition memory in F1.

Perinatal exposure to L-DE-71 compromises short-term novel object recognition memory in F1 females and F0
Having found that DE-71 exposure produces significant impairment in the SNP and SRMT, we tested the representing reduction in investigation of familiar mouse on Day 2 relative to Day 1, indicate normal social recognition memory since mean scores were no different from 0.65 for VEH/CON and H-DE-71 groups. The mean RI value for L-DE-71 was apparently greater indicating that modest reduction of investigation of familiar mouse on Day 2 may be due to chance. c When using a different novel mouse on Day 2 (Novel`), VEH/CON and H-DE-71 F1 mice showed equal preference for novel mice on Day 2 relative to Day 1, indicating normal recognition memory and that the reduction in investigation time in a and b is specific to familiar juveniles. However, L-DE-71 mice were unable to do this. d Mean RI scores are less than 1 only for L-DE-71, indicating unequal preference for Novel` (Day 2) vs Novel stimulus (Day 1) and abnormal social recognition memory. *P < .05, ****P < .0001 compared to Day 1 (a, c). ^P < .05 compared to L-DE-71. a P < .05 compared to .65 (b) or 1.0 (d). n = 5-6 litters/group. '*', stimulus mouse in insets. D day hypothesis that DE-71 exposure also interferes with nonsocial recognition memory. Using a novel object recognition memory test (NORT), Fig. 4a shows that L-DE-71 F1 did not display preferential exploration of the novel object during the Day 1 testing, as did the VEH/CON and H-DE-71 (P < 0.01) indicating that the former did not recognize objects presented 30 min earlier in the familiarization phase. This was corroborated using a discrimination index which showed that values for only VEH/ CON and H-DE-71 were > 0, indicating memory for previously encountered objects (Fig. 4b, P < 0.05). In contrast, L-DE-71 group displayed a negative mean discrimination tive raster plots indicate no significant effect of exposure on general locomotor activity quantified as cumulative distance travelled and velocity. g-l After a 24 h retention time there was no effect of exposure on investigation time of familiar and novel, discrimination index, dwell-time maps, raster plots, distance travelled, or velocity. *P < .05, **P < .01, ***P < .001 compared to familiar object (a) or VEH/ CON (b). ^^^P < .001 compared to L-DE-71 (b). a P < .05, aa P < .01, aaa P < .001 compared to 0 (h). n = 4-10 subjects/group. F and F`, familiar object; N and N`, novel object. Scale bar, 5 cm index (greater preference for familiar object, P < 0.01), a result that was significantly reduced compared to VEH/ CON (P < 0.001). Representative dwell time maps in the open field arena showed preference for novel object (right corner) for VEH/CON and H-DE-71 on Day 1 (Fig. 4c). In contrast, L-DE-71 showed less exploration of novel relative to familiar object. On Day 2 all exposure groups preferred novel over familiar object and showed similar mean values for discrimination index and dwell times after a 24 h retention time (Fig. 4g, h, i, j). There were no effects of exposure on mean distance travelled and velocity (Fig. 4e, f, k, l) as indicated by raster plots (Fig. 4d, j). Interestingly, both L-DE-71 and H-DE-71 exposed dams showed similar short-term memory deficits as L-DE-71 exposed F1 offspring ( Supplementary Fig. 5).

Abnormal social behavior in F1 produced by DE-71 exposure is not due to deficits in general olfactory processing
To examine if DE-71-induced deficits observed in social recognition ability were due to insufficient olfactory ability, we subjected female offspring to an olfactory preference test. Figure 5a shows that all mice including those treated with DE-71 displayed increased odor sniffing duration for peanut butter over water (P < 0.05-0.0001), butyric acid (P < 0.05-0.0001), and vanilla (P < 0.05-0.001). Similar results were obtained for dams (Fig. 5b). These results indicate that, like VEH/CON, DE-71 exposed offspring and dams were able to process sensory signals from different non-social odors with enough sensitivity to show preference for peanut butter over others. No group differences were noted for F0. *P < .05, ****P < .0001 compared to water. ^P < .05, ^^^P < .001 compared to vanilla; § §P < .01, § § §P < .001, § § § §P < .0001 compared to butyric acid. a P < .05, aaa P < .001 compared to VEH/CON during habituation. b P < .05, bb P < .01, bbbb P < .0001 compared to VEH/CON during dishabituation. Additional statistical results are summarized in Table 2. n = 6-15 litters/group (a), n = 11-16 subjects/group (b), n = 12-16 subjects/group (c), n = 12-16 subjects/group (d). W water, B butyric acid, P peanut butter, V vanilla

DE-71 exposure alters olfactory discrimination of social odors
We used an olfactory habituation/dishabituation test to measure olfactory discrimination. Table 2 indicates the results of the habituation/dishabituation test. The F1 VEH/CON group displayed olfactory habituation to all non-social odors and social odors (except non-social odor 2-banana) as indicated by the decline in time spent sniffing by trial 3 (Fig. 5c). F1 VEH/CON displayed olfactory dishabituation when transitioning to a new odor except from non-social 2 (banana) to social 1 (P < 0.01, P < 0.0001). Both DE-71 groups displayed deficient habituation and/or dishabituation for more than 1 odor ( Table 2). In particular, L-DE-71 showed reduced habituation to social odor 1 (from trial 1 to 2; P < 0.05; Fig. 5c). Figure 5c shows that compared to VEH/CON, L-DE-71 and H-DE-71 showed less dishabituation from social odor 1 to 2 (P < 0.01, P < 0.0001), suggesting that DE-71 produces reduced olfactory discrimination (hyposmia) especially of social odors, which requires processing via MOE and VNO (Huckins et al. 2013).
In addition, H-DE-71 also showed reduced dishabituation from non-social odor 2 to social odor 1 (P < 0.05). An apparently significant effect was also seen for L-DE-71 (P = 0.07).
These findings indicate altered social odor discrimination after perinatal exposure to DE-71 potentially associated with altered signaling through the VNO and less associated with MOE which also processes non-social odors which are normal. Olfactory discrimination of odors in F0 showed normal habituation/dishabituation profiles compared to VEH/CON (Fig. 5d, Table 2). There were no exposure group differences found for F0.

DE-71 exposure does not promote anxiety nor depressive-like behavior
Mice were evaluated for anxiety using the EPM test and time spent in closed arms relative to open arms was significantly greater in all exposure groups in F1 and F0 (P < 0.0001).
There was no effect of exposure on the number of total arm entries for F1. In contrast, the H-DE-71 F0 group exhibited significantly fewer total entries relative to VEH/CON (Supplementary Fig. 6). Using a forced swim test, depressive-like behavior was measured as time spent immobile and there was no significant effect of exposure on time spent immobile for F1 nor F0 as compared to VEH/CON ( Supplementary  Fig. 6).

Selective effects of DE-71 exposure on Suok test
Using Suok, we measured the effects of DE-71 exposure on locomotion, exploratory behavior, sensorimotor coordination and anxiety. Relative to VEH/CON, H-DE-71 (but not L-DE-71) F1 showed decreased horizontal activity, as represented by segments crossed (Fig. 6a), decreased locomotion (Fig. 6b), decreased exploratory activity (Fig. 6g), increased SAP (Fig. 6i) and decreased grooming (Fig. 6k). Falls were significantly decreased in H-DE-71, but not when normalized to segments crossed (Fig. 6c, d). In contrast to F1, F0 exposed to H-DE-71 showed decreased hind leg slips (Fig. 6e) and L-DE-71 showed increased SAP relative to VEH/CON (Fig. 6i). There were no significant differences on the other measures.

Early-life PBDE exposure does not alter locomotion on the open field test
The open field test informs about locomotion, habituation to novelty and anxiety. All F1 exposure groups showed similar reduced exploratory activity over time (habituation), measured as reduced distance traveled and velocity over the 1 h test (Fig. 7 a,b, P < 0.0001). Between-group comparisons showed no effect of exposure for F1 when compared to VEH/CON. These results helped us rule out concerns of hyper-or hypomobility in DE-71 exposed female offspring relative to VEH/ CON as reported after acute exposure to 0.8 mg/kg BDE-99 at PND 10 (Viberg et al. 2004); (Costa and Giordano 2007).
Other studies using chronic exposure of mouse dams to low doses of BDE-47 (0.1 mg/kg) or -99 (0.6 mg/kg) from gestation through third week of lactation have shown inconsistent results with both hypoactivity and no effect reported on the open field test in female offspring (Ta et al. 2011); (Koenig et al. 2012); (Branchi et al. 2002). Exploration time in center and periphery zones for all exposure groups (Fig. 7c,d) showed habituation only in the periphery (P < 0.0001). Figure 7e shows that total distance travelled in the periphery was similarly and significantly greater than in center in all exposure groups (P < 0.0001). Another measure of anxiety, number of fecal boli at 5 and 10 min into the test, indicated increased emotional reactivity in L-DE-71 F1 relative to VEH/CON (Fig. 7f, P < 0.05, P < 0.01, respectively). For F0 H-DE-71 showed greater distance travelled and velocity in the arena and exploration in the periphery zone as compared to VEH/CON (P < 0.05 and P < 0.01, respectively) and L-DE-71 (P < 0.05 and P < 0.0001). L-DE-71 mice produced more fecal boli at 60 min relative to VEH/CON (Supplementary Fig. 7).

DE-71 exposure alters prosocial gene expression in brain regions involved in social behavior
To correlate the behavioral findings with changes in gene expression of the social neuropeptide systems that are key mediators of complex social behavior, such as vasopressin (Avp), oxytocin (Oxt), PACAP (Adcyap1) and their receptors, we measured the relative expression of their genes from micropunches of discrete brain nuclei involved in social behavior: lateral septum, amygdala, bed nucleus of the stria terminalis (BNST), SON and PVN.

Perinatal exposure to DE-71 exaggerates plasma levels of arginine-8-vasopressin but not oxytocin levels in adult F1 female offspring
We next measured plasma oxytocin and Arg8-vasopressin concentrations to examine their association with social behavior phenotypes. Figure 9 shows that plasma AVP levels in L-DE-71 F1 females were significantly elevated relative to VEH/CON (P < 0.05). In contrast, there were no group differences in plasma OXT levels.

Discussion
Growing evidence suggests a positive association between early-life exposure to PBDEs and neurodevelopmental alterations (Branchi et al. 2003). Environmental factors, including xenobiotic chemical exposures, may provide a plausible explanation for the rising incidence of NDDs characterized by social deficits (Maenner et al. 2020); however, experimental evidence has not established a direct link with specific candidate chemicals. With this purpose, our study is the first to investigate the effects of the penta-PBDE mixture DE-71 on behaviors and neurochemical/endocrine profiles relevant to several core ASD symptom domains concurrently. Our experimental design exposes progeny to the full complement of congeners found in human breast milk (Lind et al. 2003). The major findings reveal that developmental DE-71 exposure produces durable deficits in social recognition, repetitive behavior and social odor discrimination in exposed female offspring. The behavioral phenotypes occurred concomitantly with changes in plasma AVP and neuromolecular markers of Oxt and Avp signaling pathways Female offspring and dams were tested on SUOK for: a, b locomotion; c-f sensorimotor coordination; g exploratory activity; h-j anxiety behaviors; and k autogrooming. Only H-DE-71 F1 showed decreased mean values in a, b, g, k and increased i whereas F0 exposed to L-DE-71 showed increased mean value in i. *P < .05, **P < .01 compared to corresponding VEH/CON. ^P < .05, ^^ P < .01 compared to corresponding L-DE-71. n for F1 (litters/group): a 10-11; b 10; c 9-11; d 10-11; e 11-12; f 10-12; g 10-12; h 10-14; i 7-8; j 9-11; k 10-12. Distance traveled in the open field arena. All F1 exposure groups showed similar reduced exploratory activity and velocity over the 1 h. c, d Exploration time in periphery and center for all groups showed habituation only in the periphery. e Exploration time in center was significantly less than in periphery for all groups, suggesting no exposure effects on anxiety. f Another measure of anxiety, number of fecal boli, indicated increased emotional reactivity in the L-DE-71 F1 relative to VEH/CON. *P < .05, **P < .01. ****P < .0001 compared to center (e) or VEH/CON (f). ^P < .05 compared to corresponding L-DE-71. a P < .0001 compared to initial time bin for corresponding treatment group. n = 19-23 subjects/group. C center zone, P periphery zone Heatmap representation (double gradient, blue-minus; redplus) of RT-qPCR analysis with the respective fold-change value (mean) of each gene studied by brain region. n = 4-17/group. *P < .05 compared to VEH/ CON. ^P < .05 compared to L-DE-71. BNST, bed nucleus of the stria terminalis; AMG amygdala, LS lateral septum, SON supraoptic nucleus, PVN paraventricular nucleus 1 3 in brain regions that coordinate complex social behaviors. Together, the behavioral, sensory and neurochemical phenotypes produced by DE-71 may provide a novel, comprehensive ASD-relevant model with high translational impact. Our results are congruent with a disrupted developmental trajectory of the social processing domain as outlined in the 2010 NIMH Research Domain Criteria (RDoC) framework (Social Processes: Workshop Proceedings, 2012); a key characteristic of ASD pathology. Our work is further strengthened by the use of the litter as the unit of statistical analysis, thus minimizing risk of bias (RoB) of individual studies (Dorman et al. 2018) and inter-individual variability. Further, DE-71 produced the common hormetic response, such that the 0.1 but not 0.4 mg/kg group exhibited most of the behavior changes even though there was a dose-dependent increase in brain accumulation of ∑PBDE congeners. Moreover, we confirmed the augmented susceptibility to developmental relative to adult exposure to PBDEs, highlighting the significance of chemical exposures during critical neurodevelopmental windows. Collectively, these data support the conclusion that environmental xenobiotics impact social behavior and related neurochemical signaling pathways in mice relevant to NDDs.

Perinatal DE-71 exposure produces deficient social recognition and increases repetitive behavior in adult female offspring
Our main finding was that in utero and lactational transfer of DE-71 from exposed mothers produces behavioral phenotypes resembling two core behavioral features of ASD DSM-V diagnosis: deficits in social reciprocity and communication and repetitive/stereotyped behaviors (American Psychiatric Association 2013). With respect to the latter, female offspring exposed to L-DE-71 showed increased activity on a marble burying test indicative of repetitive behaviors in rodent models of ASD (Angoa-Pérez et al. 2013). At 0.1 mg/kg, developmental DE-71 exposure also produced deficient short-term social memory (SNP) and long-term social recognition memory (SRM), while sociability (SOC) was not affected, ruling out a lack of the 'social motivation' component of social cognition as an underlying contributor. SRM is considered to be another distinct behavioral domain and important for the 'knowledge of self and others' component of social cognition (Bicks et al. 2020).
Though much is still unknown about the neural correlates of social behavior, the social motivation and social recognition domains have been shown to be independent of each other. For instance, deficits in SNP can occur without decrements in sociability in other models of deficient social behavior induced by high fat diet (Hayashi et al. 2020) or C-section delivery (Morais et al. 2020). The former can be restored by OXT administration. In other reports, restoring OXT content in the PVN with probiotic therapy (L. reuteri), in maternal high fat diet and valproic acid offspring, rescues SOC and SNP, but not other ASD endophenotypes (Sgritta et al. 2019); (Buffington et al. 2016). Taken together, these results suggest that different mechanisms and/or circuitry govern the various social behavior domains that can be selectively isolated by experimental contrast and susceptibility to early-life PBDEs. Specifically, perinatal DE-71 exposure significantly compromises the social recognition domain of social cognition, which is more relevant to ASD since the behaviors related to knowledge of self and others such as facial recognition, empathy and evaluation of emotion of others are disrupted in ASD patients (Bicks et al. 2020); (Bradshaw et al. 2011).
Our findings indicate deficient short-term social recognition and long-term social recognition memory in L-DE-71 F1, suggesting our results may be translational to ASD and other NDDs characterized by psycho-social deficiencies. While findings of epidemiological studies evaluating associations between PBDEs and social deficits/ASD are mixed (Vuong et al. 2018); (Braun et al. 2014); (Gibson et al. 2018), a higher risk of poor social competence has been found with increasing postnatal exposure to BDE-47 (4 yr old child serum; (Gascon et al. 2011). BDE-47 levels in cord blood have also been positively associated with poor Fig. 9 Perinatal exposure to DE-71 exaggerates plasma levels of Arg8-vasopressin but not oxytocin (OXT) in adult F1 female offspring. a Plasma Arg8-vasopressin measured using EIA using blood taken at sacrifice. L-DE-71 exposed offspring showed elevated levels. b OXT levels showed no exposure effects. *P < .05 compared to VEH/CON. n = 8-13 subjects/ group (a); n = 6-8 subjects/ group (b) social domain development in 24-month-old toddlers (Ding et al. 2015). Previous rodent studies examining the effects of environmental pollutants on social behavior have produced inconsistent results perhaps due to heterogeneity of brominated (BFR) flame retardants used, timing of exposure, sex and/or model organism used. Importantly, the only other study examining the effects of DE-71 (0.3 and 1.6 ppm) on social behavior supports our findings. Fernie and colleagues (2005) found fewer and less appropriate pair-bonding and courtship behaviors in exposed captive kestrels (Fernie et al. 2005). In contrast, female mice offspring exposed to BDE-47 perinatally via the mother showed reduced sociability relative to controls ) but no effect of BDE-47 (at 0.03 mg/kg) was detected on SNP unless administered to genetically altered mice lacking methyl-CPG binding protein 2 (Mecp2), a frontal cortical protein negatively associated with ASD (Nagarajan et al. 2006). Also in contrast to our findings, male CD-1 mice developmentally exposed to BDE-47 (0.2 mg/kg) display reduced time with conspecfics but show no effect on SNP relative to controls (Kim et al. 2015). Moreover, using rats perinatally exposed to high doses of BDE-47 (50 mg/kg) Li and others (2021) report normal preference for stranger over familiar conspecific and for social stimulus over empty corral but with reduced time spent in exploration (Li et al. 2021). Using low doses of BDE-209 (0.12 ng/mouse/day, s.c.), Chen and colleagues (2019) did not observe deficient sociability nor SNP in exposed male mice offspring (Chen et al. 2019). Therefore, it appears that PBDE effects on social behavior may be congener-and/or dose-specific. Our findings are supported by perinatal exposure studies using another flame retardant mixture Firemaster 550 (6.6 mg/kg/day) and its BFR and organophosphate components (3.3 mg/kg/d), which produce deficits in social recognition after 24 h retention in a sexand exposure-specific manner in rats (Witchey et al. 2020). Perinatal exposure to Firemaster 550 also produces abnormal partner preference in female prairie voles (1 mg/kg; Gillera et al. 2020). Other studies have reported adverse effects of PCBs (A1221) on mate preference and social behavior (Reilly et al. 2015); (Hernandez Scudder et al. 2020).

Perinatal DE-71 exposure produces deficient novel object recognition memory in dams and adult female offspring
A complex interplay between forebrain regions is responsible for normal social recognition (Bicks et al. 2020); (Tanimizu et al. 2017); (Ferguson et al. 2000) including hippocampal circuits underlying social memory formation and amygdalar circuits that process social signals such as volatile odorant pheromones that trigger social and reproductive behaviors (Kogan et al. 2000); (Noack et al. 2010). DE-71-exposed socially deficient mice also showed abnormal NOR memory suggesting abnormal function in hippocampus since the latter serves as an integration hub underlying both social and recognition memory (Hernandez et al. 2008, Pfaffl 2001, Pinson et al. 2016. Toxicological studies of developmentally administered single BDE congeners or DE-71 have not examined effects on NOR or SRM (Dorman et al. 2018). However, previous studies using peri/postnatally administered single BDE congeners such as BDE-153 (0.9 mg/kg bw) and -47 (0.03 mg/kg bw) have shown neurotoxic actions on hippocampal-dependent function related to spatial memory (Koenig et al. 2012); (Dorman et al. 2018); (Viberg et al. 2003). In support of our findings, evidence from human studies suggest that more than one environmental BDE congener may produce risk for cognitive impairments in children. For example, several PBDEs found in maternal samples 99,100,153) are associated with children's lowered IQ and cognitive scores (Herbstman et al. 2010); (Lam et al. 2017); (Azar et al. 2021), mental/ physical development (Eskenazi et al. 2013) and fine motor skills, attention and cognition (Chen et al. 2014).
It is not surprising that L-DE-71 F1 mice showed coincident deficient NOR memory and SRM. However, while deficient SRM was seen at both short and long-term retention times, NOR memory deficits were evident with short-term retention time only. Moreover, F0 showed deficits only in short-term NOR memory indicating that short-term social recognition ability and short-term novel object recognition memory are distinct constructs. Therefore, PBDEs may target different brain circuits participating in general and social memory processes and/or different neurochemical systems within each circuit. For example, hippocampal OXTRs are necessary for short-term social recognition but not novel object recognition memory in male mice (Raam et al. 2017).

Perinatal exposure to DE-71 alters social odor discrimination in adult female offspring
Recognition of conspecifics in rodents depends on proper identification, discrimination and processing of olfactory cues present in urine and secretions from skin, reproductive tract and scent glands (Kogan et al. 2000); (Noack et al. 2010). We found that the disruption of social behavior after perinatal DE-71 exposure is coincident with abnormal profiles of olfactory habituation/dishabituation to social odors. For example, socially deficient L-DE-71 mice also displayed a reduction in habituation to social odors and in dishabituation from one to another social odor. In H-DE-71 mice the deficits in olfactory discrimination were relatively less severe as was that of their social behavior (reduced sociability with normal social recognition ability and memory). In combination, these results suggest that DE-71 effects on olfactory discrimination are specific to social odors, since exposure produced no deficits in general olfactory processing, and that deficiency was coincident with social deficits in a dose-dependent manner. It is unclear why PBDEs are more neurotoxic to social odor processing but it may depend on differential targeting of CNS pathways taken by signals from neutral and social odors. Chemosensory cues are processed through two olfactory systems; neutral odors (banana and almond) are processed through the main olfactory epithelium (MOE) and social odors through both MOE and the vomeronasal organ (VNO) (Huckins et al. 2013). Signals are then processed through amygdala and hypothalamus to trigger innate social and reproductive behaviors.
There are no previous studies on PBDEs and olfactory function although  can concentrate in the epithelium of the nasal cavity (Darnerud and Risberg 2006) and developmental exposure to BDE-209 impairs subventricular zone (SVZ) neurogenesis and olfactory granule cell morphology in mice (Xu et al. 2018). However, a recent study using PCBs (Aroclor 1221, 1 mg/kg) indicates that prenatal exposure impairs mate preference behavior based on olfactory cues concomitant with impaired odor preference for mates with different hormone status in adult female offspring (Hernandez Scudder et al. 2020). Our findings that early-life PBDE exposure alters social odor discrimination may translate to autistic humans which are prone to hypo-or hyper-reactivity to sensory stimulation (American Psychiatric Association 2013); recent studies suggest that this may include atypical olfaction (Martin and Daniel 2014). Indeed, several olfactory outcomes have been reported in children with ASD, i.e., abnormal odor responses, difficulties in emotional reaction to odors, impaired detection thresholds and odor identification as well as heightened olfactory sensitivity (Martin and Daniel 2014); (Rogers et al. 2003); (Legiša et al. 2013). Further research is needed to discern the mechanisms by which PBDEs may act to alter social odor discrimination and if this underlies, in part, their social recognition deficits. Interestingly, extrahippocampal OXT and AVP systems that contribute to short-term social recognition, also modulate detection and processing of social odors (Wacker and Ludwig 2012); (Oettl et al. 2016).

Perinatal DE-71 Alters AVP and OXTergic neuromolecular phenotypes in brain regions that coordinate complex social behaviors
Our lab has previously shown that in vitro and early-life exposure to PBDEs (and PCBs) produce neuroendocrine disruption of the prosocial neuropeptide, vasopressin, under physiologically stimulated conditions (Coburn et al. 2005); (Coburn et al. 2007); (Mucio-Ramírez et al. 2017). Therefore, the observed PBDE-induced deficits in SNP and SRM may result from altered prosocial functions of AVP and/or the structurally related neuropeptide, OXT. Here we show that L-DE-71 downregulates Avp in BNST, which provides sexually dimorphic AVPergic innervation to LS (Bychowski et al. 2013). Diminished AVPergic signaling to LS may explain reduced social recognition memory in L-DE-71 F1 females, since local AVP1a receptor antagonism compromises social discrimination especially well in females (Veenema et al. 2012). DE-71-mediated upregulation in BNST Avp1ar may represent a compensatory effect to maintain signaling at normal levels. Interestingly, Avp1ar in the ventromedial nucleus (VMN) is upregulated by the PCB mixture A1221 in female rat (but not male) offspring and is not dependent on estrogenic pathways (Topper et al. 2019). The observed downregulation of Avp in SON may also impact social recognition ability indirectly via reduced AVPergic-mediated activation of BNST (Lukas et al. 2011).
At 0.1 mg/kg, DE-71 also produced elevated plasma AVP which is consistent with less inhibitory regulation over axonal secretion of AVP hormone resulting from potentially reduced levels of central AVP (Ludwig et al. 2016). Additionally, DE-71 can increase exocytosis in PC12 pheochromocytoma endocrine cells (Dingemans et al. 2008) possibly by interfering with intracellular calcium dynamics (Coburn et al. 2008) and potentially increasing secretion of stored AVP depots in axonal terminals located in the posterior pituitary and, thus, releasing AVP into the bloodstream. We also found that DE-71 alters the central OXTergic system which is also necessary for social recognition and partner preference (Ferguson et al. 2000). For example, mice with OXT gene deletion fail to remember recently encountered individuals and do not show the typical decline in preference during subsequent exposures to a familiar mouse, an effect which can be rescued by central administration of OXT (Dluzen et al. 1998). Here we show that in the BNST, L-DE-71 female F1 display significantly reduced Oxt mRNA transcripts. Assuming that there is a positive correlation between gene and peptide content and release, one interpretation of our data is that L-DE-71 exposure reduces central OXTergic signaling which is necessary for normal social discrimination (Dumais et al. 2016). Our results further indicate that BNST-originating OXT may be sufficiently important for activating BNST OXTR relative to PVN-originating OXT (Knobloch et al. 2012). Indeed, a recent report has demonstrated that OXT receptor blockade, in the extrahypothalamic population of oxytocinergic neurons of the BNST, impairs social recognition in female and male rats (Dumais et al. 2016). OXT in the BNST also drives stress-induced social vigilance and avoidance that may be at play in social behavior domains examined here 137 . Ultimately, OXTergic signaling in the BNST may not be sufficient to affect social recognition since H-DE-71 F1 females also displayed reduced BNST Oxt mRNA. Importantly, L-DE-71 also reduced Oxt transcripts in the SON. Because local release of OXT from SON dendrites that extend to MeA promotes social recognition through amygdalar OXTRs (Takayanagi 1 3 et al. 2017); (Abramova et al. 2020), downregulated Oxt in SON may underlie, in part, the associated SNP and SRM deficits. Since the promoter regions for genes of both oxytocin (Mamrut et al. 2013) and vasopressin systems (Auger et al. 2011) are susceptible to epigenetic modification (Auger et al. 2011), these genes may be altered by global DNA methylation that is possible after developmental BDE-47 exposure ; (Poston and Saha 2019). Our findings may have translational value since altered OXT and AVP mechanisms in humans have been implicated in ASD and other psychiatric and neurological disorders (Zhang et al. 2016); (Kobylinska et al. 2019); (Hendaus et al. 2019); (Oztan et al. 2018).
Growing evidence from a myriad of studies using different model organisms show that brain vasopressin and oxytocin systems may be targets of other EDCs. EDC chemicals include phytoestrogens like genestein, chlorpyrifos, diethystilbestrol, bisphenol A (BPA), organophosphate pesticides and organohalogens like PCBs, PBDEs and Firemaster 500 (Patisaul 2017). For example, we and others have shown that some of these EDCs such as Firemaster 550 (Gillera et al. 2021), PCBs such as A122 (Reilly et al. 2021) and A1254 (Coburn et al. 2005); (Coburn et al. 2007) and BPA (Witchey et al. 2020) have endocrine disruptive effects on the oxytocin and/or vasopressin systems and, in some cases, these are concomitant with their effects on social behavior (Witchey et al. 2020); (Gillera et al. 2020); (Hernandez Scudder et al. 2020); (Reilly et al. 2021).

Specificity and comprehensive profile of PBDE toxicant model of ASD
A recent meta-review has put forth recommendations to improve ASD model characterization in rodent studies such that information about reciprocal social communication and stereotyped repetitive behavior domains are characterized in the same animals (Pelch et al. 2019). To this end, we used established protocols to measure ASD relevant and other comorbid behaviors to fully characterize the DE-71-induced phenotypes (Moy et al. 2004); (Crawley 2012). We found that the effects of DE-71 were specific to social novelty preference and social recognition memory as well as repetitive behavior and olfactory discrimination of social odors. Alterations were specific to offspring exposed perinatally via maternal transfer of environmentally relevant BDE congeners; adult exposed mothers were mostly unaffected. DE-71 had little to no effects on behaviors representing the domains of anxiety, depression and locomotion indicating ASD-relevant specificity without general neurological effects. In addition, there were no indications of reduced general health, i.e., body weight in pups nor gross abnormalities in maternal nest conditions. We have recently reported that similarly exposed (L-DE-71) female offspring, and to a lesser degree, their exposed mothers, display diabetic symptomatology, effects which may relate to the present findings . Importantly, we used multiple behavioral tests to validate social and other constructs studied (locomotion and anxiety). For example, for all F1 groups, the frequency of total entries on EPM and distance travelled on OFT yielded similar results regarding locomotion. In addition, time spent in open arm on EPM, and latency to leave center on Suok was consistent with duration in center of OFT.
Our DE-71 model of ASD also shows altered prosocial peptide neurotransmitters/neurohormones that are critical to ensuring proper development of social brain networks. In particular, the vasopressin and oxytocin systems are critically involved in social cognition with mutations having socio-behavioral impact that have been implicated in core symptoms of autism (Frank and Landgraf 2008). These neurochemical systems are being actively studied as potential targets of future therapeutic interventions for ASD (Meyer-Lindenberg et al. 2011); (Bolognani et al. 2019). In light of incongruent findings reported by past rodent and human studies (Pelch et al. 2019), we believe that our findings brings us closer to understanding the risk of ASD posed by xenobiotic endocrine disrupting chemicals. Nevertheless, human and rodent studies reporting on the relationship between PBDE exposure and autistic phenotype are few in number and have yielded inconclusive results and this field would benefit from additional detailed epidemiological and animal studies on the relationship between persistent organic pollutants (POPs) and risk of ASD.

Maternal transfer of BDE congeners in DE-71
and their brain accumulation in female offspring is dose-and time-dependent BDE congener composition found in PND15 exposed brains mimics that found in humans. BDE-28, -47, -99, -100, -153 were common congeners found at ppb in both DE-71-exposed offspring groups at PND15 with three-fold greater levels in H-DE-71 than L-DE-71. ∑PBDE values for adult serum are 30-100 ng/g lipid (Costa and Giordano 2007) and 3 to ninefold higher in infants because of exposure through breastmilk and in toddlers because of exposure through house dust and the diet (Toms et al. 2008); (Rose et al. 2010); (Schecter et al. 2005); (Fischer et al. 2006). Serum ∑PBDE values can reach 482 ng/g l.w. in toddlers (California 18 month-old) (Fischer et al. 2006) but lesser values have also been reported, i.e., 127 ng/g l.w. (Vuong et al. 2017) in Ohio 2-year-olds and 100 ng/g l.w. for North Carolina 12-36 month-old toddlers (Stapleton et al. 2012). Using a divisor factor of 0.095 to convert w.w. to l.w. (unpublished observations), we estimate our mean ∑PBDE in L-DE-71 F1 at PND 15 to be 1.7-to 8.2-fold greater, suggesting ours represents a translational model of maternal PBDE transfer. The main congeners in PND15 brains, BDE-47, -85, 1 3 -99, -100, -153, and -154, accounting for 97% of the mean ∑PBDEs, also comprise the majority of congeners (96%) in DE-71 ). These and other congeners found in offspring brain samples, i.e., 28,49,138,139,140,183 and 184, have also been detected in human serum and/or breastmilk (Chao et al. 2010). Importantly, to our knowledge,  have not been previously detected in DE-71 exposed rodent brain .  have reported brain levels of PBDEs in adult offspring after exposing dams during pregnancy and lactation to high levels of DE-71 (30 mg/ kg) ). Other papers have reported brain BDE levels after maternal transfer or directly exposing postnatal or adult mice with BDE congeners (Darnerud and Risberg 2006); (Staskal et al. 2006). Our study is unique because we used environmentally relevant congeners and doses (0.1 and 0.4 mg/kg) in a maternal transfer model with translational value. We also examined penetration of congener in exposed offspring during both early postnatal development and as adults.
Although using a mixture like DE-71 closely models the PBDE contamination previously shown in human breastmilk, there are some congeners, found at low levels in breastmilk, that we did not detect in offspring brain, i.e., BDE-7, 15, 71, 77, 119, 126 (Zhao et al. 2013); (Chao et al. 2010). Of these BDE-71 and -126 are present in DE-71 (LaA Guardia et al. 2006). Little or no information is available about the penetrance and/or neuroactivity of the missing congeners. Most rodent studies, which have focused on a single PBDE congener, BDE-47, dominantly detected in humans, have not reported pervasive effects on social behavior as we do here using DE-71. We speculate that BDE-47 alone is not effective in producing deficits in social recognition ability and social memory and that, instead, several PBDE congeners may act synergistically and/or additively to generate these abnormal phenotypes, reinforcing the need for in vivo studies using PBDE formulations that mimic child exposure.
By PND 110, the BDE composition in F1 brain was limited to BDE-153. Minimal metabolism of this congener is observed in rodents due to its high lipophilicity as determined by a high octanol-water partition coefficient (Log K ow ) (Sanders et al. 2006). BDE-153 has been positively associated with lower IQ in children and can cause impaired learning and memory in animal studies (Viberg et al. 2003); (Azar et al. 2021). However, while BDE-153 (and an additional 6 congeners) are detected at ppb in postmortem brain samples from 4-71 year-old born 1940 to 2000, it seems to be significantly depleted in autistics relative to normal subjects (Mitchell et al. 2012) for unknown reasons although this findings is based on a small sample size. The relatively lower retention of BDE-47 is in line with a previous report of differential tissue accumulation and disposition of BDE congeners attributed to their toxicokinetic properties (Staskal et al. 2006). Cyp-mediated biotransformation of BDE-47 and -99 (but not 153) may contribute since these congeners contain sites with adjacent unsubstituted carbons where the metabolism occurs (Sanders et al. 2006). Our findings suggest that elevated brain levels of BDE congeners in DE-71 during the critical window of early postnatal neurodevelopment may predispose children to neurobehavioral alterations related to ASD.

Conclusion
Though the role of environmental toxicants in the etiology of NDDs is poorly understood, our data support a link between maternal toxicant exposures and abnormal social and repetitive behavior in offspring that is relevant to ASD. We have shown that early-life exposure to DE-71 leading to these phenotypes is associated with humanrelevant levels and composition of BDE congeners penetrating the postnatal offspring brain via maternal transfer. DE-71 exposure has prominent actions if it occurs during perinatal development as compared to adulthood, supporting previous studies showing the particular susceptibility of developing nervous system to neurotoxic actions of PBDEs. The abnormal social behavior phenotypes produced by DE-71 are mostly specific to social novelty preference and social recognition memory and are also associated with excessive repetitive behavior, as well as neurochemical and social odor processing correlates-suggesting that discrete brain functional systems are targeted by PBDEs to promote neurodevelopmental abnormalities. Future studies are needed to discern if DE-71 actions are sexually dimorphic and extend to exposed male offspring. We believe that our environmental toxicant mouse model has utility in future studies examining the relationship between environmental xenobiotics, neurodevelopmental reprogramming and the rising incidence of NDDs.

Limitations of the study
The results of the PCR analysis provide novel results on the effects of PBDE exposure on the expression of gene markers for 'prosocial' neuropeptides and their receptors in specific regions of the social brain network. However, these restricted regions vary in cell density and limit the RNA yield for genes of interest (GOIs), especially in the amygdala and LS. Moreover, relative expression was more variable for ROIs that have low expression of GOIs, i.e., Oxt for LS. To improve our experimental data, we followed MIQE guidelines to optimize oligonucleotide primer efficiency and target specificity. Nevertheless, since the methodological approach we outlined depends on the level and variability of gene expression and quantity of RNA collected, our results should be interpreted alongside these limitations. BDE congener analysis was performed using two mass spectrometry methods utilized by teams at different institutions. The GC/ECNI-MS method uses an ECNI ionization mode to improve sensitivity. This method provides equal sensitivity to that provided by HRGC/ HRMS that uses electron impact ionization. Therefore, the reduction in brain BDE congeners at PND110 is likely due to elimination after accumulation measured at PND 15 and not to methodological factors. The great majority of behavioral tests were analyzed using litter as the unit of statistical analysis. However, for practical reasons most other tests used individual subjects. Our findings pertain to exposed female offspring and their mothers but male offspring were omitted due to limited resources. Further research is needed to determine if the ASD phenotypes evoked using the PBDE model are sex-specific.
1 3 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:// creat iveco mmons. org/ licen ses/ by/4. 0/.