Ketamine disinhibits dendrites and enhances calcium signals in prefrontal dendritic spines

A subanesthetic dose of ketamine causes acute psychotomimetic symptoms and sustained antidepressant effects. In prefrontal cortex, the prevailing disinhibition hypothesis posits that N-methyl-d-aspartate receptor (NMDAR) antagonists such as ketamine act preferentially on GABAergic neurons. However, cortical interneurons are heterogeneous. In particular, somatostatin-expressing (SST) interneurons selectively inhibit dendrites and regulate synaptic inputs, yet their response to systemic NMDAR antagonism is unknown. Here, we report that ketamine acutely suppresses the activity of SST interneurons in the medial prefrontal cortex of the awake mouse. The deficient dendritic inhibition leads to greater synaptically evoked calcium transients in the apical dendritic spines of pyramidal neurons. By manipulating NMDAR signaling via GluN2B knockdown, we show that ketamine’s actions on the dendritic inhibitory mechanism has ramifications for frontal cortex-dependent behaviors and cortico-cortical connectivity. Collectively, these results demonstrate dendritic disinhibition and elevated calcium levels in dendritic spines as important local-circuit alterations driven by the administration of subanesthetic ketamine.

K etamine is a non-competitive N-methyl-d-aspartate receptor (NMDAR) antagonist. In humans, systemic administration of ketamine at a subanesthetic dose have both acute and sustained effects. Acutely, within about 1 h following administration, subjects experience psychotomimetic symptoms, perceptual aberrations, and cognitive impairments 1 . Because the acute effects mimic the positive and negative symptoms of schizophrenia, ketamine has found use as a tool for studying glutamatergic dysfunction and evaluating antipsychotic compounds 2,3 . On a longer timescale, starting from 2 to 4 h following administration, patients suffering from major depressive episodes report a robust antidepressant effect that remain for about a week 4,5 . Although these behavioral effects are well established and have substantial translational implications, the mechanisms by which ketamine impacts neural activity remain unclear.
In rodents, intraperitoneal or subcutaneous injection of subanesthetic ketamine activates multiple, distributed brain regions, as indicated by brain-wide mapping of metabolic activity 6 . Among the regions with elevated metabolic activity is the medial prefrontal cortex, including the ventral prelimbic and infralimbic sub-regions, as well as the dorsal cingulate and secondary motor sub-regions (Cg1/M2). This locus is consistent with earlier works demonstrating that, for NMDAR antagonists, there is a dosedependent increase in the extracellular levels of glutamate in the rat prefrontal cortex 7 , as well as an increase in the overall firing rate of prefrontal cortical neurons 8 . Since NMDAR is an ionotropic glutamate receptor that contributes to postsynaptic depolarization, why would an antagonist lead to more spiking activity? One prevailing idea is the disinhibition hypothesis, which posits that NMDAR antagonists at low dose may act preferentially on GABAergic neurons, thereby disinhibiting excitatory neurons 9,10 . In support of the disinhibition hypothesis, fast-spiking GABAergic neurons, identified based on their narrow spike waveforms, exhibit a concomitant reduction in firing in vivo at the time of prefrontal hyperactivity during NMDAR hypofunction 11 .
However, GABAergic neurons are not homogeneous. Fastspiking cells constitute a fraction of the GABAergic neurons, and another non-overlapping subpopulation is the somatostatinexpressing (SST) interneurons. SST interneurons are a major subpopulation that constitutes 20% of the GABAergic neurons in the mouse frontal cortex 12 . Unlike the fast-spiking cells which inhibit cell bodies, SST interneurons include the so-called Martinotti cells that selectively target apical dendrites and thus are poised to control synaptic inputs and dendritic integration 13 . The cell-type-specific differences suggest that if SST interneurons were altered, then the resulting circuit-level and behavioral impacts may be quite distinct from those involving perturbations of fastspiking cells. Still unknown, however, is whether and how SST interneurons and dendritic inhibition may be impacted by subanesthetic ketamine.
In this study, we use two-photon calcium imaging to characterize SST interneurons, SST axons, pyramidal neurons, and pyramidal dendritic spines in awake, head-fixed mice. For each cellular compartment, we compare spontaneous calcium transients before and within an hour after the administration of subanesthetic ketamine, relative to saline control. We find that ketamine suppresses the activity of SST interneurons in the medial prefrontal cortex. We further show that this loss of dendritic inhibition contributes to elevated synaptic calcium transients in the apical dendritic spines of pyramidal neurons. Based on these results, we conclude that the dendritic disinhibition and calcium elevations in dendritic spines are important components of ketamine's actions on the prefrontal cortex.

