Research Articles, Systems/Circuits

Perinatal Fentanyl Exposure Leads to Long-Lasting Impairments in Somatosensory Circuit Function and Behavior

Jason B. Alipio, Catherine Haga, Megan E. Fox, Keiko Arakawa, Rakshita Balaji, Nathan Cramer, Mary Kay Lobo and Asaf Keller
Journal of Neuroscience 14 April 2021, 41 (15) 3400-3417; DOI: https://doi.org/10.1523/JNEUROSCI.2470-20.2021

Abstract

One consequence of the opioid epidemic are lasting neurodevelopmental sequelae afflicting adolescents exposed to opioids in the womb. A translationally relevant and developmentally accurate preclinical model is needed to understand the behavioral, circuit, network, and molecular abnormalities resulting from this exposure. By employing a novel preclinical model of perinatal fentanyl exposure, our data reveal that fentanyl has several dose-dependent, developmental consequences to somatosensory function and behavior. Newborn male and female mice exhibit signs of withdrawal and sensory-related deficits that extend at least to adolescence. As fentanyl exposure does not affect dams' health or maternal behavior, these effects result from the direct actions of perinatal fentanyl on the pups' developing brain. At adolescence, exposed mice exhibit reduced adaptation to sensory stimuli, and a corresponding impairment in primary somatosensory (S1) function. In vitro electrophysiology demonstrates a long-lasting reduction in S1 synaptic excitation, evidenced by decreases in release probability, NMDA receptor-mediated postsynaptic currents, and frequency of miniature excitatory postsynaptic currents (mEPSCs), as well as increased frequency of miniature inhibitory postsynaptic currents (mIPSCs). In contrast, anterior cingulate cortical neurons exhibit an opposite phenotype, with increased synaptic excitation. Consistent with these changes, electrocorticograms (ECoGs) reveal suppressed ketamine-evoked γ oscillations. Morphologic analysis of S1 pyramidal neurons indicate reduced dendritic complexity, dendritic length, and soma size. Further, exposed mice exhibited abnormal cortical mRNA expression of key receptors involved in synaptic transmission and neuronal growth and development, changes that were consistent with the electrophysiological and morphologic changes. These findings demonstrate the lasting sequelae of perinatal fentanyl exposure on sensory processing and function.

SIGNIFICANCE STATEMENT This is the first study to show that exposure to fentanyl in the womb results in behavioral, circuitry, and synaptic effects that last at least to adolescence. We also show, for the first time, that this exposure has different, lasting effects on synapses in different cortical areas.

Introduction

Opioid use dependence and addiction have increased to epidemic proportions, leading to substantial financial and societal health burdens in the United States (Ryan, 2018). While opioid use dependence is highest in the United States, globally, the prevalence of opioid use disorders has increased by 47% from 1990 to 2016, with women representing one-third of the increase (GBD 2016 Alcohol and Drug Use Collaborators, 2018). In the United States, >20% of pregnant women enrolled in Medicaid are prescribed opioids. In recent years, there was a 14-fold increase in the proportion of pregnant women self-reporting opioid use. As a result of these trends, the incidence of infants born to opioid-using mothers has increased >400% between 2000 and 2012 (Haight et al., 2018; Winkelman et al., 2018; Honein et al., 2019). Not only does opioid use increase the risk of miscarriage, premature birth, and stillbirth (Whiteman et al., 2014), newborns exposed to opioids in utero may develop neonatal opioid withdrawal syndrome (NOWS), lower birth weight, smaller head circumference, and a higher risk of sudden infant death syndrome (Kandall et al., 1976; McPherson et al., 2015). Current pharmacologic treatment options for newborns with NOWS include opioid formulations of morphine, methadone, and buprenorphine (Sutter et al., 2014). While opioid therapy for newborns with NOWS improves acute withdrawal symptoms, such treatment further increases their exposure to these substances. Of particular concern is the fact that neurodevelopmental deficits may be permanent, as some have been found to persist at least to early adolescence (Lee et al., 2020).

There is a higher risk for sensory-related deficits in adolescents who were exposed to opioids prenatally (Ornoy et al., 2001; Kivisto et al., 2015). Perinatal opioid exposure also increases the risk for other disorders that are characterized by sensory processing deficits, such as attention deficit and autism spectrum disorders (Ayres, 1964; Rogers and Ozonoff, 2005; Robertson and Baron-Cohen, 2017; Balasco et al., 2019; Kilroy et al., 2019), which are also increased by perinatal opioid exposure (Ornoy, 2003; Hunt et al., 2008; Rubenstein et al., 2019).

