Multigenerational effects of bisphenol A or ethinyl estradiol exposure on F2 California mice (Peromyscus californicus) pup vocalizations
Multigenerational effects of bisphenol A or ethinyl estradiol exposure on F2 California mice (Peromyscus californicus) pup vocalizations
Sarah A. Johnson 0 1 2 3 4
Michelle J. Farrington 1 2 3
Claire R. Murphy 1 2 3
Paul D. Caldo 1 2 3
Leif A. McAllister 1 2 3
Sarabjit Kaur 1 2 3
Catherine Chun 1 2 3
Madison T. Ortega 1 2 3
Brittney L. Marshall 1 2 3
Frauke Hoffmann 3
Mark R. Ellersieck 3
A. Katrin Schenk 3
Cheryl S. Rosenfeld 1 2 3
0 Department of Animal Sciences, University of Missouri , Columbia, Missouri , United States of America
1 Department of Biomedical Sciences, University of Missouri , Columbia, Missouri , United States of America
2 Bond Life Sciences Center, University of Missouri , Columbia, Missouri , United States of America
3 Editor: Jon T Sakata, McGill University , CANADA
4 Department of Gastroenterology, School of Medicine, University of Missouri , Columbia , Missouri, United States of America, 5 Department of Chemicals and Product Safety, The German Federal Institute for Risk Assessment , Berlin, Germany , 6 Department of Agriculture Experimental Station-Statistics, University of Missouri , Columbia , Missouri, United States of America, 7 Department of Physics, Randolph College, Lynchburg, Virginia, United States of America, 8 Genetics Area Program, University of Missouri , Columbia , Missouri, United States of America, 9 Thompson Center for Autism and Neurobehavioral Disorders, University of Missouri , Columbia, Missouri , United States of America
Rodent pups use vocalizations to communicate with one or both parents in biparental species, such as California mice (Peromyscus californicus). Previous studies have shown California mice developmentally exposed to endocrine disrupting chemicals, bisphenol A (BPA) or ethinyl estradiol (EE), demonstrate later compromised parental behaviors. Reductions in F1 parental behaviors might also be due to decreased emissions of F2 pup vocalizations. Thus, vocalizations of F2 male and female California mice pups born to F1 parents developmentally exposed to BPA, EE, or controls were examined. Postnatal days (PND) 2±4 were considered early postnatal period, PND 7 and 14 were defined as mid-postnatal period, and PND 21 and 28 were classified as late postnatal period. EE pups showed increased latency to emit the first syllable compared to controls. BPA female pups had decreased syllable duration compared to control and EE female pups during the early postnatal period but enhanced responses compared to controls at late postnatal period; whereas, male BPA and EE pups showed greater syllable duration compared to controls during early postnatal
Data Availability Statement: All relevant data are
within the paper and its Supporting Information
Funding: This work was supported by NIEHS Grant
5R21ES023150 (to C.S.R.). The funders had no
role in study design, data collection and analysis,
decision to publish, or preparation of the
Competing interests: The authors have declared
that no competing interests exist.
period. In mid-postnatal period, F2 BPA and EE pups emitted greater number of phrases
than F2 control pups. Results indicate aspects of vocalizations were disrupted in F2 pups
born to F1 parents developmentally exposed to BPA or EE, but their responses were not
always identical, suggesting BPA might not activate estrogen receptors to the same extent
as EE. Changes in vocalization patterns by F2 pups may be due to multigenerational
exposure to BPA or EE and/or reduced parental care received.
Diverse rodent species primarily communicate by ultrasonic vocalizations (USVs, sounds
above 20 kHz). We, and others, have characterized USV production by neonatal California
mice (Peromyscus californicus) and other rodent species [1±7]. Pup vocalizations may be
essential in triggering parental care [
2, 5, 8
]. This may be especially important in those that are
unable to regulate their own body temperature until later in life, which includes neotropical
singing mice (Scotinomys spp.)  and California mice [9±11]. The latter species are
monogamous and biparental [10, 12±14], and while the dam is foraging, her male partner will huddle
over the pups to prevent a decline in pup body temperature [
]. We have previously shown
that California mice pup vocalizations begin to decline around the time, postnatal day (PND)
15, they can regulate their own body temperature.
In our past study, we have shown that male and female California mice developmentally
exposed to the endocrine disrupting chemicals (EDCs), bisphenol A (BPA), which can exhibit
weak estrogenic activity or ethinyl estradiol (EE, estrogen present in birth control pills and
estrogen receptor positive control for the current studies) demonstrate reduced biparental
care, including duration of time female spent nursing the pups, time spent in the nest for both
parents, and time spent grooming the pups [
]. While the reduced parental care by the F1
parents developmentally exposed to BPA or EE could be due to neurobehavioral disruptions
in these individuals, it could also be attributed to decreased vocalizations emitted by their F2
offspring exposed multigenerationally to these chemicals. Developmental exposure of
Sprague-Dawley rats to polychlorinated biphenyls (PCBs, Aroclor 1221 PCR mixture) disrupt the
number of vocalizations produced by sexually mature female rats [
]. It has also recently
been reported that maternal exposure of C57BL/6J mice to 20 μg BPA per day via a
chocolateflavored treat one week prior to breeding and throughout gestation increased the duration and
median frequency of USVs emitted by F1 pups at PND 8 during maternal separation [
Several EDCs, vinclozolin, EE, p,p'-Dichlordiphenyldichloroethylene (p,p'-DDE), Levonorgesterel
(synthetic progesterone) disrupt mating calls produced by male South African clawed frogs
(Xenopus laevis) [18±22]. Zebra finch (Taeniopygia guttata) male and female progeny derived
from females exposed to PCB, Aroclor 1248, prior to egg-laying have significant reduction in
the song control nuclei robustus arcopallialis (RA) or now called the robust nucleus of the
arcopallium, which could compromise their song production ability [
]. Similarly, treatment
of adult male starlings (Sturnus vulgaris) with estrogenic EDCs was shown to result in louder,
longer and more complex songs compared to vocalizations of control males [
]. These effects
were suggested to be due to an increase in volume of the principle nucleus in the songbird
brain, the HVC(used as proper name) [
Based on these above studies showing direct exposure to various EDCs in a variety of taxa
or exposure of F1 rodent pups to the EDCs, BPA or Aroclor 1221 PCR mixture, alters
vocalization patterns [16±24], we hypothesized that audible and ultrasonic vocalizations might be
reduced in F2 pups multigenerationally exposed to BPA or EE, which could partially account
for the reduced parental care provided by their F1 parents. Another possibility that was
explored in the current study is whether pup vocalizations are compromised in F2 pups
derived from F1 offspring developmentally exposed to these EDCs and who demonstrate
decreased parental care to their F2 pups. In other words, we sought to examine whether F2
pups who are not directly exposed to BPA or EE, demonstrate changes in vocalization
parameters. The current studies thus are the first to examine vocalization behaviors in F2 offspring of
ancestors directly exposed to BPA or EE. However, the current studies were not aimed at
determining the underlying causes (reduced F1 parental care and/or multigenerational
exposure to BPA or EE) of potential alterations in F2 acoustic patterns. To test whether vocalization
2 / 18
patterns of F2 pups derived from F1 pairings of BPA- or EE-exposed parents were disrupted
compared to those born to control F1 parents, pups were examined from postnatal day (PND)
2 to 28. Additionally, we examined potentiation behavioral responses in the tested pups. In
this case, a potentiated behavioral response refers the to the fact that a pup who has already
been isolated may vocalize more after being returned to their parents and then isolated again,
in essence this behavioral response is even greater after being isolated and tested a second time
[7, 25±27]. Mid postnatal age rat pups returned to a passive (anesthetized dam) show increased
vocalization responses when isolated again, which is considered a potentiation response in
that re-isolation after being placed back with one or both parents causes an enhanced response
relative to the vocalization responses when initially isolated from one or both parents [
No maternal potentiation effects were observed on PND 8 in mice developmentally exposed to
BPA . Passive (anesthetized) and active maternal potentiation but not paternal
potentiation has been reported in prairie voles (Microtus ochrogaster), which are also monogamous
and biparental [
]. In the current studies, California mice pups from all treatment groups
were returned to their dam and sire, and thus, the current studies tested whether a
combination of maternal and paternal potentiation exists and whether it can be affected by
multigenerational exposure to BPA or EE in F2 pups.
Materials and methods
Animal husbandry and treatments
Founder outbred adult (60±90 days of age) male and female California mice that were
pathogen-free, were purchased from the Peromyscus Genetic Stock Center (PGSC) at the University
of South Carolina (Columbia, SC). Upon arrival, individuals were placed in quarantine at the
University of Missouri for a minimum of 8 weeks to verify their disease-free status. All
experiments were approved by University of Missouri Animal Care and Use Committee (Protocol
#8643) and performed in accordance with the recommendations in the Guide for the Care and
Use of Laboratory Animals of the National Institutes of Health. Two weeks prior to breeding,
virgin P0 females, 8 to 12 wks of age were randomly assigned to receive one of three diets: 1) a
low phytoestrogen AIN 93G diet supplemented with 7% by weight corn oil to minimize
potential phytoestrogenic contamination that would otherwise be present with inclusion of soybean
oil in the diet (control), 2) this diet supplemented with 50 mg BPA/kg feed weight, which we
have documented to lead to internal serum concentrations close to those measured in
pregnant women unknowingly exposed to this chemical [
], and 3) AIN93G diet
supplemented with 0.1 parts per billion of EE, as the FDA required positive control for BPA studies
based on the fact that BPA is considered a weak estrogen . Additionally, EE-exposed
individuals might provide insights into whether BPA effects might be due to BPA binding and
activating estrogen receptors (ESRs) to the same extent as EE. We have previously validated
similar effects of the chosen BPA and EE doses in deer mice (P. maniculatus bairdii) [
and California mice [15, 33±35]. The P0 dams remained on the diet throughout gestation and
lactation, as described previously [29, 32±36]. The F1 generation sons and daughters were
weaned at 30 days of age, singly housed and placed on the AIN93G (control) diet. Food and
water were provided ad libitum.
When F1 animals reached adulthood (~90 days of age), males and females from each group
were paired with respective partners who were developmentally exposed to the same diet, such
that there were three F1 breeding pairs: BPA-exposed for males and females, EE-exposed for
males and female, and AIN93G control for males and females. When individuals were not in
the Phenotyper system (Noldus Technologies, Leesburg, VA), they were housed in white
polypropylene cages (27.8 x 7.5 x 13 cm) and maintained on a 12:12 h light:dark cycle (lights on at
3 / 18
7:00 A.M. CST, lights off at 7:00 P.M. CST). As Peromyscus species do not form a visible plug,
F1 females were weighed weekly to provide some idea of when they might give birth. The day
of birth was considered PND 0.
F2 pup ultrasonic vocalization measurements and potentiation experiments
Each individual pup was tested in isolation and away from their home-cage, as we have
described previously [
]. The total number of individual pups and litters tested are indicated
in Table 1. This number of pups and litters is similar to previous studies that have examined
the effects of developmental exposure of other rodents to the EDCs, Aroclor 1221 PCR mixture
and BPA [
]. Pup USVs and potentiation experiments were initiated on PND 2. Each pup
was demarcated with a distinguishing tattoo mark on either the front or back paw (Fine
Science Tools, Foster City, CA). California female mice can only nurse up to four pups with an
average litter size of 2.2, and pups are weaned at 30 days of age [
7, 10, 15
]. All F2 pups in the
litter were recorded throughout the entire postnatal period (PND 28), assuming they survived
up until this time-point. Pup recordings were performed at ~10.00 hrs (3 hrs after lights turned
on in the room) and 14.00 hrs (7 hrs after lights turned on in the room). Before testing, the
experimenter noted whether they were or were not latched on to one of the dam's four teats
(S1 Table). The pups were placed individually on a laminated sheet in a polypropylene box (22
inches by 27 inches) surrounded on three sides with a 20 inch by 20 inch black polypropylene
box lined with two inches of convoluted acoustic foam panel (SoundproofCow,
Chambersburg, PA), as depicted in our previous study . The sheet demarcated pup and litter
information, date, and time of day when testing was done and was disinfected with 70% ethanol
between trials. No substrate material, including shavings or bedding material, was present in
the recording box. The box was equipped with an Avisoft Bioacoustics CM16/CMPA40-5V
microphone (Glienicke, Germany) plugged into an National Instruments USB 6351 data
collection board hooked up to a Dell computer (Roundrock, TX). USV data were collected using
a custom LabView data collection code written by the laboratory of Dr. Katrin Schenk at
Randolph College. A Logitech Carl Zeiss Tessar HD 1080P (Newark, CA) camera was mounted
onto a Joby- Gorilla Pod Original Tripod (Daymen US Inc., Petaluma, CA) in the recording
box to monitor isolated pups. Pups were audio-recorded for three minutes, which minimized
the time away from their parents but permitted relevant data to be collected for the time period
assessed. The pup was then placed back in the home cage with both parents for five minutes.
All pups within a litter were recorded for a given day before testing the next litter. However,
the order of testing was randomized between treatments such that pups within the same litter
were not always tested first or last on a given day.
After being returned to the home cage, the pup was isolated again and returned to the
audio-recording box and its vocalizations recorded for another three minutes as part of a
combination maternal and paternal potentiation experiment. The audio-recording box was
sprayed down between trials with 70% ethanol for disinfection and to remove any potential
4 / 18
confounding olfactory cues for subsequent pups tested. Pup vocalizations were measured on
PND (age of the pups) 2, 3, 4, 7, 14, 21, and 28. The previous study examining effects of
developmental exposure to BPA on mouse pups only recorded them at PND 8 [
]. However, as
California mice pups are hypothermic up until about PND 15 [9±11], we sought to examine
the F2 vocalization responses from the early PND to the late PND to determine how exposure
to BPA or EE might affect the full gamut of their vocalization responses from this period, as we
had previously done when characterizing the ontogenic profile of vocalization responses in
control California mice pups [
] and has been done with mandarin vole (Microtus
mandarinus) pups, where the parents are also monogamous and biparental [
]. Recordings on PND
2±4 were considered early postnatal period, those on days 7 and 14 were deemed
mid-postnatal period, and those on days 21 and 28 were considered late postnatal period. These
timepoints were chosen to span the early, mid, and late postnatal periods and key milestones in
California mice neonatal development. Based on our observations, the fur begins to emerge
and darken between two to four days of age. Eyes open around 14±15 days of age in this
species, which is also when they can begin to thermoregulate on their own [
7, 9, 10, 38, 39
California mice pups begin walking and rearing around 22 and 22.75 days of age [
]. Pups raised
by both parents begin self-grooming around 25 days of age. California mice pups are weaned
at 30 days of age and eat solid food at this time [
Vocalizations throughout the entire duration of recording (3 minutes for both trials) were
detected and measured by a MATLAB-based program developed by our collaborator and
coauthor (Dr. A. Katrin Schenk), which allows USVs up to 200 KHz to be detected and
analyzed, as described previously [
]. Briefly, for each recording, all syllables were segmented for
analysis. The experimenter was blind to treatment when denoising and segmenting the calls.
Syllables are synonymous with calls, as defined by other rodent EDC studies [
of syllables, syllable duration, syllable's mean median frequency, average syllable power (a
summation of power percent below and above 20 kHz), power percent below 20 kHz, and
power percent above 20 kHz were determined. The average power was calculated as follows: 1)
For each time bin in the short-time Fourier transform (spectrogram) of the syllable
determined the mean power in units of decibels above the noise level. 2) The mean of these values
for all time bins contained in the syllable's spectrogram was taken. The noise level used in step
one was calculated in the following way: 1) A 2 second region of the raw data spectrogram was
picked as a representative of the background noise. This region was identified by examining
the spectrogram by eye to make sure that there were no syllables or artifacts present. 2) The
spectrogram of this region was averaged over all time bins to create a baseline noise per
frequency bin. 3) Then the noise level per frequency was defined as 3 standard deviations about
the baseline noise. Thus, the average power is a measure of intensity. The median frequency
per call was averaged for each pup, hence the mean of the median frequencies. The median
frequency was chosen per Liu et al 2003  as a read out of the spectral content of the call.
While there are many different read outs that one might choose, we elected to use the median
frequency since it takes into account all components of the call (S1 Fig).
Additionally, the number of phrases, which designates a group of syllables separated from
another group of syllables by an interval in time that is statistically longer than the mean time
between individual syllables, duration of phrases, and number of syllables in a phrase were
measured. Phrase is equivalent to what has also been called a burst in a study examining the
effects of developmental exposure to BPA on mouse vocalizations [
]. The terminology of
syllables and phrases is used herein to be consistent with previous studies that have measured
vocalizations in California mice [
3, 6, 7, 41
]. Example syllables and phrases in California mice
are shown in Fig 1. Additional representative audible and ultrasonic vocalizations in neonatal
California mice pups are shown in .
5 / 18
Fig 1. Example audible and ultrasonic vocalizations in California mice pups. The top panel depicts syllables, which are the same thing
as calls, that begin in the audible range of below 20 kHz and then extend into the ultrasonic range. The bottom panel depicts USVs alone
that are above 20 kHz. A phrase or burst is a group of syllables separated from another group of syllables by an interval in time that is
statistically longer than the mean time between individual syllables. The ear lobe at the left side of each panel represents the fact that
human hearing only extends up to 20 kHz.
Data were analyzed by using SAS version 9.4 Software (SAS Institute, Cary, NC). Partial eta
squared values (Ƞ2p) were calculated as measure of effect size. A split plot in space and time
(an ANOVA-based method) was used to analyze the data, as detailed by [
]. The main plot
consisted of F1 maternal/paternal treatment, and the subplot was F2 pup sex and F1 parental
treatment X F2 pup sex. S2 Fig provides a diagram of how the split plot in space and time
controls for litter effects (split plot in space) and repeated measure analyses (time). For pups
where the sex was not determined, they were only considered in main plot. Additionally, the
sub-sub plots consisted of postnatal period (early, mid, or late), time of day (AM or PM), and
trial (Trial 1 or Trial 2) and all possible interactions with these variables and the main plot and
sub-plot effects. We sought to analyze the F2 pup vocalization responses over the entire course
of the post-natal period and in the AM and PM to determine whether the maternal/paternal
treatments altered the pattern of the vocalization responses as the pups matured and to
examine for potential circadian differences, as it relates to maternal/paternal treatment. In a
previous study with F1 California mice, we found that developmental exposure to BPA altered later
voluntary physical activity predominantly during the dark but not the light cycle [
we postulated that we might detect similar circadian differences in terms of F2 pup vocalization
responses. Mean differences were determined using Fisher's Protected Least Significant
Difference (LSD), as described by [
]. In this case, the LSD is only calculated if the overall F test is
significant and is considered an acceptable method under these conditions [
we only ran pre-planned and relevant comparisons where pup sex, day, time, and/or trial were
hold constant, and the only thing varied was F1 maternal/paternal treatment. Additionally,
only the highest order of interactions with maternal/paternal treatment was considered. Both
of these approaches reduced the number of comparisons and likelihood of committing a Type
I error. We decided to use the Fisher's protected LSD, as recommended by [
], as other
methods, including Bonferroni, may mask real differences, and thereby, increase the likelihood
of committing a Type II error. The Fisher's protected LSD helps avoid committing a type I
error but at the same time reducing the potential of committing a type II error. With the split
6 / 18
plot design, change in the standard error of the differences reflects more than one error term
(as detailed by [
]), which is taken into account with Fisher's protected LSD. Additionally, to
test for normality, we analyzed ranked data, which are the only ones presented in the Results
section. All data are presented though as actual means (x) rather averaged rank data. The error
bars for all figures and reported data represent the standard error of the mean (SEM).
The first variable considered was latency to produce the first syllable (otherwise considered a
call in [
]). A two-way interaction existed between maternal/paternal treatment (main effect)
and day (Fig 2A, η2p = 0.95, p = 0.0002). Latency to first syllable in the early postnatal days was
greater in EE compared to control pups (p = 0.003). Likewise, the percentage of syllable power
below and above 20 kHz showed a two-way interaction with maternal/paternal treatment
(main effect) and day (Fig 2B, η2p = 0.94, p 0.0001). During the early postnatal period,
syllable power above 20 kHz was greater in EE compared to control and BPA pups (p = 0.002 and
0.02, respectively, Fig 2B). During those postnatal days, the opposite was the case for syllable
power below 20 kHz with control and BPA pups showing greater syllable power than EE
individuals (p = 0.002 and 0.02, respectively).
Syllable duration demonstrated a three-way interaction with maternal/paternal treatment
(main effect), day, and pup sex (Fig 2C & 2D, η2p = 0.98, p 0.0001). During the early
postnatal period, female control and EE pups had longer syllable duration than female BPA pups
(p = 0.01 and p = 0.006, respectively Fig 2C). However, by the late postnatal period, BPA
females exhibited longer syllable duration than control females (p = 0.03, respectively). In the
early period, BPA and EE males had greater syllable duration than control males (p = 0.02 for
both, Fig 2D). In the mid-period, EE males continued to exhibit increased syllable duration
relative to control and BPA males (p = 0.006 and 0.0001, respectively, Fig 2D).
The mean of the median frequencies of the syllables was significantly affected by the
interaction of maternal/paternal treatment (main effect) and day (Fig 3A, η2p = 0.96, p 0.0001).
Fig 2. Syllable parameter results in F2 pups. A) Latency to first syllable. B) Syllable Power Percent Above 20 KHz. C)
Average syllable for female pups. D) Average syllable duration for male pups. Error bars represent SEM.
7 / 18
Fig 3. Other syllable parameter results in F2 pups. A) Median syllable frequency. B) Circadian effects on total
number of syllables. Error bars represent SEM.
In the early postnatal period the median syllable frequency was greater in EE compared to
control and BPA pups (p = 0.007 and 0.02, respectively, Fig 3A). Total number of syllables was
altered by maternal/paternal treatment (main effect) and time of day (Fig 3B, η2p = 0.96,
p = 0.03). In the morning, BPA pups called more than controls (p = 0.04, Fig 3B). None of the
other syllable variables showed significant maternal/paternal treatment differences.
The average number of phrases (also called a burst in [
]) was influenced by the interactions
of: 1) maternal/paternal treatment (main effect) and day and 2) maternal/paternal treatment
(main effect) and day (η2p = 0.96, p = 0.002 and η2p = 0.96, p = 0.04, respectively). In the
midpostnatal period, BPA and EE pups had greater average number of phrases than control
individuals (p = 0.003 and 0.0009, respectively, Fig 4). Circadian differences were also noted
between these groups with this category of vocalizations being greater in the morning for BPA
and EE pups compared to control pups (p = 0.01 and 0.03, respectively, S3 Fig). None of the
other phrase variables showed significant maternal/paternal treatment differences.
Examination of possible maternal and paternal potentiation effects
Potentiation experiments were also performed to test whether the pup responses were greater
after being isolated and tested in Trial 1, then placed back with both parents, and then
reFig 4. Average number of phrases in F2 pups. Error bars represent SEM.
8 / 18
isolated again to measure their vocalization responses for Trial 2. No differences in any of the
vocalizations parameters measured though showed potentiation differences, which is similar
to data observed when comparing F1 BPA-exposed and control mice [
The current studies tested the hypothesis that F2 California mice pups born to F1 parents
developmentally exposed to BPA or EE might emit reduced vocalizations. Potential reductions in F2
pup vocalizations might partially account for reduced biparental care previously identified in
F1 California mice developmentally exposed to BPA or EE [
]. However, vocalization
parameters that differed in BPA and EE pups were generally enhanced relative to control
counterparts, but temporal and sex differences were noted. While BPA female pups had increased
syllable duration compared to control pups during the late postnatal period, they showed
reduced response during the early postnatal period. In the mid-postnatal period, BPA and EE
pups emitted greater number of phrases than respective control pups. F2 BPA and EE pups
emitted greater number of phrases in the morning compared to F2 control pups. These
changes suggest that the reduced parental care provided early on by F1 BPA and EE parents
resulted in them emitting greater vocalizations, including syllables and phrases at later
The primary vocalization differences detected in the current studies between F2 BPA or EE
pups compared to control pups are summarized in Table 2. These F2 communication
responses might be due to direct effects of multigenerational exposure to these EDCs as the
primordial germ cells that gave rise to the F2 generation would have been simultaneously
exposed when the P0 parents consumed the diets and the F1 generation was developing (as
shown in Figure 11.1 in [
] and other relevant reviews [47±50]). Other studies have shown
that multigenerational exposure to BPA can lead to neurobehavioral, reproductive, and
cardio-metabolic changes in F2 rodent offspring [51±58].
Another possible but not mutually exclusive explanation is that in response to the reduced
parental care provided by their F1 parents, the F2 pups had augmented vocalization responses
in attempt to attract their parents to increase their parental investment. The enhanced
vocalization responses even up into the late postnatal period, when the pups can thermoregulate and
· Early postnatal period: EE pups> control pups
· Early postnatal period: control and EE females > BPA females.
· Late postnatal period: BPA females > control females
· Early postnatal period: BPA and EE males >control males.
· Mid postnatal period: EE males > control and BPA males
· Early postnatal period: EE pups > control and BPA pups
· AM: BPA pups > controls pups
· Early postnatal period: EE pups > control and BPA pups for
those above 20 kHz
· Early postnatal period: control and BPA pups > EE pups for
those below 20 kHz
· Mid postnatal period: BPA and EE pups > control pups
· AM: BPA and EE pups >control pups
9 / 18
begin to feed on their own might also suggest increased anxiogenic behaviors. Pup body
weight might also affect acoustic parameters with heavier pups presumably being able to
generate increased magnitude of vocalization responses. While we did not measure F2 pup weight
in the current studies, we had done so in previous studies [
]. We found that F2 pups of F1
parents developmentally exposed to BPA had no difference in body weight gain compared to
F2 pups of F1 control parents. However, F2 pups of F1 parents developmentally exposed to EE
showed reduce weight gain from PND 13±32. Yet, in the current studies F2 EE females and
males showed increased magnitude of vocalization responses, suggesting body weight was not
a contributory factor to treatment differences observed in this group. In mandarin voles,
which are also monogamous and biparental, there is surprisingly a positive relationship
between pup USV emissions and level of parental care received by pups [
]. Follow up studies
are thus needed to determine whether various behavioral responses, including vocalization,
anxiogenic, cognitive, and sociosexual behaviors continue to be affected as the F2 pups mature.
To tease apart whether the postnatal vocalization responses in the F2 pups are
multigenerational in origin or due to poor parental care provided by F1 parents, cross-fostering approaches
might be useful. This procedure was not done in these initial studies, as we were concerned
whether the foster parents would successfully rear the pups. Additionally, we wished to use the
same approach where reduced parental care was identified in F1 BPA and EE exposed parents
]. Such future cross-fostering studies would entail placing control F2 pups with F1 BPA or
EE exposed parents and F2 BPA/EE pups with F1 control parents. If the elevated vocalization
responses in F2 BPA and EE pups were multigenerational in origin, the prediction would be
that they would continue when placed with control parents. However, if they originate due to
receiving inadequate parental care by their F1 exposed parents, then F2 BPA and EE pups
should have reduced vocalizations when placed with control F1 parents who will presumably
be more attentive.
The expectation at the outset was that similar vocalization responses would be observed in
F2 BPA or EE pups. While such was the case for select parameters measured, F2 EE pups
differed from both control and BPA counterparts for various vocalization aspects, such as latency
to emit first syllable, median call frequency, percentage of syllable power below and above
20KHz. The more pronounced effects in the EE group might be attributed to differences in
estrogen receptor binding and activation, as BPA is considered a weak-estrogen [
might also suggest that multigenerational exposure to BPA and EE can result in differing
effects. In our study examining effects of developmental exposure to BPA or EE on biparental
behaviors in California mice, both chemicals led to similar disruptions for nursing behaviors
in the case of the dams and time spent in the nest for both parents . However, those
exposed to EE groomed their F2 pups more relative to control parents with control F2 pups (as
determined based on measuring the frequency and duration of grooming from archived videos
taken with an infra-red camera/ Phenotyper (Noldus Technologies, Leesburg, VA), that
recorded at several time-points throughout the light and dark cycle and on several days
spanning from PND 2 to PND 30), but such responses were not observed in those exposed to BPA.
In other studies examining neurobehavioral disruptions induced by developmental exposure
to BPA and EE in California mice, deer mice (P. maniculatus bairdii), Sprague-Dawley rats,
and painted turtles (Chrysemys picta), we have found that these EDCs do not always induce
similar gene expression changes in the brain or behavioral alterations [29, 32±36, 61, 62].
Thus, the current data suggest that EE may not be the ideal positive control for BPA
developmental exposure studies and other steroidogenic compounds, such as androgens, should
possibly be tested as well.
Another prediction at the outset was that F2 BPA and EE exposed male and female pups
would should similar vocalization changes. However, sexually dimorphic differences were
10 / 18
observed for syllable duration. While BPA females showed reductions in this category during
the early postnatal period relative to control and EE females, BPA and EE males showed
enhanced response relative to controls. BPA females demonstrated increased syllable duration
compared to controls during the late postnatal period, but no differences were detected at this
time period between BPA and control pups. Other studies indicate that wild-type and
transgenic rodents demonstrate sex differences in vocalization patterns at various postnatal ages.
For instance, male rat pups placed in isolation elicit substantially more USV calls that were
typified by a significantly lower frequency and amplitude compared to isolate female rat pups
]. This same study also showed that dams are more likely to retrieve and place male pups
back into the nest than females. Testing of vocalization patterns of fragile X mental retardation
1 (Fmr1) knockout male and female mice from postnatal days 9±14 reveals that males and
females diverge in the total number of calls and call duration depending on the postnatal day
], which is similar to the results identified in the current study with F2 BPA female and
To date, two other studies have examined the effects of exposure to EDCs on rodent
vocalizations. Developmental exposure to the PCB mixture, Aroclor 122, disrupts the number of
vocalizations produced later on by sexually mature adult female Sprague-Dawley rats [
Maternal exposure of C57BL/6J mice to 20 μg BPA per day one week prior to breeding and
throughout gestation increased the duration and median frequency of USVs, as well as the
total number of ªburstsº (which is the same as phrases) emitted by F1 pups during maternal
separation when tested at PND 8 [
]. However, in this study, the number of calls was not
altered by in utero exposure to BPA. The initial study testing the effects of PCBs on rat
vocalizations was done prior to some of the recent advancements in detecting and analyzing rodent
vocalization, and thus, EDC-induced disturbances on other facets of both audible and
ultrasonic vocalizations might have been missed. The latter study with F1 mice pups used an
identical microphone and software program as that used herein. Notably, BPA exposure of F1 pups
in the Harris et al study [
] could have induced direct effects on brain development in these
pups that might have in turn affected their vocalization responses. However, the developing
brain and other organs of F2 pups were not directly exposed to BPA but only
multigenerationally exposed to BPA via the primordial germ cells contained within the developing F1 offspring
(as shown in Figure 11.1 in [
] and other relevant reviews [47±50]). Thus, any differences in
responses might be attributed to generation examined or developmental vs. multigenerational
exposure (F1 vs. F2, respectively), species (Mus musculus vs. P. californicus), and original dose
of maternal BPA exposure. In the current study, F2 California mice pups were analyzed
throughout the neonatal period, whereas in the BPA study above, F1 mouse pups were only
assessed at PND 8 [
], which represents approximately the mid-postnatal period in mice that
are generally weaned PND 21. This original study also sought to examine for other behavioral
differences due to maternal vs. paternal exposure [
]. This above study tested laboratory mice
where the male's only involvement in the offspring is at the time of fertilization in terms of
contribution of DNA and possible epigenetic factors [
]. These males are not involved in
parental care for this species. In contrast, California mice pups rely on biparental care, as
detailed in past studies [9±11]. Thus, in this latter species, it is difficult to dissect the paternal
effects due to germ cells vs. that due to disrupted paternal responses. One potential way to
address this issue would be to use in vitro fertilization and embryo transfer procedures. To
date, such procedures, however, have not been performed with California mice. However,
blastocysts have been derived in culture for their related cousins, P. maniculatus [
A previous study has shown that ultrasonic vocalizations can differ between males of
different mice species [
]. It is interesting to note that the study testing effects of BPA on F1
vocalizations also noted increased magnitude of responses [
], which might also be attributed to
11 / 18
reduced parental care provided by P0 dams directly exposed to BPA and/or developmental
effects on brain regions regulating vocalization responses.
As multigenerational studies require several rounds of breeding and California mice, which
produce relatively small litters, we limited ourselves to testing F2 pups derived from F1 males
and females exposed to a single dose of BPA that had previously been shown to disrupt F1
parental care in this species [
] and leads to internal serum concentrations comparable to
those detected in humans unknowingly exposed to this chemical [
]. Future studies
though should consider testing other doses of BPA to examine for potential non-monotonic
dose responses .
Additional studies thus should also examine how developmental and multigenerational
exposure to BPA affects molecular responses, in the brain regions, including the ventral
medial hypothalamus, thalamus, central amygdala, cortex, striatum, cerebellar vermis, and
nucleus ambiguous, governing pup audible and ultrasonic vocalizations [63, 69±71].
Potential hormones/genes that regulate rodent pup vocalizations and might be altered by
developmental or multigenerational exposure to BPA or EE include serotonin (5-HT), forkhead
box P1 (Foxp1), forkhead box 2 (Foxp2), dopamine receptor D2 (Drd2), neuroligin 2 (Nlgn2)
[63, 70, 72±74]. Exposure of zebra finch male and female progeny derived from females
exposed to the PCB Aroclor 1248, prior to egg-laying have significant reduction in the song
control nuclei RA, which would likely compromise their song production ability in
]. The increased song complexity in adult male starlings exposed to estrogenic EDCs
is likely due to an increase in volume of the principle nucleus in the songbird brain, the high
vocal center (HVC) [
As part of the current studies, we also examined whether maternal and paternal
potentiation responses might occur after the pups are placed back with their parents and then their
vocalizations responses tested again after re-isolation. In our original studies testing solely
control California mice pups, we did not notice any such differences [
], although maternal
potentiation responses have been reported in prairie voles, which are also monogamous and
]. When all pups are considered in the current study, the opposite effects of
decreased number of syllables, suggestive of potential habituation, was observed in the
second trial. However, the opposite was the case for the number of phrases that increased in
Trial 2. EE pups tended to communicate more in the audible range in Trial 1 but in the
ultrasonic range in Trial 2 in mid-postnatal period. Taken together, it is difficult to ascertain
whether California mice pups demonstrate definitive biparental potentiation responses.
Additional studies possibly testing longer periods of pup isolation might verify one way or
In summary, the current studies demonstrate that during specific postnatal periods, BPA
and EE F2 pups demonstrate augmented vocalization responses. However, temporal and
sexspecific differences are observed. Additionally, the responses of F2 EE and BPA pups were not
always identical. The vocalization alterations in these groups might be due to
multigenerational EDC exposure. On the other hand, increased vocalization responses by these F2 pups
might be their attempts at signaling to their inattentive F1 parents to provide additional
parental care. The current studies though suggest that the previously observed reduced parental care
provided by F1 parents developmentally exposed to BPA or EE is likely not due to
communication deficits by their F2 pups. Additional studies are needed to sort out the underlying
mechanisms driving the BPA and EE F2 pup vocalization changes and whether it is attributed to
intrinsic (e.g. neural in origin due to multigenerational exposure to these EDCs) and/or
extrinsic factors (i.e. reduced F1 parental care).
12 / 18
S1 Fig. Median frequency calculation. A) Spectrogram showing one syllable. B) Points in the
spectrogram that are above the noise level (background noise + 3 standard deviations). C)
Histogram of frequencies are shown in the middle figure. The median of this distribution is the
syllable's ªmedian frequencyº and is demarcated by an "X".
S2 Fig. Diagram of the general experimental design and split plot in space analyses. The
ªsplit plot in spaceº refers to the fact that the vocalization data for each pup are nested together
under the dam and sire, which controls for potential litter effects. The ªtimeº portion considers
the fact that the pup vocalization responses were repeatedly measured (time of day, trial, and
over the postnatal period).
S3 Fig. Circadian effects on total number of phrases by F2 pups. Error bars represent SEM.
S1 Table. Number of pups in parentheses latched onto dam prior to audio recordings in
the morning and afternoon on each trial day. Total number of pups for which data were
collected are listed in parentheses.
We are grateful to Angela Javurek, Emily Howald, and Paige Yeager who assisted with the
animal husbandry and experiments. Additionally, we appreciate Wayne Shoemaker, Leon
Toebben, and Danny Patterson for constructing the soundproof boxes.
Conceptualization: Sarah A. Johnson, Frauke Hoffmann, A. Katrin Schenk, Cheryl S.
Formal analysis: Sarah A. Johnson, Michelle J. Farrington, Claire R. Murphy, Paul D. Caldo,
Leif A. McAllister, Sarabjit Kaur, Catherine Chun, Madison T. Ortega, Brittney L. Marshall,
Frauke Hoffmann, Mark R. Ellersieck, A. Katrin Schenk, Cheryl S. Rosenfeld.
Funding acquisition: Cheryl S. Rosenfeld.
Investigation: Sarah A. Johnson, Michelle J. Farrington, Claire R. Murphy, Paul D. Caldo, Leif
A. McAllister, Sarabjit Kaur, Catherine Chun, Madison T. Ortega, Brittney L. Marshall,
Cheryl S. Rosenfeld.
Methodology: Sarah A. Johnson.
Project administration: Cheryl S. Rosenfeld.
Software: A. Katrin Schenk.
Supervision: Cheryl S. Rosenfeld.
Writing ± original draft: Sarah A. Johnson, Michelle J. Farrington, A. Katrin Schenk, Cheryl
13 / 18
Writing ± review & editing: Sarah A. Johnson, Michelle J. Farrington, Claire R. Murphy, Paul
D. Caldo, Leif A. McAllister, Sarabjit Kaur, Catherine Chun, Madison T. Ortega, Brittney L.
Marshall, Frauke Hoffmann, Mark R. Ellersieck, A. Katrin Schenk, Cheryl S. Rosenfeld.
14 / 18
15 / 18
californicus): A monogamous animal model. PloS one. 2013; 8(2):e55698. Epub 2013/02/14. https://
doi.org/10.1371/journal.pone.0055698 PMID: 23405200
16 / 18
17 / 18
1. Vieira M , Brown R . Ultrasonic vocalizations and ontogenetic development in California mice (Peromyscus californicus) . Behavioural processes . 2002 ; 59 ( 3 ): 147 . Epub 2002/09/25. PMID: 12270517 .
Wright SL , Brown RE . Sex differences in ultrasonic vocalizations and coordinated movement in the California mouse (Peromyscus californicus) . Behavioural processes . 2004 ; 65 ( 2 ): 155 ± 62 . Epub 2004/06/ 30. https://doi.org/10.1016/j.beproc. 2003 . 09 .004 PMID: 15222963 .
3. Pultorak JD , Fuxjager MJ , Kalcounis-Rueppell MC , Marler CA . Male fidelity expressed through rapid testosterone suppression of ultrasonic vocalizations to novel females in the monogamous California mouse . Hormones and behavior . 2015 ; 70 : 47 ± 56 . Epub 2015/03/01. https://doi.org/10.1016/j.yhbeh. 2015 . 02 .003 PMID: 25725427 .
4. Campbell P , Pasch B , Warren AL , Phelps SM . Vocal ontogeny in neotropical singing mice (Scotinomys) . PloS one . 2014 ; 9 ( 12 ): e113628 . Epub 2014 /12/04. https://doi.org/10.1371/journal.pone. 0113628 PMID: 25469986
5. Portfors CV . Types and functions of ultrasonic vocalizations in laboratory rats and mice . J Am Assoc Lab Anim Sci . 2007 ; 46 ( 1 ): 28 ± 34 . Epub 2007/01/06. PMID: 17203913 .
6. Briggs JR , Kalcounis-Rueppell MC . Similar acoustic structure and behavioural context of vocalizations produced by male and female California mice in the wild . Animal behaviour . 2011 ; 82 ( 6 ): 1263 ± 73 . http://dx.doi.org/10.1016/j.anbehav. 2011 . 09 .003.
7. Johnson SA , Painter MS , Javurek AB , Murphy CR , Howald EC , Khan ZZ , et al. Characterization of vocalizations emitted in isolation by California mouse (Peromyscus californicus) pups throughout the postnatal period . Journal of comparative psychology (Washington, DC: 1983 ). 2017 ; 131 ( 1 ): 30 ± 9 . Epub 2017/02/10. https://doi.org/10.1037/com0000057 PMID: 28182483
8. Geyer LA . Ontogeny of ultrasonic and locomotor responses to nest odors . American Zoologist . 1981 ; 21 : 117 ± 28 .
9. Gubernick DJ , Alberts JR . The biparental care system of the California mouse, Peromyscus californicus . Journal of Comparative Psychology . 1987 ; 101 ( 2 ): 169 ± 77 . Epub 1987/06/01. PMID: 3608423 .
10. Rosenfeld CS , Johnson S.A. Ellersieck M.R. , and Roberts R.M. Interactions between parents and parents and pups in the monogamous California mouse (Peromyscus californicus) PloS one . 2013 ; 8 e75725. https://doi.org/10.1371/journal.pone. 0075725 PMID: 24069441
11. Gubernick DJ , Teferi T . Adaptive significance of male parental care in a monogamous mammal . Proceedings Biological sciences. 2000 ; 267 ( 1439 ): 147 ± 50 . Epub 2000/02/25. https://doi.org/10.1098/rspb. 2000 .0979 PMID: 10687819
12. Chandler T. Population biology of coastal chaparral rodents .: University of California, Los Angeles; 1979 .
13. Eisenberg J . Studies on the behavior of Peromyscus maniculatus gambelii and Peromyscus californicus parasiticus . Behaviour . 1962 ; 19 ( 3 ): 177 ± 207 .
14. Eisenberg J. The intraspecific social behavior of some Cricetine rodents of the genus Peromyscus . American Midland Naturalist. 1963 ; 69 : 240 ± 6 .
15. Johnson SA , Javurek AB , Painter MS , Peritore MP , Ellersieck MR , Roberts RM , et al. Disruption of parenting behaviors in California Mice, a monogamous rodent species, by endocrine disrupting chemicals . PloS one . 2015 ; 10 ( 6 ): e0126284 . Epub 2015 /06/04. https://doi.org/10.1371/journal.pone.0126284 PMID: 26039462 .
16. Steinberg RM , Juenger TE , Gore AC . The effects of prenatal PCBs on adult female paced mating reproductive behaviors in rats . Hormones and behavior . 2007 ; 51 ( 3 ): 364 ± 72 . Epub 2007/02/06. https://doi. org/10.1016/j.yhbeh. 2006 . 12 .004 PMID: 17274994
17. Harris EP , Allardice HA , Schenk AK , Rissman EF . Effects of maternal or paternal bisphenol A exposure on offspring behavior . Hormones and behavior . 2017. Epub 2017 /10/02. https://doi.org/10.1016/j. yhbeh. 2017 . 09 .017 PMID: 28964733 .
18. Hoffmann F , Kloas W. An environmentally relevant endocrine-disrupting antiandrogen, vinclozolin, affects calling behavior of male Xenopus laevis . Hormones and behavior . 2010 ; 58 ( 4 ): 653 ± 9 . Epub 2010/07/06. https://doi.org/10.1016/j.yhbeh. 2010 . 06 .008 PMID: 20600051 .
19. Hoffmann F , Kloas W. Estrogens can disrupt amphibian mating behavior . PloS one . 2012 ; 7 ( 2 ): e32097 . Epub 2012 /02/23. https://doi.org/10.1371/journal.pone. 0032097 PMID: 22355410
20. Hoffmann F , Kloas W. The antiestrogens tamoxifen and fulvestrant abolish estrogenic impacts of 17alpha-ethinylestradiol on male calling behavior of Xenopus laevis . PloS one . 2012 ; 7 ( 9 ): e44715 . Epub 2012 /10/03. https://doi.org/10.1371/journal.pone. 0044715 PMID: 23028589
21. Hoffmann F , Kloas W. The synthetic progestogen, Levonorgestrel, but not natural progesterone, affects male mate calling behavior of Xenopus laevis . General and comparative endocrinology . 2012 ; 176 ( 3 ): 385 ± 90 . Epub 2012/03/07. https://doi.org/10.1016/j.ygcen. 2012 . 02 .009 PMID: 22391239 .
22. Hoffmann F , Kloas W. p,p' -Dichlordiphenyldichloroethylene (p,p'-DDE) can elicit antiandrogenic and estrogenic modes of action in the amphibian Xenopus laevis . Physiology & behavior . 2016 ; 167 : 172 ± 8 . Epub 2016/10/25. https://doi.org/10.1016/j.physbeh. 2016 . 09 .012 PMID: 27640133 .
23. Hoogesteijn AL , Kollias GV , Quimby FW , De Caprio AP , Winkler DW , DeVoogd TJ. Development of a brain nucleus involved in song production in zebra finches (Taeniopygia guttata) is disrupted by Aroclor 1248 . Environmental toxicology and chemistry . 2008 ; 27 ( 10 ): 2071 ± 5 . Epub 2008/04/03. https://doi.org/ 10. 1897 /07- 484 .1 PMID: 18380522 .
24. Markman S , Leitner S , Catchpole C , Barnsley S , Muller CT , Pascoe D , et al. Pollutants increase song complexity and the volume of the brain area HVC in a songbird . PloS one . 2008 ; 3 ( 2 ): e1674 . Epub 2008 /02/28. https://doi.org/10.1371/journal.pone. 0001674 PMID: 18301751
25. Shair HN , Brunelli SA , Masmela JR , Boone E , Hofer MA . Social, thermal, and temporal influences on isolation-induced and maternally potentiated ultrasonic vocalizations of rat pups . Dev Psychobiol . 2003 ; 42 ( 2 ): 206 ± 22 . Epub 2003/01/30. https://doi.org/10.1002/dev.10087 PMID: 12555284 .
26. Hofer MA , Masmela JR , Brunelli SA , Shair HN . Behavioral mechanisms for active maternal potentiation of isolation calling in rat pups . Behavioral neuroscience . 1999 ; 113 ( 1 ): 51 ± 61 . Epub 1999/04/10. PMID: 10197906 .
27. Myers MM , Ali N , Weller A , Brunelli SA , Tu AY , Hofer MA , et al. Brief maternal interaction increases number, amplitude, and bout size of isolation-induced ultrasonic vocalizations in infant rats (Rattus norvegicus) . Journal of comparative psychology (Washington, DC: 1983 ). 2004 ; 118 ( 1 ): 95 ± 102 . Epub 2004/03/11. https://doi.org/10.1037/ 0735 - 7036 . 118 .1.95 PMID: 15008677 .
28. Robison WT , Myers MM , Hofer MA , Shair HN , Welch MG . Prairie vole pups show potentiated isolationinduced vocalizations following isolation from their mother, but not their father . Dev Psychobiol . 2016 ; 58 ( 6 ): 687 ± 99 . Epub 2016/03/19. https://doi.org/10.1002/dev.21408 PMID: 26990108 .
29. Jasarevic E , Sieli PT , Twellman EE , Welsh TH Jr., Schachtman TR , Roberts RM , et al. Disruption of adult expression of sexually selected traits by developmental exposure to bisphenol A . Proceedings of the National Academy of Sciences of the United States of America . 2011 ; 108 ( 28 ): 11715 ± 20 . Epub 2011/06/29. https://doi.org/10.1073/pnas.1107958108 PMID: 21709224
30. Sieli PT , Jasarevic E , Warzak DA , Mao J , Ellersieck MR , Liao C , et al. Comparison of serum bisphenol A concentrations in mice exposed to bisphenol A through the diet versus oral bolus exposure . Environmental health perspectives . 2011 ; 119 ( 9 ): 1260 ± 5 . Epub 2011/06/07. https://doi.org/10.1289/ehp. 1003385 PMID: 21642047
31. vom Saal FS , Richter CA , Ruhlen RR , Nagel SC , Timms BG , Welshons WV . The importance of appropriate controls, animal feed, and animal models in interpreting results from low-dose studies of bisphenol A . Birth Defects Res A Clin Mol Teratol . 2005 ; 73 ( 3 ): 140 ±5. https://doi.org/10.1002/bdra.20120 PMID: 15751043 .
32. Jasarevic E , Williams SA , Vandas GM , Ellersieck MR , Liao C , Kannan K , et al. Sex and dose-dependent effects of developmental exposure to bisphenol A on anxiety and spatial learning in deer mice (Peromyscus maniculatus bairdii) offspring. Hormones and behavior . 2013 ; 63 : 180 ± 9 . Epub 2012/10/12. https://doi.org/10.1016/j.yhbeh. 2012 . 09 .009 PMID: 23051835 .
33. Johnson SA , Javurek AB , Painter MS , Ellersieck MR , Welsh TH Jr., Camacho L , et al. Effects of developmental exposure to bisphenol A on spatial navigational learning and memory in rats: A CLARITYBPA study . Hormones and behavior . 2015. Epub 2015 /10/06. https://doi.org/10.1016/j.yhbeh. 2015 . 09 . 005 PMID: 26436835 .
34. Johnson SA , Spollen WG , Manshack LK , Bivens NJ , Givan SA , Rosenfeld CS . Hypothalamic transcriptomic alterations in male and female California mice (Peromyscus californicus) developmentally exposed to bisphenol A or ethinyl estradiol . Physiological reports . 2017 ; 5 ( 3 ). Epub 2017 /02/16. https:// doi.org/10.14814/phy2.13133 PMID: 28196854
35. Johnson SA , Painter MS , Javurek AB , Ellersieck MR , Wiedmeyer CE , Thyfault JP , et al. Sex-dependent effects of developmental exposure to bisphenol A and ethinyl estradiol on metabolic parameters and voluntary physical activity . Journal of developmental origins of health and disease . 2015 ; 6 ( 6 ): 539 ± 52 . Epub 2015/09/18. https://doi.org/10.1017/S2040174415001488 PMID: 26378919 36 .
37. Yu P , Wang J , Tai F , Broders H , An S , Zhang X , et al. The effects of repeated early deprivation on ultrasonic vocalizations and ontogenetic development in mandarin vole pups . Behavioural processes . 2011 ; 88 ( 3 ): 162 ± 7 . Epub 2011/09/20. https://doi.org/10.1016/j.beproc. 2011 . 09 .001 PMID: 21925576 .
38. Vieira ML , Brown RE . Effects of the presence of the father on pup development in California mice (Peromyscus californicus) . Developmental Psychobiology . 2003 ; 42 ( 3 ): 246 ± 51 . Epub 2003/03/07. https:// doi.org/10.1002/dev.10097 PMID: 12621650 .
39. Harris BN , de Jong TR , Yang V , Saltzman W. Chronic variable stress in fathers alters paternal and social behavior but not pup development in the biparental California mouse (Peromyscus californicus) . Hormones and behavior . 2013 ; 64 ( 5 ): 799 ± 811 . http://dx.doi.org/10.1016/j.yhbeh. 2013 . 10 .007. PMID: 24157379
40. Liu RC , Miller KD , Merzenich MM , Schreiner CE . Acoustic variability and distinguishability among mouse ultrasound vocalizations . The Journal of the Acoustical Society of America . 2003 ; 114 ( 6 Pt 1 ): 3412 ± 22 . Epub 2004/01/13. PMID: 14714820 .
41. Kalcounis-Rueppell MC , Metheny JD , Vonhof MJ . Production of ultrasonic vocalizations by Peromyscus mice in the wild . Front Zool . 2006 ; 3:3 . Epub 2006/03/02. https://doi.org/10.1186/ 1742 -9994-3-3 PMID: 16507093
42. Steel RG . Principles and Procedures of Statistics: A Biometrical Approach . 3rd ed: McGraw-Hill Higher Education; 1997 .
43. Chew V , United S. Comparisons among treatment means in an analysis of variance . [ Washington]: Dept. of Agriculture, Agricultural Research Service; 1977 . iv, 64 p. p.
44. Saville DJ . Multiple Comparison Procedures: The Practical Solution . The American Statistician . 1990 ; 44 ( 2 ): 174 ± 80 . https://doi.org/10.2307/2684163
45. Cochran WG , Cox G.M. Experimental Designs . 2nd Edition . Wiley, editor1992 .
46. Rosenfeld CS . Chapter 11ÐAnimal Models of Transgenerational Epigenetic Effects A2ÐTollefsbol, Trygve . Transgenerational Epigenetics . Oxford: Academic Press; 2014 . p. 123 ± 45 .
47. Grandjean P , Barouki R , Bellinger DC , Casteleyn L , Chadwick LH , Cordier S , et al. Life-Long Implications of Developmental Exposure to Environmental Stressors: New Perspectives . Endocrinology. 2015 ; 156 ( 10 ): 3408 ± 15 . Epub 2015/08/05. https://doi.org/10.1210/EN.2015-1350 PMID: 26241067
48. Leon-Olea M , Martyniuk CJ , Orlando EF , Ottinger MA , Rosenfeld C , Wolstenholme J , et al. Current concepts in neuroendocrine disruption . General and comparative endocrinology . 2014 ; 203 : 158 ± 73 . Epub 2014/02/18. https://doi.org/10.1016/j.ygcen. 2014 . 02 .005 PMID: 24530523
49. Rosenfeld CS . Animal models to study environmental epigenetics . Biology of reproduction . 2010 ; 82 ( 3 ): 473 ± 88 . Epub 2009/10/16. https://doi.org/10.1095/biolreprod.109.080952 PMID: 19828779 .
50. Rosenfeld CS . Bisphenol A and phthalate endocrine disruption of parental and social behaviors . Frontiers in neuroscience. 2015 ; 9 : 57 . Epub 2015/03/19. https://doi.org/10.3389/fnins. 2015 .00057 PMID: 25784850
51. Berger A , Ziv-Gal A , Cudiamat J , Wang W , Zhou C , Flaws JA . The effects of in utero bisphenol A exposure on the ovaries in multiple generations of mice. Reproductive toxicology (Elmsford , NY). 2016 ; 60 : 39 ± 52 . Epub 2016/01/10. https://doi.org/10.1016/j.reprotox. 2015 . 12 .004 PMID: 26746108
52. Jang YJ , Park HR , Kim TH , Yang WJ , Lee JJ , Choi SY , et al. High dose bisphenol A impairs hippocampal neurogenesis in female mice across generations . Toxicology . 2012 ; 296 ( 1 ±3): 73 ± 82 . Epub 2012/ 04/10. https://doi.org/10.1016/j.tox. 2012 . 03 .007 PMID: 22484357 .
53. Lombo M , Fernandez-Diez C , Gonzalez-Rojo S , Navarro C , Robles V , Herraez MP . Transgenerational inheritance of heart disorders caused by paternal bisphenol A exposure. Environmental pollution (Barking , Essex: 1987 ). 2015 ; 206 : 667 ± 78 . Epub 2015/09/01. https://doi.org/10.1016/j.envpol. 2015 . 08 . 016 PMID: 26322593 .
Wolstenholme JT , Edwards M , Shetty SR , Gatewood JD , Taylor JA , Rissman EF , et al. Gestational exposure to bisphenol a produces transgenerational changes in behaviors and gene expression . Endocrinology . 2012 ; 153 ( 8 ): 3828 ± 38 . Epub 2012/06/19. https://doi.org/10.1210/en.2012-1195 PMID: 22707478
55. Ziv-Gal A , Wang W , Zhou C , Flaws JA . The effects of in utero bisphenol A exposure on reproductive capacity in several generations of mice . Toxicology and applied pharmacology . 2015 ; 284 ( 3 ): 354 ± 62 . Epub 2015/03/17. https://doi.org/10.1016/j.taap. 2015 . 03 .003 PMID: 25771130
56. Li G , Chang H , Xia W , Mao Z , Li Y , Xu S. F0 maternal BPA exposure induced glucose intolerance of F2 generation through DNA methylation change in Gck . Toxicology letters . 2014 ; 228 ( 3 ): 192 ± 9 . Epub 2014/05/06. https://doi.org/10.1016/j.toxlet. 2014 . 04 .012 PMID: 24793715 .
57. Mahalingam S , Ther L , Gao L , Wang W , Ziv-Gal A , Flaws JA . The effects of in utero bisphenol A exposure on ovarian follicle numbers and steroidogenesis in the F1 and F2 generations of mice. Reproductive toxicology (Elmsford , NY). 2017 ; 74 : 150 ± 7 . Epub 2017/10/04. https://doi.org/10.1016/j.reprotox. 2017 . 09 .013 PMID: 28970132 .
58. Susiarjo M , Xin F , Bansal A , Stefaniak M , Li C , Simmons RA , et al. Bisphenol a exposure disrupts metabolic health across multiple generations in the mouse . Endocrinology . 2015 ; 156 ( 6 ): 2049 ± 58 . Epub 2015/03/26. https://doi.org/10.1210/en.2014-2027 PMID: 25807043
59. Vandenberg LN , Colborn T , Hayes TB , Heindel JJ , Jacobs DR Jr., Lee DH , et al. Hormones and Endocrine-Disrupting Chemicals : Low-Dose Effects and Nonmonotonic Dose Responses . Endocrine Reviews . 2012 . Epub 2012/03/16. https://doi.org/10.1210/er.2011-1050 PMID: 22419778 .
60. Vandenberg LN , Maffini MV , Sonnenschein C , Rubin BS , Soto AM . Bisphenol-A and the great divide: a review of controversies in the field of endocrine disruption . Endocr Rev . 2009 ; 30 ( 1 ): 75 ± 95 . Epub 2008/ 12/17. https://doi.org/10.1210/er.2008-0021 PMID: 19074586
61. Manshack LK , Conard CM , Bryan SJ , Deem SL , Holliday DK , Bivens NJ , et al. Transcriptomic alterations in the brain of painted turtles (Chrysemys picta) developmentally exposed to bisphenol A or ethinyl estradiol . Physiological genomics . 2017 ; 49 ( 4 ): 201 ± 15 . https://doi.org/10.1152/physiolgenomics. 00103. 2016 PMID: 28159858 .
62. Manshack LK , Conard CM , Johnson SA , Alex JM , Bryan SJ , Deem SL , et al. Effects of Developmental Exposure to Bisphenol A and Ethinyl Estradiol on Spatial Navigational Learning and Memory in Painted Turtles (Chrysemys picta) . Hormones and behavior . 2016. Epub 2016 /08/02. https://doi.org/10.1016/j. yhbeh. 2016 . 07 .009 PMID: 27476434 .
63. Bowers JM , Perez-Pouchoulen M , Edwards NS , McCarthy MM . Foxp2 mediates sex differences in ultrasonic vocalization by rat pups and directs order of maternal retrieval . The Journal of neuroscience: the official journal of the Society for Neuroscience . 2013 ; 33 ( 8 ): 3276 ± 83 . Epub 2013/02/22. https://doi. org/10.1523/jneurosci.0425- 12 . 2013 PMID: 23426656
64. Reynolds CD , Nolan SO , Jefferson T , Lugo JN . Sex-specific and genotype-specific differences in vocalization development in FMR1 knockout mice . Neuroreport . 2016 ; 27 ( 18 ): 1331 ± 5 . Epub 2016/11/09. https://doi.org/10.1097/WNR.0000000000000701 PMID: 27824730
65. Rando OJ , Simmons RA. I' m eating for two: parental dietary effects on offspring metabolism . Cell . 2015 ; 161 ( 1 ): 93 ± 105 . Epub 2015/03/31. https://doi.org/10.1016/j.cell. 2015 . 02 .021 PMID: 25815988 .
66. Veres M , Duselis AR , Graft A , Pryor W , Crossland J , Vrana PB , et al. The biology and methodology of assisted reproduction in deer mice (Peromyscus maniculatus) . Theriogenology . 2012 ; 77 ( 2 ): 311 ± 9 . Epub 2011/09/20. https://doi.org/10.1016/j.theriogenology. 2011 . 07 .044 PMID: 21924468
67. Musolf K , Meindl S , Larsen AL , Kalcounis-Rueppell MC , Penn DJ . Ultrasonic Vocalizations of Male Mice Differ among Species and Females Show Assortative Preferences for Male Calls . PLoS One . 2015 ; 10 ( 8 ): e0134123 . Epub 2015 /08/27. https://doi.org/10.1371/journal.pone. 0134123 PMID: 26309246
68. Vandenberg LN . Non-monotonic dose responses in studies of endocrine disrupting chemicals: bisphenol a as a case study. Dose-response: a publication of International Hormesis Society . 2014 ; 12 ( 2 ): 259 ± 76 . Epub 2014/06/10. https://doi.org/10.2203/dose-response. 13 - 020 .Vandenberg PMID: 24910584
69. Cox ET , Hodge CW , Sheikh MJ , Abramowitz AC , Jones GF , Jamieson-Drake AW , et al. Delayed developmental changes in neonatal vocalizations correlates with variations in ventral medial hypothalamus and central amygdala development in the rodent infant: effects of prenatal cocaine . Behavioural brain research . 2012 ; 235 ( 2 ): 166 ± 75 . Epub 2012/08/08. https://doi.org/10.1016/j.bbr. 2012 . 07 .035 PMID: 22867871
70. Bowers JM , Perez-Pouchoulen M , Roby CR , Ryan TE , McCarthy MM . Androgen modulation of Foxp1 and Foxp2 in the developing rat brain: impact on sex specific vocalization . Endocrinology . 2014 ; 155 ( 12 ): 4881 ± 94 . Epub 2014/09/24. https://doi.org/10.1210/en.2014-1486 PMID: 25247470
71. Wetzel DM , Kelley DB , Campbell BA . Central control of ultrasonic vocalizations in neonatal rats: I. Brain stem motor nuclei . Journal of comparative and physiological psychology . 1980 ; 94 ( 4 ): 596 ± 605 . Epub 1980/08/01. PMID: 7410625 .
72. Curry T , Egeto P , Wang H , Podnos A , Wasserman D , Yeomans J . Dopamine receptor D2 deficiency reduces mouse pup ultrasonic vocalizations and maternal responsiveness. Genes, brain, and behavior . 2013 ; 12 ( 4 ): 397 ± 404 . Epub 2013/03/26. https://doi.org/10.1111/gbb.12037 PMID: 23521753 .
73. Mosienko V , Beis D , Alenina N , Wohr M. Reduced isolation-induced pup ultrasonic communication in mouse pups lacking brain serotonin . Molecular autism . 2015 ; 6 : 13 . Epub 2015/04/23. https://doi.org/10. 1186/s13229-015 -0003-6 PMID: 25901271