Abstract
Psychedelics are a broad class of drugs defined by their ability to induce an altered state of consciousness1,2. These drugs have been used for millennia in both spiritual and medicinal contexts, and a number of recent clinical successes have spurred a renewed interest in developing psychedelic therapies3,4,5,6,7,8,9. Nevertheless, a unifying mechanism that can account for these shared phenomenological and therapeutic properties remains unknown. Here we demonstrate in mice that the ability to reopen the social reward learning critical period is a shared property across psychedelic drugs. Notably, the time course of critical period reopening is proportional to the duration of acute subjective effects reported in humans. Furthermore, the ability to reinstate social reward learning in adulthood is paralleled by metaplastic restoration of oxytocin-mediated long-term depression in the nucleus accumbens. Finally, identification of differentially expressed genes in the ‘open state’ versus the ‘closed state’ provides evidence that reorganization of the extracellular matrix is a common downstream mechanism underlying psychedelic drug-mediated critical period reopening. Together these results have important implications for the implementation of psychedelics in clinical practice, as well as the design of novel compounds for the treatment of neuropsychiatric disease.
Classically, psychedelics have been defined to include drugs such as lysergic acid diethylamide (LSD), mescaline, phenylcyclohexyl piperidine (PCP), ibogaine, 3,4-methylenedioxymethamphetamine (MDMA), psylocibin and ketamine, because each of these compounds produces alterations to sensory, self, time and space perception that are “so alien to everyday experience that they shed new light on the workings of these everyday mental functions”. Although more recent attempts have been made to subcategorize psychedelics on the basis of the subjective character of the altered state that they induce (for example, hallucinogenic, empathogenic, oneirogenic or dissociative), their chemical structure (for example, tryptamines, phenethylamines or arylcyclohexamines), or their principal binding target (for example, serotonin receptor 2A (5-HT2AR), monoamine transporter, κ-opioid receptor (KOR) or N-methyl-D-aspartate receptor (NMDAR)), the importance of these categories for therapeutic applications remains unclear, since psychedelics that span the diversity of classification systems have shown remarkable promise for the treatment of addiction, post-traumatic stress disorder (PTSD) and depression. Thus, identification of a common neurobiological mechanism that can account for the shared therapeutic effects of psychedelics is an obvious priority for translational neuroscience.
During specific periods of brain development, the nervous system exhibits heightened sensitivity to ethologically relevant stimuli, as well as increased malleability for synaptic, circuit and behavioural modifications. These mechanistically constrained windows of time are called critical periods and neuroscientists have long sought methods to reopen them for therapeutic benefit. Recently, we have discovered a novel critical period for social reward learning and shown that the empathogenic psychedelic MDMA is able to reopen this critical period. This mechanism shares a number of features with the therapeutic effects of MDMA-assisted psychotherapy for the treatment of PTSD, including rapid onset, durability and context dependence. At the same time, cocaine does not reopen the social reward learning critical period, and since cocaine does not share the psychedelics’ therapeutic profile, these results lend further support for the view that the reinstatement of social reward learning in adulthood underlies the therapeutic efficacy of MDMA.
Whether the ability of MDMA to reopen the critical period for social reward learning generalizes across psychedelics remains an open question. MDMA is classified as an ‘empathogen’ because its acute subjective effects are distinctly prosocial in quality. The fact that this quality is not shared by hallucinogenic psychedelics such as psilocybin and LSD, dissociative psychedelics such as ketamine, or oneirogenic psychedelics such as ibogaine challenges the idea that these drugs could reopen the social reward learning critical period. However, the psychotropic effects of MDMA include an altered state of consciousness shared by all psychedelics, and if it is this characteristic rather than its prosocial properties that embodies the subjective experience of reopening critical periods, then the ability to reinstate social reward learning in adulthood might generalize across psychedelics.
Critical period reopening is a shared property
To test whether the ability of MDMA to reopen the social reward learning critical period generalizes across psychedelics, we began by examining the effect of psilocybin pretreatment on the magnitude of social reward learning in adulthood using the social reward conditioned place preference (sCPP) assay (Extended Data Fig. 1). We administered a single intraperitoneal (i.p.) dose of psilocybin(0.3 mg kg−1) to adult male mice (at postnatal day 96 (P96)) and 48 h later (at P98), we assessed the magnitude of sCPP (Fig. 1a). Mice pretreated with psilocybin, but not saline, exhibited a significant sCPP at P98 (Fig. 1b–d). To formally designate ‘open’ and ‘closed’ states of this critical period, we next generated a natural spline regression model to previously published data with knots at P35 and P98 (P = 1.003 × 10−6; root mean square error (r.m.s.e.) = 0.19; R2 = 0.11), as shown in Extended Data Fig. 2. When compared with this derived curve, the magnitude of sCPP in saline-treated mice did not deviate significantly from the closed state (P = 0.72), whereas the fit derived from psilocybin-treated mice demonstrated a significant mean shift (P = 1.12 × 10−6) in range of the open state (Fig. 1e). Similarly, pretreatment with LSD (i.p. 1 µg kg−1) but not saline, also reopened the critical period for social reward learning (saline P = 0.90, LSD P = 1.76 × 10−9) (Fig. 1f–i). Next, we examined the effects of ketamine (i.p. 3 mg kg−1) and ibogaine (i.p. 40 mg kg−1). Mice pretreated with either drug also exhibited sCPP in adulthood (P = 8.78 × 10−4 and P = 3.17 × 10−5, respectively) (Fig. 1k–m). As with MDMA, these effects were dose-dependent (Extended Data Fig. 3). In juveniles, MDMA (i.p. 10 mg kg−1) pretreatment did not lead to a further increase the magnitude of social reward learning (Extended Data Fig. 3). In contrast to its effects on social reward learning behaviour, pretreatment with psychedelics had no effect on the magnitude of two addiction-like behaviours: cocaine reward learning and amphetamine-induced locomotor sensitization (Extended Data Fig. 4). Together, these studies demonstrate that as with empathogenic psychedelic, hallucinogenic, oneirogenic and dissociative psychedelics are able to reopen the critical period for social reward learning.'