Nearly all studies of perinatal exposure have focused on the naturally occurring opioid, morphine. However, the number of encounters with synthetic opioids, primarily fentanyl, has increased 300% since 2014 (O'Donnell et al., 2017; Jannetto et al., 2019). Fentanyl is 50–100 times more potent than morphine and is responsible for the majority of opioid-related overdose deaths (Volpe et al., 2011; Spencer et al., 2019). Understanding the effects that synthetic opioids have on the developmental trajectory of offspring is critically needed to develop effective treatment and prevention plans. To investigate the persistent neurobiological behavioral, circuit, network, and molecular abnormalities following perinatal fentanyl exposure, we developed a mouse model of such exposure, in which pregnant mouse dams are exposed to fentanyl (Alipio et al., 2020). We showed that newborn pups exhibit signs of spontaneous somatic withdrawal. During adolescence, fentanyl-exposed mice displayed abnormal affective-like behavior. By adulthood, fentanyl-exposed mice exhibited deficits in auditory processing.

In the present study, we further validate this model by determining the effects of different doses of fentanyl on maternal care and health. We then use this model to study the consequences of perinatal fentanyl exposure on sensory behaviors during adolescence, as well as the synaptic mechanisms that subserve them. We focus on the primary somatosensory cortex (S1) and anterior cingulate cortex (ACC), as they play key roles in sensory processing and integration (Mountcastle, 1998).

Materials and Methods

Animals

All procedures adhered to the Guide for the Care and Use of Laboratory Animals and approved by the Institutional Animal Care and Use Committee at the University of Maryland School of Medicine. Male and female C57BL/6J mice were used and bred in our temperature and humidity-controlled vivarium. Separate cohorts of animals were used for each behavioral test to prevent possible crossover effects between tests. When copulatory plugs were identified, we removed the sires, and added fentanyl or vehicle control to the water hydration pouches. Offspring were weaned at postnatal day (PD)21 and housed two to five per cage in single-sex groups. Food and water were available ad libitum, and lights were maintained on a 12/12 h light/dark cycle.

Statistical analyses

Statistical tests were conducted using Prism 8 (GraphPad) and sample size was determined using G*Power software suite (Heinrich-Heine, Universität Düsseldorf). Detailed statistical data are listed in Table 1. If there was no statistically significant sex difference or exposure/sex interaction, we grouped animals and analyzed according to exposure conditions, per NIH recommendations for studying sex differences. Parametric tests were used when appropriate assumptions were met, otherwise, nonparametric tests were used. Partial η2 were used to calculate effect size, and Cohen's d, Glass' δ, or Hedges' g were used for two group comparisons and selected based on whether parametric assumptions were met. All experimenters were blind to treatment conditions throughout data collection, scoring, and analysis.

View this table:
Table 1.

mRNA expression of receptors involved in synaptic transmission

Fentanyl citrate

We used 1, 10, or 100 µg/ml fentanyl citrate (calculated as free base) in 2% (w/v) saccharin, or 2% saccharin (vehicle control) in the water hydration pouches, which were replenished weekly until litters were weaned.

Maternal care behavior

We assessed maternal behavior as previously described (Alipio et al., 2020). We conducted daily observation sessions of pregnant dams (n = 9–10 dams per group) in the vivarium from PD1–PD7, between 9 and 10:30 A.M., and between 2 and 3 P.M. We scan sampled each dam in real time, once a minute for 30 min, yielding 30 scans per session; two sessions per day for 7 d yielded a total of 420 scans per dam. We used ethological parameters for assessing in and out of nest care-related behaviors, which comprised of: “licking/grooming,” “active nursing,” “passive nursing,” and “nest building.” Additionally, “self-grooming” and “eating/drinking” were taken as measures of dam self-maintenance whereas “pups out of nest” and “climbing/digging” were taken as measures of neglectful behavior.

Pup retrieval test

On PD7, we briefly removed the dam from the home cage and disturbed the nesting material, distributing it throughout the cage (n = 9–11 dams per group, with a minimum of 4 pups in the litter). We then placed two pups in the corners away from the nest end of the home cage. Then, we reintroduced the dam and measured the latency to sniff a pup, retrieve each of the pups, start nest building, and crouch over pups. We terminated the test if a dam did not complete the task within 15 min, resulting in a latency of 900 s for any behaviors not observed during that time frame.

Spontaneous somatic withdrawal behavior

We tested mice 24 h after cessation of drug effect on PD22 (n = 8–11 mice per group). We first habituated mice to the testing room for 1 h and scored behaviors in real time in 5-min time bins for a 15 min observation period. Spontaneous somatic withdrawal signs were scored using a modified rating scale (Schulteis et al., 1999). Counted signs included escape jumps and wet dog shakes. Escape jumps equaling 2–4 were assigned 1 point, 5–9 were assigned 2 points, and 10 or more were assigned 3 points. Wet dog shakes equaling 1–2 were assigned 2 points, and 3 or more were assigned 4 points. All other presence signs were assigned 2 points and consisted of 10 distinct behaviors, including persistent trembling, abnormal postures, abnormal gait, paw tremors, teeth chattering, penile erection/genital grooming, excessive eye blinking, ptosis (orbital tightening), swallowing movements, and diarrhea. Counted and presence scores were summed from all three of the 5-min time bins to compute the global withdrawal score.

Sensory threshold and adaptation

To assess tactile sensitivity, we habituated mice (n = 10–12 mice per group) to an elevated clear Plexiglas box with a mesh bottom for 10 min. We applied von Frey filaments of increasing forces to the plantar surface of the hind paw. The filament was applied to the same paw throughout the test. We used the up-down method to determine withdrawal threshold (Dixon, 1965; Chaplan et al., 1994; Deuis et al., 2017). To assess sensory adaptation, we applied a von Frey filament above threshold to the plantar surface of the hind paw, opposite to the hindpaw used during tactile sensitivity testing. The filament was applied once every 30 s until the animal stopped responding.

In vitro slice electrophysiology

We used a modified slice collection method (Ting et al., 2014). We anesthetized adolescent mice (PD35–PD45; see Results for sample size used in each experiment) with ketamine/xylazine, removed their brains, and prepared coronal slices (300 µm thick) containing the S1 or ACC. For recordings, we placed slices in a submersion chamber and continually perfused (2 ml/min) with artificial cerebrospinal fluid (ACSF) containing the following: 119 mm NaCl, 2.5 mm KCl, 1.2 mm NaH2PO4, 2.4 mm NaHCO3, 12.5 mm glucose, 2 mm MgSO4·7H2O, and 2 mm CaCl2·2H2O. Solutions were saturated with carbogen (95% O2 and 95% CO2) throughout use.

We obtained whole-cell patch-clamp recordings, in voltage-clamp mode (−70 mV), through pipettes containing the following: 130 mm cesium methanesulfonate, 10 mm HEPES, 1 mm magnesium chloride, 2.5 mm ATP-Mg, 0.5 mm EGTA, 0.2 mm GTP-Tris, 5 mm QX-314, and 2% biocytin. In S1 and ACC, we recorded from layer 5 neurons. For miniature postsynaptic current recordings, tetrodotoxin (1 μm) was included in the ACSF. All electrically evoked current responses were recorded at 1.5 times threshold, the minimum stimulation required to produce a current response. We positioned the stimulating electrode below the recorded neuron, in layer 6. To isolate EPSCs, gabazine (1 μm) was included in the ACSF. To isolate IPSCs, 6-cyano-7-nitroquinoxaline-2,3-dione (CNQX; 20 μm) and dl-2-amino-5-phosphonopentanoic acid (APV; 50 μm) was included in the ACSF. The N-methyl-D-aspartate receptor (NMDAR)-mediated component of the EPSC was pharmacologically isolated by including CNQX (20 μm) in the ACSF and by voltage clamping at 40 mV. Impedance of patch electrodes was 4–6 MΩ. Series resistance <40 MΩ was monitored throughout the recording, and recordings were discarded if series resistance changed by >20%. All recordings were obtained at room temperature.

In vivo electrocorticogram (ECoG) recordings

We anesthetized mice (PD30; see Results for sample size used in each experiment) with ketamine/xylazine and subcutaneously implanted a radiotelemetry transmitter (ETA-F10, Data Sciences International) with its leads implanted over the dura above S1 (−1.5 mm from bregma) and the cerebellum (−6.4 mm from bregma). Mice were allowed 10 d to recover from surgery. Mice were acclimated to the behavior testing room for 1 h before ECoG recordings. ECoGs were recorded with Dataquest A.R.T. acquisition system (Data Sciences International) with cortical ECoG recordings referenced to the cerebellum. Baseline 10-min ECoG recordings were followed by an intraperitoneal injection of vehicle saline or 10 mg/kg ketamine, and 50 min of postinjections recordings. Each mouse was injected with either ketamine or saline on consecutive days, with the order counterbalanced between animals. We analyzed ECoG recordings with custom-written MATLAB scripts (version 2019a, MathWorks) and the mtspecgramc routine in the Chronux Toolbox (http://chronux.org; Bokil et al., 2010). Oscillation power in the γ-bandwidth (30–80 Hz) was computed in 10-s time bins from spectrograms for each animal and averaged in 10-min bins.

Morphology

We filled cells with biocytin (0.1%) during slice electrophysiology recordings and immediately fixed slices in 4% paraformaldehyde for 24–48 h following recording (see Results for sample size used in each experiment). We washed slices three times for 10 min each with 1× PBS. We incubated slices overnight in streptavidin conjugated to Cy3 (Jackson ImmunoResearch; #016-160-084). We washed slices three times for 10 min each with 1× PBS and mounted the next day. To assess dendritic arbors, Z-stacks were reconstructed using Neurolucida software (MBF Bioscience). Sholl analysis was performed by counting basal and apical dendrite intersections across 10-µm concentric circles.

RNA isolation and RT-qPCR

We collected S1 and ACC tissue punches from adolescent (PD40) mice which were stored at −80°C until processing (n = 12–14 mice per group). RNA was isolated and extracted from these samples using TRIzol (Invitrogen) and the MicroElute total RNA kit (Omega Bio-tek, Inc) with a DNase step. RNA concentrations were measured on a NanoDrop spectrophotometer, and 400-ng cDNA was then synthesized using an iScript cDNA synthesis kit (Bio-Rad Laboratories). mRNA expression changes were measured by RT-qPCR with Perfecta SYBR Green FastMix (Quanta) and quantified using a CFX384 system (Bio-Rad Laboratories) using the ΔΔCT method as described previously (Fox et al., 2020) using GAPDH as a housekeeping gene. Primer sequences can be found in Table 2.

View this table:
Table 2.

RT-qPCR Primer sequences

Results

Fentanyl exposure has no effect on dam health or maternal care behavior

As we previously described (Alipio et al., 2020), we model perinatal fentanyl exposure by administering fentanyl citrate (0, 1, 10, or 100 µg/ml) in the drinking water of pregnant mouse dams throughout their pregnancy and until their litters were weaned at PD21. Consistent with our previous study, in which we also administered 10 µg/ml fentanyl in the drinking water, we predicted that fentanyl administration would not influence the dam's general health (n = 9–10 dams per group; Fig. 1). We compared body weight between different concentrations of fentanyl and vehicle control dams throughout pregnancy until weaning (Fig. 1B). There was no interaction between daily weight and fentanyl exposure groups [two-way repeated measures (RM) ANOVA, F(93,1023) = 0.88, p = 0.76], nor a main effect of fentanyl exposure on dam weights (two-way RM ANOVA, F(3,33) = 2.05, p = 0.12). We also assessed whether fentanyl exposure influences dams' liquid consumption (Fig. 1C). We found no interaction between daily liquid consumption and fentanyl exposure groups (two-way RM ANOVA, F(93,1023) = 0.87, p = 0.79), nor did fentanyl exposure affect liquid consumption between groups across days (two-way RM ANOVA, F(3,33) = 0.82, p = 0.49). We observed similar results with dams' food consumption (Fig. 1D). There was an interaction, with a medium effect size, between daily food consumption and fentanyl exposure (two-way RM ANOVA, F(93,1023) = 1.49, p < 0.002, partial η2 = 0.11); however, a post hoc multiple comparison analysis revealed no differences between exposure groups across each day (Tukey's post hoc test, p > 0.05). These data suggest that there is no difference in daily food consumption or effect of drug exposure between groups. Together, these data suggest that fentanyl exposure during pregnancy, at any of the concentrations tested, does not influence dams' total body weight, or liquid and food consumption.

Figure 1.
Figure 1.

Perinatal fentanyl exposure does not influence dam health or maternal care. A, Timeline depicting fentanyl exposure on dams, maternal care behavior, and pup retrieval test. There was no difference in dam weight (B), liquid consumption (C), or food consumption (D) across all concentrations of fentanyl tested. E, Fentanyl exposure did not adversely affect passive maternal care behaviors during PD1–PD7. F–J, There were no differences between exposure groups in latencies to sniff, retrieve pups, start nest building, or crouch over pups. Data depict means for parametric or medians for non-parametric comparisons with 95% confidence intervals.

Changes in maternal care during the early postnatal period can influence the development of neural systems, and their associated neurobiological and behavioral outputs (Meaney, 2001; Weaver et al., 2004). To assess maternal care behavior, we compared the time dams spent performing maternal care behaviors during 30-min sessions, twice a day, from PD1 to PD7 (n = 10 dams per group; Fig. 1E). There was no interaction between fentanyl exposure and maternal care behaviors (two-way ANOVA, F(24,324) = 0.89, p = 0.60), nor a main effect of fentanyl exposure (two-way ANOVA, F(3,324) = 0.18, p = 0.90), which suggests that fentanyl exposure does not influence maternal care behavior.

To assess a more active form of maternal care, we performed a pup retrieval test on PD7 (n = 9–11 dams per group; Fig. 1F–J). We compared latencies to sniff a pup, retrieve each of two pups, start nest building, and crouch over pups between fentanyl-exposed and control dams. We found no differences between groups in the latency to sniff a pup (ANOVA, F(3,35) = 1.04, p = 0.38; Fig. 1H), retrieve the first pup (Kruskal–Wallis test, H = 3.26, p = 0.35; Fig. 1G), retrieve the second pup (ANOVA, F(3,37) = 2.36, p = 0.08; Fig. 1H), nest build (ANOVA, F(3,36) = 0.91, p = 0.44; Fig. 1I), or crouch over their pups (ANOVA, F(3,37) = 1.43, p = 0.24; Fig. 1J). Taken together, these data suggest that fentanyl exposure during pregnancy does not adversely affect maternal care behaviors at any of the concentrations tested. We recognize, however, that we did not assess maternal care during the full postnatal period (PD8–PD21) and cannot exclude possible effects of maternal care during that period.

Fentanyl exposure results in smaller litters and higher mortality

Women who use opioids during pregnancy have increased incidence of sudden infant death syndrome and miscarriages (Kandall et al., 1975; Kahila et al., 2010; Brogly et al., 2018). We predicted that this increased risk will be observed in litters perinatally exposed to opioids (n = 12 litters per group Fig. 2B,C). We compared the number of pups born from dams exposed to different concentrations of fentanyl or vehicle (Fig. 2B). There was a large effect of fentanyl exposure on the number of pups per litter (ANOVA, F(3,44) = 14.89, p < 10−4, partial η2 = 0.50), with a large effect decrease in the number of pups per litter at birth that were exposed to each concentration of fentanyl tested, when compared with vehicle controls (Tukey's post hoc test, 1 µg/ml: p < 10−4, Cohen's d = 2.51, 10 µg/ml: p = 0.001, Cohen's d = 2.39, 100 µg/ml: p < 10−4, Cohen's d = 2.30). There were no differences in the number of pups per litter between fentanyl exposure groups (p > 0.05). We assessed the litter mortality rate at weaning on PD21 (Fig. 2C). There was a large effect of fentanyl exposure on litter mortality rate (Kruskal–Wallis test, H = 13.92, p = 0.003, partial η2 = 0.21). Post hoc analyses indicate a higher mortality rate with large effect sizes for litters that were exposed to each concentration of fentanyl tested, when compared with vehicle controls (Dunn's post hoc test, 1 µg/ml: p = 0.01, Glass' δ = 3.79, 10 µg/ml: p = 0.01, Glass' δ = 2.92, 100 µg/ml: p = 0.01, Glass' δ = 3.11). There were no differences in litter mortality rate between fentanyl exposure groups (p > 0.05). These data indicate that dams exposed to fentanyl during pregnancy, regardless of concentration, have smaller litters and a higher litter mortality rate compared with controls.

Figure 2.
Figure 2.

Perinatal fentanyl exposure results in aberrant effects at birth, withdrawal behavior, and impaired sensory function. A, Timeline depicting withdrawal behavior 24 h after weaning and sensory behavior test during adolescence. Exposure resulted in smaller litter size (B) and a higher litter mortality rate (C). D, Mice exposed to fentanyl perinatally exhibited abnormal weight during early development. E, At weaning (PD21), exposed male mice weighed less than controls. F, At adolescence (PD35), exposed male mice weighed more than controls and males weighed more than females at each concentration of fentanyl tested. G, By adulthood (PD55), males weighed more than females, but there were no differences between fentanyl exposure groups. H, Perinatal fentanyl exposure induces spontaneous somatic withdrawal behavior 24 h after cessation, at all concentrations of fentanyl tested. I, Perinatal fentanyl exposure impaired sensory adaptation to continuous application of tactile stimuli (I) but did not influence paw withdrawal threshold (J). Data depict means with 95% confidence intervals. *p < 0.05 compared to vehicle control, †p < 0.05 compared to opposite sex condition.

Abnormal weight across early development

Children prenatally exposed to opioids have lower weights and smaller head circumference compared with age matched controls (Wilson et al., 1979; McPherson et al., 2015). We predicted that pups perinatally exposed to fentanyl would have abnormal body weight across development (PD21–PD55). To test this prediction, we measured body weight across development of pups perinatally exposed to different concentrations of fentanyl and compared them to controls (n = 8–11 mice per group; Fig. 2D–G). There was an interaction, with a large effect size, between perinatal fentanyl exposure and weight (two-way RM ANOVA, F(238,2380) = 21.91, p < 10−4, partial η2 = 0.65), with a large main effect of fentanyl exposure (two-way RM ANOVA, F(7,70) = 30.41, p < 10−4, partial η2 = 0.83). These data suggest perinatal fentanyl exposure influences body weight across development.

To assess body weight at specific developmental periods, we tested whether pups perinatally exposed to fentanyl had abnormal body weight at weaning (PD21; Fig. 2E). There was an interaction, with a large effect size, between perinatal fentanyl exposure and sex (two-way ANOVA, F(3,71) = 6.45, p = 0.001, partial η2 = 0.23). Exposed males, but not females, weighed less than vehicle controls, with a large effect size (Tukey's post hoc test, 1 µg/ml: p < 10−4, Cohen's d = 2.82, 10 µg/ml: p = 0.004, Cohen's d = 1.96, 100 µg/ml: p < 10−4, Cohen's d = 2.51). There were no differences in weight between fentanyl exposure groups (p > 0.05) and, as expected, vehicle males weighed more than vehicle females (Tukey's post hoc test, p = 0.01, Cohen's d = 1.59).

At adolescence (PD35; Fig. 2F) there was no interaction between perinatal fentanyl exposure and sex (two-way ANOVA, F(3,70) = 2.49, p = 0.06). There were main effects of drug exposure (two-way ANOVA, F(3,70) = 4.55, p = 0.005, partial η2 = 0.16), and sex (two-way ANOVA, F(1,70) = 161.2, p < 10−4, partial η2 = 0.64). Exposed males, but not females, weighed more than vehicle controls, with a large effect size (Tukey's post hoc test, 1 µg/ml: p = 0.03, Cohen's d = 2.96, 10 µg/ml: p = 0.003, Cohen's d = 1.65, 100 µg/ml: p = 0.001, Cohen's d = 1.54). There were no differences in weight between fentanyl exposure groups (p > 0.05). Males weighed more than females at all concentrations tested (Tukey's post hoc test, VEH: p = 0.01, Cohen's d = 2.82, 1 µg/ml: p < 10−4, Cohen's d = 5.39, 10 µg/ml: p < 10−4, Cohen's d = 2.99, 100 µg/ml: p < 10−4, Cohen's d = 2.69).

By adulthood (PD55; Fig. 2G) there was no interaction between perinatal fentanyl exposure and sex (two-way ANOVA, F(3,70) = 0.20, p = 0.89). There were large main effects of drug exposure (two-way ANOVA, F(3,70) = 7.64, p = 0.0002, partial η2 = 0.22), and sex (two-way ANOVA, F(1,70) = 351.9, p < 10−4, partial η2 = 0.77). In male and female groups, there were no differences in weight when comparing across all fentanyl concentrations tested (Tukey's post hoc test, p > 0.05). Males weighed more than females at all concentrations tested (Tukey's post hoc test, VEH: p < 10−4, Cohen's d = 5.94, 1 µg/ml: p < 10−4, Cohen's d = 6.00, 10 µg/ml: p < 10−4, Cohen's d = 3.78, 100 µg/ml: p < 10−4, Cohen's d = 3.83).

Taken together, these results indicate that perinatal fentanyl exposure only affects the body weight of males at weaning. At all other age groups examined, neither males nor females that were treated perinatally differed in weight from age or sex-matched controls.

Fentanyl withdrawal behavior

Newborns exposed to opioids in utero may develop NOWS. Twenty-four hours after weaning (PD21), at a time when fentanyl was expected to have cleared in treated mice (Hug and Murphy, 1981), we tested whether perinatally exposed mice exhibited spontaneous somatic withdrawal signs (n = 8–11 mice per group; Fig. 2H). There was no interaction between sex and exposure (two-way ANOVA, F(3,73) = 1.81, p = 0.15). There were large main effects of sex (two-way ANOVA, F(1,73) = 12.69, p = 0.0007, partial η2 = 0.14) and fentanyl exposure (two-way ANOVA, F(3,73) = 189.5, p < 10−4, partial η2 = 0.88). Exposed males had higher withdrawal scores compared with controls, with a large effect size, at all concentrations tested (Tukey's post hoc test, 1 µg/ml: p < 10−4, Cohen's d = 6.90, 10 µg/ml: p < 10−4, Cohen's d = 3.29, 100 µg/ml: p < 10−4, Cohen's d = 8.25). Males exposed to 100 µg/ml had higher withdrawal scores compared with males exposed to 1 µg/ml (Tukey's post hoc test, p < 10−4, Cohen's d = 2.08) or 10 µg/ml (Tukey's post hoc test, p < 10−4, Cohen's d = 3.46). Males exposed to 1 µg/ml had higher withdrawal scores compared with males exposed to 10 µg/ml (Tukey's post hoc test, p = 0.0009, Cohen's d = 1.75).

Exposed females also had higher withdrawal scores compared with controls, with a large effect size (Tukey's post hoc test, 1 µg/ml: p < 10−4, Cohen's d = 4.95, 10 µg/ml: p = 0.009, Cohen's d = 2.47, 100 µg/ml p < 10−4, Cohen's d = 7.00). Females exposed to 100 µg/ml had higher withdrawal scores compared with females exposed to 1 µg/ml (Tukey's post hoc test, p < 10−4, Cohen's d = 1.99) or 10 µg/ml (Tukey's post hoc test, p < 10−4, Cohen's d = 4.53). Males exposed to 10 µg/ml had higher withdrawal scores compared with females exposed to 10 µg/ml (Tukey's post hoc test, p = 0.01, Cohen's d = 1.42). The expression of withdrawal signs supports the conclusion that our model recapitulates opioid exposure during pregnancy in humans.

Impaired sensory adaptation in adolescent mice

Adolescent children prenatally exposed to opioids may exhibit lasting somatosensory deficits (Kivisto et al., 2015) that may be independent of the expression of NOWS (Bakhireva et al., 2019). Adaptation, the cessation of a paw withdrawal response, to repeated stimuli is a key measure of sensory processing. We predicted that adolescent mice perinatally exposed to different concentrations of fentanyl would have impaired adaptation to repeated tactile stimuli (n = 10–12 mice per group; Fig. 2I). We compared the number of paw withdrawal responses to repeated application of a punctate stimulus that evoked withdrawal responses at threshold forces (see Materials and Methods). Compared with control mice, fentanyl-exposed mice failed to adapt, that is, they continued to respond to the stimuli. There was an interaction, with a medium effect size, between sex and fentanyl exposure (two-way ANOVA, F(3,78) = 2.736, p = 0.04, partial η2 = 0.09). Fentanyl-exposed males had a higher number of paw withdrawal responses compared with controls, with a large effect size at all concentrations tested (Tukey's post hoc test, 1 µg/ml: p = 0.01, Cohen's d = 1.23, 10 µg/ml: p < 10−4, Cohen's d = 3.06, 100 µg/ml: p < 10−4, Cohen's d = 4.28). Likewise, exposed females had a higher number of paw withdrawal responses compared with controls, with a large effect size at all concentrations tested (Tukey's post hoc test, 1 µg/ml: p = 0.0002, Cohen's d = 2.30, 10 µg/ml: p = 0.02, Cohen's d = 1.57, 100 µg/ml: p = 0.0007, Cohen's d = 2.80).

The failure to adapt was not because of reduced sensitivity to stimuli, because the responses to punctate stimuli elicited similar responses in both groups. We used von Frey filaments to test paw withdrawal thresholds (n = 10–12 mice per group; Fig. 2J). There was no interaction between sex and fentanyl exposure (two-way ANOVA, F(3,75) = 0.38, p = 0.76). There was a medium main effect of sex (two-way ANOVA, F(1,75) = 11.39, p = 0.001, partial η2 = 0.13) and a large effect of fentanyl exposure (two-way ANOVA, F(1,75) = 4.15, p = 0.008, partial η2 = 0.14). However, there were no differences in paw withdrawal thresholds within sex or exposure groups compared with controls (Tukey's post hoc test, p > 0.05). Taken together, these findings suggest that perinatal fentanyl exposure impairs adaptation to repeated tactile stimulation but does not influence mice's ability to sense tactile stimuli.

Decreased excitatory synaptic transmission in S1

Behavioral deficits in tactile sensory perception and processing are associated with changes in synaptic transmission in S1 (Fox et al., 2000). Therefore, we predicted that adolescent mice perinatally exposed to different concentrations of fentanyl would have impaired excitatory synaptic transmission in S1 (Fig. 3). We recorded miniature EPSCs (mEPSCs) in S1 layer 5 neurons from adolescent mice, in the presence of gabazine, a GABAA receptor antagonist (N = 4–5 mice per group, n = 2–5 neurons per mouse). Figure 3B depicts mEPSCs recorded from adolescent male and female mice perinatally exposed to different concentrations of fentanyl or vehicle control. We assessed the inter-event intervals, the reciprocal of their frequency, from 3-min recordings (Fig. 3C). There was a rightward shift in the cumulative frequency plot of the inter-event interval in fentanyl-exposed mice compared with controls (Kruskal–Wallis test, H = 3244, p < 10−4, partial η2 = 0.50). The cumulative frequency was shifted to the right in fentanyl-exposed male mice compared with controls, with a large effect size at all concentrations tested (Dunn's post hoc test, 1 µg/ml: p < 10−4, Hedges' g = 4.72, 10 µg/ml: p < 10−4, Hedges' g = 2.23, 100 µg/ml: p < 10−4, Hedges' g = 5.16). Similarly, there was a rightward shift in the cumulative frequency of the inter-event intervals in fentanyl exposed females compared with controls, with a large effect size (Dunn's post hoc test, 1 µg/ml: p < 10−4, Hedges' g = 4.51, 10 µg/ml: p < 10−4, Hedges' g = 1.09, 100 µg/ml: p < 10−4, Hedges' g = 3.52).

Figure 3.
Figure 3.

Perinatal fentanyl exposure impairs synaptic transmission in somatosensory cortical neurons. A, Timeline depicting slice electrophysiology recordings in S1 layer 5 neurons of adolescent mice. B, Example traces of mEPSCs. C, Cumulative frequency plot of the interevent intervals of mEPSCs are shifted to the right in fentanyl-exposed mice compared with controls. D, Grouped data reflect decreased mEPSC frequency in fentanyl-exposed mice. There were no differences in the cumulative frequency of the mEPSC amplitude (E) nor in the grouped data for mEPSC amplitude (F). G, Example traces of evoked paired pulse and NMDAR-mediated response. Perinatal fentanyl exposure results in increased paired pulse ratio in fentanyl-exposed mice (H). There were no differences in AMPAR-mediated response amplitude (I). J, There was decreased NMDAR-mediated response amplitude, which is also reflected in the NMDA/AMPA ratio (K). L, Example traces of mIPSCs. M, Cumulative frequency plot of the interevent intervals of mIPSCs are shifted to the left in fentanyl-exposed mice compared with controls. N, Grouped data reflect increased mIPSC frequency in fentanyl-exposed mice. O, Cumulative frequency plot of the mIPSC amplitude was shifted to the right. P, There were no differences in grouped data of mIPSC amplitude. Data depict means for parametric or medians for non-parametric comparisons with 95% confidence intervals.

The frequency of mEPSCs in these samples were lower in fentanyl-exposed mice, compared with controls. We averaged data from all neurons recorded from each animal and considered each animal as a single sample (Fig. 3D). This analysis revealed that mEPSCs in fentanyl-exposed mice were 3-fold less frequent than those in controls. There was no interaction between sex and fentanyl exposure (two-way ANOVA, F(3,29) = 0.57, p = 0.63). There was a main effect of fentanyl exposure with a large effect size (two-way ANOVA, F(3,29) = 30.58, p < 10−4, partial η2 = 0.75) and no effect of sex (two-way ANOVA, F(1,29) = 0.54, p = 0.46). Fentanyl-exposed males had a lower frequency of mEPSCs compared with controls, with a large effect size at all concentrations tested (Tukey's post hoc test, 1 µg/ml: p < 10−4, Cohen's d = 3.66, 10 µg/ml: p = 0.0003, Cohen's d = 2.66, 100 µg/ml: p < 10−4, Cohen's d = 3.93). Similarly, fentanyl exposed females also had a lower frequency of mEPSCs compared with controls, with a large effect size at all concentrations tested (Tukey's post hoc test, 1 µg/ml: p = 0.0009, Cohen's d = 3.02, 10 µg/ml: p = 0.002, Cohen's d = 2.04, 100 µg/ml: p = 0.001, Cohen's d = 2.96). There were no differences in the frequency of mEPSCs between male or female exposure groups regardless of fentanyl concentration (Tukey's post hoc test, p > 0.05) These data indicate that perinatal fentanyl exposure results in lower frequency of mEPSCs in S1 neurons of adolescent mice.

There were shifts in the cumulative frequency of mEPSC amplitudes between fentanyl-exposed mice compared with controls, with a small effect size (Kruskal–Wallis test, H = 543.9, p < 10−4, partial η2 = 0.21; Fig. 3E). The cumulative frequency of the amplitude was shifted to the left in fentanyl-exposed male mice compared with controls at all concentrations tested (Dunn's post hoc test, 1 µg/ml: p < 10−4, Hedges' g = 0.35, 10 µg/ml: p < 10−4, Hedges' g = 0.29, 100 µg/ml: p < 10−4, Hedges' g = 0.57). In female mice, there were no differences between vehicle control and the 1 µg/ml fentanyl group (Dunn's post hoc test, p = 0.17), a rightward shift in the 10 µg/ml fentanyl group, with a small effect size (p < 10−4, Hedges' g = 0.24), and a leftward shift in the 100 µg/ml fentanyl group, with a small effect size (p = 0.0008, Hedges' g = 0.21).

When data from all cells from a single animal are averaged, and comparisons are made between treated and control animals, we find no differences in amplitudes of mEPSCs in S1 neurons between sex or exposure groups (Fig. 3F). There was no interaction between sex and fentanyl exposure (two-way ANOVA, F(3,29) = 0.78, p = 0.51). There was no main effect of fentanyl exposure (two-way ANOVA, F(1,29) = 2.70, p = 0.06) nor a main effect of sex (two-way ANOVA, F(1,29) = 3.21, p = 0.08).

These data suggest that perinatal fentanyl exposure does not impact postsynaptic current responses in adolescent mice.

Since there were no significant interactions between sex and fentanyl exposure, nor a main effect of sex in any of the mEPSC analyses, we grouped and analyzed male and female mice according to exposure conditions. We performed subsequent experiments in mice exposed to 10 µg/ml fentanyl, since there were no robust concentration-dependent difference in the dependent variables tested above. Additionally, this concentration of fentanyl elicits affective and sensory deficits, and is the optimal concentration that mice will readily self-administer (Wade et al., 2008; Alipio et al., 2020).

Evoked responses

To determine whether fentanyl exposure affects evoked responses, we recorded electrically evoked EPSCs from S1 layer 5 neurons in ACSF containing gabazine to block inhibitory activity. Stimulating electrodes were placed below the recording neuron in layer 6. Figure 3G depicts paired pulse responses recorded from mice perinatally exposed to fentanyl or vehicle control. The paired pulse ratio is calculated by dividing the amplitude of the second pulse by the amplitude of the first pulse. Note that the ratio is higher in the mouse perinatally exposed to fentanyl, compared with the control. Group analysis revealed that perinatal fentanyl exposure led to an increased paired pulse ratio of evoked EPSCs recorded from S1 neurons, which was of a large effect size (N = 13–18 mice per group, n = 2–5 neurons per mouse, unpaired t test, t(29) = 3.14, p = 0.003, Cohen's d = 1.14; Fig. 3H). Because changes in paired pulse ratios primarily reflect changes in presynaptic mechanisms (Zucker and Regehr, 2002), these findings suggest that perinatal fentanyl exposure results in reductions in presynaptic glutamate release in adolescent mice.

NMDA/AMPA receptor ratios

Figure 3G depicts examples of evoked responses in which the AMPAR and NMDAR components of the response were isolated, as described in Materials and Methods. Note the smaller NMDA amplitude in fentanyl-exposed mouse, compared with the control. Group analysis revealed a reduction in the NMDAR/AMPAR ratio, with a medium effect size (N = 13–17 mice per group, n = 2–5 neurons per mouse, unpaired t test, t(31) = 2.23, p = 0.03, Cohen's d = 0.77; Fig. 3K). This is reflected by a large effect decrease in the NMDAR-mediated EPSC amplitude (Mann–Whitney test, U = 19, p = 0.0002, Glass' δ = 1.04; Fig. 3J), and no difference in AMPAR-mediated EPSC amplitude (unpaired t test, t(27) = 0.79, p = 0.43; Fig. 3I).

Together, these data indicate that perinatal fentanyl exposure impairs spontaneous and evoked excitatory synaptic transmission in S1 neurons in adolescent mice. And that this impairment is mediated both by a reduction in presynaptic glutamate release and a reduction in postsynaptic, NMDAR-mediated responses.

Increased inhibitory synaptic transmission in S1

Cortical activity involves both excitatory and inhibitory synaptic transmission. To determine whether perinatal fentanyl exposure influences inhibitory synaptic transmission in S1, we compared miniature IPSCs (mIPSCs) in S1 layer 5 neurons from adolescent mice perinatally exposed to fentanyl or vehicle control (N = 7–10 mice per group, n = 2–5 neurons per mouse). Figure 3L depicts mIPSCs recorded from a mouse perinatally exposed to fentanyl or a vehicle control. The frequency of mIPSCs in these samples is higher in the fentanyl-exposed mouse. The cumulative probability of the interevent interval was shifted to the left in fentanyl-exposed mice compared with controls, with a small effect size (Kolmogorov–Smirnov test, D = 0.26, p < 10−4, Hedges' g = 0.20; Fig. 3M). As a group, fentanyl-exposed mice exhibited a large effect increase in mIPSC frequency (unpaired t test, t(16) = 2.48, p = 0.02, Cohen's d = 1.38; Fig. 3N). These data suggest that perinatal fentanyl exposure increases the frequency of inhibitory spontaneous vesicle release events in S1 neurons from adolescent mice.

There was a rightward shift in the cumulative frequency of mIPSC amplitudes in fentanyl-exposed mice, which is of a small effect size (Kolmogorov–Smirnov test, D = 0.21, p < 10−4, Hedges' g = 0.29; Fig. 3O). This difference is driven by a larger number of events at lower amplitudes. However, when data from all cells from a single animal are averaged, and comparisons are made between treated and control animals, we find no difference in mIPSC amplitude (unpaired t test, t(18) = 1.08, p = 0.29; Fig. 3P).

These data suggest that perinatal fentanyl exposure increases spontaneous inhibitory synaptic transmission in S1 neurons from adolescent mice, and that these changes may be expressed through a presynaptic mechanism.

Increased excitatory synaptic transmission in ACC

The ACC has been implicated in complex sensory processing behaviors (Kennerley et al., 2006) and behavioral adaptation to environmental stimuli (Brockett et al., 2020). Since perinatal fentanyl exposure impaired sensory adaptation in adolescent mice, we predicted that synaptic transmission in ACC from exposed mice would also be impaired (Fig. 4). Surprisingly, the changes in the adolescent ACC of perinatally treated mice were opposite to those we identified in S1.

Figure 4.
Figure 4.

Perinatal fentanyl exposure impairs synaptic transmission in anterior cingulate cortical neurons. A, Timeline depicting slice electrophysiology recordings in ACC layer 5 neurons of adolescent mice. B, Example traces of mEPSCs. C, Cumulative frequency plot of the interevent intervals of mEPSCs are shifted to the left in fentanyl-exposed mice compared with controls. D, Grouped data reflect increased mEPSC frequency in fentanyl-exposed mice. E, Cumulative frequency of the mEPSC amplitude was shifted to the right in fentanyl-exposed mice. F, There were no differences in grouped data of mEPSC amplitudes. G, Example traces of evoked paired pulse and NMDAR-mediated responses. H, Perinatal fentanyl exposure results in decreased paired pulse ratio in fentanyl-exposed mice. There were no differences in AMPAR-mediated response amplitudes (I), NMDAR-mediated response amplitudes (J), or the NMDA/AMPA ratio (K). L, Example traces of mIPSCs. M, Cumulative frequency plot of the interevent intervals of mIPSCs are shifted to the left in fentanyl-exposed mice compared with controls. N, There were no differences in grouped data of mIPSC frequency. O, Cumulative frequency plot of the mIPSC amplitude was shifted to the right. P, There were no differences in grouped data of mIPSC amplitude. Data depict means for parametric or medians for non-parametric comparisons with 95% confidence intervals.

We recorded mEPSCs in ACC layer 5 neurons from adolescent mice perinatally exposed to fentanyl or vehicle control (N = 8–10 mice per group, n = 2–5 neurons per mouse). Figure 4B includes sample recordings depicting increased frequency and amplitude of mEPSCs recorded in ACC layer 5 neurons from an exposed mouse or a vehicle control. The cumulative frequency plot of the interevent intervals was shifted to the left in fentanyl-exposed mice compared with controls, with a large effect size (Kolmogorov–Smirnov test, D = 0.58, p < 10−4, Hedges' g = 1.79; Fig. 4C). As a group, fentanyl-exposed mice exhibited a large effect increase of mEPSC frequency in ACC neurons (unpaired t test, t(18) = 3.43, p = 0.003, Cohen's d = 1.53; Fig. 4D).

The cumulative frequency of the amplitude was shifted to the right in fentanyl-exposed mice, with a medium effect size (Kolmogorov–Smirnov test, D = 0.38, p < 10−4, Hedges' g = 0.59; Fig. 4E). However, as a group, there was no difference in the average mEPSC amplitude (unpaired t test, t(17) = 0.75, p = 0.46; Fig. 4F). As discussed above, this discrepancy might reflect a specific effect on a subset of mEPSCs.

These data suggest that perinatal fentanyl exposure increases the frequency of spontaneous excitatory synaptic transmission in ACC neurons from adolescent mice.

Evoked responses

Similar to S1, we investigated whether perinatal fentanyl exposure impairs evoked excitatory synaptic transmission in ACC layer 5 neurons. Stimulating electrodes were placed below the recording neuron in layer 6. Figure 4G depicts the paired pulse responses recorded from the ACC of a mouse perinatally exposed to fentanyl or control. Note the decreased paired pulse ratio in the fentanyl-exposed mouse. Group analysis revealed that fentanyl-exposed mice exhibit a large effect decrease in the paired pulse ratio of AMPAR-mediated EPSCs recorded from ACC neurons (N = 9–10 mice per group, n = 2–5 neurons per mouse, unpaired t test, t(18) = 5.73, p < 10−4, Cohen's d = 2.56; Fig. 4H).

NMDA/AMPA receptor ratios

Figure 4G depicts example traces of evoked AMPAR and NMDAR-mediated currents. There was no difference in the AMPAR-mediated EPSC amplitude (N = 8–11 mice per group, n = 2–5 neurons per mouse, unpaired t test, t(19) = 0.32, p = 0.74; Fig. 4I), the NMDAR-mediated EPSC amplitude (Mann–Whitney test, U = 16, p = 0.05; Fig. 4J), or in the NMDAR/AMPAR ratio (Mann–Whitney test, U = 37, p = 0.14; Fig. 4K). These data suggest that perinatal fentanyl exposure increases evoked excitatory transmission to ACC neurons, likely through presynaptic mechanisms.

No changes in inhibitory synaptic transmission in ACC

Figure 4L depicts mIPSCs recorded from an adolescent mouse perinatally exposed to fentanyl and a control. There was a leftward shift in the cumulative frequency of the interevent interval of fentanyl-exposed mice compared with controls, with a large effect size (N = 14–18 mice per group, n = 2–5 neurons per mouse, Kolmogorov–Smirnov test, D = 0.43, p < 10−4, Hedges' g = 1.16; Fig. 4M). There was no difference in the average mIPSC frequency between groups (unpaired t test, t(30) = 0.68, p = 0.49; Fig. 4N).

Likewise, the cumulative frequency of the amplitude was shifted to the left in fentanyl-exposed mice compared with controls, with a small effect size (N = 7–12 mice per group, n = 2–5 neurons per mouse, Kolmogorov–Smirnov test, D = 0.08, p < 10−4, Hedges' g = 0.08; Fig. 4O). There was no difference in average mIPSC amplitude between groups (unpaired t test, t(17) = 0.21, p = 0.83; Fig. 4P). These data suggest that perinatal fentanyl exposure does not influence spontaneous inhibitory synaptic transmission to ACC neurons of adolescent mice.

Taken together, these findings suggest that, in adolescent mice, the largest effects of perinatal fentanyl exposure is decreased presynaptic release of glutamate in S1, and an increase of this release in ACC.

Reduced cortical oscillatory response

The changes in synaptic efficacies described above suggest that network activity in mice perinatally exposed to fentanyl may be affected. Using electrocorticography (ECoG), we assessed synchronous network activity from large neuronal populations. We and others have previously shown that reproducible ECoG oscillations, particularly in the γ (30 to 80 Hz) range, can be evoked by injecting a subanesthetic dose of ketamine (10 mg/kg), an NMDAR antagonist (Raver et al., 2013; Raver and Keller, 2014). We used this approach to compare γ power change in fentanyl exposed and control mice (n = 7–9 mice per group; Fig. 5). Figure 5B depicts ECoGs recorded from adolescent mice treated perinatally with either fentanyl or vehicle, along with the spectrograms (Fig. 5C,D) of these recordings. Note that the increase in γ-band activity in the treated mouse is smaller than that in the control. Group analysis revealed no interaction between sex and exposure across all time points (two-way RM ANOVA, F(4,56) = 1.41, p = 0.24; Fig. 5E), therefore, mice were grouped according to fentanyl exposure. There was a large main effect of fentanyl exposure (F(1,14) = 36.74, p < 10−4, partial η2 = 0.41). Exposed mice exhibited smaller γ-band activity at both the 10–20 min (Bonferroni's post hoc test, p = 0.005) and 20–30 min (p = 0.01) postketamine time points. These data suggest that perinatal fentanyl exposure reduces ketamine-evoked γ-band cortical activity in awake adolescent mice.

Figure 5.
Figure 5.

Perinatal fentanyl exposure reduces ketamine-evoked cortical oscillations. A, Timeline depicting transmitter implant and in vivo ECoG recordings. B, Example ECoG traces of cortical oscillations recorded from fentanyl or vehicle control adolescent mice, before and after intraperitoneal injection of subanesthetic dose of ketamine (10 µg/ml). Normalized power of cortical oscillations in control (C) and perinatal fentanyl exposed (D) adolescent mice. E, Grouped data of γ power across 10-min time bins after ketamine injection. Data depict means with 95% confidence intervals. *p < 0.05.

Reduced branching of pyramidal neurons

Perinatal morphine exposure can delay dendritic development in cortical pyramidal neurons (Ricalde and Hammer, 1990). To determine whether perinatal fentanyl exposure influences dendritic morphology we analyzed the morphologies of layer 5 neurons in S1 and ACC (Fig. 6). Note the reduced dendritic complexity in pyramidal neurons from fentanyl-exposed mice (Fig. 6B).

Figure 6.
Figure 6.

Perinatal fentanyl exposure reduces morphology of basal dendrites of pyramidal neurons in S1. A, Timeline depicting morphologic assays of pyramidal neurons. B, Example of biocytin-filled layer 5 pyramidal neurons in S1 and ACC. Sholl analysis reveals that perinatal fentanyl exposure results in reduced branching of S1 basal (C) but not apical dendrites (D). There were no differences in branching of ACC basal (E) or apical dendrites (F). Fentanyl-exposed mice exhibited, in S1, decreased soma diameter (G) and decreased total dendritic length (H). Fentanyl-exposed mice also had decreased mRNA expression of TrkB receptors (K), with no difference in expression of BDNF (L). In ACC, there were no differences in soma diameter (I), total dendritic length (J), mRNA expression of TrkB receptors (M), or BDNF (N). Data depict means for parametric or medians for non-parametric comparisons with 95% confidence intervals. *p < 0.05.

S1 basal dendrites from fentanyl-exposed mice had less complex branching, evidenced as a smaller number of Sholl intersections (n = 5–7 mice per group; Fig. 6C). There was a large effect interaction between fentanyl exposure and the number of intersections across dendritic branch distance from the soma (two-way RM ANOVA, F(23,207) = 3.21, p < 10−4, partial η2 = 0.26), with a smaller number of intersections between 70–120 µm distance from the soma (Bonferroni's post hoc test; p < 0.05). Sholl analyses of S1 apical dendrites also revealed a main effect of fentanyl exposure (two-way RM ANOVA, F(1,9) = 5.99, p = 0.03, partial η2 = 0.13; Fig. 6D), and smaller number of intersections (two-way RM ANOVA, F(3,28) = 10.89, p < 10−4, partial η2 = 0.54).

Fentanyl exposure was associated with a large effect reduction in somata diameter (unpaired t test, t(10) = 2.95, p = 0.01, Cohen's d = 1.75; Fig. 6G), as well as a reduction in total dendritic length (unpaired t test, t(10) = 3.40, p = 0.006, Cohen's d = 2.03; Fig. 6H) of S1 pyramidal neurons.

As in S1, ACC basal dendrites from fentanyl-exposed mice had smaller number of Sholl intersections, with a small effect size (n = 12–19 mice per group; two-way RM ANOVA, F(31,899) = 1.50, p = 0.03, partial η2 = 0.04; Fig. 6E). There was no difference in the number of intersections of ACC apical dendrites (Fig. 6F). There was no interaction between fentanyl exposure and the number of intersections across dendritic branch distance from the soma (two-way RM ANOVA, F(41,1189) = 0.64, p = 0.96), nor a main effect of drug exposure (two-way RM ANOVA, F(1,29) = 0.17, p = 0.67).

There were no differences in somata diameter (unpaired t test, t(28) = 0.41, p = 0.68; Fig. 6I), nor in total dendritic length (unpaired t test, t(28) = 0.94, p = 0.35; Fig. 6J) in ACC pyramidal neurons.

Together, these data suggest perinatal fentanyl exposure decreases the dendritic arbor, length, and soma diameter of S1 pyramidal neurons. The dendritic arbor of basal dendrites in ACC were similarly decreased.

Decreased mRNA expression of tropomyosin receptor kinase B (TrkB) in S1

Brain-derived neurotrophic factor (BDNF) and its receptor, TrkB, are involved in growth, development, and maturation of neurons in S1 (Lush et al., 2005). Therefore, the changes in dendritic morphology in pyramidal neurons suggest a corresponding change in the expression of BDNF and TrkB in S1 of adolescent mice (n = 13–14 mice per group). Indeed, there was a large effect decrease in the expression of TrkB in S1 (unpaired t test, t(25) = 2.48, p = 0.01, Cohen's d = 0.95; Fig. 6K). There was no difference in the expression of BDNF (Mann–Whitney test, U = 69, p = 0.30; Fig. 6L). There were no differences in expression of TrkB (unpaired t test, t(25) = 1.74, p = 0.09; Fig. 6M) or BDNF in the ACC (unpaired t test, t(24) = 0.48, p = 0.63; Fig. 6N). RT-qPCR primers for BDNF target transcript variants 1–12 and primers for TrkB target full-length and short-length isoforms.

These results indicate that decreased mRNA expression of TrkB in S1 corroborates with decreased dendritic arbor, length, and soma size of S1 layer 5 pyramidal neurons. As predicted, there were no differences in expression of BDNF or TrkB in the ACC since there were small effect differences in the morphologic analysis of the dendritic arbor.

Increased cannabinoid receptor 1 (CB1R) mRNA expression in S1, decreased in ACC

CB1R is a Gi/o coupled heteroreceptor and are presynaptically expressed on both inhibitory and excitatory neurons in S1 and ACC (Lomazzo et al., 2017; Yeh et al., 2017). Activation of CB1Rs suppresses presynaptic release. We tested whether perinatal fentanyl exposure affects mRNA expression of CB1Rs in S1 and ACC of adolescent mice. In S1, there was a large effect increased expression of CB1Rs (unpaired t test, t(25) = 2.41, p = 0.02, Cohen's d = 0.92; Table 1), whereas in ACC there was a large effect decreased expression (unpaired t test, t(25) = 2.75, p = 0.01, Cohen's d = 1.05). These findings demonstrate that changes in CB1R expression parallel the changes in presynaptic release occurring after perinatal fentanyl exposure, suggesting that these changes may be causally related.

Increased GABAB receptor mRNA expression in ACC

A third mechanism of presynaptic transmitter regulation in neocortex involves GABAB receptors. In ACC, where there was an increase in mEPSC frequency, there was a large effect decrease mRNA expression of GABAB1 (unpaired t test, t(24) = 2.54, p = 0.01, Cohen's d = 0.99; Table 1). There was no difference in expression of GABAB2 receptors (unpaired t test, t(25) = 0.28, p = 0.77). In S1, where there was a decrease in mEPSC frequency, there was no difference in expression of GABAB1 (unpaired t test, t(25) = 0.27, p = 0.78), nor in GABAB2 (Mann–Whitney test, U = 64, p = 0.46). Thus, in ACC, but not in S1, changes in GABAB receptor mRNA expression were correlated with the electrophysiological findings.

Metabotropic glutamate receptor (mGluR) mRNA expression

In addition to CB1R, mGluRs regulate synaptic transmission. We compared the mRNA expression of several mGluR subtypes in adolescent mice perinatally exposed to fentanyl, with controls. In S1, there was a large effect decrease in expression of mGluR8 (Mann–Whitney test, U = 42, p = 0.01, Glass' δ = 0.78; Table 1). There was no difference in expression of the other group-3 mGluRs: mGluR4 (unpaired t test, t(25) = 0.53, p = 0.59) or mGluR7 (unpaired t test, t(25) = 0.31, p = 0.75). There were no differences in group-1 mGluRs: mGluR1 (unpaired t test, t(24) = 0.46, p = 0.13) and mGluR5 (unpaired t test, t(25) = 0.14, p = 0.88), nor in group-2 mGluRs: mGluR2 (unpaired t test, t(25) = 0.13, p = 0.89) and mGluR3 (unpaired t test, t(25) = 0.77, p = 0.44). Group-3 mGluR8 receptors inhibit transmitter release and are expressed on both glutamatergic and GABAergic terminals (Shigemoto et al., 1997; Ferraguti et al., 2005). Given the 3-fold decrease of S1 excitatory transmission, the decreased expression of mGluR8 was unexpected. These results suggest that decreased S1 excitatory transmission does not involve changes to mGluR8 mRNA expression.

In ACC, where we found a large decrease in presynaptic glutamate release, there was a large effect decrease mRNA expression of mGluR1 (unpaired t test, t(25) = 5.82, p < 10−4, Cohen's d = 2.26), but no difference in expression of mGluR5 (unpaired t test, t(25) = 1.23, p = 0.22). Among the group-2 mGluRs, there was a large effect decrease expression of mGluR2 (unpaired t test, t(25) = 2.55, p = 0.01, Cohen's d = 0.98), and mGluR3 (unpaired t test, t(25) = 2.13, p = 0.04, Cohen's d = 0.81). There were no differences in group-3 mGluRs: mGluR4 (unpaired t test, t(25) = 1.45, p = 0.15), mGluR7 (unpaired t test, t(25) = 0.82, p = 0.41), or mGluR8 (unpaired t test, t(25) = 2.01, p = 0.06). The decreases in group-1 and group-2 mGluRs are consistent with a reduction in presynaptic inhibition of glutamate release in ACC, and with our finding of increased excitatory transmission in this cortical area.

AMPAR and NMDAR mRNA expression

In S1, despite an absence of change in the amplitude of AMPAR currents, there was a large effect decrease in expression of GluR1 (unpaired t test, t(25) = 2.09, p = 0.04, Cohen's d = 0.79; Table 1), and GluR3 (Mann–Whitney test, U = 48, p = 0.03, Glass' δ = 1.06). There was no difference in expression of GluR2 (Mann–Whitney test, U = 62, p = 0.16), or GluR4 (unpaired t test, t(25) = 1.80, p = 0.08). We found no change in the amplitude of AMAR currents in ACC. Consistent with this, there was no difference in expression of GluR1 (unpaired t test, t(25) = 1.05, p = 0.30), GluR2 (unpaired t test, t(25) = 0.05, p = 0.95), GluR3 (unpaired t test, t(25) = 1.96, p = 0.06), or GluR4 (unpaired t test, t(25) = 0.09, p = 0.92) in ACC of mice perinatally exposed to fentanyl compared with controls. That changes in the expression of GluR1 and GluR3 did not result in changes in AMPA responses in S1 might relate to the complex interactions and functions among these subunits (Diering and Huganir, 2018).

The reduction in NMDAR-mediated currents in S1 neurons suggests that expression of NMDARs was affected by perinatal fentanyl exposure. However, we found no difference in expression of GluN2A (unpaired t test, t(25) = 0.73, p = 0.46; Table 1), GluN2B (Mann–Whitney test, U = 84, p = 0.75), or GluN2D (unpaired t test, t(25) = 0.67, p = 0.50) in mice perinatally exposed to fentanyl compared with controls. We did find a large effect increase expression of GluN2C (Mann–Whitney test, U= 36, p = 0.01, Glass' δ = 2.04); however these subunits are primarily expressed in S1 astrocytes (Palygin et al., 2011; Ravikrishnan et al., 2018).

In ACC, where there was no change in NMDA currents, there was a large effect decreased expression of GluN2A (unpaired t test, t(25) = 2.52, p = 0.01, Cohen's d = 0.96) and GluN2D (unpaired t test, t(24) = 6.17, p < 10−4, Cohen's d = 2.42). There was no difference in expression of GluN2B (unpaired t test, t(22) = 1.72, p = 0.09) or GluN2C (unpaired t test, t(25) = 0.83, p = 0.41).

These results may reflect the limited sensitivity of our mRNA expression assay. It is important to note that other molecular mechanisms may be affected by perinatal fentanyl exposure including receptor trafficking, expression of the receptor on the membrane, receptor affinity, or posttranslational changes in these subunits. Our expression data are consistent with previous reports showing that prenatal morphine exposure alters kinetic properties of NMDARs (Yang et al., 2000, 2003).

Discussion

We report that perinatal fentanyl exposure results in neurobiological deficits that last at least until adolescence. This exposure suppresses adaptation to sensory stimuli, impairs synaptic transmission in S1 and ACC, suppresses cortical oscillations, results in abnormal dendritic morphology of cortical pyramidal neurons, and alters mRNA expression of gene targets related to synaptic transmission and dendritic morphology.

The model

We confirm and expand on our previous description of a model of perinatal fentanyl exposure, in which dams are administered fentanyl in their drinking water (Alipio et al., 2020). The treatment extended from conception through weaning, to mimic the embryonic development period in humans (Chen et al., 2017). We expand on our previous findings by testing the effects of a range of fentanyl concentrations. Here, we find that this exposure has no effect on dam's health or maternal care behaviors, suggesting that the lasting effects we observe in perinatal fentanyl exposed adolescents reflects the direct actions of fentanyl during early development. Consistent with findings in humans, perinatal fentanyl exposure in mice results in smaller litters, higher newborn mortality rates, and signs of withdrawal shortly after birth. Thus, this model of perinatal opioid exposure has face validity for NOWS and for the lasting effects of this exposure.

Cortical excitatory tone

In S1 layer 5 neurons of adolescent mice, perinatal fentanyl exposure resulted in a reduction in synaptic excitation, and an increase in synaptic inhibition. The frequency of mEPSCs was increased and that of mIPSCs were decreased, consistent with changes in presynaptic release (Zucker and Regehr, 2002). The increase in the paired-pulse ratio of evoked EPSCs is consistent with reduced vesicle release from excitatory synapses. The amplitude of evoked excitatory synapses was decreased, due, in part, to a reduction in their NMDAR-mediated component.

In ACC layer 5 neurons of adolescent mice, perinatal fentanyl exposure enhanced excitatory synaptic transmission. Spontaneous excitatory vesicle release was enhanced, likely through a presynaptic mechanism. This presynaptic effect was mirrored in electrically evoked excitatory synaptic transmission. These findings are demonstrated by increases in the frequency of mEPSCs and depression of evoked paired pulse EPSCs.

Our findings are consistent with previous reports of long-lasting changes in neocortical and hippocampal synaptic activity and plasticity after prenatal exposure to morphine or heroin (Slotkin et al., 2003; Velísek et al., 2003; Yanai et al., 2003; Niu et al., 2009; Tan et al., 2015).

The changes in synaptic transmission in S1 and ACC were large and affected both male and female mice. For example, the frequency of mEPSC were 3-fold smaller in S1 and 2-fold larger in ACC than in controls. The NMDA/AMPA ratio was 2-fold smaller in S1. These findings indicate that the changes resulting from perinatal fentanyl exposure are physiologically salient.

The contrasting changes in excitatory synaptic transmission in S1 and ACC may be because of differential expression of opioid receptors during cortical development. Opioid receptors are detectable in the central nervous system of mice as early as embryonic day (ED)11.5 and develop in the cortical subplate at ED15.5 (Zhu et al., 1998). This is around the same time of migration and differentiation of neurons in layer IV of S1 (Angevine and Sidman, 1961). ACC development does not begin until the last few days of gestation, with full lamination occurring around PD10 (Van Eden and Uylings, 1985). Activity at these opioid receptors during this critical period of development constrains growth and development of cortical neurons (Zagon and McLaughlin, 1987; Hammer et al., 1989; Ricalde and Hammer, 1990). Thus, the contrasting effects of perinatal fentanyl exposure on S1 and ACC might reflect the different maturational state of opioid systems in these areas at the time of fentanyl exposure.

Morphologic changes

The changes in synaptic functions were accompanied by abnormalities in basal dendrites of S1 layer 5 pyramidal neurons. Consistent with this finding, there was a decrease in mRNA expression of TrkB, which is involved in neuronal growth, development, and maturation (Klein et al., 1989). These findings are consistent with prior studies reporting that perinatal morphine exposure reduces dendritic growth and development of S1 pyramidal neurons by acting on opioid receptors during early development (Ricalde and Hammer, 1990; Maharajan et al., 2000; Mei et al., 2009). More recent studies reveal that chronic morphine administration reduces TrkB mRNA expression in D1-type medium spiny neurons within the nucleus accumbens (Koo et al., 2014).

Neural oscillations

Cortical oscillations, particularly in the γ range, depend critically on the balance of inhibitory and excitatory synaptic activity (Buzsáki and Wang, 2012). Thus, we predicted that cortical oscillations will become disorganized by perinatal fentanyl exposure. Indeed, treated mice exhibited suppressed oscillations. These oscillations are related to large scale brain network activity and to cognitive phenomena, such as sensory processing, working memory, and attention (Lutz et al., 2004; Steriade, 2006; Buzsáki and Tingley, 2018). Perinatal exposure resulted also in decreased expression of the CB1R mRNA, consistent with reports of such decreases after morphine exposure (Rubino et al., 1997). A balance of excitatory and inhibitory transmission is reflected in γ oscillations, that are regulated by CB1R activity (Holderith et al., 2011; Buzsáki and Wang, 2012), suggesting that these synaptic anomalies may be causally related to the oscillation abnormalities. The important role γ oscillations in normal sensory processing and other cognitive tasks suggest that the abnormalities we observed may lead to the lasting deficits in sensory adaptation in mice treated perinatally with fentanyl (Mably and Colgin, 2018; Adaikkan and Tsai, 2020)

Sensory adaptation

A consequence of the constellation of circuit, network, and morphologic effects, reported here, of perinatal fentanyl exposure is compromised sensory adaptation. Here, we show that this exposure results in compromised adaptation to repeated tactile stimuli. Longitudinal clinical studies of children and adolescents exposed perinatally to opioids report cognitive, behavioral, and sensory deficits, as well as increased vulnerability to future stressors (Logan et al., 2013; Kivisto et al., 2015; Maguire et al., 2016; Tobon et al., 2019). Impaired sensory processing is associated with attention deficit disorders, autism spectrum, schizophrenia, and synesthesia (Freedman et al., 1987; Lou et al., 1989; Leekam et al., 2007; van Leeuwen et al., 2020).

The synaptic and circuit sequelae of perinatal fentanyl exposure reported here suggest that these lasting behavioral deficits reflect specific and contrasting lasting changes in different cortical regions. Further studies of these changes may direct interventions to ameliorate or prevent these neuropsychiatric deficits.

Footnotes

  • This work was supported by the Opioid Use Disorders Initiative, MPowering The State, from the State of Maryland. Support was provided also by NIH/NIDA R01 DA038613 (to M.K.L.) and K99 DA050575 (to M.E.F.). We thank Lace M. Riggs for providing valuable comments on the manuscript.

  • The authors declare no competing financial interests.

  • Correspondence should be addressed to Asaf Keller at akeller{at}som.umaryland.edu

References

  1. ↵
    2. Adaikkan C,
    3. Tsai LH
    (2020) Gamma entrainment: impact on neurocircuits, glia, and therapeutic opportunities. Trends Neurosci 43:24–41. doi:10.1016/j.tins.2019.11.001 pmid:31836315
    OpenUrlCrossRefPubMed
  2. ↵
    2. Alipio JB,
    3. Brockett AT,
    4. Fox ME,
    5. Tennyson SS,
    6. deBettencourt CA,
    7. El-Metwally D,
    8. Francis NA,
    9. Kanold PO,
    10. Lobo MK,
    11. Roesch MR,
    12. Keller A
    (2020) Enduring consequences of perinatal fentanyl exposure in mice. Addict Biol. Advance online publication. Retrieved March 18, 2020. doi: 10.1111/adb.12895.
    OpenUrlCrossRef
  3. ↵
    2. Angevine JB,
    3. Sidman RL
    (1961) Autoradiographic study of cell migration during histogenesis of cerebral cortex in the mouse. Nature 192:766–768. doi:10.1038/192766b0 pmid:17533671
    OpenUrlCrossRefPubMed
  4. ↵
    2. Ayres AJ
    (1964) Tactile functions. Their relation to hyperactive and perceptual motor behavior. Am J Occup Ther 18:6–11. pmid:14116444
    OpenUrlPubMed
  5. ↵
    2. Bakhireva LN,
    3. Holbrook BD,
    4. Shrestha S,
    5. Leyva Y,
    6. Ashley M,
    7. Cano S,
    8. Lowe J,
    9. Stephen JM,
    10. Leeman L
    (2019) Association between prenatal opioid exposure, neonatal opioid withdrawal syndrome, and neurodevelopmental and behavioral outcomes at 5-8 months of age. Early Hum Dev 128:69–76. doi:10.1016/j.earlhumdev.2018.10.010 pmid:30554024
    OpenUrlCrossRefPubMed
  6. ↵
    2. Balasco L,
    3. Provenzano G,
    4. Bozzi Y
    (2019) Sensory abnormalities in autism spectrum disorders: a focus on the tactile domain, from genetic mouse models to the clinic. Front Psychiatry 10:1016. doi:10.3389/fpsyt.2019.01016 pmid:32047448
    OpenUrlCrossRefPubMed
  7. ↵
    2. Bokil H,
    3. Andrews P,
    4. Kulkarni JE,
    5. Mehta S,
    6. Mitra PP
    (2010) Chronux: a platform for analyzing neural signals. J Neurosci Methods 192:146–151. doi:10.1016/j.jneumeth.2010.06.020 pmid:20637804
    OpenUrlCrossRefPubMed
  8. ↵
    2. Brockett AT,
    3. Tennyson SS,
    4. deBettencourt CA,
    5. Gaye F,
    6. Roesch MR
    (2020) Anterior cingulate cortex is necessary for adaptation of action plans. Proc Natl Acad Sci USA 117:6196–6204. doi:10.1073/pnas.1919303117 pmid:32132213
    OpenUrlAbstract/FREE Full Text
  9. ↵
    2. Brogly SB,
    3. Saia KE,
    4. Werler MM,
    5. Regan E,
    6. Hernández-Díaz S
    (2018) Prenatal treatment and outcomes of women with opioid use disorder. Obstet Gynecol 132:916–922. doi:10.1097/AOG.0000000000002881 pmid:30204704
    OpenUrlCrossRefPubMed
  10. ↵
    2. Buzsáki G,
    3. Wang XJ
    (2012) Mechanisms of gamma oscillations. Annu Rev Neurosci 35:203–225. doi:10.1146/annurev-neuro-062111-150444 pmid:22443509
    OpenUrlCrossRefPubMed
  11. ↵
    2. Buzsáki G,
    3. Tingley D
    (2018) Space and time: the hippocampus as a sequence generator. Trends Cogn Sci 22:853–869. doi:10.1016/j.tics.2018.07.006 pmid:30266146
    OpenUrlCrossRefPubMed
  12. ↵
    2. Chaplan SR,
    3. Bach FW,
    4. Pogrel JW,
    5. Chung JM,
    6. Yaksh TL
    (1994) Quantitative assessment of tactile allodynia in the rat paw. J Neurosci Methods 53:55–63. doi:10.1016/0165-0270(94)90144-9 pmid:7990513
    OpenUrlCrossRefPubMed
  13. ↵
    2. Chen VS,
    3. Morrison JP,
    4. Southwell MF,
    5. Foley JF,
    6. Bolon B,
    7. Elmore SA
    (2017) Histology atlas of the developing prenatal and postnatal mouse central nervous system, with emphasis on prenatal days E7.5 to E18.5. Toxicol Pathol 45:705–744. doi:10.1177/0192623317728134 pmid:28891434
    OpenUrlCrossRefPubMed
  14. ↵
    2. Deuis JR,
    3. Dvorakova LS,
    4. Vetter I
    (2017) Methods used to evaluate pain behaviors in rodents. Front Mol Neurosci 10:284. doi:10.3389/fnmol.2017.00284 pmid:28932184
    OpenUrlCrossRefPubMed
  15. ↵
    2. Diering GH,
    3. Huganir RL
    (2018) The AMPA receptor code of synaptic plasticity. Neuron 100:314–329. doi:10.1016/j.neuron.2018.10.018 pmid:30359599
    OpenUrlCrossRefPubMed
  16. ↵
    2. Dixon WJ
    (1965) The up-and-down method for small samples. J Am Stat Assoc 60:967–978. doi:10.1080/01621459.1965.10480843
    OpenUrlCrossRef
  17. ↵
    2. Ferraguti F,
    3. Klausberger T,
    4. Cobden P,
    5. Baude A,
    6. Roberts JD,
    7. Szucs P,
    8. Kinoshita A,
    9. Shigemoto R,
    10. Somogyi P,
    11. Dalezios Y
    (2005) Metabotropic glutamate receptor 8-expressing nerve terminals target subsets of GABAergic neurons in the hippocampus. J Neurosci 25:10520–10536. doi:10.1523/JNEUROSCI.2547-05.2005 pmid:16280590
    OpenUrlAbstract/FREE Full Text
  18. ↵
    2. Fox K,
    3. Glazewski S,
    4. Schulze S
    (2000) Plasticity and stability of somatosensory maps in thalamus and cortex. Curr Opin Neurobiol 10:494–497. doi:10.1016/s0959-4388(00)00112-4 pmid:10981619
    OpenUrlCrossRefPubMed
  19. ↵
    2. Fox ME,
    3. Chandra R,
    4. Menken MS,
    5. Larkin EJ,
    6. Nam H,
    7. Engeln M,
    8. Francis TC,
    9. Lobo MK
    (2020) Dendritic remodeling of D1 neurons by RhoA/Rho-kinase mediates depression-like behavior. Mol Psychiatry 25:1022–1034. doi:10.1038/s41380-018-0211-5 pmid:30120419
    OpenUrlCrossRefPubMed
  20. ↵
    2. Freedman R,
    3. Adler LE,
    4. Gerhardt GA,
    5. Waldo M,
    6. Baker N,
    7. Rose GM,
    8. Drebing C,
    9. Nagamoto H,
    10. Bickford-Wimer P,
    11. Franks R
    (1987) Neurobiological studies of sensory gating in schizophrenia. Schizophr Bull 13:669–678. doi:10.1093/schbul/13.4.669 pmid:2894074
    OpenUrlCrossRefPubMed
  21. ↵
    GBD 2016 Alcohol and Drug Use Collaborators (2018) The global burden of disease attributable to alcohol and drug use in 195 countries and territories, 1990-2016: a systematic analysis for the Global Burden of Disease Study 2016. Lancet Psychiatry 5:987–1012.
    OpenUrl
  22. ↵
    2. Haight SC,
    3. Ko JY,
    4. Tong VT,
    5. Bohm MK,
    6. Callaghan WM
    (2018) Opioid use disorder documented at delivery hospitalization - United States, 1999-2014. MMWR Morb Mortal Wkly Rep 67:845–849. doi:10.15585/mmwr.mm6731a1
    OpenUrlCrossRefPubMed
  23. ↵
    2. Hammer RP,
    3. Ricalde AA,
    4. Seatriz JV
    (1989) Effects of opiates on brain development. Neurotoxicology 10:475–483. pmid:2696899
    OpenUrlPubMed
  24. ↵
    2. Holderith N,
    3. Németh B,
    4. Papp OI,
    5. Veres JM,
    6. Nagy GA,
    7. Hájos N
    (2011) Cannabinoids attenuate hippocampal γ oscillations by suppressing excitatory synaptic input onto CA3 pyramidal neurons and fast spiking basket cells. J Physiol 589:4921–4934. doi:10.1113/jphysiol.2011.216259 pmid:21859823
    OpenUrlCrossRefPubMed
  25. ↵
    2. Honein MA,
    3. Boyle C,
    4. Redfield RR
    (2019) Public health surveillance of prenatal opioid exposure in mothers and infants. Pediatrics 143:e20183801. doi:10.1542/peds.2018-3801
    OpenUrlFREE Full Text
  26. ↵
    2. Hug CC,
    3. Murphy MR
    (1981) Tissue redistribution of fentanyl and termination of its effects in rats. Anesthesiology 55:369–375. doi:10.1097/00000542-198110000-00006 pmid:7294371
    OpenUrlCrossRefPubMed
  27. ↵
    2. Hunt RW,
    3. Tzioumi D,
    4. Collins E,
    5. Jeffery HE
    (2008) Adverse neurodevelopmental outcome of infants exposed to opiate in-utero. Early Hum Dev 84:29–35. doi:10.1016/j.earlhumdev.2007.01.013 pmid:17728081
    OpenUrlCrossRefPubMed
  28. ↵
    2. Jannetto PJ,
    3. Helander A,
    4. Garg U,
    5. Janis GC,
    6. Goldberger B,
    7. Ketha H
    (2019) The fentanyl epidemic and evolution of fentanyl analogs in the United States and the European Union. Clin Chem 65:242–253. doi:10.1373/clinchem.2017.281626 pmid:30305277
    OpenUrlAbstract/FREE Full Text
  29. ↵
    2. Kahila H,
    3. Gissler M,
    4. Sarkola T,
    5. Autti-Rämö I,
    6. Halmesmäki E
    (2010) Maternal welfare, morbidity and mortality 6-15 years after a pregnancy complicated by alcohol and substance abuse: a register-based case-control follow-up study of 524 women. Drug Alcohol Depend 111:215–221. doi:10.1016/j.drugalcdep.2010.04.014
    OpenUrlCrossRefPubMed
  30. ↵
    2. Kandall SR,
    3. Albin S,
    4. Dreyer E,
    5. Comstock M,
    6. Lowinson J
    (1975) Differential effects of heroin and methadone on birth weights. Addict Dis 2:347–355. pmid:1163371
    OpenUrlPubMed
  31. ↵
    2. Kandall SR,
    3. Albin S,
    4. Lowinson J,
    5. Berle B,
    6. Eidelman AI,
    7. Gartner LM
    (1976) Differential effects of maternal heroin and methadone use on birthweight. Pediatrics 58:681–685. pmid:980601
    OpenUrlAbstract/FREE Full Text
  32. ↵
    2. Kennerley SW,
    3. Walton ME,
    4. Behrens TEJ,
    5. Buckley MJ,
    6. Rushworth MFS
    (2006) Optimal decision making and the anterior cingulate cortex. Nat Neurosci 9:940–947. doi:10.1038/nn1724 pmid:16783368
    OpenUrlCrossRefPubMed
  33. ↵
    2. Kilroy E,
    3. Aziz-Zadeh L,
    4. Cermak S
    (2019) Ayres theories of autism and sensory integration revisited: what contemporary neuroscience has to say. Brain Sci 9:68.
    OpenUrl
  34. ↵
    2. Kivisto K,
    3. Nevalainen P,
    4. Lauronen L,
    5. Tupola S,
    6. Pihko E,
    7. Kivitie-Kallio S
    (2015) Somatosensory and auditory processing in opioid-exposed newborns with neonatal abstinence syndrome: a magnetoencephalographic approach. J Matern Fetal Neonatal Med 28:2015–2019.
    OpenUrl
  35. ↵
    2. Klein R,
    3. Parada LF,
    4. Coulier F,
    5. Barbacid M
    (1989) trkB, a novel tyrosine protein kinase receptor expressed during mouse neural development. EMBO J 8:3701–3709. pmid:2555172
    OpenUrlPubMed
  36. ↵
    2. Koo JW,
    3. Lobo MK,
    4. Chaudhury D,
    5. Labonté B,
    6. Friedman A,
    7. Heller E,
    8. Peña CJ,
    9. Han MH,
    10. Nestler EJ
    (2014) Loss of BDNF signaling in D1R-expressing NAc neurons enhances morphine reward by reducing GABA inhibition. Neuropsychopharmacology 39:2646–2653. doi:10.1038/npp.2014.118 pmid:24853771
    OpenUrlCrossRefPubMed
  37. ↵
    2. Lee SJ,
    3. Bora S,
    4. Austin NC,
    5. Westerman A,
    6. Henderson JMT
    (2020) Neurodevelopmental outcomes of children born to opioid-dependent mothers: a systematic review and meta-analysis. Acad Pediatr 20:308–318. doi:10.1016/j.acap.2019.11.005 pmid:31734383
    OpenUrlCrossRefPubMed
  38. ↵
    2. Leekam SR,
    3. Nieto C,
    4. Libby SJ,
    5. Wing L,
    6. Gould J
    (2007) Describing the sensory abnormalities of children and adults with autism. J Autism Dev Disord 37:894–910. doi:10.1007/s10803-006-0218-7 pmid:17016677
    OpenUrlCrossRefPubMed
  39. ↵
    2. Logan BA,
    3. Brown MS,
    4. Hayes MJ
    (2013) Neonatal abstinence syndrome: treatment and pediatric outcomes. Clin Obstet Gynecol 56:186–192. doi:10.1097/GRF.0b013e31827feea4 pmid:23314720
    OpenUrlCrossRefPubMed
  40. ↵
    2. Lomazzo E,
    3. König F,
    4. Abassi L,
    5. Jelinek R,
    6. Lutz B
    (2017) Chronic stress leads to epigenetic dysregulation in the neuropeptide-Y and cannabinoid CB1 receptor genes in the mouse cingulate cortex. Neuropharmacology 113:301–313. doi:10.1016/j.neuropharm.2016.10.008 pmid:27737789
    OpenUrlCrossRefPubMed
  41. ↵
    2. Lou HC,
    3. Henriksen L,
    4. Bruhn P,
    5. Børner H,
    6. Nielsen JB
    (1989) Striatal dysfunction in attention deficit and hyperkinetic disorder. Arch Neurol 46:48–52. doi:10.1001/archneur.1989.00520370050018 pmid:2783366
    OpenUrlCrossRefPubMed
  42. ↵
    2. Lush ME,
    3. Ma L,
    4. Parada LF
    (2005) TrkB signaling regulates the developmental maturation of the somatosensory cortex. Int J Dev Neurosci 23:523–536. doi:10.1016/j.ijdevneu.2005.04.003 pmid:16009525
    OpenUrlCrossRefPubMed
  43. ↵
    2. Lutz A,
    3. Greischar LL,
    4. Rawlings NB,
    5. Ricard M,
    6. Davidson RJ
    (2004) Long-term meditators self-induce high-amplitude gamma synchrony during mental practice. Proc Natl Acad Sci USA 101:16369–16373. doi:10.1073/pnas.0407401101 pmid:15534199
    OpenUrlAbstract/FREE Full Text
  44. ↵
    2. Mably AJ,
    3. Colgin LL
    (2018) Gamma oscillations in cognitive disorders. Curr Opin Neurobiol 52:182–187. doi:10.1016/j.conb.2018.07.009 pmid:30121451
    OpenUrlCrossRefPubMed
  45. ↵
    2. Maguire DJ,
    3. Taylor S,
    4. Armstrong K,
    5. Shaffer-Hudkins E,
    6. Germain AM,
    7. Brooks SS,
    8. Cline GJ,
    9. Clark L
    (2016) Long-term outcomes of infants with neonatal abstinence syndrome. Neonatal Netw 35:277–286. doi:10.1891/0730-0832.35.5.277 pmid:27636691
    OpenUrlAbstract/FREE Full Text
  46. ↵
    2. Maharajan P,
    3. Prencipe R,
    4. Di Francesco P,
    5. Paino G,
    6. Ravagnan G,
    7. Maharajan V
    (2000) Maternal morphine alters parvalbumin immunoreactivity patterns in neonatal mouse brain. Synapse 35:265–271. doi:10.1002/(SICI)1098-2396(20000315)35:4<265::AID-SYN4>3.0.CO;2-6
    OpenUrlCrossRefPubMed
  47. ↵
    2. McPherson C,
    3. Haslam M,
    4. Pineda R,
    5. Rogers C,
    6. Neil JJ,
    7. Inder TE
    (2015) Brain injury and development in preterm infants exposed to fentanyl. Ann Pharmacother 49:1291–1297. doi:10.1177/1060028015606732 pmid:26369570
    OpenUrlCrossRefPubMed
  48. ↵
    2. Meaney MJ
    (2001) Maternal care, gene expression, and the transmission of individual differences in stress reactivity across generations. Annu Rev Neurosci 24:1161–1192. doi:10.1146/annurev.neuro.24.1.1161 pmid:11520931
    OpenUrlCrossRefPubMed
  49. ↵
    2. Mei B,
    3. Niu L,
    4. Cao B,
    5. Huang D,
    6. Zhou Y
    (2009) Prenatal morphine exposure alters the layer II/III pyramidal neurons morphology in lateral secondary visual cortex of juvenile rats. Synapse 63:1154–1161. doi:10.1002/syn.20694 pmid:19670310
    OpenUrlCrossRefPubMed
  50. ↵
    2. Mountcastle VB
    (1998) Perceptual neuroscience. Cambridge: Harvard University Press.
  51. ↵
    2. Niu L,
    3. Cao B,
    4. Zhu H,
    5. Mei B,
    6. Wang M,
    7. Yang Y,
    8. Zhou Y
    (2009) Impaired in vivo synaptic plasticity in dentate gyrus and spatial memory in juvenile rats induced by prenatal morphine exposure. Hippocampus 19:649–657. doi:10.1002/hipo.20540 pmid:19115391
    OpenUrlCrossRefPubMed
  52. ↵
    2. O'Donnell JK,
    3. Halpin J,
    4. Mattson CL,
    5. Goldberger BA,
    6. Gladden RM
    (2017) Deaths involving fentanyl, fentanyl analogs, and U-47700 - 10 states, July-December 2016. MMWR Morb Mortal Wkly Rep 66:1197–1202. doi:10.15585/mmwr.mm6643e1 pmid:29095804
    OpenUrlCrossRefPubMed
  53. ↵
    2. Ornoy A
    (2003) The impact of intrauterine exposure versus postnatal environment in neurodevelopmental toxicity: long-term neurobehavioral studies in children at risk for developmental disorders. Toxicol Lett 140–141:171–181. doi:10.1016/S0378-4274(02)00505-2 pmid:12676464
    OpenUrlCrossRefPubMed
  54. ↵
    2. Ornoy A,
    3. Segal J,
    4. Bar-Hamburger R,
    5. Greenbaum C
    (2001) Developmental outcome of school-age children born to mothers with heroin dependency: importance of environmental factors. Dev Med Child Neurol 43:668–675. doi:10.1017/s0012162201001219 pmid:11665823
    OpenUrlCrossRefPubMed
  55. ↵
    2. Palygin O,
    3. Lalo U,
    4. Pankratov Y
    (2011) Distinct pharmacological and functional properties of NMDA receptors in mouse cortical astrocytes. Br J Pharmacol 163:1755–1766. doi:10.1111/j.1476-5381.2011.01374.x pmid:21449975
    OpenUrlCrossRefPubMed
  56. ↵
    2. Raver SM,
    3. Keller A
    (2014) Permanent suppression of cortical oscillations in mice after adolescent exposure to cannabinoids: receptor mechanisms. Neuropharmacology 86:161–173. doi:10.1016/j.neuropharm.2014.07.006 pmid:25036610
    OpenUrlCrossRefPubMed
  57. ↵
    2. Raver SM,
    3. Haughwout SP,
    4. Keller A
    (2013) Adolescent cannabinoid exposure permanently suppresses cortical oscillations in adult mice. Neuropsychopharmacology 38:2338–2347. doi:10.1038/npp.2013.164 pmid:23822952
    OpenUrlCrossRefPubMed
  58. ↵
    2. Ravikrishnan A,
    3. Gandhi PJ,
    4. Shelkar GP,
    5. Liu J,
    6. Pavuluri R,
    7. Dravid SM
    (2018) Region-specific expression of NMDA receptor GluN2C subunit in parvalbumin-positive neurons and astrocytes: analysis of GluN2C expression using a novel reporter model. Neuroscience 380:49–62. doi:10.1016/j.neuroscience.2018.03.011 pmid:29559384
    OpenUrlCrossRefPubMed
  59. ↵
    2. Ricalde AA,
    3. Hammer RP
    (1990) Perinatal opiate treatment delays growth of cortical dendrites. Neurosci Lett 115:137–143. doi:10.1016/0304-3940(90)90444-e pmid:2172870
    OpenUrlCrossRefPubMed
  60. ↵
    2. Robertson CE,
    3. Baron-Cohen S
    (2017) Sensory perception in autism. Nat Rev Neurosci 18:671–684. doi:10.1038/nrn.2017.112 pmid:28951611
    OpenUrlCrossRefPubMed
  61. ↵
    2. Rogers SJ,
    3. Ozonoff S
    (2005) Annotation: what do we know about sensory dysfunction in autism? A critical review of the empirical evidence. J Child Psychol Psychiatry 46:1255–1268. doi:10.1111/j.1469-7610.2005.01431.x pmid:16313426
    OpenUrlCrossRefPubMed
  62. ↵
    2. Rubenstein E,
    3. Young JC,
    4. Croen LA,
    5. DiGuiseppi C,
    6. Dowling NF,
    7. Lee LC,
    8. Schieve L,
    9. Wiggins LD,
    10. Daniels J
    (2019) Brief report: maternal opioid prescription from preconception through pregnancy and the odds of autism spectrum disorder and autism features in children. J Autism Dev Disord 49:376–382. doi:10.1007/s10803-018-3721-8 pmid:30132098
    OpenUrlCrossRefPubMed
  63. ↵
    2. Rubino T,
    3. Tizzoni L,
    4. Viganò D,
    5. Massi P,
    6. Parolaro D
    (1997) Modulation of rat brain cannabinoid receptors after chronic morphine treatment. Neuroreport 8:3219–3223. doi:10.1097/00001756-199710200-00007 pmid:9351646
    OpenUrlCrossRefPubMed
  64. ↵
    2. Ryan SA
    (2018) Calculating the real costs of the opioid crisis. Pediatrics 141:e20174129. doi:10.1542/peds.2017-4129
    OpenUrlFREE Full Text
  65. ↵
    2. Schulteis G,
    3. Heyser CJ,
    4. Koob GF
    (1999) Differential expression of response-disruptive and somatic indices of opiate withdrawal during the initiation and development of opiate dependence. Behav Pharmacol 10:235–242. doi:10.1097/00008877-199905000-00001 pmid:10780239
    OpenUrlCrossRefPubMed
  66. ↵
    2. Shigemoto R,
    3. Kinoshita A,
    4. Wada E,
    5. Nomura S,
    6. Ohishi H,
    7. Takada M,
    8. Flor PJ,
    9. Neki A,
    10. Abe T,
    11. Nakanishi S,
    12. Mizuno N
    (1997) Differential presynaptic localization of metabotropic glutamate receptor subtypes in the rat hippocampus. J Neurosci 17:7503–7522. pmid:9295396
    OpenUrlAbstract/FREE Full Text
  67. ↵
    2. Slotkin TA,
    3. Seidler FJ,
    4. Yanai J
    (2003) Heroin neuroteratogenicity: delayed-onset deficits in catecholaminergic synaptic activity. Brain Res 984:189–197. doi:10.1016/s0006-8993(03)03152-4 pmid:12932853
    OpenUrlCrossRefPubMed
  68. ↵
    2. Spencer MR,
    3. Warner M,
    4. Bastian BA,
    5. Trinidad JP,
    6. Hedegaard H
    (2019) Drug overdose deaths involving fentanyl, 2011-2016. Natl Vital Stat Rep 68:1–19.
    OpenUrlPubMed
  69. ↵
    2. Steriade M
    (2006) Grouping of brain rhythms in corticothalamic systems. Neuroscience 137:1087–1106. doi:10.1016/j.neuroscience.2005.10.029 pmid:16343791
    OpenUrlCrossRefPubMed
  70. ↵
    2. Sutter MB,
    3. Leeman L,
    4. Hsi A
    (2014) Neonatal opioid withdrawal syndrome. Obstet Gynecol Clin North Am 41:317–334. doi:10.1016/j.ogc.2014.02.010 pmid:24845493
    OpenUrlCrossRefPubMed
  71. ↵
    2. Tan JW,
    3. Duan TT,
    4. Zhou QX,
    5. Ding ZY,
    6. Jing L,
    7. Cao J,
    8. Wang LP,
    9. Mao RR,
    10. Xu L
    (2015) Impaired contextual fear extinction and hippocampal synaptic plasticity in adult rats induced by prenatal morphine exposure. Addict Biol 20:652–662. doi:10.1111/adb.12158 pmid:24903743
    OpenUrlCrossRefPubMed
  72. ↵
    2. Ting JT,
    3. Daigle TL,
    4. Chen Q,
    5. Feng G
    (2014) Acute brain slice methods for adult and aging animals: application of targeted patch clamp analysis and optogenetics. Methods Mol Biol 1183:221–242. doi:10.1007/978-1-4939-1096-0_14 pmid:25023312
    OpenUrlCrossRefPubMed
  73. ↵
    2. Tobon AL,
    3. Habecker E,
    4. Forray A
    (2019) Opioid use in pregnancy. Curr Psychiatry Rep 21:118. doi:10.1007/s11920-019-1110-4 pmid:31734808
    OpenUrlCrossRefPubMed
  74. ↵
    2. Van Eden CG,
    3. Uylings HB
    (1985) Postnatal volumetric development of the prefrontal cortex in the rat. J Comp Neurol 241:268–274. doi:10.1002/cne.902410303 pmid:4086657
    OpenUrlCrossRefPubMed
  75. ↵
    2. van Leeuwen TM,
    3. Neufeld J,
    4. Hughes J,
    5. Ward J
    (2020) Synaesthesia and autism: different developmental outcomes from overlapping mechanisms. Cogn Neuropsychol 37:433–449. doi:10.1080/02643294.2020.1808455 pmid:32845799
    OpenUrlCrossRefPubMed
  76. ↵
    2. Velísek L,
    3. Slamberová R,
    4. Vathy I
    (2003) Prenatal morphine exposure suppresses mineralocorticoid receptor-dependent basal synaptic transmission and synaptic plasticity in the lateral perforant path in adult male rats. Brain Res Bull 61:571–576. doi:10.1016/s0361-9230(03)00194-1 pmid:14519453
    OpenUrlCrossRefPubMed
  77. ↵
    2. Volpe DA,
    3. McMahon Tobin GA,
    4. Mellon RD,
    5. Katki AG,
    6. Parker RJ,
    7. Colatsky T,
    8. Kropp TJ,
    9. Verbois SL
    (2011) Uniform assessment and ranking of opioid μ receptor binding constants for selected opioid drugs. Regul Toxicol Pharmacol 59:385–390. doi:10.1016/j.yrtph.2010.12.007
    OpenUrlCrossRefPubMed
  78. ↵
    2. Wade CL,
    3. Schuster DJ,
    4. Domingo KM,
    5. Kitto KF,
    6. Fairbanks CA
    (2008) Supraspinally-administered agmatine attenuates the development of oral fentanyl self-administration. Eur J Pharmacol 587:135–140. doi:10.1016/j.ejphar.2008.04.007 pmid:18495108
    OpenUrlCrossRefPubMed
  79. ↵
    2. Weaver ICG,
    3. Cervoni N,
    4. Champagne FA,
    5. D'Alessio AC,
    6. Sharma S,
    7. Seckl JR,
    8. Dymov S,
    9. Szyf M,
    10. Meaney MJ
    (2004) Epigenetic programming by maternal behavior. Nat Neurosci 7:847–854. doi:10.1038/nn1276 pmid:15220929
    OpenUrlCrossRefPubMed
  80. ↵
    2. Whiteman VE,
    3. Salemi JL,
    4. Mogos MF,
    5. Cain MA,
    6. Aliyu MH,
    7. Salihu HM
    (2014) Maternal opioid drug use during pregnancy and its impact on perinatal morbidity, mortality, and the costs of medical care in the United States. J Pregnancy 2014:906723. doi:10.1155/2014/906723 pmid:25254116
    OpenUrlCrossRefPubMed
  81. ↵
    2. Wilson GS,
    3. McCreary R,
    4. Kean J,
    5. Baxter JC
    (1979) The development of preschool children of heroin-addicted mothers: a controlled study. Pediatrics 63:135–141. pmid:86983
    OpenUrlAbstract/FREE Full Text
  82. ↵
    2. Winkelman TNA,
    3. Villapiano N,
    4. Kozhimannil KB,
    5. Davis MM,
    6. Patrick SW
    (2018) Incidence and costs of neonatal abstinence syndrome among infants with Medicaid: 2004-2014. Pediatrics 141:e20173520. doi:10.1542/peds.2017-3520
    OpenUrlAbstract/FREE Full Text
  83. ↵
    2. Yanai J,
    3. Huleihel R,
    4. Izrael M,
    5. Metsuyanim S,
    6. Shahak H,
    7. Vatury O,
    8. Yaniv SP
    (2003) Functional changes after prenatal opiate exposure related to opiate receptors' regulated alterations in cholinergic innervation. Int J Neuropsychopharmacol 6:253–265. doi:10.1017/S1461145703003523 pmid:12974992
    OpenUrlCrossRefPubMed
  84. ↵
    2. Yang SN,
    3. Yang JM,
    4. Wu JN,
    5. Kao YH,
    6. Hsieh WY,
    7. Chao CC,
    8. Tao PL
    (2000) Prenatal exposure to morphine alters kinetic properties of NMDA receptor-mediated synaptic currents in the hippocampus of rat offspring. Hippocampus 10:654–662. doi:10.1002/1098-1063(2000)10:6<654::AID-HIPO1003>3.0.CO;2-T
    OpenUrlCrossRefPubMed
  85. ↵
    2. Yang SN,
    3. Huang LT,
    4. Wang CL,
    5. Chen WF,
    6. Yang CH,
    7. Lin SZ,
    8. Lai MC,
    9. Chen SJ,
    10. Tao PL
    (2003) Prenatal administration of morphine decreases CREBSerine-133 phosphorylation and synaptic plasticity range mediated by glutamatergic transmission in the hippocampal CA1 area of cognitive-deficient rat offspring. Hippocampus 13:915–921. doi:10.1002/hipo.10137 pmid:14750654
    OpenUrlCrossRefPubMed
  86. ↵
    2. Yeh ML,
    3. Selvam R,
    4. Levine ES
    (2017) BDNF-induced endocannabinoid release modulates neocortical glutamatergic neurotransmission. Synapse 71:e21962. doi:10.1002/syn.21962
    OpenUrlCrossRef
  87. ↵
    2. Zagon IS,
    3. McLaughlin PJ
    (1987) Endogenous opioid systems regulate cell proliferation in the developing rat brain. Brain Res 412:68–72. doi:10.1016/0006-8993(87)91440-5 pmid:3607463
    OpenUrlCrossRefPubMed
  88. ↵
    2. Zhu Y,
    3. Hsu MS,
    4. Pintar JE
    (1998) Developmental expression of the mu, kappa, and delta opioid receptor mRNAs in mouse. J Neurosci 18:2538–2549. pmid:9502813
    OpenUrlAbstract/FREE Full Text
  89. ↵
    2. Zucker RS,
    3. Regehr WG
    (2002) Short-term synaptic plasticity. Annu Rev Physiol 64:355–405. doi:10.1146/annurev.physiol.64.092501.114547 pmid:11826273
    OpenUrlCrossRefPubMed
Back to top

In this issue

Email
Print
View Full Page PDF
Citation Tools
Respond to this article
Perinatal Fentanyl Exposure Leads to Long-Lasting Impairments in Somatosensory Circuit Function and Behavior
Jason B. Alipio, Catherine Haga, Megan E. Fox, Keiko Arakawa, Rakshita Balaji, Nathan Cramer, Mary Kay Lobo, Asaf Keller
Journal of Neuroscience 14 April 2021, 41 (15) 3400-3417; DOI: 10.1523/JNEUROSCI.2470-20.2021
del.icio.us logo Digg logo Reddit logo Twitter logo Facebook logo Google logo Mendeley logo

Jump to section

Keywords

Responses to this article

Respond to this article

Jump to comment:

No eLetters have been published for this article.