Results
Cell-type-specific responses to subanesthetic ketamine. We performed two-photon microscopy on awake, head-fixed mice (Fig. 1a), targeting the Cg1/M2 sub-regions of the medial prefrontal cortex (Fig. 1b). Initially, while mice were head-fixed under the two-photon microscope, we recorded body motion using an infrared camera. This was because systemic administration of subanesthetic ketamine induces hyperlocomotion in rodents 14 , and we wanted to avoid movement as a confound in our imaging experiments. We observed that ketamine (10 mg/kg, s.c.) increased body motion but only transiently, and therefore limited all of our data collection to 30-60 min post-injection (Fig. 1c, d). For calcium imaging, we used AAV1-CamKII-GCaMP6f-WPRE-SV40 to express the calcium-sensitive fluorescent protein GCaMP6f in pyramidal neurons in Cg1/M2, and imaged spontaneous fluorescence transients from the awake mouse (Fig. 1e). From the fluorescence transients, we detected calcium events using a peeling method based on template matching 15 . Ketamine (10 mg/kg) increased the rate of spontaneous calcium events in pyramidal neuron cell bodies in layer 2/3 (200-400 µm from the dura) (ketamine: 23.7 ± 2.1%, saline: 9.4 ± 1.9%, relative to pre-injection, mean ± s.e.m.; P = 3 × 10 −8 , twosample t-test; Fig. 1f, g; Supplementary Fig. 1a). The elevation in calcium event rates did not correlate with animal movement ( Supplementary Fig. 1b) and was consistent across individual mice ( Supplementary Fig. 2a). Ketamine's effect lasted up tõ 50 min post-injection ( Supplementary Fig. 3a). Somatic calcium transients have been shown to directly relate to the firing rate of cortical neurons 16,17 , therefore the elevated calcium event rates reflect hyperactivity of pyramidal neurons, consistent with previous reports 11,18 . To characterize the effect of ketamine on the SST subpopulation of GABAergic neurons, we injected AAV1-Syn-DIO-GCaMP6s-WPRE-SV40 into Cg1/M2 of SST-IRES-Cre mice to express GCaMP6s selectively in SST interneurons (Fig. 1h). Following ketamine, the rate of spontaneous calcium events was reduced for the cell bodies of layer 2/3 SST interneurons (ketamine: −13 ± 3%, saline: 13 ± 6%, relative to preinjection, mean ± s.e.m.; P = 1 × 10 −4 , two-sample t-test; Fig. 1i, j, Supplementary Fig. 1c, d; Supplementary Fig. 2b; Supplementary  Fig. 3b). These results indicate that ketamine increases the activity of pyramidal neurons while suppressing the activity of SST GABAergic neurons.
Ketamine disinhibits dendrites and raises calcium in spines. SST interneurons inhibit selectively dendrites including the dendritic spines of pyramidal neurons 13 . The ketamine-induced reduction of SST interneuron activity is therefore expected to diminish dendritic inhibition. To determine if this may apply to the apical dendritic compartment, we imaged the inhibitory inputs-SST axonal boutons in the superficial layer of Cg1/M2 (<200 µm from the dura) (Fig. 2a), which showed fewer calcium events after ketamine (10 mg/kg) ( Fig. 2b; Supplementary Fig. 2c; Supplementary Fig. 3c). We also imaged directly the postsynaptic calcium signals in the apical dendrite tufts of pyramidal neurons. We focused on the dendritic spines (Fig. 2c), because subthreshold synaptic activation can be characterized by their accompanying calcium elevations in the spine compartment with minimal outflow to the dendritic shaft 19 . To restrict our analyses to these localized fluorescence transients, we used a standard procedure to regress out the contribution from the dendritic shaft 17,19 , following which we applied the same calcium event detection procedure. Ketamine (10 mg/kg) induced a 43.42 ± 0.01% increase in spontaneous calcium event rate for dendritic spines in Cg1/M2, significantly higher than saline controls (4.34 ± 0.01%, P = 0.02, two-sample t-test; Fig. 2c Supplementary Fig. 3d). To characterize the temporal dynamics of the aberrant transients, we identified calcium events that co-occurred in each image frame, thereby separating amplitude (number of events in each frame) from frequency (number of frames with at least one event). Both the amplitude and frequency parameters were elevated following ketamine (Fig. 2e). The results did not depend on the specifics of the template in the event detection procedure, because an alternative analysis based on threshold crossing yielded similar conclusions ( Supplementary Fig. 4). Interestingly, the ketamine-induced  Elevated synaptic calcium responses to long-range inputs. By focusing on fluorescence transients that are localized to dendritic spines, we are characterizing calcium signals related to synaptic activation. It follows then that a more stringent and explicit test for a synaptically activated signal to directly stimulate the afferent inputs. To test this, we take advantage of the knowledge that apical dendritic spines receive long-range cortico-cortical inputs 20 , and Cg1/M2 receives an abundance of monosynaptic excitatory inputs from the retrosplenial cortex (RSC). We implanted electrodes to stimulate RSC (biphasic pulse with duration of 10 ms and a peak amplitude of ±150 µA), and simultaneously imaged the evoked calcium responses in Cg1/M2 (Fig. 4a) Fig. 4b). Ketamine (10 mg/kg), but not saline, increased evoked responses in Cg1/M2 dendritic spines (ketamine: F(6,1146) = 6.6, P = 7 × 10 −7 , saline: F(6,1218) = 20.0, P = 0.06, two-way interaction between time (pre vs. post) and stimulation strength, two-way mixed ANOVA after three-way mixed ANOVA; Fig. 4c, d). Across all imaged spines, the increase in evoked calcium response was more pronounced at higher stimulation strengths ( Fig. 4d; Supplementary  Fig. 5f-g). By directly controlling the activation of long-range inputs, these results provide evidence that ketamine elevates the synaptically evoked component of calcium transients in prefrontal dendritic spines.
Ketamine's effect on dendrites rely on GluN2B in SST neurons. We wanted to determine the causal relations between ketamine, SST interneuron activity, and synaptic calcium signals. We predict that if ketamine exerts its effects through SST interneurons, then knocking down NMDAR signaling in SST interneurons should reproduce the local-circuit alterations and render additional ketamine ineffective. To test this prediction, we created AAV1-CMV-dsRed-pSico-GluN2BshRNA to induce Credependent expression of short hairpin RNA (shRNA) against GluN2B ( Fig. 5a) 21,22 . We targeted the GluN2B subunit because neocortical SST interneurons contain abundant GluN2B transcripts 23,24 , and the subunit contributes a sizable fraction of NMDAR-mediated current in the related low-threshold spiking interneurons in the rodent medial prefrontal cortex 25 . This virus was confirmed to reduce GluN2B expression in a Cre-dependent manner in vitro (−42 ± 3% relative to controls; Supplementary  Fig. 6) and to decrease NMDAR-mediated currents in brain slices 26 . Injected into Cg1/M2 of SST-IRES-Cre animals (Fig. 5b), viral-mediated knockdown of GluN2B in SST interneurons (GluN2B-SST KD) reduced SST interneuron activity at baseline (GluN2B-SST KD: 1.5 ± 0.2 Hz, controls: 1.9 ± 0.1 Hz, mean ± s.e. m.; P = 0.04, two-sample t-test; Fig. 5c), mimicking the effects of ketamine. Notably, GluN2B-SST KD blocked further ketamineinduced (10 mg/kg) influence on SST interneuron activity (Fig. 5d). In a separate group of animals, we imaged dendritic spines, finding that GluN2B-SST KD increased the rate of spontaneous calcium events at baseline (GluN2B-SST KD: 0.85 ± 0.03 Hz, controls: 0.73 ± 0.04 Hz; P = 0.01, two-sample t-test; Fig. 5e), and additional ketamine (10 mg/kg) had no detectable effect (Fig. 5f). The causal manipulations thus confirm the predicted role of SST interneurons and highlight the involvement of GluN2B in mediating ketamine's effect on the dendritic spines.
Behavioral consequence of downregulating SST neuron activity. Ketamine exerts numerous acute behavioral effects, so next we asked to what extent dendritic disinhibition in the medial prefrontal cortex may be involved. We tested animals on three behavioral assays: trace fear conditioning, pre-pulse inhibition and locomotor activity, and determined the effect of ketamine with or without bilateral GluN2B-SST KD (Fig. 6a). In trace fear conditioning, an auditory stimulus is paired with a footshock with an intervening trace period. This assay requires associative learning across a temporal gap, which relies on the medial prefrontal cortex 27 . SST-IRES-Cre animals received ketamine (10 mg/kg) 15 min before conditioning, and subsequently showed impaired freezing 24 h later in the 15-s trace but not in the delay (no trace) paradigm (Fig. 6b), indicating an impairment specific to temporal association but not learning in general. After bilateral GluN2B-SST KD in Cg1/M2, ketamine no longer had detectable effects on trace fear learning (Fig. 6b). Pre-pulse inhibition, a cross-species measure of sensorimotor gating, is regulated by prefrontal cortical activity 28 . Ketamine (40 mg/kg) reduced prepulse inhibition in control SST-IRES-Cre animals (main effects of b Coronal histological section, showing the extent of AAV-mediated expression of GCaMP6f in the mouse medial prefrontal cortex. Cg1, cingulate cortex. M2, secondary motor cortex. PrL, prelimbic cortex. c Example motion traces from a head-fixed animal during a two-photon imaging session before and after saline injection (black), and another session before and after ketamine (10 mg/kg) injection (red). Gray shading denotes the 5 min during which the injection was made and imaging was done to verify that the field of view had not shifted. d Time-averaged motion (mean ± s.e.m.) for epochs including preinjection, 5-30 min post-injection, and 30-60 min post-injection for saline (left) and ketamine (right). For saline injection (black), there were no detectable differences between epochs (pre-injection vs. 5-30 min post-injection: P = 0.6; pre-injection vs. 30-60 min post-injection: P = 0.1; Wilcoxon signed rank test). For ketamine (red), hyperlocomotion was detected transiently following injection (pre-injection vs. 5-30 min post-injection: P = 0.06; pre-injection vs. 30-60 min post-injection: P = 1; Wilcoxon signed rank test). Each point is an imaging session. n = 5 animals each for saline and ketamine. e Schematic of imaging location, and an in vivo two-photon image of GCaMP6s-expressing pyramidal neurons in Cg1/M2. Inset, magnified view of neuronal cell bodies. f The normalized difference in the rate of spontaneous calcium events of pyramidal neurons. Normalized difference was calculated as post-injection minus pre-injection values normalized by the pre-injection value (ketamine (10 mg/kg): 23.7 ± 2.1%, saline: 9.4 ± 1.9%, mean ± s.e.m.; P = 3 × 10 −8 , two-sample t-test). For ketamine, n = 613 cells from 5 animals. For saline, n = 681 cells from 5 animals. g Each row shows time-lapse fluorescence transients from the same pyramidal cell in the pre-injectiion (left) and post-injection (right) periods. Two example cells were plotted for saline injection (black), and two other examples were plotted for ketamine (10 mg/kg) injection (pre-injection: black; post-injection: red). h-j Same as (e-g) for GCaMP6s-expressing SST interneurons in Cg1/M2 of SST-IRES-Cre animals (ketamine (10 mg/kg): −12 ± 3%, saline: 13 ± 6%, mean ± s.e.m.; P = 1 × 10 −4 , two-sample t-test). For ketamine, n = 198 cells from 5 animals. For saline, n = 179 cells from 5 animals. *P < 0.05; **P < 0.01; ***P < 0.001; n.s., not significant.
manipulations, an acute downregulation of prefrontal SST interneuron activity fully occludes the behavioral effects of ketamine on trace fear learning and pre-pulse inhibition.

Discussion
The goal of this study was to delineate the impact of subanesthetic ketamine on the prefrontal cortical circuit. There were two major findings. One, ketamine reduces the activity of dendrite-targeting SST interneurons. Two, as a consequence of the diminished dendritic inhibition, ketamine increases the amplitude of the synaptically evoked calcium transients in the apical dendritic spines of pyramidal neurons. As we will elaborate below, these local-circuit alterations provide insights into the neural mechanisms for the acute psychotomimetic and sustained antidepressant effects of ketamine.
The novelty of our results is to highlight SST interneurons as a key target of systemic NMDAR antagonism. Although in principle the disinhibition hypothesis is general and need not be tied to particular GABAergic neuron subtypes, so far empirical studies of NMDAR hypofunction have focused on fast-spiking cells for technical reasons. For example, recordings of inhibitory postsynaptic currents have provided clear evidence for disinhibition 32,33 , but the method cannot distinguish the cell types that provide the GABAergic inputs, and an electrode placed at the soma tends to favor the detection of perisomatic inputs. In vivo electrophysiology has also supported disinhibition by showing reduced activity for fast-spiking cells 11 , but SST interneurons are likely not examined in this earlier study because their spike waveforms have overlapping features with excitatory neurons and thus cannot be easily isolated 16 . By contrast, using a GABAergic Credriver mouse line, our fluorescence signals would originate predominantly from SST interneurons. The ability to characterize responses from genetically identified SST interneurons distinguishes our results from prior studies of ketamine and neural activity. Freezing behavior for control SST-IRES-Cre animals in delay fear conditioning with no trace period (saline: 47 ± 6%, ketamine (10 mg/kg): 50 ± 5%, mean ± s.e.m.; P = 0.5, Wilcoxon rank-sum test; n = 8 and 12 animals for saline and ketamine, respectively), in trace fear conditioning (saline: 41 ± 7%, ketamine (10 mg/kg): 2 ± 4%, mean ± s.e.m.; P = 6 × 10 −6 , Wilcoxon rank-sum test; n = 8 animals each for saline and ketamine), and for GluN2B-SST KD animals in trace fear conditioning (saline: 33 ± 7%, ketamine (10 mg/kg): 43 ± 6%, mean ± s.e.m.; P = 0.2, Wilcoxon rank-sum test; n = 8 animals each for saline and ketamine). c Pre-pulse inhibition as a measure of sensorimotor gating. For control SST-IRES-Cre animals, a two-way within-subjects ANOVA with prepulse intensity (3, 6, and 9 dB) and drug (saline, ketamine (40 mg/kg)) as within-subjects factors found significant main effects of pre-pulse intensity (F(2,12) = 27.5, P = 4 × 10 −6 ) and drug (F(1,24) = 2.0, P = 2 × 10 −4 ), but a non-significant interaction (F(2,24) 1.8, P = 0.18). Post-hoc Tukey-Kramer's tests for saline vs. ketamine were significant at 3 dB (40 ± 7%, 5 ± 6%, mean ± s.e.m., P = 0.01), at 6 dB (59 ± 5%, 17 ± 8%, mean ± s.e.m., P = 9 × 10 −4 ), at 9 dB (67 ± 4%, 38 ± 6%, mean ± s.e.m., P = 0.003). n = 13 animals each for saline and ketamine. For GluN2B-SST KD animals, two-way within-subjects ANOVA revealed a significant main effect of pre-pulse intensity (F(2,24) = 6.3, P = 1.7 × 10 −4 ), but non-significant drug effect (F(1,12) = 4.0, P = 0.93) or two-way interaction (F(2,24) = 1.4, P = 0.26). Saline vs. ketamine yielded no difference at 3 dB (32 ± 9%, 39 ± 6%, mean ± s.e.m.), at 6 dB (58 ± 6%, 56 ± 6%, mean ± s.e.m.), and at 9 dB (67 ± 6%, 63 ± 6%, mean ± s.e.m.). n = 13 animals each for saline and ketamine. Perisomatic and dendritic disinhibition may represent separate, parallel ways by which ketamine acts on the prefrontal cortex. This is probable because when we manipulated dendritic inhibition, the behavioral and electrophysiological alterations are largely orthogonal to those previously reported for NMDAR manipulations in fast-spiking cells. More specifically, when the obligatory GluN1 subunit is deleted from fast-spiking cells, behavioral deficits were reported for delay fear conditioning, but not for pre-pulse inhibition 34 . This is opposite to what we saw. In another study with postnatal GluN1 removal in many fast-spiking cells, although there was a deficit in pre-pulse inhibition similar to our animals, the NMDAR antagonist-induced locomotion was reduced 35 , unlike our manipulation of SST interneurons which had no effect on the ketamine-induced hyperlocomotion. Finally, enhanced baseline gamma-band oscillations were observed following NMDAR deletion in fast-spiking cells 34,36 and may even be produced by removing NMDARs in pyramidal neurons 37 , whereas we show that dendritic inhibition was involved in controlling the envelope of oscillatory signals on a much slower timescale. This comparison comes with the caveat that we are interpreting behavioral outcomes from studies involving NMDAR manipulations to different extent and at different developmental time points. Notwithstanding the caveat, it seems possible that ketamine's suppressive effects on perisomatic and dendritic inhibition may be responsible for non-overlapping sets of behavioral and electrophysiological phenotypes.
NATURE COMMUNICATIONS | https://doi.org/10.1038/s41467-019-13809-8 ARTICLE the inhibitory influence of SST interneurons is expected to be greater, as the cells have high spontaneous firing rates in vivo 40 , setting the stage for the suppression by ketamine. The positive skew in the ketamine-induced increase of synaptic calcium transients (Fig. 2d) suggests that only a subset of dendritic spines were impacted by ketamine. This agrees with anatomical studies, showing that only~20-50% of cortical dendritic spines are innervated by SST axonal boutons 38,41,42 . By contrast, for SST interneuron soma (Fig. 1i) and axonal boutons (Fig. 2b), ketamine appears to induce a subtractive shift of the entire activity distribution, which is consistent with most SST neurons having NMDAR-mediated currents 25 . One intriguing observation is that the synaptic calcium elevations were observed in the medial prefrontal cortex, but not the primary motor cortex. This may be because although the systemic administration is expected to deliver ketamine to all brain regions, the density of cortical SST interneurons is highest in the frontal associations regions 43 . The gradient of cortical interneuron density may explain the heterogeneous, region-specific responses to a global perturbation 44 . How might the dendritic disinhibition relate to the behavioral effects of subanesthetic ketamine? Acutely, ketamine produces psychotomimetic symptoms, perceptual aberrations, and cognitive impairments. In the mouse, the Cg1/M2 sub-regions of the medial prefrontal cortex are involved in linking antecedent events to current decisions 45 . Dendritic inhibition is thought to be crucial for the context-dependent gating of inputs 46 . Dysregulation of dendritic inhibition is therefore expected to allow aberrant sensory inputs into the prefrontal cortex, thereby disorganizing the task-specific activity patterns of prefrontal cortical neurons 47 . Experimentally, this matches with our finding of ketamine-induced prefrontal cortical hyperconnectivity. Our mouse experiments were motivated by human neuroimaging studies that have found abnormally high functional connectivity for the prefrontal cortex in humans with ketamine or early-course schizophrenia 48,49 . Crucially, here we have extended to test the potential neural mechanism by showing that the prefrontal cortical hyperconnectivity in the mouse depends on the ability of ketamine to perturb dendritic inhibition.
What about the antidepressant effects of subanesthetic ketamine? One striking consequence of ketamine is the synaptic adaptations in the medial prefrontal cortex that can be detected 24 h later, including increased dendritic spine density and strengthened synaptic responses [50][51][52] . It is thought that these positive influences on synapse function might counteract the impaired synaptic transmission observed in rodents after chronic stress 53  and depressed subjects 54,55 . One possibility is that the heightened calcium transients in dendritic spines within the first hour of systemic administration is important for synaptic plasticity. A potential link between the acute calcium elevations and subsequent synaptic strengthening is supported by studies that have shown postsynaptic calcium influx to be required and sufficient to induce activity-dependent synaptic plasticity 56,57 . It agrees with the observation that the antidepressant effect is blunted if ketamine is administered while excitatory transmission is blocked 50,58 , and that manipulations of SST interneurons can tune depressive-like behaviors in rodents 21,26,59 . SST interneurons are also known to connect densely with pyramidal networks 60 and thereby regulate their correlated firing 61 , which may relate to a very recent observation that ketamine restores coordinated activity in prefrontal ensembles 52 . Therefore, there are reasons to suspect that dendritic disinhibition might be relevant for ketamine's antidepressant effects. Additional studies will be required to more firmly establish the potential relationship. Surgery. For all stereotaxic surgical procedures, the mouse was first anesthetized with isoflurane in oxygen (3-4% during induction, 1-1.5% for the remainder of the surgery) while placed in a stereotaxic apparatus (David Kopf Instruments). The animal lay on top of a water-circulating pad (Gaymar Stryker) set at a constant temperature of 38°C. Eyes were lubricated with ophthalmic ointment. Carprofen (5 mg/kg, s.c.; 024751, Henry Schein Animal Health) and dexamethasone (3 mg/kg, i.m.; 002459, Henry Schein Animal Health) was used pre-operatively. The skin above the skull was then sterilized by swabs of ethanol and betadine before an incision was made to expose the skull. We made a craniotomy (small or large depending on the experiment as detailed below) using a dental drill. Immediately after surgery and on each day for the following 3 days, carprofen (5 mg/kg, s.c.) was given to the animal. All animals had at least 1 week to recover before starting experiments unless stated otherwise. For imaging, we made a small skin incision and craniotomy (~0.5 mm diameter) in the right hemisphere centered at 1.5 mm anterior of bregma (anteriorposterior, AP) and 0.3 mm lateral of midline (medial-lateral, ML) to target the medial prefrontal cortex (Cg1 and M2 sub-regions) 63 . A glass pipette (3-000-203-G/X, Drummond Scientific) was pulled to a fine tip (P-97 Flaming/Brown Micropipette Puller, Sutter Instrument), front-filled with the relevant viruses (see below) and then slowly lowered into the brain using the stereotaxic apparatus at 4 sites corresponding to vertices of a 0.2-mm wide square centered at the above coordinates. At depth of 0.3 to 0.5 mm below the dura, we ejected the virus using a microinjector (Nanoject II, Drummond Scientific). To target the retrosplenial cortex (RSC) and primary motor cortex (M1), we used the same approach centering at AP = −1.4 mm/ML = 0.3 mm, and AP = 1.8 mm/ML = 1.6 mm, respectively. We injected 9.2 or 18.2 nl per injection pulse, with~30 s in between each injection pulse. To reduce backflow of the virus, we waited at least 5 min after finishing injection at one site before retracting the pipette to move on to the next site. For experiments that required more than one virus, we injected each virus sequentially. The brain was kept moist with artificial cerebrospinal fluid (ACSF, in mM: 5 KCl, 5 HEPES, 135 NaCl, 1 MgCl2, 1.8 CaCl2; pH 7.3). After completing all injections, the craniotomy was covered with silicone elastomer (0318, Smooth-On, Inc.), and the skin sutured (Henry Schein). After~3 weeks, the animal underwent a second surgery for cranial window implant. After the initial steps as outlined above, an incision was made to remove all of the skin above the skull, which was cleaned to remove connective tissues. We made a 3-mm diameter craniotomy around the previously targeted location. A glass window was put together by gluing two 3-mmdiameter, #1-thickness circular glass (640720, Warner Instruments), using UVsensitive optical adhesive (NOA61, Norland Products). The glass window was then placed on the exposed brain with slight downward pressure while applying highviscosity adhesive (Loctite 454) around the perimeter of the glass window. After the adhesive cured, a custom-made stainless steel headplate (eMachineShop.com) was affixed to the skull using quick adhesive cement (C&B Metabond, Parkell).

Methods
For behavior, we followed the surgical procedures for imaging, with the following exceptions. There was no second surgery for cranial window and head plate implant. All injections were bilateral injections, in which we used one of the coordinates above, plus the mirrored coordinates for the other hemisphere. We injected at all the locations in one hemisphere first, before front-filling more virus (if necessary), and moving onto the other hemisphere. We waited at least 4 weeks before behavioral testing.
For local field potential recordings, we fabricated bundles of stainless steel wire electrodes (2-3 wires per bundle) (790500, A-M Systems). The diameter of an electrode was 114.3 µm and 50.8 µm with and without coating, respectively. There was~350 µm of exposed tip for each electrode which was separated by~250 µm from other electrodes in the bundle. Impedance was typically less than 100 kΩ when tested at 1 kHz. Following the procedures above for surgery and small craniotomy, a bundle of electrodes was each lowered into the right Cg1/M2 and right RSC, centered at the stereotaxic coordinates as described above, to a depth of 0.5 mm below the dura. We used a miniature stainless steel screw (#0000-160) placed over the right cerebellum as the reference electrode. All electrodes were then soldered to a miniature connector (A79002-001 or A79000-001, Omnetics). Care was taken to avoid shorted connections. All pieces of the implant along with a stainless steel headplate were affixed to the skull using cyanoacrylate (Loctite 454) and quick adhesive cement (C&B Metabond, Parkell). For experiments combining viral injections and electrophysiology, we first injected the virus, and then waited 3 weeks before proceeding with a second surgery for electrode implant.
For electrical stimulation, we first injected viruses using the procedures above. Three weeks later, we performed a second surgery to implant stimulating electrodes. A bundle of stainless steel wires (3-4 wires per bundle) (790500, A-M Systems) was inserted into RSC. The wire implant was soldered to a miniature connector (A79002-001 or A79000-001, Omnetics). A glass-covered cranial window was made over Cg1/M2 using procedures described above. All pieces of the microstimulation implant and a stainless steel headplate were affixed to the skull using cyanoacrylate (Loctite 454) and quick adhesive cement (C&B Metabond, Parkell).
Histology. Mice were transcardially perfused with paraformaldehyde solution (4% (v/v) in phosphate-buffered saline). The brains stayed in the fixative for at least 48 h, and then were sectioned coronally with a vibratome and mounted on slides with coverslips. We imaged the sections with an upright fluorescence microscope (Zeiss Axio Imager M2).
Drugs. We used saline (0.9% Sodium chloride, Hospira). Ketamine working solution (1 mg/ml in saline) was prepared from stock solution (McKesson), approximately every 30 days. All ketamine and saline injections were subcutaneous (s.c.). All imaging, locomotor activity and fear conditioning experiments used 10 mg/kg. Only the PPI experiments used a higher dose of 40 mg/kg. Initial testing at 10 mg/kg in control mice did not yield reliable effects on PPI performance. We therefore used a higher dosage, following a prior study 65 . Clozapine N-oxide (CNO) (Sigma Millipore) working solution (1 mg/ml) was prepared fresh daily before each experiment from a stock solution dissolved in DMSO (10 mg/ml) stored at 4°C. Vehicle (1:10 DMSO in saline) was used in control injections. All CNO and vehicle injections were intraperitoneal (i.p.).
Two-photon calcium imaging. We used a laser-scanning two-photon microscope (Movable Objective Microscope, Sutter Instrument) controlled by the ScanImage software 66 . A tunable Ti:Sapphire femotosecond laser (Chameleon Ultra II, Coherent) was the excitation source which was focused onto the brain using a water immersion objective (XLUMPLFLN, ×20/0.95 N.A., Olympus). The timeaveraged laser power was typically less than 100 mW during experiments. Filters had the following center excitation (ex) and emission (em) wavelengths (Semrock or Chroma): GCaMP6s/f, 920 nm (ex) and 525 nm (em); dsRed, 1040 nm (ex) and 605 nm (em). Emitted photons were collected by GaAsP photomultiplier tubes. We acquired images that were 256 × 256 pixels at 0.68 or 0.54 µm per pixel resolution for dendritic spines and axonal boutons or 1.35 µm per pixel resolution for cell bodies. Frame rate was 3.62 Hz using bidirectional scanning on a set of galvanometer-based scanners. A software program (Presentation, NeuroBehavioral Systems) sent TTL pulse via a data acquisition device (USB-201, Measurement Computing) to the two-photon microscope at 30-s intervals, as well as to trigger all the other equipment. The Presentation software would log timestamps so that data from imaging and other equipment (e.g., microstimulation and body motion) can be synchronized in later analyses. Mice were habituated to head fixation under the microscope for 3-5 days of increasing durations before the commencement of any data collection. While head fixed, the mouse would sit in an acrylic tube, which permits postural adjustments but limits gross body movements. Subcutaneous injections were given through a drilled hole in the tube so as to allow the mouse to remain head fixed. To examine the acute effects of ketamine, we would image a field of view for 10-15 min to obtain 'pre-injection' data. The imaging was then interrupted for ketamine (10 mg/kg) or saline (10 ml/kg) injection (s.c.). At 30 min after the injection, we would verify that the same field of view was still in focus, and then obtain imaging data for another 10-15 min. Imaging 30 min after the injection ensured that we bypass the transient hyperlocomotion effects of ketamine (Fig. 1c,  d). In a subset of animals, we imaged continuously for 60 min after injection, but analyzed only the data after the 30-min post-injection mark. Each subject was imaged for ketamine and saline injections, with the order counterbalanced across subjects. We imaged at 0-200 µm below dura for dendritic spines, and 200-500 µm below dura for cell bodies. We targeted the medial portion of Cg1/M2, within 0-500 µm of the midline as visualized by the dark band of the midsagittal sinus.
For imaging experiments with cell-type-specific knockdown experiments, during live imaging, we searched for and imaged fields of view with many cells doublelabeled with dsRed and GCaMP6s.
Analysis of calcium imaging data. We first concatenated all of the .tiff image files from the same experiment and corrected for lateral motion using either TurboReg implemented as a plugin in ImageJ or NoRMCorre 67 in MATLAB. A custom graphical user interface (GUI) in MATLAB was used for the manual selection of region of interests (ROI). The rest of the analysis depended on experiment type as described below.
For dendritic spines, we manually scrolled through the imaging frames to look for ROIs that likely represented dendritic spines and dendritic shaft-neurite segments with multiple protrusions that display correlated patterns of fluorescence transients. The dendritic shaft ROIs were selected to be within 20 µm of the corresponding dendritic spine ROIs. We averaged values of pixels within a ROI to generate F ROI (t). For each ROI, we quantified the contribution of neuropil to the fluorescence signal. Specifically, we would calculate a radius, r, treating the ROI area as the area of a circle, and then creating an annulus-shaped neuropil ROI with inner and outer radii of 2r and 3r, respectively. We excluded pixels that belonged to the ROIs of other spines and shaft. To also exclude pixels that may belong to unselected dendritic structures, we calculated the time-averaged signal for each pixel, and then determined the median for all pixels within the neuropil ROI. We then excluded pixels if their time-averaged signal was higher than the median. We finally averaged across the non-excluded pixels in the neuropil ROI to generate F neuropil (t). To subtract the neuropil signal, we used the formula: where c = 0.4 is the neuropil correction factor. The fractional change in fluorescence, ΔF/F(t), for each spine and shaft was calculated as follows: where F 0 (t) is the 10th percentile of F(t) within a 2-min sliding window. Because we were presently interested in calcium signals due to subthreshold synaptic inputs, we regressed out fluorescence signals from non-local sources 17 . For each spine with ΔF/F spine (t), we computed ΔF/F synaptic (t) by subtracting out the scaled version of fluorescence transients in the corresponding dendritic shaft, ΔF/F shaft (t) , using the following equation: where α was determined by a linear regression forced through the origin of ΔF/F spine (t) vs. ΔF/F shaft (t).
For axonal boutons, we drew ROIs that are likely single axonal boutonsvaricosities on a relatively thin neurite branch with no protrusions. This manual identification is supported by the known neuroanatomy. For SST interneurons, they send axons but very few dendrites, if any, to layer 1. For retrosplenial cortical neurons, they send long-range axons to layer 1 of Cg1/M2. To generate F ROI (t), values of pixels within an ROI were averaged. For each ROI, we would estimate the contribution of neuropil to the fluorescence signal, by calculating and subtracting F neuropil (t) to obtain F candidate_bouton (t), using the procedures described above. Some boutons may belong to the same neuron. To limit duplicates, we used a modified a published procedure 68 . Briefly, we first analyzed a subset of imaging experiments to identify pairs of boutons on the same branch, based on distance (within 10 µm of each other) and visual confirmation of a connecting axonal segment. From this subset of bouton pairs (n = 46 bouton pairs, 2 animals), we used their F candidate_bouton (t) values to compute a correlation coefficient for each pair. The 5th percentile value for all the correlation coefficients in this subset was 0.59, which was used as the cut-off threshold value. Then, for each axonal imaging experiment, we performed an iterative procedure to select unique boutons. In each iteration, a correlation matrix using F candidate_bouton (t) between all pairs of candidate boutons in a field of view was computed. The pair of candidate boutons with the highest correlation coefficient was identified. If the value exceeded the threshold of 0.59, we would select one out of that pair of candidate boutons randomly, and looked for other boutons whose correlation coefficient with the selected bouton exceeded the threshold of 0.59. The original identified pair and boutons from subsequent search then would form one cluster. Entries in the correlation matrix corresponding to this cluster would be excluded for subsequent iterations. The process continued again to identify another cluster until no pair of candidate boutons in the correlation matrix had value above the threshold of 0.59. At the end of the iterative procedure, each of the remaining boutons was assigned to be its own cluster. We randomly selected one bouton as the representative for each cluster. The fluorescence signals of these representative boutons, F bouton (t), were then converted to fractional changes in fluorescence, ΔF/F bouton (t), as described above and used for subsequent analyses.
For cell bodies, we selected ROIs that are likely cell bodies-fluorescent objects with a round boundary and a lateral extent of about 10 µm. Values of pixels within a ROI were averaged to generate F ROI (t). For each ROI, we would remove the contribution of neuropil to the fluorescence signal, by calculating and subtracting F neuropil (t) to obtain F soma (t), and then calculating the fractional changes in fluorescence ΔF/F soma (t), using procedures outlined in the previous sections. For experiments with cell-type-specific GluN2B manipulation, we selected only GCaMP6s-expressing cells that were co-labeled with dsRed.
After obtaining ΔF/F synaptic (t), ΔF/F bouton (t), or ΔF/F soma (t) for a ROI, we did the following analysis to detect calcium events. We employed a validated "peeling" algorithm 15 , which used an iterative, template-matching procedure to decompose ΔF/F(t) into a series of elementary calcium events. We used a template for an elementary calcium event with instantaneous onset, an amplitude of 0.3, and a decay time constant of 1 s for a single-exponential function. The algorithm would search for a match to this specified template in ΔF/F(t). When a match occurred, a calcium event was recorded, and the template subtracted from ΔF/F(t) (i.e., peeling). The remaining ΔF/F trace was then searched for another match, in an iterative manner until no further matches were found. The times of recorded calcium events formed the output of the algorithm. The temporal resolution of the output is limited by the imaging frame rate. There can be multiple events associated with the same event time. In the paper describing the peeling algorithm, extensive numerical simulations suggested excellent detection performance, even in worstcase scenarios with low signal-to-noise ratios and low sampling rate 15 . For each imaging session, we calculated calcium event rate by dividing the number of calcium events by the duration of the imaging period. For Figs. 2e, 3c, and Supplementary Fig. 1a, c, we characterized the calcium events in terms of amplitude and frequency. We identified calcium events that co-occurred in the same imaging frame (same event time), distinguishing amplitude (mean number of calcium events per frame, for frames with at least one event) from frequency (number of frames with at least one event divided by the duration of the imaging session). We imaged the same compartment (spine, bouton, or cell body) before and after injection, and therefore were able to quantify the change in calcium event rate, amplitude, and frequency for each compartment by calculating the postinjection minus pre-injection values normalized by the pre-injection value.
To determine if our results depend on the event detection procedure, we used an alternative method to quantify calcium transients in dendritic spines ( Supplementary Fig. 4). We started from ΔF/F synaptic (t). For each ROI and for each condition (i.e., pre-injection or post-injection), we would calculate a threshold corresponding to 3 times the median absolute deviation of ΔF/F synaptic (t). We noted imaging frames in which ΔF/F synaptic (t) was above the threshold. From this subset of frames, we determined the mean amplitude above threshold (ΔF/F synaptic (t) minus the threshold), and the fraction of time spent above the threshold (the number of above-threshold frames divided by the number of all frames) (see Supplementary Fig. 4a for an illustration). We quantified the changes in these measures-mean amplitude above threshold and fraction of time above threshold -by calculating the post-injection minus pre-injection values normalized by the pre-injection value.
Recording animal movements during two-photon imaging. To assess the mouse's movements during two-photon imaging (Fig. 1c, d), we illuminated the animal with infrared light and recorded videos (640 × 480 pixels, 20 Hz) of the animal's face and body using an infrared camera (See3CAM_12CUNIR, e-con Systems). Although the animal was head-fixed, it could still move parts of its body such as the face, trunk, and limbs. To synchronize the infrared video with imaging, the infrared illuminator was triggered to turn off for 100 ms by the same TTL pulse used to trigger the two-photon microscope. During analysis, the interval between dark frames therefore indicated a duration of 30 s. We computed motion m in 1-s bins: where N is the number of video frames in the 1-s bin, X and Y are the width and height of the frame, p x,y (f + 1) and p x,y (f) are values of the pixel at coordinates x and y in frames f + 1 and f, respectively.
Simultaneous electrical stimulation and imaging. Imaging in Cg1/M2 was performed as described above. For microstimulation, we injected electrical currents into the implanted stimulation electrodes in RSC, using an analog stimulus isolator (Model 2200, A-M Systems). The timing and strength of the stimulation were specified by a script written in a software program (Presentation, NeuroBehavioral Systems) via command voltages out of a multifunctional input/output device (USB-6001, National Instruments). During each experiment, to calibrate the stimulation current (which could vary due to depleting battery for the stimulus isolator or high impedance of the implanted electrode), an electronic circuit was placed in-series, allowing us to measure the stimulation current online and tune the electrical current to within 5% of the desired level. For the stimulation pattern, each pulse had a biphasic waveform, lasted 5 ms in each phase, and had peak amplitudes of ±150 µA. We used 1, 2, 4, 8, 16, 32, or 64 pulses as the stimulation levels. If more than 1 pulse was used, the pulses were applied one after the other with no interpulse interval. A block consisted of the 7 stimulation levels, presented in random order with a random inter-trial period of 10-14 s. For each experiment, we presented 10 blocks successively before injection of ketamine or saline, and then presented another 10 blocks successively 30-60 min after injection. Prior to the post-injection stimulation, we imaged to confirm that the field of view has not changed. Using procedures described in earlier sections, we determined ΔF/F in each dendritic spine and axonal bouton. The log files generated by the Presentation software were used to align the timestamps of the imaging data to the stimulation events. For each dendritic spine or axonal bouton, we determined the ΔF/F measured within 1 s following the onset of stimulation, and reported the averaged values across the 10 blocks.
Electrophysiology. Animals were head-fixed and placed in a Faraday cage. An interface cable was mated to the connector on the animal's head in order to amplify local field potential (LFP) signals using a differential amplifier (RHD2216, Intan Technology), high-pass filtered at 0.5 Hz and digitized at 1 kHz (RHD2000, Intan Technology). Before actual electrophysiological recordings, mice were habituated to head-fixation for 3-5 days of increasing durations. We recorded for 30 min preinjection, and then for 60 min post-injection. We discarded data from electrodes that were noisy (impedance >1 GΩ). One electrode from each site (Cg1/M2 or RSC) was chosen randomly for analysis. For analysis, we compared the 30-min segment recorded pre-injection to the 30-min segment recorded from 30-60 min post-injection. We focused on the integrated gamma band power, because it was previously found to correlate with hemodynamic signals underlying bloodoxygenation-level-dependent resting-state connectivity 31 . To determine the integrated gamma band power, spectrograms were computed from the local field potential recordings using a sliding window with duration of 2 s and step size of 0.4 s. At each window, spectral power within the gamma band from 30 to 80 Hz was summed, and this yielded the integrated gamma band power. Correlation between integrated gamma band power in Cg1/M2 and RSC was calculated based on Pearson correlation coefficient at zero lag. Coherence between integrated gamma band power in Cg1/M2 and RSC was computed using the Chronux 69 package in MATLAB, and the magnitude of coherence was reported.
Behavior. Fear learning was assessed using a conditioning box utilizing an infrared video camera (320 × 240 pixels, 30 Hz; Med Associates, Inc.). To enable more flexible protocols involving multiple stimuli and repetitions, we modified the box to be controlled by a behavioral control software (Presentation, NeuroBehavioral Systems). Each of the conditioned stimuli (CS) was a 20-s long series of auditory pips (500 ms on, 500 ms off). The pip was either 2.5 kHz (85 dB, calibrated by A-weighting at 15 cm from speaker) or 11 kHz (75 dB). For delay conditioning experiments, one CS (CS+) co-terminated with the unconditioned stimulus of footshock (US; 0.65 mA, 0.5 s), and the other CS (CS−) was not associated with any event. For trace conditioning experiments, the end of CS+ and US was separated by 15 s of no stimulus. In the acquisition phase, animals were injected with either ketamine (10 mg/kg, s.c.) or saline (10 ml/kg, s.c.), returned to the home cage, and then placed in the conditioning box 15 min later. For DREADD experiments, CNO (5 mg/kg, i.p.) or vehicle (10 mg/ml, i.p.) was injected 30 min before trace fear conditioning followed by ketamine or saline 15 min before commencement of conditioning. After 6 min of habituation, conditioning started during which 5 CS− and 5 CS+ were presented in a random order with a random inter-trial period of 100-110 s. We counter-balanced across subjects the auditory characteristics of CS+ and CS−, as well as the ketamine and saline injections. In the recall phase, which occurred about 24 h later, we tested the animals in a different context. Context was altered using different odors, wall inserts, light, and ambient noise. After 6 min of habituation, testing commenced, during which 4 CS− were presented followed by 4 CS+ (with no US) with a random inter-trial period of 100-110 s. Motion was quantified from the infrared video data using automated procedures in the Video Freeze software (Med Associates, Inc.). Briefly, according to the publication related to the software 70 , motion was quantified by calculating the sum of pixel-by-pixel value changes across successive video frames. Time periods when the motion values fell below the software-preset threshold for at least 1 s were considered freezing periods. For 35 s after the onset of each CS+ or CS− (i.e., the 20-s long CS presentation plus the 15-s long trace period), we determined the fraction of time in which freezing was detected. We reported the differential value, i.e., the fraction of time spent freezing in response to CS+, subtracted by the fraction of time spent freezing in response to CS−, as measured in the recall phase.
To measure sensorimotor gating, we used a commercial system (SR-Lab, San Diego Instruments). An animal would be injected with ketamine (40 mg/kg, s.c.) or saline (10 ml/kg, s.c.) and then immediately placed in a cylinder on a platform with an accelerometer for recording motion, all of which lie inside a soundproof chamber. For DREADD experiments, CNO (5 mg/kg, i.p.) or vehicle (10 mg/ml, i.p.) was given 30 min before injection of ketamine or saline and immediate commencement of PPI. The animal was habituated for 5 min to background noise (66 dB, calibrated by A-weighting at 25 cm from speaker) and startle stimuli (120 dB noise, 50 ms). For testing, there were 5 types of trials: startle, pre-pulse trials during which the onset of a pre-pulse stimulus (3, 6, or 9 dB noise above background, 50 ms) preceded a startle stimulus by 100 ms, and catch trials with no stimulus. A block consisted of 5 trials, 1 of each type, in a random order with a random inter-trial period of 8-23 s. For each animal, we would test for 10 blocks (i.e., 50 trials in total). The pre-pulse inhibition (PPI) is calculated as follows: where S startle and S prepulse were the rectified motion signals recorded by the accelerometer during the 50-ms long startle stimulus averaged across trial repeats, for the startle and prepulse trials, respectively. To assess locomotor activity, we placed animals in a conditioning box equipped with an infrared video camera (320 × 240 pixels, 30 Hz; Med Associates, Inc.). After 30 min of habituation, they were injected with either ketamine (10 mg/kg, s.c.) or saline (10 ml/kg, s.c.) and returned to the box for another 60 min. Motion was assessed from the infrared video data quantified using automated procedures in the Video Freeze software (Med Associates, Inc.). To plot the results for Fig. 6d, we smoothed the motion signals with a 5-min long moving window at 1 s steps, for visualization purposes. For statistical comparisons, we used the timeaveraged motion for each of the epochs: pre-injection, 5-30 min post-injection and 30-60 min post-injection.
Statistics. Statistical tests as described in figure legends were done in MATLAB. For the imaging experiments, all planned two-group comparisons were done with parametric t-tests (two-tailed, α = 0.05). For experiments with multiple parameters and conditions, we employed ANOVA tests, with groupings as indicated in the figure legend, followed by post-hoc tests using Tukey-Kramer for statistically significant main or interaction effects. The sampling distribution of the mean was assumed to be normal, but this was not formally tested. To test the robustness of our results, we also ran all the statistical tests in Figs. 1-3 using the two-sample Kolmogorov-Smirnov (K-S) test. This nonparametric test allows us to compare the two cumulative fraction plots (saline vs. ketamine) directly without any assumption of normality. None of the conclusions of whether to reject the null hypothesis at α = 0.05 were affected by a change in statistical tests (data not shown). For behavior and electrophysiology experiments, we used non-parametric tests (twotailed, α = 0.05). In the figures, P-values were represented as n.s., not significant; P > 0.05; *P < 0.05; **P < 0.01; ***P < 0.001. See figure legends for exact P-values.
Reporting summary. Further information on research design is available in the Nature Research Reporting Summary linked to this article.

Data availability
The data that support the findings of this study are available from the corresponding author upon reasonable request.