Duration of the psychedelic open state
The duration of acute subjective effects and the durability of the therapeutic response vary considerably across psychedelics. For example, in humans, the acute subjective effects of ketamine last 30–120 min, whereas its antidepressant effectslast for 1 week. By contrast, the subjective effects of psilocybin and MDMA last for 3–6 h, whereas the acute effects of LSD and ibogaine persist for 8–10 h and 36–72 h, respectively; these long-lasting subjective effects correspond to highly durable therapeutic effects that last months to years. Previously, we showed that MDMA-induced critical period reopening lasts for two weeks, but returns to the closed state by four weeks. Here, to further probe the time course of the critical period open state induced by psychedelics, we examined the duration of critical period reopening following treatment with ketamine, psilocybin, LSD and ibogaine (Fig. 2a). One week following psychedelic treatment, psilocybin-treated mice, but not those treated with ketamine, exhibited significant social reward learning (Fig. 2b–e). Two weeks following psychedelic treatment, the social reward learning critical period remained open for both psilocybin- and LSD-treated mice (Fig. 2f–i). At three weeks, LSD-treated mice, but not those treated with psilocybin, exhibited significant social reward learning (Fig. 2j–m), whereas at four weeks, the social reward learning critical period remained open for mice treated with ibogaine but not those treated with LSD (Fig. 2n–q). For each psychedelic, we examined at least three time points; increasing the LSD dose to 50 µg kg−1 did not extend the duration of the open state (Extended Data Fig. 5). As shown in Fig. 3, the progressively longer-lasting open states induced by ketamine (Figs. 1f–i and 2b–e and Extended Data Fig. 5), followed by psilocybin (Fig. 2b–i), MDMA11 (Extended Data Fig. 5), LSD (Fig. 2j–q) and ibogaine (Fig. 2n–q and Extended Data Fig. 5) are proportional to the duration of the acute subjective effects of these drugs in humans. These results provide a mechanistic explanation for the importance of the post-treatment integration period for clinical implementation of psychedelics, and inform the design of novel compounds for clinical applications.


