
==== Front
bioRxiv
BIORXIV
bioRxiv
2692-8205
Cold Spring Harbor Laboratory

39005445
10.1101/2024.07.02.601766
preprint
2
Article
The ultrasonic vocalization (USV) syllable profile during neonatal opioid withdrawal and a kappa opioid receptor component to increased USV emissions in female mice
http://orcid.org/0000-0002-9387-9063
Wingfield Kelly K. 12
Misic Teodora 1
Jain Kaahini 1
McDermott Carly S. 1
Abney Nalia M. 1
Richardson Kayla T. 13
Rubman Mia B. 4
http://orcid.org/0000-0001-6517-7614
Beierle Jacob A. 25
Miracle Sophia A. 16
Sandago Emma J. 1
Baskin Britahny M. 18
Lynch William B. 156
Borrelli Kristyn N. 1256
Yao Emily J. 1
Wachman Elisha M. 7
Bryant Camron D. 1*
1 Laboratory of Addiction Genetics, Center for Drug Discovery, Department of Pharmaceutical Sciences, Northeastern University, Boston, MA USA
2 T32 Biomolecular Pharmacology Training Program, Boston University Chobanian & Avedisian School of Medicine, Boston, MA USA
3 Post-Baccalaureate Research Education Program, Boston University Chobanian & Avedisian School of Medicine
4 NIH/NIDA Summer Undergraduate Fellowship Program
5 Transformative Training Program in Addiction Science, Boston University Chobanian & Avedisian School of Medicine
6 Graduate Program for Neuroscience, Boston University, Boston, MA USA
7 Department of Pediatrics, Boston University Chobanian & Avedisian School of Medicine, Boston Medical Center, Boston MA USA
8 T32 Training Program on Development of Medications for Substance Use Disorders Fellowship, Center for Drug Discovery, Northeastern University
Author Contributions

K.K.W.: Data collection, analysis, and writing of manuscript; T.M.: Data collection; K.J.: Data collection; C.S.M.: Data collection; N.M.A.; Data collection; K.T.R.: Data collection; M.B.R.; Data collection; J.A.B.: Data collection and RNAseq analysis; S.A.M.: Data collection; E.J.S.: Data collection; B.M.B.: Data collection; W.B.L: Gene pathway enrichment and analysis; K.N.B.; CFW Data collection; E.J.Y.; CFW Data collection; E.M.W.: Writing of manuscript; C.D.B.: Experimental design, analysis input, writing of manuscript.

* Corresponding Author: Camron D. Bryant, Ph.D., Laboratory of Addiction Genetics, Department for Pharmaceutical Sciences, Center for Drug Discovery, Northeastern University, 140 The Fenway, X138, Boston, MA USA 02115, c.bryant@northeastern.edu
01 9 2024
2024.07.02.601766https://creativecommons.org/licenses/by-nc-nd/4.0/ This work is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which allows reusers to copy and distribute the material in any medium or format in unadapted form only, for noncommercial purposes only, and only so long as attribution is given to the creator.
nihpp-2024.07.02.601766.pdf
Rationale

Opioid use during pregnancy can lead to negative infant health outcomes, including neonatal opioid withdrawal syndrome (NOWS). NOWS comprises gastrointestinal, autonomic nervous system, and neurological dysfunction that manifest during spontaneous withdrawal. Variability in NOWS severity necessitates a more individualized treatment approach. Ultrasonic vocalizations (USVs) in neonatal mice are emitted in isolation as a stress response and are increased during opioid withdrawal, thus modeling a negative affective state that can be utilized to test new treatments.

Objectives

We sought to identify the behavioral and USV profile, brainstem transcriptomic adaptations, and role of kappa opioid receptors in USVs during neonatal opioid withdrawal.

Methods

We employed a third trimester-approximate opioid exposure model, where neonatal inbred FVB/NJ pups were injected twice-daily with morphine (10mg/kg, s.c.) or saline (0.9%, 20 ul/g, s.c.) from postnatal day(P) 1 to P14. This protocol induces reduced weight gain, hypothermia, thermal hyperalgesia, and increased USVs during spontaneous morphine withdrawal.

Results

On P14, there were increased USV emissions and altered USV syllables during withdrawal, including an increase in Complex 3 syllables in FVB/NJ females (but not males). Brainstem bulk mRNA sequencing revealed an upregulation of the kappa opioid receptor (Oprk1), which contributes to withdrawal-induced dysphoria. The kappa opioid receptor (KOR) antagonist, nor-BNI (30 mg/kg, s.c.), significantly reduced USVs in FVB/NJ females, but not males during spontaneous morphine withdrawal. Furthermore, the KOR agonist, U50,488h (0.625 mg/kg, s.c.), was sufficient to increase USVs on P10 (both sexes) and P14 (females only) in FVB/NJ mice.

Conclusions

We identified an elevated USV syllable, Complex 3, and a female-specific recruitment of the dynorphin/KOR system in increased USVs associated with neonatal opioid withdrawal severity.

neonatal opioid withdrawal syndrome
kappa opioid receptor
ultrasonic vocalizations
sex differences
spectrotemporal profile
morphine
emotional-affective withdrawal
RNA-seq
transcriptome
brainstem
This study was funded by the National Institute on Drug Abuse (NIDA) grants U01DA050243 (C.D.B.), U01DA55299 (C.D.B.), T32DA055553 (B.M.B.)
==== Body
pmcINTRODUCTION

Opioid use during pregnancy is a major public health concern (Goetz et al., 2021; Haight, 2018; Hales et al., 2020; Ko et al., 2020; VanHouten et al., 2019). Exposure to opioids during gestation can lead to negative health outcomes for infants, such as preterm birth, growth restriction, and neonatal opioid withdrawal syndrome (NOWS). Approximately 59 infants in the United States are diagnosed with NOWS daily and are at an increased risk for admission to the neonatal intensive care unit (NICU), longer hospital stays, and greater hospital costs compared to other infants (Haight, 2018; Milliren et al., 2018; Strahan et al., 2020; Tolia et al., 2015; Wachman & Werler, 2019; Winkelman et al., 2018). NOWS consists of symptoms induced by spontaneous cessation of in-utero opioid exposure and is characterized by gastrointestinal irregularities (weight loss, poor feeding, diarrhea, vomiting), disrupted autonomous system (sweating, respiratory difficulties, sneezing) and neurological dysfunction (sleep disruptions, irritability, and excessive, high-pitched crying) (Jansson et al., 2009; Ko et al., 2017; Strahan et al., 2020). Current treatments for NOWS involve non-pharmacologic interventions, such as rooming-in and breastfeeding, as well as mu-opioid receptor agonist therapy (Abdel-Latif et al., 2006; Howard et al., 2017; Kraft et al., 2017; Wachman et al., 2018; Wachman & Werler, 2019). NOWS severity is highly variable, which can be attributed to multiple factors, including genetics, biological sex, type of opioid exposure, maternal polysubstance use, gestational age of exposure, and hospital care model (Liu et al., 2010; O’Connor et al., 2013; Seligman et al., 2008; Wachman et al., 2013, 2015). Thus, it is challenging to investigate the impact of individual factors on NOWS severity in a clinical setting. However, understanding the factors contributing to symptom severity can inform personalized treatment strategies to improve infant health outcomes.

Rodent models for NOWS traits are used to study factors influencing withdrawal symptom onset, duration, and severity (Byrnes & Vassoler, 2018; Ferrante & Blendy, 2024; Richardson et al., 2006). To effectively model NOWS traits in mice, we employ a third trimester-approximate opioid exposure paradigm, which is sufficient to induce a withdrawal state in pups (Borrelli et al., 2021; Dunn et al., 2023; Ferrante et al., 2022; Richardson et al., 2006; S. A. Robinson et al., 2020). P1–P14 in mice largely mirrors the neurodevelopmental events that occur during the third trimester of human pregnancy, including oligodendrocyte maturation, myelination, synaptogenesis, and establishment of the blood-brain barrier (Byrnes & Vassoler, 2018; Craig et al., 2003; Rice & Barone, 2000; Saunders et al., 2012; Semple et al., 2013). Other models include prenatal opioid exposure to the pups and dam; however, these models create challenges for accomplishing precise delivery of doses to individual pups, which is necessary for the study of factors (e.g., genetic) underlying individual differences in neonatal opioid withdrawal presence and severity. Prenatal models also necessitate maternal drug exposure, which, while valid for studying the human condition, complicates the understanding of innate individual differences underlying neonatal opioid withdrawal severity.

Here, we implemented a third trimester-approximate NOWS model (Borrelli et al., 2021; S. A. Robinson et al., 2020)in inbred FVB/NJ and outbred Swiss Webster Carworth Farms White (CFW) mouse pups to investigate the behavioral and transcriptional adaptations of chronic morphine exposure and withdrawal. We first used supervised machine learning to classify ultrasonic vocalization (USV) syllable profiles during spontaneous morphine withdrawal. Neonatal USVs are emitted exclusively in isolation to communicate discomfort and promote maternal attention (Branchi et al., 2001; Caruso et al., 2022; D’Amato et al., 2005; Scattoni et al., 2008). Thus, USVs can serve to model the negative affective state associated with morphine withdrawal. Next, we examined the transcriptomic changes in the brainstem (medulla and pons), a brain region containing multiple nuclei (e.g., locus coeruleus, rostroventral medulla) that contribute to neurobehavioral adaptations during opioid withdrawal, including enhanced nociception, dysphoria, and anhedonia (Basinger & Hogg, 2023; Bruchas et al., 2010; Downs & McElligott, 2022; Land et al., 2008; Shippenberg et al., 2007). We then investigated whether kappa opioid receptors (KOR), which are known to contribute to dysphoria during opioid withdrawal, contribute to USV emissions and syllable profiles during spontaneous morphine withdrawal. Together, our results identify a specific USV syllable, Complex 3, as a potential marker for the severity of withdrawal-associated dysphoria in mice, as well as a female-specific KOR component to increased USVs during spontaneous morphine withdrawal. These results support the extensive body of literature documenting sex differences in the dynorphin/KOR system in responses comprising negative affective-motivational states and emphasize the importance of biological sex in considering treatment strategies for opioid withdrawal (Becker & Chartoff, 2019; Conway et al., 2019; Fox & Sinha, 2009; Madurai et al., 2024).

METHODS

Mice