Metaplasticity, not hyperplasticity
Dynamic regulation of the extent to which synaptic plasticity can be induced is called ‘metaplasticity’ and is thought to be one of the mechanisms underlying the establishment of critical periods. Previously, we showed that oxytocin induces a novel form of presynaptically expressed long-term depression, and implicated this plasticity in encoding social reward learning. Here, to determine whether the ability to induce metaplastic upregulation of oxytocin plasticity generalizes across psychedelics, we pretreated adult mice with either saline, cocaine or psychedelics. Forty-eight hours or two weeks later we prepared ex vivo acute slices containing the nucleus accumbens (NAc) and conducted whole-cell voltage-clamp recordings from medium spiny neurons (MSNs) (Fig. 4a–c). A 10-min bath application of oxytocin induced a significant decrease in the frequency (Fig. 4d–k) but not the amplitude (Fig. 4l–s) of miniature excitatory post-synaptic currents (mEPSCs) following pretreatment with MDMA, LSD, psilocybin, ketamine and ibogaine, but not with saline or cocaine, at 48 h; this metaplasticity persisted for 2 weeks in the LSD pretreatment group, but not in the ketamine pretreatment groups. We did not observe significant changes in baseline mEPSC amplitude or frequency following pretreatment with psychedelics in the NAc or in layer 5 of the medial prefrontal cortex (mPFC) (Extended Data Fig. 6). Together, these results provide evidence that psychedelics induce metaplasticity rather than hyperplasticity, a distinction that is especially important for designing biomarkers to test therapeutic profiles and abuse liability of novel compounds.