All experiments involving mice were conducted in accordance with the National Institutes of Health Guide for the Care and Use of Laboratory Animals and were approved by the Institutional Animal Care and Use Committees at Boston University and Northeastern University. FVB/NJ inbred mice were purchased at 8 weeks old from The Jackson Laboratory (Strain #001800). The FVB/NJ strain was chosen based on their excellent fecundity and the unexplored prospect of potentially exploiting closely related FVB substrains for mapping the genetic basis of NOWS model traits (Bryant et al., 2020). Outbred Cartworth Farm Webster (Swiss Webster) mice (8 weeks old) were purchased from Charles River Laboratories, and their historical dataset (Borrelli et al., 2021) was used to test for generalization of our findings, as they represent a commonly used and a relatively more diverse outbred genetic background. Mice were provided ad libitum laboratory breeding chow (Teklad 18% Protein Diet, Envigo) and tap water and maintained on a 12h light/dark cycle (lights on at 0630h). Breeders were paired in-house after one week of acclimation to the vivarium. All breeder cages contained nestlets. Sires were removed seven days after breeder pairing to strictly control for cage environment and avoid untimed pregnancies. Phenotyping occurred during spontaneous withdrawal (16h) between 0900h and 1100h and was performed by female experimenters to control for the effect of experimenter sex on rodent behavior (Driel & Talling, 2005; Sorge et al., 2014).

Tattooing of mice

Pup tails were tattooed (ATS-3 General Rodent Tattoo System, AIMS) on postnatal day (P) 7 following behavioral testing and morning injections for identification and were returned to their home cage.

Morphine administration in FVB/NJ pups from P1 to P15

We employed a third-trimester approximate mouse model of opioid exposure where injections of either morphine sulfate pentahydrate (10mg/kg, 20ul/g, s.c.; Sigma-Aldrich) or saline (0.9%, 20ul/g, s.c.) were administered twice-daily at 0900h and 1700h from P1 – P14, which is functionally equivalent to the third trimester of human pregnancy (Barr et al., 2011; Byrnes & Vassoler, 2018; Ferrante & Blendy, 2024; Richardson et al., 2006). Additional details on the rationale and procedure are provided in Supplementary Information. On P7 and P14, morning injections were administered following phenotyping during spontaneous morphine withdrawal at approximately 1100h.

Morphine administration in CFW pups from P1 to P14 (Borrelli et al., 2021)

Data analyzed in CFW mice was from a previous study that used a similar morphine treatment regimen, except a higher dose (15 mg/kg, s.c.). In the current study, we used a lower dose in FVB/NJ mice to minimize lethality. See Supplementary Information for additional details.

Recording of ultrasonic vocalizations (USVs)

Individual pups were each placed into a Plexiglass box (43 cm length × 20 cm width × 45 cm height; Lafayette Instruments) within a sound-attenuating chamber (Med Associates). USVs were recorded using the Ultrasound Recording Interface (Avisoft Bioacoustics UltrasoundGate 816H) for 10min (P7) or 15min (P14). Locomotor activity during all USV testing sessions was recorded using infrared cameras and tracked with ANY-maze software (Stoelting).

Thermal nociception testing on P7 and P14 in FVB/NJ pups

After USV recordings, each pup was removed from the sound-attenuating chamber and placed in a Plexiglass cylinder (diameter, 15 cm; height, 33 cm) on a 52.5°C hot plate (IITC Life Science). On P7, the nociceptive response was defined as the latency for the pup to roll onto its back (Borrelli et al., 2021). On P14, the nociceptive response was defined as the latency to jump, attempt to jump, hind paw lick or attempt to hind paw lick. Pups were removed from the hot plate immediately after observing a nociceptive response or after the 30s cut-off (P7) or the 60s cut-off (P14) if no response was observed. Upon completion of nociceptive testing, each pup was weighed and administered their morning injection and returned to their home cage.

Supervised ultrasonic vocalization classification

Labeled syllables and accompanying acoustic features were used to train a custom random forest classifier in Python. Additional classifier information is provided in the Supplementary Information.

Bulk RNA-seq of brainstem from morphine-withdrawn pups

The mice used for RNA-seq also underwent phenotyping on P8 and P15 while under the influence of morphine in the a.m. to assess the effect of alleviation of opioid withdrawal on USVs. Additional details are provided in Supplementary Information. Maintenance doses resumed for the p.m. injections on P8 and P15. During spontaneous withdrawal (16 h post-morphine) on P16, mice were sacrificed by live, rapid decapitation, and brains were removed from the skull. Brainstem tissue, including the pons, medulla and part of the spinal cord, was immediately collected and stored in RNAlater at 4°C for 5 days, then blotted dry and stored at −80°C for later RNA extractions. Further details on regional brain dissections, RNA extractions, RNA-seq, and analysis are provided in the Supplementary Information.

Real-time quantitative PCR (RT-qPCR) of Oprk1, Pdyn, and Slc6a3 in the brainstem and midbrain on P16

Brainstem and midbrain tissue for RT-qPCR was collected as described in the Supplementary Information. cDNA was synthesized using oligo-DT primers from total RNA using a cDNA reverse transcription kit (ThermoFisher Scientific, Cat#4368814). RT-qPCR was performed using PowerUP SYBR Green (ThermoFisher Scientific, Cat#A25741). The primer sequences are listed in the Supplemental Information.

Brainstem gene pathway enrichment analysis

Pathway enrichment analysis was performed on the lists of detectable genes for FVB/NJ and CFW mice. Pathways were obtained from a curated list across multiple different sources (Merico et al., 2010). Background genes were set to only include genes within our data sets. Gene were ranked by the absolute log10 p-values representing differential expression significance multiplied by the sign of the log2FC expression of each gene to preserve directionality. Pathway enrichment analysis was performed on subsequent gene rankings using the fgsea R package with an adaptive multilevel splitting approach and a minimum pathway size set to 15 genes and maximum set to 500 genes (Cérou et al., 2019; Korotkevich et al., 2016). Pathways were separated by those containing upregulated vs. downregulated genes. The collapsePathways command in fgsea was used to remove redundant/similar pathways that were significant (p < 0.05). The top 10 enriched pathways are listed by gene ratio (leading-edge genes/total genes in pathway), with circle size indicating number of leading-edge genes within each pathway and color indicating the pathway p value. Pathway clustering was performed using the treeplot function in enrichplot, using a ward.D clustering algorithm with the number of clusters set at 5 (Ward, 1963). The pathway circle size indicates the number of leading-edge genes within each pathway, and colors indicate the pathway p value. Figures were created in R studio (https://www.r-project.org/).

Effect of the kappa opioid receptor antagonist nor-BNI on ultrasonic vocalizations

A separate cohort of FVB/NJ mice was used to test the effects of the KOR antagonist nor-BNI (nor-Binaltorphimine Dihydrochloride; Sigma-Aldrich) on USVs. From P1–P14, injections of either morphine (10mg/kg, s.c) or saline (0.9%, 20 ul/g, s.c.) were administered twice-daily at 0900 h and 1700 h. Nor-BNI has an extremely slow onset of KOR antagonism and a prolonged duration of action at KOR in rodents (Kishioka et al., 2013; Marchette et al., 2021; Munro et al., 2012). Given the absence of pharmacokinetic studies of nor-BNI in neonatal mice, we chose the highest dose (30 mg/kg, s.c.) noted in adult rodent studies (CITE REF) to ensure a complete and prolonged KOR antagonism. Thus nor-BNI (30 mg/kg, s.c.) or saline 0.9%, (20 ul/g, s.c.) was administered at 1300 h on P14, approximately 4 h after the morning morphine or saline injections. Other studies waited between 5 h – 48 h before testing the effects of nor-BNI on rodent behavior, thus we allowed 20h before phenotyping on P15 to ensure we were capturing the KOR-specific effects of nor-BNI on USV emissions (Kelsey et al., 2015; Schlosburg et al., 2013; Walker et al., 2011). USVs were recorded at 0900h for 15 min.

Effect of the kappa opioid receptor agonist U50,488h on ultrasonic vocalizations

A separate cohort of morphine-naïve FVB/NJ mice was used to determine whether administration of the KOR agonist U50,488h ((±)-U-50488h hydrochloride; Fisher Scientific) would be sufficient to induce USVs mirroring the opioid withdrawal profile. From P1–P14, pups were injected with saline (0.9%, 20 ul/g, s.c.) twice daily at 0900 h and 1700 h to replicate the NOWS injection regimen. On P10 and P14, pups were injected with either U50,488h (0.625 mg/kg, s.c.) or saline (s.c.) and then returned to their home cage. We allowed 10 min before recording to ensure the effects of U50,488h on USV emission within previously reported time frames (Carden et al., 1994; Huang et al., 2020). USVs were recorded at 0900 h for 15 min. Pups continued to receive saline injections following P10 testing until testing for the behavioral effects of U50,488h again on P14, which is when we normally record USVs. See Supplemental Information for U50,488h dose-response pilot and P10 justification.

Statistical analysis

Analysis was performed in R (https://www.r-project.org/). All data are presented as the mean ± standard error of the mean (SEM), and p<0.05 was considered significant. Figures were created in GraphPad Prism (https://graphpad.com/). Body weight, temperature, USV locomotion, and USVs over time were analyzed using linear mixed models with Saline Treatment and Female Sex as the reference variables and Pup as a random effect (repeated measures). Sex was removed from the model if there were no interactions. Significant interactions of interest were followed up with least-square means (predicted marginal means) and Tukey’s Honestly Significant Difference (HSD) tests. All other data were analyzed using linear models with Saline Treatment and Female Sex as the reference variables. For nor-BNI experiments, we were mostly interested in the effect of nor-BNI on USV emissions in morphine-treated females, thus Morphine Treatment (P1 – P14 treatment), Saline Treatment (pre-treatment) and Female sex were used as reference variables. Based on the RNA-seq results, we had a specific hypothesis for the direction of fold-change in gene expression, therefore, for RT-qPCR analysis, one-tailed t-tests were performed. For Slc6a3, one sample was identified as an outlier (ROUT (Q = 1%)) and was removed from the analysis.

RESULTS

Morphine exposure from P1–P14 is sufficient to induce multiple signs of spontaneous opioid withdrawal in neonatal FVB/NJ pups

An experimental timeline is provided in Fig.1a. Reduced body weight and poor body temperature regulation are commonly observed in infants diagnosed with NOWS (Finnegan et al., 1975; Hudak et al., 2012; Jansson & Velez, 2012; Patrick et al., 2020). Twice-daily morphine injections (10 mg/kg, s.c.) from P1–P14 significantly reduced body weight (Fig.1b) and temperature during withdrawal (Fig.1c). Hyperalgesia is common in opioid withdrawal (Corder et al., 2018). Morphine-withdrawn pups displayed thermal hyperalgesia on P7, as indicated by a reduced latency to elicit a nociceptive response on the hot plate (Fig.1d). In examining hot plate velocity on P7 during the time interval of exposure, there was no effect of morphine treatment, indicating that reduced locomotor activity did not prolong hot plate exposure and thus did not confound the group difference (Fig.1e). On P14, morphine-withdrawn pups displayed reduced hot plate latency (Fig.1f) and increased hot plate velocity (Fig.1g). Because increased velocity decreases paw exposure to the hot plate, we can again conclude that increased locomotor activity did not confound hyperalgesia in morphine-withdrawn pups.

USV locomotor activity of FVB/NJ pups during spontaneous morphine withdrawal on P7 and P14

On P7, morphine-withdrawn pups traveled a greater distance during the first 2–5min of the USV recording session compared to saline pups (Fig.2a). For the entire 10 min, morphine-withdrawn pups traveled a greater distance versus saline-treated pups (Fig.2b). There was no effect of morphine withdrawal on P7 average velocity (Fig.2c). On P14, there was no significant difference in total distance traveled over 15 min (Fig.2d) nor average velocity between treatment groups (Fig.2e).

USV syllable profile of FVB/NJ pups during spontaneous morphine withdrawal on P14

Neonatal USVs communicate negative internal states. Thus, we measured USVs as an indicator of enhanced distress during spontaneous morphine withdrawal (Branchi et al., 2001; D’Amato et al., 2005; Ehret, 2005). On P7, there was a significant Morphine Treatment × Time interaction (Fig.S1a); however, there was no overall effect of Morphine Treatment on the total number of USVs over 10 min (Fig.S1b). On P14, there was a time-dependent increase in vocalizations in morphine-withdrawn female pups (Fig.3a). Overall, Morphine-withdrawn females emitted significantly more USVs compared to saline-treated females, whereas there was no Treatment effect in males (Fig.3b).

USVs are classified into syllable types based on their acoustic features (frequency, duration) (Grimsley et al., 2011; Heckman et al., 2016; Portfors, 2007). We used a random forest model to classify USV syllables based on these known classifications (Table.S1) during neonatal morphine withdrawal. On P7, morphine-withdrawn pups emitted fewer Complex 3 syllables than saline pups (Fig.S1c). On P14, morphine-withdrawn pups showed the opposite, namely a robust increase in the percentage of Complex 3 syllables (Fig.3c). The increase in Complex 3 syllables was completely driven by morphine-withdrawn females, as there was no treatment difference in males (Fig.3d,e). Morphine-withdrawn females also showed a decrease in the percentage of Flat, Short, and Upward syllables, whereas males only showed an increase in the percentage of Complex syllables (Fig.3d,e).

The identification of increased Complex 3 syllables during spontaneous opioid withdrawal was a robust and novel, female-selective observation. To determine whether increased Complex 3 emissions in females generalize to another genetic background, we also examined USV syllable subtypes in morphine-withdrawn outbred CFW stock mice from a historical dataset (Borrelli et al., 2021) comprising a nearly identical morphine regimen, except with a higher twice daily dose (15 mg/kg, s.c.; see Supplementary Information). On P7, treatment had no effect on the USV syllable profile (Fig.S2a). On P14, the USV profile in morphine-withdrawn CFW pups showed a similar increase in the percentage of Complex 3 syllables (Fig.S2b), except in this case, both sexes showed this effect (perhaps due to the higher historical morphine dose that was used in the CFW study (Borrelli et al., 2021)). Additional similarities between morphine-withdrawn CFW and FVB pups included a decrease in the percentage of Flat and Upward syllables on P14 (Fig.S2b).

Brainstem transcriptomic analysis identified an upregulation of the kappa opioid receptor

In a subset of FVB/NJ mice, we administered an a.m. maintenance dose of morphine (10 mg/kg) prior to USV recording on P8 and P15 to evaluate the semi-acute effects of morphine exposure on the USV syllable profile (Fig.S3). We collected brainstem tissue during spontaneous withdrawal on P16 (16 h post-final morphine injection) to evaluate the effects of chronic morphine exposure and withdrawal on gene expression (Fig.4a). A schematic of regional brain dissections is provided in Fig.S4. Transcriptomic analysis via RNA-seq revealed a significant upregulation of Oprk1 (log2FC=0.36, **p=0.0049), which codes for the kappa opioid receptor (KOR) and Slc6a3 (dopamine transporter (DAT); log2FC=3.69, **p=0.0067) during spontaneous morphine withdrawal (Table.S2). We also observed a downregulation of genes associated with myelination, such as Mbp (myelin basic protein; log2FC= 0.41, ****p<0.0001) and Plp1 (proteolipid protein 1; log2FC=−0.38, ****p<0.00032) (Table.S3).

We performed RT-qPCR on brainstem and midbrain tissue to validate these mRNA transcripts with a larger sample size. Consistent with our brainstem RNA-seq results, we observed withdrawal-induced upregulation of Oprk1 (Fig.4b), no effect on Pdyn expression (Fig.4c), and upregulation of Slc6a3 (Fig.4d). Interestingly, for the midbrain, we did not observe a withdrawal effect on Oprk1 expression (Fig.4e). However, we found upregulation of Pdyn (Fig.4f). There was no effect of withdrawal on Slc6a3 (Fig.4g) expression.

Enrichment analysis of upregulated genes in the brainstem revealed processes associated with synaptic organization, structure, and activity (Fig.5a), and enriched pathways of downregulated genes associated with extracellular matrix organization and glial development (Fig.5b).

The kappa opioid receptor antagonist nor-BNI reduced the number of USVs in females during spontaneous morphine withdrawal

KOR activation contributes to the dysphoric aspects of withdrawal, such as hyperirritability, anhedonia, and anxiety (Land et al., 2008; Margolis & Karkhanis, 2019; Shippenberg et al., 2007). Furthermore, sex differences in the dynorphin/KOR system on motivated behaviors in the context of drug addiction are well-documented in rodents and humans (Chartoff & Mavrikaki, 2015). Therefore, we aimed to determine whether the dynorphin/KOR system contributes to increased USVs and Complex 3 emissions during opioid withdrawal. On P15, 20 h post-administration of the KOR antagonist, nor-BNI, there was a significant decrease in USVs in morphine-withdrawn females but no effect of nor-BNI on USVs in males (Fig.6a). Furthermore, there was a reduction in USVs over time in morphine-treated females, but not males, pre-treated with nor-BNI compared to morphine controls (Fig.6b). These results identify a KOR contribution to the female-specific enhancement of USVs in morphine-withdrawn pups. However, KOR activation was not specific for a particular syllable type, as nor-BNI had no effect on the overall USV syllable profile (Fig.6c).

The kappa opioid receptor agonist U50,488h is sufficient to induce an increase in USVs

Previous studies reported that the administration of KOR agonists, such as U50,488h, is aversive and can induce anhedonia and decreased reward sensitivity (Conway et al., 2019; Huang et al., 2020; Vonvoigtlander et al., 1983). Here, we wished to determine whether KOR activation with the KOR agonist U50,488h was sufficient to increase USVs and induce a syllable profile similar to morphine-withdrawn pups, particularly in females. To conduct this experiment, we considered a postnatal day where saline-treated controls would show sufficient baseline USVs to detect either an increase or decrease in USVs. In neonates, the rate of USVs peaks around P7–P10 and declines during the second postnatal week (Elwood & Keeling, 1982). Thus, we tested the KOR agonist U50,488h on P10 and then again on P14, when we normally assess USVs (Fig.3). Based on a pilot study, the U50,488h dose of 0.625mg/kg was chosen because it induced an increase in USVs without affecting locomotor activity (Fig.S5a,b). On P10, U50,488h increased USVs relative to saline controls during the first 10 min of the recording session (Fig.7a). Additionally, in considering the entire 15 min, U50,488h increased overall USVs in both sexes (Fig.7b). In contrast to our original hypothesis (although in line with the null effect of the KOR antagonist in morphine-withdrawn mice), there was no effect of KOR activation on the percentage of Complex 3 emissions on P10 in either sex (Fig.7c). In considering other syllable types, there was a significant decrease in the percentage of Downward and Short emissions (Fig.7c). On P14, there was a time-dependent effect of U50,488h, whereby KOR activation increased USVs during the 5–9 min of the recording session (Fig.7d). Furthermore, on P14, U50,488h increased USVs in females, but not in males (Fig.7e). In contrast to our original hypothesis, U50,488h decreased the percentage of Complex 3 syllables in both sexes (Fig.7f). U50,488h also decreased the emission of other multi-component syllables, including Complex 2 and Complex 4 in both sexes (Fig.7f).

DISCUSSION

Consistent with previous studies, twice-daily injections of morphine from P1 to P14 was sufficient to induce withdrawal-associated model traits in FVB/NJ pups on P14, including low body weight, hypothermia, hyperalgesia, and increased USVs on P14 (Borrelli et al., 2021; Carden et al., 1991; S. A. Robinson et al., 2020). Neonatal mice emit USVs exclusively in isolation as a distress signal to communicate negative internal states and promote maternal attention (Branchi et al., 2001; D’Amato et al., 2005; Dirks et al., 2002; Ehret, 2005). Rodent studies observed increased USVs during opioid withdrawal in neonatal pups (mice and rats) and adult rats (Barr & Wang, 1992; Borrelli et al., 2021; S. A. Robinson et al., 2020; Vivian & Miczek, 1991). On P14, we observed increased USVs, specifically in FVB/NJ female pups during spontaneous morphine withdrawal.

The severity of NOWS and the need for treatment are determined by observer-rated scales using the Finnegan Neonatal Abstinence Scoring System (FNASS) (Finnegan et al., 1975). Thus, NOWS assessment is subjective and may contribute to unnecessary treatment interventions and lengthy hospital stays (Jansson et al., 2009; Milliren et al., 2018; Timpson et al., 2018). Therefore, unbiased methods for assessing withdrawal symptoms may lead to better prediction of optimal treatment plans, improved infant outcomes, and reduced hospital costs. NOWS infants present with excessive, high-pitched cries that can be differentiated from non-NOWS cries (Blinick et al., 1971). Again, the severity of NOWS cries is subjective and can be interpreted differently between listeners. A recent study revealed that infants with NOWS emit cries with distinct features unrecognizable by humans, such as altered frequency formants, utterances, and amplitude (Manigault et al., 2022). Thus, alterations in acoustic features of cries may be associated with the severity of NOWS, which could be useful for predicting the need for pharmacological treatment compared to subjective evaluation.

A major goal of this study was to identify the USV spectrotemporal syllable profile associated with the negative affective state of neonatal morphine withdrawal. For P14, we identified a robust increase in Complex 3 syllables in mice from two genetic backgrounds, implicating Complex 3 syllables as a potential preclinical biobehavioral marker for symptomatic severity of the negative affective state associated with morphine withdrawal. Furthermore, morphine-withdrawn FVB/NJ females emitted a greater proportion of Complex 3 syllables compared to morphine-withdrawn males, suggesting female pups experience a more severe affective withdrawal state. However, sex differences in Complex 3 emissions during neonatal opioid withdrawal may depend on genetic background as both female and male outbred CFW pups emitted more Complex 3 syllables compared to saline-treated pups. This female-specific effect potentially depends on genetic background, as morphine-withdrawn CFW mice of both sexes showed increased USVs. Alternatively, it should be noted that the lack of sex difference in CFW mice could potentially be explained by the FVB/NJ pups receiving a lower dose of morphine (10mg/kg) rather than genetic background. For example, there could be impaired morphine metabolism and clearance at 10 mg/kg that could underlie the female-specific increase in Complex 3 syllables, whereas the higher dose in CFWs (15 mg/kg) could override this pharmacokinetic effect on behavior.

We analyzed the transcriptome of the brainstem and midbrain tissue due to its role in morphine withdrawal and the production of vocalizations in rodents. On P16 during morphine withdrawal, there was an upregulation of Oprk in the brainstem, which codes for kappa opioid receptors (KOR), and increased expression of Pdyn, which codes for dynorphin, in the midbrain. KORs and their endogenous ligand, dynorphin (Mansour et al., 1995), mediate withdrawal-associated dysphoria in rodents and humans (Chartoff & Mavrikaki, 2015; Land et al., 2008; Leconte et al., 2022; Shippenberg et al., 2007). Studies have shown that projections from the periaqueductal grey (PAG) in the midbrain to the nucleus retroambiguus (RAm) of the brainstem are necessary for USV production in adult mice (Jürgens, 2002, 2002; Tschida et al., 2019). Therefore, increased dynorphin release in the midbrain may activate KOR in the brainstem and contribute to enhanced withdrawal-induced dysphoria and USVs. We assessed KOR antagonism as a potential target for reducing negative affect in the morphine-withdrawn state. Pre-treatment with the KOR antagonist, nor-BNI on P14 (a known postnatal day for morphine withdrawal-induced increase in USVs), decreased USVs exclusively in morphine-treated females during withdrawal on P15 (20 h pre-treatment) but had no effect on the composition of syllable emissions. This result suggests a female-specific component of KOR in mediating USV emissions and by extension, the severity of the aversive state, but it does not account for the increase in Complex 3 syllables observed during withdrawal. Studies in adult rats have shown that the KOR antagonist nor-BNI reduces withdrawal symptoms (wet dog shakes, jumps, weight loss, and conditioned place aversion) in morphine-dependent rats (Kelsey et al., 2015), attenuates ethanol withdrawal-induced decreases in immobility during the Forced Swim Test (Jarman et al., 2018), and blocks the decrease in open arm explorations in the elevated plus maze during ethanol withdrawal (Valdez & Harshberger, 2012). Most importantly, nor-BNI can block the increase in aversive 22-kHz USVs in rats during ethanol withdrawal and following administration of U50,488h (Berger et al., 2013). Administration of KOR agonist U50,488h was sufficient to increase USVs on P10 in both sexes and on P14, specifically in females. Rat studies have observed enhanced USVs in neonates following U50,488h administration (Carden et al., 1991, 1994; Nazarian et al., 1999). However, few, if any, of these studies have described the spectrotemporal syllable profile. In contrast to our prediction, U50,488h induced a decrease (rather than an increase) in Complex 2, −3, and −4 syllables on P14. It is likely that exogenously administered U50,488h is activating additional circuitry beyond that which is driving increased USVs during morphine withdrawal (e.g., sensory, affective, and/or motor circuitry) that underlies the reduction in Complex 3 emissions. Together, these data are consistent with a female-specific role of KOR activation on P14–P15 in mediating overall USVs but not the type of syllables emitted, indicating that changes in syllable profiles during morphine withdrawal involve additional mechanisms.

Interestingly, there was a significant upregulation of Slc6a3 (dopamine (DA) transporter, DAT) in the brainstem. Activation of KOR at DAergic terminals inhibits DA release and can promote negative affective states that drive drug-seeking behavior (Shippenberg et al., 1993; Zachry et al., 2021). Accordingly, increased DAT expression would be expected to increase synaptic DA uptake and further support a negative affective state. In neonatal rats, KOR activation with U50,488h increased USV emissions (Carden et al., 1991, 1994; Nazarian et al., 1999), which can be attenuated by DA receptor agonists (Nazarian et al., 1999). Although DAergic neurons in the midbrain do not project to lower parts of the the brainstem (Reyes et al., 2009), activation of KOR on DAergic projections in other regions, such as the ventral tegmental area (VTA), to forebrain structures, including nucleus accumbens (Margolis et al., 2003, 2006; Speciale et al., 1993), may contribute to reduced brainstem DA levels. Thus, increased Oprk1 and Slc6a3 expression in the brainstem together could reflect changes in expression in, e.g., DAergic forebrain terminals and together contribute to the aversive state during withdrawal and mediate enhanced USV emission.

Sex differences in reward and stress-related pathways mediated by the dynorphin/KOR system have been observed in rodents and humans (Becker & Chartoff, 2019; Conway et al., 2019; Russell et al., 2014; Vijay et al., 2016). In adult rats, U50,488h-mediated KOR activation decreases reward motivation and induces greater thermal analgesic effects in females compared to males (Bartok & Craft, 1997; Conway et al., 2019; Russell et al., 2014). Human studies have also observed increased analgesic efficacy of KOR agonists in females versus males (Gear et al., 1996; Mogil et al., 2003). Female rats are less sensitive to the DA-reducing and depressive-like effects of KOR agonists (Conway et al., 2019; Laman-Maharg et al., 2018; Russell et al., 2014). These discrepant results in adult rodents versus neonatal females in our study could be due to differences in the immature and rapidly evolving circuitry in neonates versus adults. Rodent studies have shown spatiotemporal variations in KOR expression during development in rats and mice (Georges et al., 1998; Tan et al., 2018; Zhu et al., 1998). Additionally, adult human studies have shown sex differences in KOR availability and expression (Vijay et al., 2016). Time-dependent sex differences in the dynorphin/KOR circuitry might also explain why U50,488h increased USVs in both sexes on P10 but only in females on P14.

There are several limitations to this study. First, the third trimester-approximate model involves neonatal injections, potentially affecting maternal care quality between morphine- and saline-treated pups. There is conflicting evidence regarding the effects of opioid exposure (both to the dam and pups) on maternal care behavior, although evidence indicates that it is dependent on maternal opioid exposure rather than pup exposure (Alipio et al., 2021; Grecco et al., 2021; Smith et al., 2022). In our model, dams were not exposed to morphine – thus, behavioral differences between morphine- and saline-treated pups are most likely due to morphine withdrawal and not maternal care deficits. An additional limitation is that we used the KOR antagonist, nor-BNI, which has a delayed onset of effect (Kishioka et al., 2013; Munro et al., 2012), requiring modifications to our NOWS paradigm. Nor-BNI also has early mu-opioid receptor antagonistic effects, which could induce withdrawal shortly after administration. Nonetheless, future studies could implement faster-acting KOR antagonists, such as CERC-501(LY2456302) (Rorick-Kehn et al., 2014) to assess USVs on P14 rather than P15. However, CERC-501 also binds to mu-opioid receptors, has a lower binding affinity for KOR than nor-BNI and has long-acting effects at KOR (Rorick-Kehn et al., 2014). Additionally, using constitutive or Cre-dependent KOR-knockout mouse models could further elucidate the relationship between the dynorphin/KOR system and USVs during neonatal morphine withdrawal.

We observed robust sex differences in USV profiles in a third trimester-approximate treatment regimen for NOWS model traits. Notably, we identified a unique USV profile, as characterized by an increase in the proportion of Complex 3 syllables, associated with morphine withdrawal. We propose enhanced Complex 3 emissions as a novel marker for withdrawal-induced dysphoria, which can serve as a new phenotype for testing the efficacy of novel therapeutics in mice. We identified a KOR component of mediating USV emissions in morphine-withdrawn FVB/NJ, suggesting that the KOR system in females may contribute to greater withdrawal distress and dysphoria compared to males. These results contribute to our understanding of the biological mechanisms underlying sex differences in NOWS symptom severity.

More broadly, we provide a quantitative, machine learning-based approach to studying negative internal states associated with neonatal opioid withdrawal in mice, identifying specific USV profiles as indicators of opioid withdrawal. Examining adaptations in spectrotemporal features of neonatal mouse vocalizations unique to neonatal opioid withdrawal provides a translational approach for modeling NOWS severity beyond excessive crying in the clinic and assessing potential experimental therapeutics in reducing emotional symptom severity.

Supplementary Material

Supplement 1

ACKNOWLEDGMENTS

We thank the Boston University Genome Science Institute and Microarray and Sequencing Resource Core for the Illumina mRNA-seq library prep kit and flow cell awarded to Jacob A. Beierle.

Funding

This study was funded by the National Institute on Drug Abuse (NIDA) grants U01DA050243 (C.D.B.), U01DA55299 (C.D.B.), T32DA055553 (B.M.B.)

Fig. 1. Morphine exposure from P1–P14 is sufficient to induce opioid withdrawal traits in neonatal FVB/NJ pups.

Data are plotted as the mean ± SEM. Saline = blue lines/bars; Morphine = orange lines/bars. Closed circles = Females; Open circles = males. (a) Experimental Timeline. (b) Body Weight: The effect of Morphine Treatment on body weight was dependent on Postnatal day (β = −0.30, SE = 0.0055, t(611) = −54.38, ****p < 0.0001), where morphine-treated pups weighed significantly less than saline-treated pups from P3 (**p = 0.0013) – P14 (P4 – P14: ****p < 0.0001). (c) Temperature: The effect of Morphine Treatment on body temperature was dependent on Postnatal day (β = −0.087, SE = 0.041, t(611) = −2.12, *p = 0.035), where morphine-treated pups displayed hypothermia compared to saline-treated pups from P3–P11 and P13 – P14 (all *p ≤ 0.026). SAL, n = 25 (11F, 14M); MOR, n = 22 (11F, 11M). (d) P7 Hot Plate Latency: There was no effect of Sex (β = 0.49, SE = 2.64, t(46) = 0.19, p = 0.85) or a Morphine Treatment × Sex interaction (β = 0.26, SE = 3.81, t(46) = 0.068, p = 0.95). The simplified model revealed that Morphine Treatment was associated with a decrease in hot plate latency compared to saline controls (β = −12.88, SE = 1.86, t(48) = −6.91, ****p < 0.0001). (e) P7 Hot Plate Velocity: There was no effect of Sex (β = 0.00037, SE = 0.005, t(39) = 0.076, p = 0.94) or a Morphine Treatment × Sex interaction (β =0.0030, SE = 0.0069, t(39) = 0.43, p = 0.67). The simplified model revealed that Morphine Treatment had no effect on hot plate velocity (β = 0.001, SE = 0.0034, t(41) = 0.27, p = 0.79). SAL, n = 23–26 (11–12F, 12–14M); MOR, n = 22–24 (12–13F, 10–11M). (f) P14 Hot Plate Latency: There was no effect of Sex (β = 1.089, SE = 6.085, t(35) = 0.18, p = 0.86) or a Morphine Treatment × Sex interaction (β = 7.33, SE = 8.91, t(35) = 0.82, p = 0.67). The simplified model revealed that Morphine Treatment was associated with decreased hot plate latency compared to saline controls (β = −15.50, SE = 4.42, t(37) = −3.51, **p = 0.0012). (g) P14 Hot Plate Velocity: There was no effect of Sex (β = −0.0047, SE = 0.0041, t(33) = −1.15, p = 0.23) or a Morphine Treatment × Sex interaction (β = 0.0094, SE = 0.0059, t(33) = 1.58, p = 0.12). The simplified model revealed that Morphine Treatment was associated with increased hot plate velocity compared to saline controls (β = 0.011, SE = 0.0030, t(35) = 3.67, ***p = 0.00081). SAL, n = 20–21 (8–9F, 12M); MOR, n = 17–18 (9F, 8–9M)

Fig. 2. Locomotor activity during USV assessment of FVB/NJ pups undergoing spontaneous morphine withdrawal on P7 and P14.

Data are plotted as the mean ± SEM. Saline = blue lines/bars; Morphine = orange lines/bars. Closed circles = Females; Open circles = males. (a) P7 Distance: The effect of Morphine Treatment on USV distance was dependent on Time (β = −0.030, SE = 0.0056, t(378) = −5.34, ****p < 0.0001), where morphine-treated pups traveled a greater distance than saline controls from 2 – 5 min of the recording session (all *padj ≤ 0.031). (b) P7 Total Distance: There was no effect of Sex (β = −0.11, SE = 0.58, t(38) = −0.20, p = 0.85) or a Morphine Treatment × Sex interaction (β = −0.73, SE = 0.80, t(38) = −0.92, p = 0.36). The simplified model revealed that Morphine-treated pups traveled a greater distance during USV recordings compared to saline-treated pups (β = 1.24, SE = 0.40, t(40) = 3.12, **p = 0.0033). (c) P7 Average Velocity: There was no effect of Sex (β = 0.0022, SE = 0.0040, t(38) = 0.59, p = 0.56) or a Morphine Treatment × Sex interaction (β = −0.0020, SE = 0.0055, t(38) = −0.37, p = 0.72). The simplified model revealed that Morphine Treatment had no effect on USV velocity (β = −0.0031, SE = 0.0027, t(40) = −1.17, p = 0.25). SAL, n = 21 (9F, 12M); MOR, n = 21 (11F, 10M). (d) P14 Distance: There was no Morphine Treatment × Time interaction (β = 0.0029, SE = 0.012, t(518) = 0.24, p = 0.81), so data was collapsed across the 15 min recording session. There was no effect of Sex (β = −0.93, SE = 3.84, t(33) = −2.42, p = 0.81) or a Morphine Treatment × Sex interaction (β = −0.26, SE = 5.83, t(33) = −0.045, p = 0.96). The simplified model revealed that Morphine Treatment had no effect on the total distance traveled (β = 5.20, SE = 2.81, t(35) = 1.85, p = 0.073). (e) P14 Average Velocity: There was no effect of Sex (β = 0.019, SE = 0.016, t(33) = 1.20, p = 0.24) or a Morphine Treatment × Sex interaction (β =−0.020, SE = 0.024, t(33) = −0.82, p = 0.42). The simplified model revealed that Morphine Treatment had no effect on USV velocity (β = 0.00020, SE = 0.012, t(35) = 0.017, p = 0.99). SAL, n = 21 (9F, 12M); MOR, n = 13 (7F, 9M)

Fig. 3. USV profiles of FVB/NJ pups during spontaneous morphine withdrawal on P14.

Data are plotted as the mean ± SEM. Saline = blue lines/bars; Morphine = orange lines/bars. Closed circles = Females; Open circles = males. (a) P14 USV Emissions: There was a Morphine Treatment × Sex × Time interaction (β = −3.60, SE = 1.03, t(518) = −3.49, ***p = 0.0052), where morphine females vocalized more than saline females over time, specifically during the 4 – 15 min intervals (all *padj ≤ 0.035). (b) P14 Total USV Emissions: The effect of Morphine Treatment on USV emission was dependent on Sex (β = 919.6, SE = 180.7, t(33) = 5.088, ****p < 0.0001). Specifically, morphine-treated females vocalized more than saline females (β = 937.0, SE = 181, t(33) = 5.19, ****p < 0.0001). (c) P14 Syllable Profile: Morphine Treatment was associated with an increase in the percentage of Complex 3 syllables emitted (β = 0.14, SE = 0.033, t(35) = 4.19, ***p = 0.00018), and a decrease in the percentage of Flat (β = −0.051, SE = 0.017, t(35) = −3.0, **p = 0.0051), and Short (β = −0.030, SE = 0.014, t(35) = −2.08, *p = 0.045) syllables emitted compared to saline controls. There were Morphine Treatment × Sex interactions for Complex 3 (β = 0.15, SE = 0.060, t(33) = 2.54, *p = 0.016) and Upward (β = 0.14, SE = 0.044, t(33) = 3.13, **p = 0.0036) syllables, so we investigated the syllable profile for each sex separately. (d) Female P14 Syllable Profile: Morphine Treatment was associated with an increase in the percentage of Complex 3 syllables (β = 0.21, SE = 0.046, t(15) = 4.8, ***p = 0.00030) and a decrease in Flat (β = −0.051, SE = 0.020, t(15) = −2.55, *p = 0.022), Short (β = −0.055, SE = 0.020, t(15) = −2.69, *p = 0.017), and Upward syllables (β = −0.097, SE = 0.031, t(15) = −3.13, **p = 0.0069). (e) Male P14 Syllable Profile: Morphine Treatment was associated with an increase in the percentage of Complex syllables (β = 0.017, SE = 0.0072, t(18) = 2.41, *p = 0.027). SAL, n = 20 (8F, 12M); MOR, n = 17 (9F, 8M)

Fig. 4. Differentially expressed genes in the brainstem and midbrain during spontaneous morphine withdrawal on P16.

(a) Brainstem RNAseq: Data reflects the effect of morphine exposure on gene expression relative to saline pups (sex-collapsed; Saline, n = 4 (2F, 2M); Morphine, n = 4 (2F, 2M). Log2 Fold Change (FC) represents gene expression, and Log10 p values (unadjusted) reflect significance. To avoid crowding, we labeled genes with an absolute (abs) log2FC ≥ 1.5 and/or Log10P ≥ 4.0. Dots are individual genes. Grey = absolute (abs) log2FC ≥ 0.26; non-significant (unadjusted p < 0.05; NS). Green = abslog2FC ≥ 0.26, non-significant; Orange = abslog2FC ≤ 0.26, significant (unadjusted p < 0.05); Pink = abslog2FC ≥ 0.26, significant (unadjusted p < 0.05). RT-qPCR: Data are plotted as the mean ± SEM. Saline = blue bars; Morphine = orange bars. Closed circles = Females (SAL, n = 3–4; MOR, n = 3–4); Open circles = males (SAL, n = 6–8; MOR, n = 5–7). (b) Brainstem Oprk1: Morphine Treatment was associated with upregulation of Oprk1 (t(17) = 2.06, *p = 0.0028). (c) Brainstem Pdyn: There was no effect of Morphine Treatment on Pdyn expression (t(17) = 1.35, p = 0.097). (d) Brainstem Slc6a3: Morphine Treatment was associated with upregulation of Sl6a3 (t(16) = 2.00, *p = 0.031). SAL, n = 10 (3F, 6 – 7M); MOR, n = 9 (4F, 5M). (e) Midbrain Oprk1: Morphine Treatment had no effect on Oprk1 expression (t(20) = 0.066, p = 0.47). (f) Midbrain Pdyn: Morphine Treatment was associated with upregulation of Pdyn (t(20) = 3.47, *p = 0.0012). (g) Midbrain Slc6a3: There was no effect of Morphine Treatment on Slc6a3 expression (t(20) = 0.37, p = 0.64)

Fig. 5. Pathway enrichment analysis of differentially expressed genes in the brainstem during spontaneous morphine withdrawal on P16. Colors (pink, yellow, teal, green, blue) correspond to pathway names. The size of the circles reflects the number of genes in the pathway (the larger the circle, the more genes are associated with the pathway). The color of the circles reflects the significance (decreasing p value (more significant) from blue to red).

(a) Upregulated Pathways: clustering of the top 20 enriched pathways consisting of upregulated genes in FVB/NJ mice. (b) Downregulated Pathways: Clustering of the top 20 enriched pathways consisting of downregulated genes in FVB/NJ mice. Clustering was determined by calculating Jaccard similarity coefficients between pathways, with the number of clusters set to 5

Fig. 6. The kappa opioid receptor antagonist nor-BNI reduces USV emissions in FVB/NJ female pups during spontaneous withdrawal on P15.

Data are plotted as the mean ± SEM. Saline animals pre-treated with saline = blue lines/bars; Saline animals pre-treated with nor-BNI = yellow lines/bars; Morphine animals pre-treated with saline = orange lines/bars; Morphine animals pre-treated with nor-BNI = purple lines/bars. Closed circles = Females; Open circles = males. (a) P15 Total USV Emissions: There was a Morphine Treatment × nor-BNI Treatment × Sex interaction (β = −874.5, SE = 421.2, t(143) = −2.076, *p = 0.040), where Morphine-treated females pre-treated with nor-BNI (MOR/norBNI) emitted fewer USVs compared to morphine-treated females pre-treated with saline (MOR/SAL) (β = −272.35, SE = 127.11, t(143) = −2,143 *p = 0.034). (b) Female P15 USV Emissions: Given the three-way interaction, we analyzed females and males separately across time to simplify the linear mixed model. There was a time-dependent decrease in USVs in morphine-treated females pre-treated with nor-BNI (MOR/norBNI) compared to morphine-treated females pre-treated with saline (MOR/SAL) (β = −1.42, SE = 0.52, t(1148) = −2.72, **p = 0.0068), specifically during the 4–8 and 10 min intervals (all *padj ≤ 0.037). (c) P15 Syllable Profile: There was no effect of Sex on any Syllable emitted (all p ≥ 0.20), so Sex was removed from the models. Morphine pups pre-treated with nor-BNI emitted fewer Chevron syllables compared to morphine-treated pups pre-treated with saline (β = −0.036, SE = 0.017, t(147) = −2.065 *p = 0.041). SAL/SAL, n = 40 (22F, 18M); SAL/norBNI, n = 33 (17F, 16M); MOR/SAL, n = 39 (21F, 18M); MOR/norBNI, n = 39 (22F, 17M)

Fig. 7. The kappa-opioid receptor agonist U50,488h, increases USV emission in morphine-naïve FVB/NJ pups on P10 and P14.

Data are plotted as the mean ± SEM. Saline = blue lines/bars; U50,488h = green lines/bars. Closed circles = Females; Open circles = males. (a) P10 USV Emission: The effect of U50,488h Treatment depended on Time (β = −5.89, SE = 0.52, t(812) = −11.42, ****p < 0.0001). U50,488h-treated pups vocalized more than saline pups during the 2–9min intervals (all *padj ≤ 0.020), and less than saline pups during the 14–15min intervals (both *padj ≤ 0.028). (b) P10 Total USV Emission: There was no effect of Sex (β = −167.71, SE = 205.43, t(54) = −1.34, p = 0.19) or a U50,488h Treatment × Sex interaction (β =−34.18, SE = 288.38, t(54) = −0.12, p = 0.91). The simplified model revealed that U50,488h Treatment was associated with increased USV emission compared to saline controls (β = 406.2, SE= 148.9, t(56) = 2.73, **p = 0.0085). (c) P10 Syllable Profile: There was no effect of Sex on the percentage of syllables emitted (all p ≥ 0.42) or a U50,488h Treatment × Sex interaction (all p ≥ 0.42). The simplified model revealed that U50,488h Treatment was associated with a decrease in the percentage of Downward syllables (β = −0.021, SE = 0.0070, t(56) = −3.04, **p = 0.0036) and Short syllables (β = −0.015, SE = 0.0072, t(56) = −2.11, *p = 0.039) compared to saline controls. (d) P14 USV Emission: The effect of U50,488h on USV emission depended on Time (β = −0.56, SE = 0.025, t(966) = −2.25, *p = 0.025), where U50,488h-treated pups vocalized more than saline pups during the 5–9min intervals (all **padj ≤ 0.0087). (e) P14 Total USV Emission: There was a U50,488h Treatment × Sex interaction (β = 202.10, SE = 97.93, t(65) = 2.06, *p = 0.043), where U50,488h Treatment was associated with increased USV emission in females compared to saline-treated females (β = 164.85, SE = 67.34, t(65) = 2.45, *p = 0.017). (f) P14 Syllable Profile: There was no effect of Sex on the percentage of syllables emitted (all p ≥ 0.12) or a U50,488h Treatment × Sex interaction (all p ≥ 0.15). The simplified model revealed that U50,488h Treatment was associated with a decrease in Complex 2 (β = −0.021, SE = 0.0082, t(67) = −2.57, *p = 0.012), Complex 3 (β = −0.078, SE = 0.018, t(67) = −4.43, ****p < 0.0001), and Complex 4 (β = −0.0080, SE = 0.0028, t(67) = −2.91, **p = 0.0049) syllables. SAL, n = 29–32 (18–19F, 11–13M); U50,488h, n = 31–37 (17–21F, 12–16M)

CONFLICTS OF INTERESTS

The authors have nothing to disclose.
==== Refs
REFERENCES

Abdel-Latif M. E. , Pinner J. , Clews S. , Cooke F. , Lui K. , & Oei J. (2006). Effects of Breast Milk on the Severity and Outcome of Neonatal Abstinence Syndrome Among Infants of Drug-Dependent Mothers. Pediatrics, 117 (6 ), e1163–e1169. 10.1542/peds.2005-1561 16740817
Alipio J. B. , Brockett A. T. , Fox M. E. , Tennyson S. S. , deBettencourt C. A. , El-Metwally D. , Francis N. A. , Kanold P. O. , Lobo M. K. , Roesch M. R. , & Keller A. (2021). Enduring consequences of perinatal fentanyl exposure in mice. Addiction Biology, 26 (2 ), e12895. 10.1111/adb.12895 32187805
Barr G. A. , McPhie-Lalmansingh A. , Perez J. , & Riley M. (2011). Changing Mechanisms of Opiate Tolerance and Withdrawal during Early Development: Animal Models of the Human Experience. ILAR Journal, 52 (3 ), 329–341. 10.1093/ilar.52.3.329 23382147
Barr G. A. , & Wang S. (1992). Tolerance and withdrawal to chronic morphine treatment in the week-old rat pup. European Journal of Pharmacology, 215 (1 ), 35–42. 10.1016/0014-2999(92)90605-4 1516648
Bartok R. E. , & Craft R. M. (1997). Sex Differences in Opioid Antinociception. Journal of Pharmacology and Experimental Therapeutics, 282 (2 ), 769–778.9262340
Basinger H. , & Hogg J. P. (2023). Neuroanatomy, Brainstem. In StatPearls. StatPearls Publishing. http://www.ncbi.nlm.nih.gov/books/NBK544297/
Becker J. B. , & Chartoff E. (2019). Sex differences in neural mechanisms mediating reward and addiction. Neuropsychopharmacology, 44 (1 ), 166–183. 10.1038/s41386-018-0125-6 29946108
Berger A. L. , Williams A. M. , McGinnis M. M. , & Walker B. M. (2013). Affective Cue-Induced Escalation of Alcohol Self-Administration and Increased 22-kHz Ultrasonic Vocalizations during Alcohol Withdrawal: Role of Kappa-Opioid Receptors. Neuropsychopharmacology, 38 (4 ), 647–654. 10.1038/npp.2012.229 23212453
Blinick G. , Tavolga W. N. , & Antopol W. (1971). Variations in birth cries of newborn infants from narcotic-addicted and normal mothers. American Journal of Obstetrics and Gynecology, 110 (7 ), 948–958. 10.1016/0002-9378(71)90549-7 5567416
Borrelli K. N. , Yao E. J. , Yen W. W. , Phadke R. A. , Ruan Q. T. , Chen M. M. , Kelliher J. C. , Langan C. R. , Scotellaro J. L. , Babbs R. K. , Beierle J. C. , Logan R. W. , Johnson W. E. , Wachman E. M. , Cruz-Martín A. , & Bryant C. D. (2021). Sex Differences in Behavioral and Brainstem Transcriptomic Neuroadaptations following Neonatal Opioid Exposure in Outbred Mice. eNeuro, 8 (5 ). 10.1523/ENEURO.0143-21.2021
Branchi I. , Santucci D. , & Alleva E. (2001). Ultrasonic vocalisation emitted by infant rodents: A tool for assessment of neurobehavioural development. Behavioural Brain Research, 125 (1 ), 49–56. 10.1016/S0166-4328(01)00277-7 11682093
Bruchas M. R. , Land B. B. , & Chavkin C. (2010). The dynorphin/kappa opioid system as a modulator of stress-induced and pro-addictive behaviors. Brain Research, 1314 , 44–55. 10.1016/j.brainres.2009.08.062 19716811
Bryant C. D. , Smith D. J. , Kantak K. M. , Nowak T. S. , Williams R. W. , Damaj M. I. , Redei E. E. , Chen H. , & Mulligan M. K. (2020). Facilitating Complex Trait Analysis via Reduced Complexity Crosses. Trends in Genetics, 36 (8 ), 549–562. 10.1016/j.tig.2020.05.003 32482413
Bush N. E. , & Ramirez J.-M. (2024). Latent neural population dynamics underlying breathing, opioid-induced respiratory depression and gasping. Nature Neuroscience, 27 (2 ), 259–271. 10.1038/s41593-023-01520-3 38182835
Byrnes E. M. , & Vassoler F. M. (2018). Modeling prenatal opioid exposure in animals: Current findings and future directions. Frontiers in Neuroendocrinology, 51 , 1–13. 10.1016/j.yfrne.2017.09.001 28965857
Carden S. E. , Barr G. A. , & Hofer M. A. (1991). Differential effects of specific opioid receptor agonists on rat pup isolation calls. Developmental Brain Research, 62 (1 ), 17–22. 10.1016/0165-3806(91)90185-L 1662121
Carden S. E. , Davachi L. , & Hofer M. A. (1994). U50, 488 increases ultrasonic vocalizations in 3-, 10-, and 18-day-old rat pups in isolation and the home cage. Developmental Psychobiology, 27 (1 ), 65–83. 10.1002/dev.420270107 8112489
Caruso A. , Marconi M. A. , Scattoni M. L. , & Ricceri L. (2022). Ultrasonic vocalizations in laboratory mice: Strain, age, and sex differences. Genes, Brain and Behavior, 21 (5 ), e12815. 10.1111/gbb.12815 35689354
Cérou F. , Guyader A. , & Rousset M. (2019). Adaptive multilevel splitting: Historical perspective and recent results. Chaos: An Interdisciplinary Journal of Nonlinear Science, 29 (4 ), 043108. 10.1063/1.5082247
Chartoff E. H. , & Mavrikaki M. (2015). Sex Differences in Kappa Opioid Receptor Function and Their Potential Impact on Addiction. Frontiers in Neuroscience, 9 , 466. 10.3389/fnins.2015.00466 26733781
Coffey K. R. , Marx R. E. , & Neumaier J. F. (2019). DeepSqueak: A deep learning-based system for detection and analysis of ultrasonic vocalizations. Neuropsychopharmacology, 44 (5 ), Article 5. 10.1038/s41386-018-0303-6
Conway S. M. , Puttick D. , Russell S. , Potter D. , Roitman M. F. , & Chartoff E. H. (2019). Females are less sensitive than males to the motivational- and dopamine-suppressing effects of kappa opioid receptor activation. Neuropharmacology, 146 , 231–241. 10.1016/j.neuropharm.2018.12.002 30528327
Corder G. , Castro D. C. , Bruchas M. R. , & Scherrer G. (2018). Endogenous and Exogenous Opioids in Pain. Annual Review of Neuroscience, 41 (1 ), 453–473. 10.1146/annurev-neuro-080317-061522
Craig A. , Ling Luo N. , Beardsley D. J. , Wingate-Pearse N. , Walker D. W. , Hohimer A. R. , & Back S. A. (2003). Quantitative analysis of perinatal rodent oligodendrocyte lineage progression and its correlation with human. Experimental Neurology, 181 (2 ), 231–240. 10.1016/S0014-4886(03)00032-3 12781996
D’Amato F. R. , Scalera E. , Sarli C. , & Moles A. (2005). Pups Call, Mothers Rush: Does Maternal Responsiveness Affect the Amount of Ultrasonic Vocalizations in Mouse Pups? Behavior Genetics, 35 (1 ), 103–112. 10.1007/s10519-004-0860-9 15674537
Dirks A. , Fish E. W. , Kikusui T. , van der Gugten J. , Groenink L. , Olivier B. , & Miczek K. A. (2002). Effects of corticotropin-releasing hormone on distress vocalizations and locomotion in maternally separated mouse pups. Pharmacology Biochemistry and Behavior, 72 (4 ), 993–999. 10.1016/S0091-3057(02)00809-2 12062591
Downs A. M. , & McElligott Z. A. (2022). Noradrenergic circuits and signaling in substance use disorders. Neuropharmacology, 208 . 10.1016/j.neuropharm.2022.108997
Driel K. S. van , & Talling J. C. (2005). Familiarity increases consistency in animal tests. Behavioural Brain Research, 159 (2 ), 243–245. 10.1016/j.bbr.2004.11.005 15817187
Dunn A. D. , Robinson S. A. , Nwokafor C. , Estill M. , Ferrante J. , Shen L. , Lemchi C. O. , Creus-Muncunill J. , Ramirez A. , Mengaziol J. , Brynildsen J. K. , Leggas M. , Horn J. , Ehrlich M. E. , & Blendy J. A. (2023). Molecular and long-term behavioral consequences of neonatal opioid exposure and withdrawal in mice. Frontiers in Behavioral Neuroscience, 17 , 1202099. 10.3389/fnbeh.2023.1202099 37424750
Ehret G. (2005). Infant Rodent Ultrasounds ? A Gate to the Understanding of Sound Communication. Behavior Genetics, 35 (1 ), 19–29. 10.1007/s10519-004-0853-8 15674530
Elwood R. W. , & Keeling F. (1982). Temporal organization of ultrasonic vocalizations in infant mice. Developmental Psychobiology, 15 (3 ), 221–227. 10.1002/dev.420150306 7095288
Ferrante J. R. , & Blendy J. A. (2024). Advances in animal models of prenatal opioid exposure. Trends in Neurosciences, 0 (0 ). 10.1016/j.tins.2024.03.005
Ferrante J. R. , Dunn A. D. , Mengaziol J. , Noreck J. , Lemchi C. O. , Creus-Muncunill J. , Ehrlich M. E. , & Blendy J. A. (2022). Neonatal Opioid Withdrawal Syndrome in Mice: Acute Molecular and Long-Term Behavioral Consequences. The FASEB Journal, 36 . 10.1096/fasebj.2022.36.S1.R4002
Finnegan L. P. , Connaughton J. F. , Kron R. E. , & Emich J. P. (1975). Neonatal abstinence syndrome: Assessment and management. Addictive Diseases, 2 (1–2 ), 141–158.1163358
Fox H. C. , & Sinha R. (2009). Sex Differences in Drug-Related Stress-System Changes: Implications for Treatment in Substance-Abusing Women. Harvard Review of Psychiatry, 17 (2 ), 103. 10.1080/10673220902899680 19373619
Gear R. W. , Miaskowski C. , Gordon N. C. , Paul S. M. , Heller P. H. , & Levine J. D. (1996). Kappa–opioids produce significantly greater analgesia in women than in men. Nature Medicine, 2 (11 ), 1248–1250. 10.1038/nm1196-1248
Georges F. , Normand E. , Bloch B. , & Le Moine C. (1998). Opioid receptor gene expression in the rat brain during ontogeny, with special reference to the mesostriatal system: An in situ hybridization study. Developmental Brain Research, 109 (2 ), 187–199. 10.1016/S0165-3806(98)00082-0 9729374
Geva R. , Sopher K. , Kurtzman L. , Galili G. , Feldman R. , & Kuint J. (2013). Neonatal brainstem dysfunction risks infant social engagement. Social Cognitive and Affective Neuroscience, 8 (2 ), 158–164. 10.1093/scan/nsr082 22146141
Goetz T. G. , Becker J. B. , & Mazure C. M. (2021). Women, opioid use and addiction. FASEB Journal: Official Publication of the Federation of American Societies for Experimental Biology, 35 (2 ), e21303. 10.1096/fj.202002125R 33433026
Grecco G. G. , Mork B. E. , Huang J.-Y. , Metzger C. E. , Haggerty D. L. , Reeves K. C. , Gao Y. , Hoffman H. , Katner S. N. , Masters A. R. , Morris C. W. , Newell E. A. , Engleman E. A. , Baucum A. J. , Kim J. , Yamamoto B. K. , Allen M. R. , Wu Y.-C. , Lu H.-C. , … Atwood B. K. (2021). Prenatal methadone exposure disrupts behavioral development and alters motor neuron intrinsic properties and local circuitry. eLife, 10 , e66230. 10.7554/eLife.66230 33724184
Grimsley J. M. S. , Monaghan J. J. M. , & Wenstrup J. J. (2011). Development of Social Vocalizations in Mice. PLOS ONE, 6 (3 ), e17460. 10.1371/journal.pone.0017460 21408007
Haight S. C. (2018). Opioid Use Disorder Documented at Delivery Hospitalization—United States, 1999–2014. MMWR. Morbidity and Mortality Weekly Report, 67 . 10.15585/mmwr.mm6731a1
Hales C. M. , Martin C. B. , & Gu Q. (2020). Prevalence of Prescription Pain Medication Use Among Adults: United States, 2015–2018. NCHS Data Brief, 369 , 1–8.
Heckman J. , McGuinness B. , Celikel T. , & Englitz B. (2016). Determinants of the mouse ultrasonic vocal structure and repertoire. Neuroscience & Biobehavioral Reviews, 65 , 313–325. 10.1016/j.neubiorev.2016.03.029 27060755
Howard M. B. , Schiff D. M. , Penwill N. , Si W. , Rai A. , Wolfgang T. , Moses J. M. , & Wachman E. M. (2017). Impact of Parental Presence at Infants’ Bedside on Neonatal Abstinence Syndrome. Hospital Pediatrics, 7 (2 ), 63–69. 10.1542/hpeds.2016-0147 28137920
Huang P. , Gentile T. A. , Muschamp J. W. , & Liu-Chen L.-Y. (2020). The kappa opioid receptor agonist U50,488H did not affect brain-stimulation reward while it elicited conditioned place aversion in mice. BMC Research Notes, 13 , 384. 10.1186/s13104-020-05227-7 32799930
Hudak M. L. , Tan R. C. , THE COMMITTEE ON DRUGS, THE COMMITTEE ON FETUS AND NEWBORN, Frattarelli D. A. C. , Galinkin J. L. , Green T. P. , Neville K. A. , Paul I. M. , Van Den Anker J. N. , Papile L.-A. , Baley J. E. , Bhutani V. K. , Carlo W. A. , Cummings J. , Kumar P. , Polin R. A. , Wang K. S. , & Watterberg K. L. (2012). Neonatal Drug Withdrawal. Pediatrics, 129 (2 ), e540–e560. 10.1542/peds.2011-3212 22291123
Jansson L. M. , & Velez M. (2012). Neonatal abstinence syndrome. Current Opinion in Pediatrics, 24 (2 ), 252–258. 10.1097/MOP.0b013e32834fdc3a 22227786
Jansson L. M. , Velez M. , & Harrow C. (2009). The Opioid Exposed Newborn: Assessment and Pharmacologic Management. Journal of Opioid Management, 5 (1 ), 47–55.19344048
Jarman S. K. , Haney A. M. , & Valdez G. R. (2018). Kappa opioid regulation of depressive-like behavior during acute withdrawal and protracted abstinence from ethanol. PLOS ONE, 13 (9 ), e0205016. 10.1371/journal.pone.0205016 30265734
Jürgens U. (2002). Neural pathways underlying vocal control. Neuroscience & Biobehavioral Reviews, 26 (2 ), 235–258. 10.1016/S0149-7634(01)00068-9 11856561
Kelsey J. E. , Verhaak A. M. S. , & Schierberl K. C. (2015). The kappa-opioid receptor antagonist, nor-binaltorphimine (nor-BNI), decreases morphine withdrawal and the consequent conditioned place aversion in rats. Behavioural Brain Research, 283 , 16–21. 10.1016/j.bbr.2015.01.008 25591478
Kishioka S. , Kiguchi N. , Kobayashi Y. , Yamamoto C. , Saika F. , Wakida N. , Ko M.-C. , & Woods J. H. (2013). Pharmacokinetic evidence for the long-lasting effect of nor-binaltorphimine, a potent kappa opioid receptor antagonist, in mice. Neuroscience Letters, 552 , 98–102. 10.1016/j.neulet.2013.07.040 23933210
Ko J. Y. , D’Angelo D. V. , Haight S. C. , Morrow B. , Cox S. , Essen B. S. von , Strahan A. E. , Harrison L. , Tevendale H. D. , Warner L. , Kroelinger C. D. , & Barfield W. D. (2020). Vital Signs: Prescription Opioid Pain Reliever Use During Pregnancy — 34 U.S. Jurisdictions, 2019. Morbidity and Mortality Weekly Report, 69 (28 ), 897. 10.15585/mmwr.mm6928a1 32673301
Ko J. Y. , Wolicki S. , Barfield W. D. , Patrick S. W. , Broussard C. S. , Yonkers K. A. , Naimon R. , & Iskander J. (2017). CDC Grand Rounds: Public Health Strategies to Prevent Neonatal Abstinence Syndrome. Morbidity and Mortality Weekly Report, 66 (9 ), 242–245. 10.15585/mmwr.mm6609a2 28278146
Korotkevich G. , Sukhov V. , Budin N. , Shpak B. , Artyomov M. N. , & Sergushichev A. (2016). Fast gene set enrichment analysis. 10.1101/060012
Kraft W. K. , Adeniyi-Jones S. C. , Chervoneva I. , Greenspan J. S. , Abatemarco D. , Kaltenbach K. , & Ehrlich M. E. (2017). Buprenorphine for the Treatment of the Neonatal Abstinence Syndrome. The New England Journal of Medicine, 376 (24 ), 2341–2348. 10.1056/NEJMoa1614835 28468518
LaGasse L. L. , Neal A. R. , & Lester B. M. (2005). Assessment of infant cry: Acoustic cry analysis and parental perception. Mental Retardation and Developmental Disabilities Research Reviews, 11 (1 ), 83–93. 10.1002/mrdd.20050 15856439
Laman-Maharg A. , Williams A. V. , Zufelt M. D. , Minie V. A. , Hao R. , Carroll F. I. , & Trainor B. C. (2018). Sex Differences in the Effects of a Kappa Opioid Receptor Antagonist in the Forced Swim Test. Frontiers in Pharmacology, 9 . 10.3389/fphar.2018.00093
Land B. B. , Bruchas M. R. , Lemos J. C. , Xu M. , Melief E. J. , & Chavkin C. (2008). The Dysphoric Component of Stress Is Encoded by Activation of the Dynorphin κ-Opioid System. The Journal of Neuroscience, 28 (2 ), 407–414. 10.1523/JNEUROSCI.4458-07.2008 18184783
Leconte C. , Mongeau R. , & Noble F. (2022). Traumatic Stress-Induced Vulnerability to Addiction: Critical Role of the Dynorphin/Kappa Opioid Receptor System. Frontiers in Pharmacology, 13 . https://www.frontiersin.org/articles/10.3389/fphar.2022.856672
Liu A. J. W. , Jones M. P. , Murray H. , Cook C.-M. , & Nanan R. (2010). Perinatal risk factors for the neonatal abstinence syndrome in infants born to women on methadone maintenance therapy. The Australian & New Zealand Journal of Obstetrics & Gynaecology, 50 (3 ). 10.1111/j.1479-828X.2010.01168.x
Madurai N. K. , Jantzie L. L. , & Yen E. (2024). Sex differences in neonatal outcomes following prenatal opioid exposure. Frontiers in Pediatrics, 12 . 10.3389/fped.2024.1357970
Manigault A. W. , Sheinkopf S. J. , Silverman H. F. , & Lester B. M. (2022). Newborn Cry Acoustics in the Assessment of Neonatal Opioid Withdrawal Syndrome Using Machine Learning. JAMA Network Open, 5 (10 ), e2238783. 10.1001/jamanetworkopen.2022.38783 36301544
Mansour A. , Hoversten M. T. , Taylor L. P. , Watson S. J. , & Akil H. (1995). The cloned μ, δ and κ receptors and their endogenous ligands: Evidence for two opioid peptide recognition cores. Brain Research, 700 (1 ), 89–98. 10.1016/0006-8993(95)00928-J 8624732
Marchette R. C. N. , Gregory-Flores A. , Tunstall B. J. , Carlson E. R. , Jackson S. N. , Sulima A. , Rice K. C. , Koob G. F. , & Vendruscolo L. F. (2021). κ-Opioid receptor antagonism reverses heroin withdrawal-induced hyperalgesia in male and female rats. Neurobiology of Stress, 14 , 100325. 10.1016/j.ynstr.2021.100325 33997152
Margolis E. B. , Hjelmstad G. O. , Bonci A. , & Fields H. L. (2003). κ-Opioid Agonists Directly Inhibit Midbrain Dopaminergic Neurons. The Journal of Neuroscience, 23 (31 ), 9981. 10.1523/JNEUROSCI.23-31-09981.2003 14602811
Margolis E. B. , & Karkhanis A. N. (2019). Dopaminergic Cellular and Circuit Contributions to Kappa Opioid Receptor Mediated Aversion. Neurochemistry International, 129 , 104504. 10.1016/j.neuint.2019.104504 31301327
Margolis E. B. , Lock H. , Chefer V. I. , Shippenberg T. S. , Hjelmstad G. O. , & Fields H. L. (2006). κ opioids selectively control dopaminergic neurons projecting to the prefrontal cortex. Proceedings of the National Academy of Sciences, 103 (8 ), 2938–2942. 10.1073/pnas.0511159103
Merico D. , Isserlin R. , Stueker O. , Emili A. , & Bader G. D. (2010). Enrichment Map: A Network-Based Method for Gene-Set Enrichment Visualization and Interpretation. PLoS ONE, 5 (11 ), e13984. 10.1371/journal.pone.0013984 21085593
Milliren C. E. , Gupta M. , Graham D. A. , Melvin P. , Jorina M. , & Ozonoff A. (2018). Hospital Variation in Neonatal Abstinence Syndrome Incidence, Treatment Modalities, Resource Use, and Costs Across Pediatric Hospitals in the United States, 2013 to 2016. Hospital Pediatrics, 8 (1 ), 15–20. 10.1542/hpeds.2017-0077 29263122
Mogil J. S. , Wilson S. G. , Chesler E. J. , Rankin A. L. , Nemmani K. V. S. , Lariviere W. R. , Groce M. K. , Wallace M. R. , Kaplan L. , Staud R. , Ness T. J. , Glover T. L. , Stankova M. , Mayorov A. , Hruby V. J. , Grisel J. E. , & Fillingim R. B. (2003). The melanocortin-1 receptor gene mediates female-specific mechanisms of analgesia in mice and humans. Proceedings of the National Academy of Sciences, 100 (8 ), 4867–4872. 10.1073/pnas.0730053100
Munro T. A. , Berry L. M. , Van’t Veer A. , Béguin C. , Carroll F. I. , Zhao Z. , Carlezon W. A. , & Cohen B. M. (2012). Long-acting κ opioid antagonists nor-BNI, GNTI and JDTic: Pharmacokinetics in mice and lipophilicity. BMC Pharmacology, 12 (1 ), 5. 10.1186/1471-2210-12-5 22642416
Nazarian A. , Rodarte-Freeman A. L. , & McDougall S. A. (1999). Dopaminergic modulation of kappa opioid-mediated ultrasonic vocalization, antinociception, and locomotor activity in the preweanling rat. Behavioral Neuroscience, 113 (4 ), 816–825. 10.1037/0735-7044.113.4.816 10495089
O’Connor A. B. , O’Brien L. , & Alto W. A. (2013). Are there gender related differences in neonatal abstinence syndrome following exposure to buprenorphine during pregnancy? Journal of Perinatal Medicine, 41 (5 ), 621–623. 10.1515/jpm-2012-0288 23612625
Patrick S. W. , Barfield W. D. , Poindexter B. B. , COMMITTEE ON FETUS AND NEWBORN, C. O. S. U. A. P., Cummings J. , Hand I. , Adams-Chapman I. , Aucott S. W. , Puopolo K. M. , Goldsmith J. P. , Kaufman D. , Martin C. , Mowitz M. , Gonzalez L. , Camenga D. R. , Quigley J. , Ryan S. A. , & Walker-Harding L. (2020). Neonatal Opioid Withdrawal Syndrome. Pediatrics, 146 (5 ), e2020029074. 10.1542/peds.2020-029074
Portfors C. V. (2007). Types and Functions of Ultrasonic Vocalizations in Laboratory Rats and Mice. Journal of the American Association for Laboratory Animal Science, 46 (1 ), 28–34.17203913
Reyes B. a. s. , Chavkin C. , & van Bockstaele E. j. (2009). Subcellular targeting of kappa-opioid receptors in the rat nucleus locus coeruleus. Journal of Comparative Neurology, 512 (3 ), 419–431. 10.1002/cne.21880 19009591
Rice D. , & Barone S. (2000). Critical periods of vulnerability for the developing nervous system: Evidence from humans and animal models. Environmental Health Perspectives, 108 (Suppl 3 ), 511–533.10852851
Richardson K. A. , Yohay A.-L. J. , Gauda E. B. , & McLemore G. L. (2006). Neonatal Animal Models of Opiate Withdrawal. ILAR Journal, 47 (1 ), 39–48. 10.1093/ilar.47.1.39 16391430
Robinson M. D. , McCarthy D. J. , & Smyth G. K. (2010). edgeR: A Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics, 26 (1 ), 139. 10.1093/bioinformatics/btp616 19910308
Robinson S. A. , Jones A. D. , Brynildsen J. K. , Ehrlich M. E. , & Blendy J. A. (2020). Neurobehavioral Effects of Neonatal Opioid Exposure in Mice: Influence of the OPRM1 SNP. Addiction Biology, 25 (5 ), e12806. 10.1111/adb.12806 31267641
Rorick-Kehn L. M. , Witkin J. M. , Statnick M. A. , Eberle E. L. , McKinzie J. H. , Kahl S. D. , Forster B. M. , Wong C. J. , Li X. , Crile R. S. , Shaw D. B. , Sahr A. E. , Adams B. L. , Quimby S. J. , Diaz N. , Jimenez A. , Pedregal C. , Mitch C. H. , Knopp K. L. , … McKinzie D. L. (2014). LY2456302 is a novel, potent, orally-bioavailable small molecule kappa-selective antagonist with activity in animal models predictive of efficacy in mood and addictive disorders. Neuropharmacology, 77 , 131–144. 10.1016/j.neuropharm.2013.09.021 24071566
Russell S. E. , Rachlin A. B. , Smith K. L. , Muschamp J. , Berry L. , Zhao Z. , & Chartoff E. H. (2014). Sex Differences in Sensitivity to the Depressive-like Effects of the Kappa Opioid Receptor Agonist U-50488 in Rats. Biological Psychiatry, 76 (3 ), 213–222. 10.1016/j.biopsych.2013.07.042 24090794
Saunders N. , Liddelow S. , & Dziegielewska K. (2012). Barrier Mechanisms in the Developing Brain. Frontiers in Pharmacology, 3 . https://www.frontiersin.org/articles/10.3389/fphar.2012.00046
Scattoni M. L. , Gandhy S. U. , Ricceri L. , & Crawley J. N. (2008). Unusual Repertoire of Vocalizations in the BTBR T+tf/J Mouse Model of Autism. PLOS ONE, 3 (8 ), e3067. 10.1371/journal.pone.0003067 18728777
Schlosburg J. E. , Whitfield T. W. , Park P. E. , Crawford E. F. , George O. , Vendruscolo L. F. , & Koob G. F. (2013). Long-Term Antagonism of κ Opioid Receptors Prevents Escalation of and Increased Motivation for Heroin Intake. The Journal of Neuroscience, 33 (49 ), 19384–19392. 10.1523/JNEUROSCI.1979-13.2013 24305833
Seligman N. S. , Salva N. , Hayes E. J. , Dysart K. C. , Pequignot E. C. , & Baxter J. K. (2008). Predicting length of treatment for neonatal abstinence syndrome in methadone-exposed neonates. American Journal of Obstetrics and Gynecology, 199 (4 ), 396.e1–396.e7. 10.1016/j.ajog.2008.06.088
Semple B. D. , Blomgren K. , Gimlin K. , Ferriero D. M. , & Noble-Haeusslein L. J. (2013). Brain development in rodents and humans: Identifying benchmarks of maturation and vulnerability to injury across species. Progress in Neurobiology, 106–107 , 1–16. 10.1016/j.pneurobio.2013.04.001
Shippenberg T. S. , Bals-Kubik R. , & Herz A. (1993). Examination of the neurochemical substrates mediating the motivational effects of opioids: Role of the mesolimbic dopamine system and D-1 vs. D-2 dopamine receptors. Journal of Pharmacology and Experimental Therapeutics, 265 (1 ), 53–59.8386244
Shippenberg T. S. , Zapata A. , & Chefer V. I. (2007). Dynorphin and the pathophysiology of drug addiction. Pharmacology & Therapeutics, 116 (2 ), 306–321. 10.1016/j.pharmthera.2007.06.011 17868902
Smith C. J. , Lintz T. , Clark M. J. , Malacon K. E. , Abiad A. , Constantino N. J. , Kim V. J. , Jo Y. C. , Alonso-Caraballo Y. , Bilbo S. D. , & Chartoff E. H. (2022). Prenatal opioid exposure inhibits microglial sculpting of the dopamine system selectively in adolescent male offspring. Neuropsychopharmacology, 47 (10 ), 1755–1763. 10.1038/s41386-022-01376-4 35835992
Sorge R. E. , Martin L. J. , Isbester K. A. , Sotocinal S. G. , Rosen S. , Tuttle A. H. , Wieskopf J. S. , Acland E. L. , Dokova A. , Kadoura B. , Leger P. , Mapplebeck J. C. S. , McPhail M. , Delaney A. , Wigerblad G. , Schumann A. P. , Quinn T. , Frasnelli J. , Svensson C. I. , … Mogil J. S. (2014). Olfactory exposure to males, including men, causes stress and related analgesia in rodents. Nature Methods, 11 (6 ), 629–632. 10.1038/nmeth.2935 24776635
Speciale S. G. , Manaye K. F. , Sadeq M. , & German D. C. (1993). Opioid receptors in midbrain dopaminergic regions of the rat II. Kappa and delta receptor autoradiography. Journal of Neural Transmission / General Section JNT, 91 (1 ), 53–66. 10.1007/BF01244918
Strahan A. E. , Guy G. P. Jr , Bohm M. , Frey M. , & Ko J. Y. (2020). Neonatal Abstinence Syndrome Incidence and Health Care Costs in the United States, 2016. JAMA Pediatrics, 174 (2 ), 200–202. 10.1001/jamapediatrics.2019.4791 31841581
Tan K. Z. , Cunningham A. M. , Joshi A. , Oei J. L. , & Ward M. C. (2018). Expression of kappa opioid receptors in developing rat brain – Implications for perinatal buprenorphine exposure. Reproductive Toxicology, 78 , 81–89. 10.1016/j.reprotox.2018.04.006 29635048
Timpson W. , Killoran C. , Maranda L. , Picarillo A. , & Bloch-Salisbury E. (2018). A Quality Improvement Initiative to Increase Scoring Consistency and Accuracy of the Finnegan Tool: Challenges in Obtaining Reliable Assessments of Drug Withdrawal in Neonatal Abstinence Syndrome. Advances in Neonatal Care, 18 (1 ), 70. 10.1097/ANC.0000000000000441 29045256
Tolia V. N. , Patrick S. W. , Bennett M. M. , Murthy K. , Sousa J. , Smith P. B. , Clark R. H. , & Spitzer A. R. (2015). Increasing Incidence of the Neonatal Abstinence Syndrome in U.S. Neonatal ICUs. New England Journal of Medicine, 372 (22 ), 2118–2126. 10.1056/NEJMsa1500439 25913111
Tschida K. , Michael V. , Takatoh J. , Han B.-X. , Zhao S. , Sakurai K. , Mooney R. , & Wang F. (2019). A Specialized Neural Circuit Gates Social Vocalizations in the Mouse. Neuron, 103 (3 ), 459–472.e4. 10.1016/j.neuron.2019.05.025 31204083
Valdez G. R. , & Harshberger E. (2012). Kappa opioid regulation of anxiety-like behavior during acute ethanol withdrawal. Pharmacology Biochemistry and Behavior, 102 (1 ), 44–47. 10.1016/j.pbb.2012.03.019 22487769
VanHouten J. P. , Rudd R. A. , Ballesteros M. F. , & Mack K. A. (2019). Drug Overdose Deaths Among Women Aged 30–64 Years—United States, 1999–2017. Morbidity and Mortality Weekly Report, 68 (1 ), 1–5. 10.15585/mmwr.mm6801a1 30629574
Vijay A. , Wang S. , Worhunsky P. , Zheng M.-Q. , Nabulsi N. , Ropchan J. , Krishnan-Sarin S. , Huang Y. , & Morris E. D. (2016). PET imaging reveals sex differences in kappa opioid receptor availability in humans, in vivo. American Journal of Nuclear Medicine and Molecular Imaging, 6 (4 ), 205–214.27648372
Vivian J. A. , & Miczek K. A. (1991). Ultrasounds during morphine withdrawal in rats. Psychopharmacology, 104 (2 ), 187–193. 10.1007/BF02244177 1876663
Vonvoigtlander P. F. , Lahti R. A. , & Ludens J. H. (1983). U-50,488: A selective and structurally novel non-Mu (kappa) opioid agonist. Journal of Pharmacology and Experimental Therapeutics, 224 (1 ), 7–12.6129321
Wachman E. M. , Hayes M. J. , Brown M. S. , Paul J. , Harvey-Wilkes K. , Terrin N. , Huggins G. S. , Aranda J. V. , & Davis J. M. (2013). Association of OPRM1 and COMT Single-Nucleotide Polymorphisms With Hospital Length of Stay and Treatment of Neonatal Abstinence Syndrome. JAMA, 309 (17 ), 1821–1827. 10.1001/jama.2013.3411 23632726
Wachman E. M. , Hayes M. J. , Sherva R. , Brown M. S. , Davis J. M. , Farrer L. A. , & Nielsen D. A. (2015). Variations in opioid receptor genes in neonatal abstinence syndrome. Drug and Alcohol Dependence, 155 , 253–259. 10.1016/j.drugalcdep.2015.07.001 26233486
Wachman E. M. , Schiff D. M. , & Silverstein M. (2018). Neonatal Abstinence Syndrome: Advances in Diagnosis and Treatment. JAMA, 319 (13 ), 1362–1374. 10.1001/jama.2018.2640 29614184
Wachman E. M. , & Werler M. M. (2019). Pharmacologic Treatment for Neonatal Abstinence Syndrome: Which Medication Is Best? JAMA Pediatrics, 173 (3 ), 221–223. 10.1001/jamapediatrics.2018.5029 30667482
Walker B. M. , Zorrilla E. P. , & Koob G. F. (2011). Systemic κ-Opioid Receptor Antagonism by Nor-binaltorphimine Reduces Dependence-Induced Excessive Alcohol Self-Administration in Rats. Addiction Biology, 16 (1 ), 116–119. 10.1111/j.1369-1600.2010.00226.x 20579007
Wang Z. , Hu J. , Johnson W. E. , & Campbell J. D. (2019). scruff: An R/Bioconductor package for preprocessing single-cell RNA-sequencing data. BMC Bioinformatics, 20 (1 ), 222. 10.1186/s12859-019-2797-2 31046658
Ward J. H. (1963). Hierarchical Grouping to Optimize an Objective Function. Journal of the American Statistical Association, 58 (301 ), 236–244. 10.1080/01621459.1963.10500845
Wei X. P. , Collie M. , Dempsey B. , Fortin G. , & Yackle K. (2022). A novel reticular node in the brainstem synchronizes neonatal mouse crying with breathing. Neuron, 110 (4 ), 644–657.e6. 10.1016/j.neuron.2021.12.014 34998469
Winkelman T. N. A. , Villapiano N. , Kozhimannil K. B. , Davis M. M. , & Patrick S. W. (2018). Incidence and Costs of Neonatal Abstinence Syndrome Among Infants With Medicaid: 2004–2014. Pediatrics, 141 (4 ), e20173520. 10.1542/peds.2017-3520
Zachry J. E. , Nolan S. O. , Brady L. J. , Kelly S. J. , Siciliano C. A. , & Calipari E. S. (2021). Sex differences in dopamine release regulation in the striatum. Neuropsychopharmacology, 46 (3 ), 491–499. 10.1038/s41386-020-00915-1 33318634
Zhu Y. , Hsu M.-S. , & Pintar J. E. (1998). Developmental Expression of the μ, κ, and δ Opioid Receptor mRNAs in Mouse. The Journal of Neuroscience, 18 (7 ), 2538–2549. 10.1523/JNEUROSCI.18-07-02538.1998 9502813