5-HT2AR is not the universal mechanism
The serotonin receptor 5-HT2AR, first identified by its binding to LSD, mediates alterations of perception and cognition induced by ‘serotonergic psychedelics’ such as LSD and psilocybin. Furthermore, MDMA is thought to trigger synaptic efflux of serotonin through its binding at the serotonin transporter SERT, and some of the effects of ketamine are reportedly mediated by 5-HT2AR. Thus, we sought to determine the role of 5-HT2AR in reopening the social reward learning critical period with LSD, psilocybin, MDMA and ketamine. We administered psychedelics intraperitoneally in P96 adult mice either alone or in combination with ketanserin (HTR-A, 0.1 mg kg−1)—the 5-HT2AR antagonist used in human studies—which we injected 30 min before the psychedelic (Extended Data Fig. 7). Pre-treatment with either LSD or psilocybin induced reinstatement of sCPP measured 48 h later, and this effect was blocked by co-administration of ketanserin (Extended Data Fig. 7). However, MDMA-induced reinstatement of sCPP persisted in the presence of ketanserin (Extended Data Fig. 7). Similarly, co-administration of ketanserin did not block ketamine-induced reinstatement of social reward learning in adulthood (Extended Data Fig. 7). These results demonstrate that whereas 5-HT2ARs are required for LSD- and psilocybin-induced reopening of the social reward learning critical period (with potential contributions from serotonin 2B and 2C receptors, since ketanserin also has affinity at these serotonin receptor 2 subtypes), MDMA and ketamine reinstate social reward learning in a 5-HT2AR-independent manner. Although some have argued that psychedelics that bind 5-HT2AR (such as LSD and psilocybin) should be classified separately from those that do not (such as MDMA and ketamine), these results identify a novel property (critical period reopening) that coheres the category of psychedelics but violates the 5-HT2AR-binding boundary. Thus, combined with the data presented in Figs. 1 and 2, these results support the continued use of the established naming convention for psychedelics, rather than subclassification or renaming based on receptor binding or subjective properties.
β-arrestin-2 is not the universal mechanism
Recent studies indicate that prolonged binding at the 5-HT2AR by LSD triggers β-arrestin-2 (β-arr2)-biased signalling over canonical G-protein signalling. Moreover, the effects of MDMA and ibogaine are also thought to be mediated by metabotropic G-protein-coupled receptors (GPCRs). Although the therapeutic effects of ketamine are thought to be mediated by ionotropic NMDA receptors, the metabotropic glutamate receptor 5 has also been implicated. To test the hypothesis that β-arr2-biased signalling mediates the ability of psychedelics to reopen the social reward learning critical period, we examined their effects in commercially available β-arr2-knockout (KO) mice. We began by determining baseline sCPP in juvenile and adult β-arr2-KO mice and found that these mice exhibited the normal maturational profile of social reward learning (Extended Data Fig. 8). Next, we compared the magnitude of sCPP in adult (P98) β-arr2 wild-type and β-arr2-KO mice 48 h following administration of psychedelic drugs (Extended Data Fig. 9). LSD and MDMA reopened the social reward learning critical period in wild-type mice but did not do so in β-arr2-KO mice (Extended Data Fig. 9). Conversely, ketamine and ibogaine were able to reinstate social reward learning in both wild-type and β-arr2-KO mice (Extended Data Fig. 9). Together, these results demonstrate that whereas β-arr2 signalling is required for LSD- or MDMA-induced reopening of the social reward learning critical period, ketamine or ibogaine reinstate social reward learning in a β-arr2-independent manner.
Psychedelics induce remodelling of the ECM
Since psychedelics as a class all reopen the social reward learning critical period (Fig. 1) even though these drugs act on a diverse array of principal binding targets (Extended Data Fig. 7) and biochemical signalling pathways (Extended Data Fig. 9), we reasoned that the common mechanism that enables critical period reopening might be downstream of these cellular processes. Furthermore, given the durability of the response (Fig. 2), we hypothesized that psychedelics may modulate the expression of specific genes or pathways. To test this hypothesis, we carried out RNA sequencing of the microdissected NAc 48 h and 2 weeks after pretreatment with either saline, cocaine, ketamine, LSD or MDMA. We collected total mRNA from each sample and made strand-specific libraries for each of three replicates from each condition. Transcript-level abundances were collapsed to gene-level expression estimates for model fitting.
To directly compare treatment-related transcriptional changes specific to the shared ability of psychedelics to reopen the social reward learning critical period, we analysed the gene expression dataset between conditions in which the critical period is in the open state (48 h and 2 weeks after LSD treatment, 48 h after ketamine treatment, and 48 h after MDMA treatment) versus conditions where the critical period remains in or returns to the closed state (48 h and two weeks after saline treatment, 48 h and two weeks after cocaine treatment, and two weeks after ketamine treatment). Using this approach, we identified 65 genes that were significantly differentially expressed (likelihood ratio test; Benjamini–Hochberg-corrected q ≤ 0.1) (Fig. 5). Gene set enrichment analysis of this list identified significant enrichment of ontologies associated with endothelial development, regulation of angiogenesis, vascular development and tissue morphogenesis. Of note, many of the top scoring genes are components of the extracellular matrix (ECM) or have been implicated in its remodelling, including: Fn1), Mmp16(, Trpv4, Tinagl1, Nostrin, Cxcr4, Adgre5, Robo4 and Sema3g. Additionally, the differentially expressed gene set includes the immediate early genes (IEGs) Fos, Junb, Arc and Dusp. When we did not control for the psychedelic-specific psychoactive response (saline versus all drug conditions, including cocaine), we identified 39 differentially expressed genes (Benjamini–Hochberg-corrected q ≤ 0.15) (Extended Data Fig. 10); however, enrichment analysis identified no significant ontologies associated with this gene set, and only 6 genes (Hspa12b, Sema3g, Eng, Flt4, Cavin1 and Ube4b) overlapped with the differentially expressed genes in the open state versus closed state dataset shown in Fig. 5. These results provide evidence that the shared ability of psychedelics to reopen the social reward learning critical period converges at transcriptional regulation of the ECM. On the basis of these findings, our working model (Fig. 6) posits that psychedelics act at a diverse array of binding targets (such as SERT, 5-HT2AR, NMDA and KOR), to trigger a downstream signalling response that leads to activity-dependent (perhaps via IEG-mediated coincidence detection) degradation of the ECM, which in turn is the permissive event that enables metaplasticity. In this model, transcriptional upregulation of ECM components (for example, FN1) and downregulation of ECM proteolytic enzymes (for example, MMP-16), reflects the homeostatic response to these long-lasting cellular changes. Together, these results demonstrate novel biological effects (behavioural, temporal, electrophysiological and molecular) that—similar to therapeutic effects—are shared across psychedelics.


Conclusions
These studies provide a novel conceptual framework for understanding the therapeutic effects of psychedelics, which have shown significant promise for treating a wide range of neuropsychiatric diseases, including depression, PTSD and addiction. Although other studies have shown that psychedelics can attenuate depression-like behaviours and may also have anxiolytic, anti-inflammatory and antinociceptive properties, it is unclear how these properties directly relate to the durable and context dependent therapeutic effects of psychedelics. Furthermore, although previous in vitro studies have suggested that psychedelic effects might be mediated by their ability to induce hyperplasticity, this account does not distinguish psychedelics from addictive drugs (such as cocaine, amphetamine, opioids, nicotine and alcohol) whose capacity to induce robust, bidirectional, morphological and physiological hyperplasticity is thought to underlie their addictive properties. Moreover, our ex vivo results (Fig. 4 and Extended Data Fig. 6) are consistent with in vivo studies, which demonstrate that dendritic spine formation following administration of psychedelics is both sparse and context dependent, suggesting a metaplastic rather than a hyperplastic mechanism. Indeed, previous studies have also directly implicated metaplasticity in the mechanism of action of ketamine. At the same time, since our results show that psychedelics do not directly modify addiction-like behaviours (Extended Data Fig. 4), they provide a mechanistic clue that critical period reopening may be the neural substrate underlying the ability of psychedelics to induce psychological flexibility and cognitive reappraisal, properties that have been linked to their therapeutic efficacy in the treatment of addiction, anxiety and depression.
Although the current studies have focused on the critical period for social reward learning, critical periods have also been described for a wide variety of other behaviours, including imprinting in snow geese, song learning in finches, language learning in humans, as well as brain circuit rearrangements following sensory or motor perturbations, such as ocular dominance plasticity and post-stroke motor learning. Since the ability of psychedelics to reopen the social reward learning critical period is independent of the prosocial character of their acute subjective effects (Fig. 1), it is tempting to speculate that the altered state of consciousness shared by all psychedelics reflects the subjective experience of reopening critical periods. Consistent with this view, the time course of acute subjective effects of psychedelics parallels the duration of the open state induced across compounds (Figs. 2 and 3). Furthermore, since our results point to a shared molecular mechanism (metaplasticity and regulation of the ECM) (Figs. 4–6) that has also been implicated in the regulation of other critical periods, these results suggest that psychedelics could serve as a ‘master key’ for unlocking a broad range of critical periods. Indeed, recent evidence suggests that repeated application of ketamine is able to reopen the critical period for ocular dominance plasticity by targeting the ECM. This framework expands the scope of disorders (including autism, stroke, deafness and blindness) that might benefit from treatment with psychedelics; examining this possibility is an obvious priority for future studies.
Methods
Mice
Male wild-type mice were bred in house and weaned at 3 weeks or obtained from Jackson Laboratories (stock no. 000664). β-arr2-KO mice (stock no. 011130) were obtained from Jackson Laboratories, bred in house and weaned at 3 weeks of age. All mice were inbred to the C57BL/6J congenic ‘wild-type’ strain (as opposed to outbred ‘true wilds’, which were not used in this study). Congenic strains are generated by backcrossing for a minimum of 10 generations, a standard that is derived from the congenic interval, and the theoretical estimate that by the 10th generation, 99.99% of the congenic strain background will be from the recipient inbred. Although the β-arr2-KO mouse (Jackson Laboratories stock no. 011130), was originally derived on the 129X1/SvJ background, it was backcrossed to the C57BL/6J congenic strain at Jackson Laboratories (https://www.jax.org/strain/011130). All mice were maintained on a 12 h:12 h natural light:dark cycle, starting at 07:30 with food and water provided ad libitum. All behavioural experiments were conducted during the same circadian period (07:30–19:30) in a dedicated, sound- and odour-controlled behavioural testing room, which is separated from the vivarium, and no other experiments were conducted simultaneously in the same room. Sample size was estimated based on previous work and published literature. Experimenters were blind to the condition when subjective criteria were used as a component of data analysis, and control and test conditions were interleaved. Mice were randomly assigned to experimental and control groups. All procedures complied with the animal care standards set forth by the National Institutes of Health and were in accordance with protocols approved by the Johns Hopkins University Animal Care and Use Committee.
sCPP assay
The protocol for sCPP was adapted from previously published work. Mice were socially housed (3–5 males) in a cage containing corncob bedding (Anderson Cob, 0.25 inch cob, Animal Specialties and Provisions) until the pre-determined age for sCPP testing. Each mouse was used for only one behavioural time point. At the pre-determined age, mice were placed in an open field activity chamber (ENV-510, Med Associates) equipped with infrared beams and a software interface (Activity Monitor, Med Associates) to monitor the position of the mouse. The apparatus was partitioned into two equally sized zones using a clear Plexiglas wall, with a 5 cm diameter circular hole at the base; each zone contained one type of novel bedding (Alpha-Dri, Animal Specialties and Provisions or Kaytee Soft Granule, Petco). The amount of time spent freely exploring each zone was recorded during 30-min test sessions. For example, a score of 900 means that the mouse spent exactly 50% of its time on each of the two beddings, whereas a score of 1,800 means that it spent the full 30 min in the bedding that would be subsequently assigned as the social conditioning cue, and no time in the bedding that would be assigned as the isolation conditioning cue. After an initial pre-conditioning trial to establish baseline preference for the two sets of bedding cues, mice were assigned to receive social conditioning (with cage mates) for 24 h on one type of bedding, followed by 24 h of isolation conditioning (without cage mates) on the other bedding cue. To assure unbiased design, chamber assignments were counterbalanced for side and bedding cues. Immediately after the isolation conditioning, a 30-min post-conditioning trial was conducted to establish preference for the two conditioned cues. CPP is a learned association between a condition (for example, social) and a cue (bedding). It does not require scent from the other mice, as the bedding itself serves as the cue. Exclusion criteria for this behaviour are strictly defined as a pre-conditioning preference score of >1.5 or <0.5. Mice are never excluded based on the quality of their social interactions. Pre-conditioning versus post-conditioning social preference scores were considered significant if paired Student’s t-test P values were less than 0.05. Comparisons between experimental conditions were made using both normalized social preference scores (time spent in social zone post-treatment divided by pre-treatment) and subtracted social preference scores (time spent in social zone post minus pre); these were considered significant if unpaired Student’s t-test P values were <0.05. All experiments were performed during the mouse rest period (light cycle), since pilot experiments revealed that sCPP is most robust if assayed during this period. Prior to i.p. drug treatment experiments (MDMA, LSD, psilocybin, ketamine or ibogaine hydrochloride), mice were habituated to the injection procedure with daily saline i.p. injections in the home cage. Pharmacological delivery schedules were counterbalanced for type of drug. Unless otherwise stated (Fig. 2 and Extended Data Fig. 5), for pretreatment, experiments mice were tested 48 h after the injection to allow for complete clearance of the drug. For the experiment testing involvement of the 5-HT2AR, the 5-HT2AR antagonist ketanserin was administered i.p. 30 min prior to the injection of the drug tested.
Electrophysiology
Subjects received an i.p. injection of either LSD (1 µg kg−1), ketamine (3 mg kg−1), psilocybin (0.3 mg kg−1), MDMA (10 mg kg−1), ibogaine (40 mg kg−1) or saline. Forty-eight hours after drug treatment, either parasagittal slices containing the NAc core (250 µm thick) or coronal slices containing the PL/IL region of the mPFC (250 µm thick) were prepared from C57BL/6 mice using standard procedures. In brief, after mice were anaesthetized with isoflurane and decapitated, brains were quickly removed and placed in ice-cold low-sodium, high-sucrose dissecting solution (228 mM sucrose, 26 mM NaHCO3, 11 mM glucose, 2.5 mM KCl, 1 mM NaH2PO4, 1 mM MgSO4, 0.5 mM CaCl2). Slices were collected with a Leica VT 1200s vibrating microtome. Slices were allowed to recover for a minimum of 60 min in a submerged holding chamber (∼25 °C) containing artificial cerebrospinal fluid (ACSF) consisting of 119 mM NaCl, 2.5 mM KCl, 2.5 mM CaCl2, 1.3 mM MgCl2, 1 mM NaH2PO4, 11 mM glucose and 26.2 mM NaHCO3. For hyperplasticity recordings (Extended Data Fig. 6), slices were removed from the holding chamber and placed into the recording chamber, where they were continuously perfused with oxygenated (95% O2, 5% CO2) ACSF at 2 ml min−1 at 25 °C. For metaplasticity recordings (Fig. 4), slices were removed from the holding chamber and incubated first for 10 min in oxygenated ACSF containing picrotoxin (50 µM, Sigma), followed by 10-min incubation in oxygenated ACSF containing both picrotoxin and oxytocin (1 µM, Tocris) before being placed into the recording chamber. Whole-cell voltage-clamp recordings from MSNs or layer V pyramidal cells were obtained under visual control using a 40× objective. The NAc core was identified by the presence of the anterior commissure, and the PL/IL region of the mPFC was identified by the presence of the forceps minor of the corpus callosum. Recordings were made with electrodes (2.5–4.0 MΩ) filled with 115 mM CsMeSO4, 20 mM CsCl, 10 mM HEPES, 0.6 mM EGTA, 2.5 mM MgCl, 10 mM sodium phosphocreatine, 4 mM sodium ATP, 0.3 mM sodium GTP and 1 mM QX-314. Miniature EPSCs were collected at a holding potential of −70 mV in the presence of tetrodotoxin (0.5 μM, Tocris Biosciences) and picrotoxin (50 μM, Sigma). Two minutes after break-in, 30-s blocks of events (total of 200 events per cell) were acquired and analysed using the Recording Artist plugin in Igor Pro software with threshold parameters set at 5 pA amplitude and <3 ms rise time. All events included in the final data analysis were verified visually. Data were analysed by multivariate analysis of variance (MANOVA) with three independent variables (drug, brain area and age) and two dependent variables (frequency and amplitude). Likelihood ratio test performed comparing the full model using treatment, age, and structure to a reduced model using age and structure. All calculations were performed in either GraphPad Prism 9 or the R programming language and are available as Supplementary Code 1 and in the repository at https://github.com/genesofeve/DolenPsychedelicOpenState.
RNA extraction and sequencing
Male wild-type C57BL6/J mice were injected i.p. with LSD, ketamine, cocaine (20 mg kg−1) or saline solution either 2 weeks or 48 h before the mice were euthanized. At P98 to P112, mice were euthanized, brains were rapidly removed and ~1mm thick coronal slice (n = 3 mice per condition) containing the nucleus accumbens were sectioned using a mouse brain matrix. To microdissect the NAc, slices were placed in a petri dish containing ice-cold ACSF (125 mM NaCl, 2.5 mM KCl, 2 mM CaCl2, 1 mM MgCl2, 1.25 mM NaH2PO4, 10 mM glucose and 26 mM NaHCO3) supplemented with RNase inhibitor and oxygenated with carbogen gas (95% O2 and 5% CO2) to pH 7.3–7.4. The NAc was identified using the anterior commissure and other structural markers. Between each dissection, blades were replaced and all the instruments and the matrix were cleaned with a solution containing RNase inhibitor. Following dissection, tissue was immediately placed into 0.5 ml Trizol and subjected to a 15 s burst with a tissue homogenizer to lyse the cells. Samples were kept on ice prior to storage at −20 °C. Total RNA were extracted using the RNeasy Kit from Qiagen. The quality of purified RNA was assessed via both a nanodrop and 2100 Bioanalyzer from Agilent. Library preparation was performed using a TruSeq Stranded mRNA kit (Illumina) using the recommended protocol. Individual dual-indexed libraries were quality controlled, pooled, and sequenced on the NovaSeq 6000 platform on a single S1 flowcell to an average depth of 76,841,745 (±8,066,939.82) paired-end 100 bp reads per sample. Reads were pseudoaligned to the mouse GENCODE vM25 (ref. 71) reference transcriptome using kallisto (v0.46.2) with 100 bootstrapped samples and 6 threads. Defaults were used for all other parameters. Estimated transcript-level abundances were collapsed to gene-level expression estimates and analysed using the sleuth (v0.30.0) R/Bioconductor package. To identify genes with differential expression as a function of samples where the critical period is reopened we performed a likelihood ratio test comparing a full model which included batch, and critical period to a reduced model that only included batch. Time was not used as an explanatory variable in this model fitting. Using this test, we identified 65 genes as significantly differentially expressed at a 10% false discovery rate (Benjamini–Hochberg-corrected q ≤ 0.1). To identify genes with differential expression as a function of any drug treatment (including cocaine) versus saline we performed a likelihood ratio test comparing a full model that included batch, and ‘treated vs untreated’ to a reduced model that only included batch. Using this test, we identified 39 genes as significantly differentially expressed at a 15% false discovery rate (Benjamini–Hochberg-corrected q ≤ 0.15). Time was not used as an explanatory variable in this model fitting. Raw data will be made publicly available (Gene Expression Omnibus accession numbers: GSE230679 and GSM7231202–GSM7231228). Code to reproduce the RNA-seq analysis and associated figures is provided as Supplementary Code 2 and in the repository at https://github.com/genesofeve/DolenPsychedelicOpenState.
Statistics
All statistical details can be found in the figure legends, including the type of statistical analysis used, P values, n, degrees of freedom, t values and f values. Sample sizes were not predetermined by statistical methods; instead they were estimated based on the previously published literature. Data distributions were assumed to be normal. Homogeneity of variance was tested using Levene’s test for equality of variances. Comparisons between experimental manipulations were made using a two-tailed Students t-test (paired or unpaired, and with or without Welch’s correction as appropriate) and MANOVA for comparisons between multiple outcome measures, with P < 0.05 considered significant.
Linear, β-spline, loess smoothing and natural spline models evaluated on the previously published time course of normalized social preference scores. Loess smoothing yielded a pseudoinverse at age 41.695 and a knot point of 35 was chosen for both β-spline and natural spline models. The natural spline outperformed the β-spline (adjusted R2 of 0.1053 versus 0.5554, respectively) with fewer parameters. Residuals were plotted against fitted values and age to check model assumptions. Leave one out cross validation was also used to assess model fit. Control data from all new experiments was used as test data via the predict R function. RSME and R2values were comparable between the original model and the new data. Two-way t-tests to compare means of controls groups against matched or binned time periods was done to confirm fit to new data. The full model including coefficients for splines, experiment and condition was constructed and tested against reduced models with the final reduced model being reported. MANOVA analysis was carried out using multivariate linear models and the ANOVA function. All statistical comparisons were carried out in the R programming language and can be found in Supplementary Codes 3 and 4 as well as in the repository at https://github.com/genesofeve/DolenPsychedelicOpenState.