Classic and dissociative psychedelics induce similar hyper-synchronous states in the cognitive-limbic cortex-basal ganglia system
This rat study finds that the firing rates change differently between classical psychedelics and dissociative psychedelics. Compared to the non-psychedelic control, both increased the prevalence of brain waves in groups of brain cells (high-frequency oscillations in local fields potential).
Authors
- Barrientos, S. A.
- Brys, I.
- Halje, P.
Published
Abstract
The neurophysiological mechanisms behind the profound changes in perception and cognition induced by psychedelic drugs are not well understood. To identify neuronal activity specific to the psychedelic state, we here investigated the effects of classic psychedelics (LSD, DOI) and dissociative psychedelics (ketamine, PCP) on neuronal firing rates and local field potentials in several brain structures involved in cognitive processing in freely moving rats. The classic psychedelics had a net inhibitory effect on firing rates of putative interneurons and principal cells in all recorded regions. The dissociative psychedelics had a similar inhibitory effect on principal cells, but an opposite excitatory effect on interneurons in most regions. However, the inhibitory effect on principal cells was not specific to the psychedelic state, as similar inhibition occurred with a non-psychedelic psychotropic control (amphetamine). In contrast, both types of psychedelics dramatically increased the prevalence of high-frequency oscillations (HFOs) in local field potentials, while the non-psychedelic control did not. Further analysis revealed strong HFO phase locking between structures and very small phase differences corresponding to <1 ms delays. Such standing-wave behavior suggests local generation of HFOs in multiple regions and weak, fast coupling between structures. The observed HFO hypersynchrony is likely to have major effects on processes that rely on integration of information across neuronal systems, and it might be an important mechanism behind the changes in perception and cognition during psychedelic drug use. Potentially, similar mechanisms could induce hallucinations and delusions in psychotic disorders and would constitute promising targets for new antipsychotic treatments.
Research Summary of 'Classic and dissociative psychedelics induce similar hyper-synchronous states in the cognitive-limbic cortex-basal ganglia system'
Introduction
Brys and colleagues frame the study around two linked gaps in knowledge: first, how acute psychedelic states alter neuronal activity to produce profound changes in perception and cognition; and second, what common neural mechanisms, if any, exist across phenomenologically similar but pharmacologically distinct psychedelics. Earlier work shows classic psychedelics (for example LSD, DOI) primarily act via 5-HT2A receptors and dissociative psychedelics (ketamine, PCP) via NMDA antagonism, yet both classes have been associated with increased glutamatergic signalling in cortico-limbic networks. Prior electrophysiological evidence is sparse and partly contradictory: single-cell firing changes reported in awake animals differ between drug classes, while studies of local field potentials (LFPs) have repeatedly identified aberrant high-frequency oscillations (HFOs) after both ketamine and 5-HT2A agonists. It therefore remains uncertain which physiological signatures are specific to the psychedelic state and which reflect generic psychotropic effects such as altered locomotion or arousal. This study set out to identify neuronal activity patterns that are specific to the psychedelic state by recording single-unit firing and LFPs simultaneously from multiple cortico-limbic and basal ganglia structures in freely behaving rats. Classic psychedelics (LSD, DOI) and dissociative psychedelics (ketamine, PCP) were compared with a non-psychedelic psychostimulant control (amphetamine). The investigators aimed to determine whether (a) firing-rate changes or (b) population-level oscillatory phenomena best characterise the psychedelic state, and to map the spatial extent and inter-areal relationships of any psychedelic-specific oscillations.
Methods
Nine adult Sprague-Dawley rats were implanted with custom 128-wire microelectrode arrays and allowed at least two weeks of recovery. Electrode positions were verified post mortem in five animals using computed tomography; electrode tips received fifty unique anatomical labels that were subsequently grouped into ten broader regions by presumed functional similarity. Recordings targeted multiple cortical and subcortical areas implicated in cognitive processing, including olfactory bulb/cortex, ventral striatum, orbitofrontal cortex, medial prefrontal cortex and integrative thalamus. Animals were recorded in a round open-field arena. After about 60 minutes of baseline recording, each rat received an intraperitoneal injection of one of: LSD (0.3 mg/kg), DOI (2 mg/kg), ketamine (25–50 mg/kg), PCP (5 mg/kg) or d-amphetamine (4 mg/kg). Recordings continued for another 60–120 minutes. Baseline measures were averaged over −35 to −5 minutes relative to injection and on-drug measures over 30 to 60 minutes post-injection. Behaviour was scored from video for 1 minute every 10 minutes using ordinal prevalence ratings, and locomotion/centre occupancy were quantified from centroid tracking. Head-twitch responses (HTR) were detected with onboard accelerometers where available and by manual video scoring for recordings lacking accelerometry. Neural data acquisition used Neuralynx or OpenEphys/Intan systems to capture LFPs and single-unit activity. Spike sorting used hierarchical clustering with a 2 ms refractory-period constraint; units were classified into putative principal cells, interneurons or left unclassified when assignment probability was <0.75. To emphasise local sources, bipolar LFP time series were computed from electrode pairs within the same structure. The IRASA method separated rhythmic peaks from fractal background in power spectra, with least-squares fitting used to parametrically define peak height and frequency. Instantaneous phase and amplitude of HFOs were extracted by bandpass-filtering ±5 Hz around the median HFO frequency for each recording and applying the Hilbert transform. Amplitude auto- and cross-correlations and circular statistics (resultant vector lengths, mean phase differences) quantified co-modulation and phase relationships. Phase inversion between nearby electrodes was used as evidence for local current dipoles. Granger causality in the frequency domain was estimated from bipolar LFPs using BSMART routines with 500 ms windows and model order 5; total Granger causality between regions was summarised as the median amplitude of causality peaks. Statistical testing largely used nested ANOVA and non-parametric tests (Wilcoxon) where stated.
Results
Behavioural effects followed expected class-specific patterns. NMDA antagonists (ketamine, PCP) induced hyperlocomotion and ataxia, 5-HT2A agonists (LSD, DOI) produced head-twitch responses, and amphetamine produced hyperlocomotion with stereotypy (behavioural differences significant at p<0.001). No single overt behaviour was similarly and specifically altered by both 5-HT2A agonists and NMDA antagonists, indicating that any psychedelic-specific neural signatures are unlikely to be trivially driven by shared motor changes. Single-unit results came from 365 units recorded across 169 sites. At the global scale, 5-HT2A agonists produced widespread inhibition: 65% of cells showed significant inhibition and only 8% excitation. NMDA antagonists produced a mixed pattern (about 40% inhibited and 40% excited), while amphetamine produced 32% inhibition and 46% excitation. When units were grouped by anatomical region and cell class, 5-HT2A agonists induced net inhibitory modulation in both putative principal cells and interneurons across examined structures. By contrast, NMDA antagonists tended to inhibit principal cells while exciting interneurons in most structures. Comparing psychedelic drugs as a group with amphetamine, a psychedelic-specific firing-rate difference was observed only in integrative thalamus (psychedelics: mean z-scored change −0.3; amphetamine: +1.6; p<0.05). Cluster analyses quantified the population-response dissimilarities: the Bhattacharyya distance between 5-HT2A and NMDA response patterns was 6.6, while NMDA and amphetamine were much more similar (distance 0.43), and 5-HT2A versus amphetamine had distance 0.98. The most consistent and specific electrophysiological signature of the psychedelic state was a dramatic increase in high-frequency oscillations (HFOs) around ~150 Hz. Both 5-HT2A agonists and NMDA antagonists markedly increased HFO amplitude and prevalence across multiple structures, whereas amphetamine increased broad gamma (30–80 Hz) but did not generate HFOs. HFO detection rates shifted markedly upward in psychedelic conditions, and a classification of prevalence (persistent, prevalent, occasional, absent) produced an almost identical spatial pattern for 5-HT2A and NMDA conditions, with amphetamine resembling baseline. Analyses of HFO dynamics revealed spindle-like amplitude modulation with an autocorrelogram full-width at half-maximum of 50±18 ms. Instantaneous-amplitude cross-correlations were stronger within structures (median 0.82±0.10) than between structures (median 0.40±0.15; p<0.001), although some between-structure pairs (notably olfactory cortex, ventral striatum and orbitofrontal cortex, and mPFC–ventral striatum) showed co-modulation comparable to within-structure pairs. Simultaneous HFO frequencies across structures were tightly matched at any given moment despite variability across time and subjects. Phase analyses across 6237 electrode pairs with detectable HFOs found non-random phase relations in 86% (kappa>1), and among those 95% had absolute phase differences <π/4; in practice the mean phase deviations between the olfactory bulb and other structures ranged from 0.001 to 0.45 radians, corresponding to temporal delays of less than 1 ms. Phase inversion—interpreted as evidence of local current dipoles—was observed in about 2% of intrastructural electrode pairs and appeared in multiple regions (olfactory bulb/cortex, orbitofrontal cortex, mPFC, ventral striatum), supporting local HFO generation in several areas. Granger causality analyses showed that during the psychedelic state the spectral peak of causal influence shifted from low gamma at baseline to the HFO band. The mean Granger causality in the HFO band increased by 97% versus baseline (p<0.001). Structure-specific measures indicated the strongest directional influence was from ventral striatum to medial prefrontal cortex; among all pairwise directions only vStr→mPFC had Granger causality consistently above zero.
Discussion
Brys and colleagues interpret their findings as demonstrating a dissociation between single-cell firing-rate changes and the population-level dynamics that most specifically characterise the acute psychedelic state. Although classic and dissociative psychedelics had different effects on firing rates—5-HT2A agonists producing widespread inhibition and NMDA antagonists producing mixed inhibition/excitation with interneuron excitation—the salient common physiological signature was the emergence of widespread, phase-synchronised HFOs around 150 Hz. This hypersynchrony was absent with the non-psychedelic psychostimulant amphetamine, indicating specificity to the psychedelic state rather than to general arousal or locomotor change. The authors position the HFO hypersynchrony as a plausible mechanism by which psychedelics could disrupt information integration across brain systems, because tight phase synchronisation and near-zero phase lags across multiple cortico-limbic and basal ganglia structures would profoundly affect timing-dependent computations. They outline one mechanistic hypothesis linking receptor-level effects to oscillatory dynamics: both drug classes decrease NMDA-mediated currents and increase AMPA-mediated currents, and because AMPA receptors have much faster deactivation kinetics (order of a few milliseconds versus tens of milliseconds for NMDA), the net effect could be to speed up excitatory postsynaptic temporal dynamics and thereby facilitate the emergence of faster stable oscillatory states. At the network level, the observed near-zero phase lags and the presence of local phase inversions argue against a single dominant source (for example the olfactory bulb) simply propagating HFOs unidirectionally to other regions. Instead, the data are more consistent with multiple locally generated oscillators that are weakly and rapidly coupled, producing standing-wave-like synchronous states with very small inter-regional delays. The authors note that mathematical models of coupled oscillators can produce such outcomes even with local connectivity and delayed interactions, although the biophysical substrate for rapid long-range synchronisation remains to be clarified; gap junctions, ephaptic coupling or specialised fast pathways are discussed as possibilities but not resolved by the present data. Key limitations and uncertainties acknowledged by the investigators include the remaining puzzle of how such fast oscillations synchronise across macroscopic distances given known synaptic and axonal delays, and the variability in HFO prevalence across structures and individuals. They also emphasise that this work is an initial step toward integrating single-cell pharmacology, network dynamics and global brain states; causal links between HFO hypersynchrony and subjective aspects of the psychedelic experience, or their relevance to human psychoses, require further study. Finally, the authors suggest that HFO hypersynchrony might represent a mechanistic target for understanding hallucinations and delusions in psychotic disorders and for the development of novel antipsychotic interventions, while cautioning that more mechanistic and translational work is needed to substantiate such implications.
View full paper sections
INTRODUCTION
Converging evidence show that psychedelics are effec2ve at trea2ng several neuropsychiatric condi2ons and that their therapeu2c effect depends mainly on their ability to induce neuroplas2city. Less is known about how psychedelics alter brain ac2vity to induce the acute changes in percep2on and cogni2on that they are most known for. The acute psychedelic state is important to study in a medical context, both for its poten2al contribu2on to long-term therapeu2c effects but also as a model for psychosis. More fundamentally, psychedelicsinduced changes in brain ac2vity might reveal processes important for the study of consciousness. Psychedelic drugs are primarily classified phenomenologically based on their ability to induce a psychedelic experience. Nevertheless, it is well established that classic psychedelics like lysergic acid diethylamide (LSD) and 2,5-dimethoxy-4-iodo-amphetamine (DOI) exert their effects mainly through their agonis2c ac2on on 5-HT2A receptors, while dissocia2ve psychedelics, such as ketamine and phencyclidine (PCP), act mainly as non-compe22ve antagonists on N-methyl-D-aspartate (NMDA) glutamate receptors. Despite these differences, increased glutamatergic neurotransmission in cor2colimbic networks has been iden2fied as a common downstream effect linked to the psychedelic state. In par2cular, the psychedelic state has been linked to glutamate-dependent depolarizing membrane currents in a subpopula2on of pyramidal cells in medial prefrontal cortex (mPFC;. However, this increased glutamate signaling does not result in general network excita2on and the effect on neuronal spiking ac2vity is complex. The very few inves2ga2ons performed in awake animals to date indicate that the highly selec2ve 5-HT2A agonist DOI decreases spiking ac2vity in the orbitofrontal cortex, anterior cingulate cortex and motor cortex of rodents. In contrast, of 3 28 (which was not certified by peer review) is the author/funder. All rights reserved. No reuse allowed without permission. The copyright holder for this preprint this version posted September 28, 2022. ;doi: bioRxiv preprint the NMDA antagonists ketamine, PCP and MK801 promote net excita2on in rat mPFC, even if many individual cells also respond with inhibi2on. Inves2ga2ons of synchronized neuronal ac2vity, in the form of local field poten2als, have more consistently found overlapping effects for both classic and dissocia2ve psychedelics. Ketamine has been found to induce aberrant high-frequency oscilla2onsin several cor2colimbic structures in rodents. Similarly, 5-HT2A agonists induce HFOs in the ventral striatum and mPFC of rodents. HFOs occur to a lesser degree in normal condi2ons in some of these structures and are believed to enable the integra2on of otherwise isolated neuronal informa2on through synchronous ac2vity. The olfactory bulb has been iden2fied as an important source of psychedelic-induced HFOs. However, local infusion experiments suggest that HFO genera2on can be ini2ated independently in mul2ple structures and then spread to other regions. Here, we inves2gated spike ac2vity and HFO oscilla2ons in recordings from several brain structures in parallel in freely behaving rats to iden2fy both local and system-wide neuronal ac2vity changes specific to the psychedelic state. We compared classic psychedelics (LSD, DOI) and dissocia2ve psychedelics (ketamine, PCP) to a non-psychedelic psychoac2ve control (amphetamine) and found that aberrant HFOs were consistently and specifically present during the psychedelic state in ventral striatum and in medial prefrontal, orbital and olfactory cor2ces. Moreover, we found that HFOs were phase synchronized between brain structures in a way that could severely influence the integra2on of informa2on across neuronal systems. Hence, we propose that such hypersynchrony is a mechanism by which an altered state of consciousness can be induced, either during a psychedelic experience or as part of a psycho2c episode.
ANIMALS
Nine Sprague-Dawley rats were used (Taconic, Denmark) in this study. All procedures were approved in advance by the Malmö/Lund ethical commibee of animal experiments.
CONSTRUC+ON AND IMPLANTA+ON OF ELECTRODE ARRAYS
Microelectrode arrays with 128 wires were built and implanted as previously described. Three screws on the occipital bone (over the cerebellum) were connected to a silver wire and used as ground for the recording system. The animals were allowed to recover for at least two weeks ater implanta2on. Electrode posi2ons were verified post mortem in 5 animals using computed tomography (see Supplementary methods). A total of 50 unique anatomical labels were abributed to the electrode 2ps in this way and they were further grouped into 10 broader regions based on assumed func2onal similarity (see Tableand).
PHARMACOLOGICAL TREATMENTS
To record the behavioral and electrophysiological effects of pharmacological treatments, animals were placed in a round open field arena. Ater ~60 min of baseline recording, the animal was intraperitoneally injected with LSD (lysergic acid diethylamid, 0.3 mg/kg, Lipomed AG, Switzerland), DOI (2,5-dimethoxy-4iodoamphetamine hydrochloride, 2 mg/kg, Lipomed AB, Switzerland), ketamine (Ketaminol, 25 -50 mg/ kg, Intervet AB, Sweden), PCP (phencyclidine hydrochloride, 5 mg/kg, Lipomed AG, Switzerland) or amphetamine (d-amphetamine hydrochloride, 4 mg/kg, Tocris, UK) and recorded for another 60-120 minutes. Data was averaged over -35 to -5 minutes for baseline measurements and 30 to 60 minutes for of 5 28 (which was not certified by peer review) is the author/funder. All rights reserved. No reuse allowed without permission. on-drug measurements (rela2ve to drug injec2on). The experiment was repeated ater at least 48 hours of rest (see Figure).
SIGNAL ACQUISI+ON
Local field poten2al (LFP) and single unity ac2vity were recorded with the Neuralynx mul2channel recording system using a unity gain preamplifier (Neuralynx, MT, USA) or with the OpenEphys acquisi2on systemusing 4 Intan RHD2132 amplifier boards with on-board AD-converters and accelerometers (Intan technologies, CA, USA). See Supplementary methods for more details. Synchronized video was recorded at 25 fps with a camera placed above the arena.
BEHAVIORAL SCORING
Behavior was scored offline from the videos for 1 minute every 10 minutes. Behaviors were scored from 0 -3 depending on their prevalence (0 = not present, 1 = present for more than 5 seconds, 2 = present for more than 30 seconds, 3 = present con2nuously). See Tableand Supplementary methods for detailed defini2ons of the scored behaviors.
VIDEO TRACKING
The centroid posi2on of the animal was tracked using algorithms adapted from, and locomo2on speed and distance traveled was then calculated. A center area was defined as a circle with 2/3 the radius of the arena and the frac2on of 2me spent in the center was quan2fied.
HEAD-TWITCH RESPONSE
Head-twitch responses (HTR) were detected using the on-board accelerometers on the Intan RHD2132 amplifier boards together with custom-made Matlab code available at hbps://github.com/NRC-Lund/ of 6 28 htrdetector. See Figureand Supplementary methods for more details. Neuralynx recordings lacked accelerometer data and were analyzed manually from the videos.
SPIKE SOR+NG
Extracted spikes were clustered according to a hierarchical clustering scheme using a 2 ms refractory period. See Supplementary methods for more details. Units were classified into puta2ve cell types based on waveform features as previously described. A neuron remained unclassified if the probability of belonging to a cell type was <0.75 for both types. Figureand Tableshow the clustering and summarize the waveform features.
LFP POWER SPECTRAL DENSI+ES
To emphasize local sources of the measured electrical poten2al, bipolar LFP 2me series were computed from all unique pairs of electrodes from the same structure. Power spectra for truly rhythmic ac2vity was separated from arrhythmic "fractal" ac2vity using the IRASA method. Leastsquare fi{ng of each spectrum to the func2on allowed us to define peak height (a1) and peak frequency (a2) parametrically and to define criteria for when a peak was detected (see Figureand Supplementary methods for more details).
INSTANTANEOUS PHASE AND AMPLITUDE
To quan2fy the instantaneous phase and amplitude of HFOs, the analy2cal signal was calculated (Matlab hilbert func2on) from monopolar LFP 2me series ater bandpass filtering ±5 Hz around the median HFO frequency of each recording. Amplitude auto-and cross-correla2ons were calculated from the instantaneous amplitude 2me series. Mean phase differences and resultant vector lengths were of calculated using the CircStat toolbox. A wire pair was said to be phase inverted if and if , where is the average phase difference and is the concentra2on parameter of the von Mises distribu2on.
GRANGER CAUSALITY
Granger causality was calculated from bipolar LFP 2me series using the one_bi_ga.m func2on of the BSMART toolboxwith a 500 ms window length and a model order of 5. Spectral peaks were detected using the Matlab findpeaks.m func2on with default se{ngs. Only the highest peak was analyzed further if mul2ple peaks were detected. The total Granger causality between regions was es2mated as the median of the amplitude of the Granger causality peaks from all relevant wire pairs.
STA+S+CAL ANALYSIS
Nested ANOVAwas used unless otherwise stated. Please see Supplementary Methods for more details.
BEHAVIOR
There are known behavioral cues in rodents that predict psychedelic effects in humans. However, those cues are different for classic and dissocia2ve psychedelics; 5-HT2A agonists mainly induce head-twitch responses, while NMDA antagonists mainly induce hyperlocomo2on and ataxia. Therefore, the first step in our study was to compare the effects of classic and dissocia2ve psychedelics to the non-psychedelic control on a broad set of behaviors and look for psychedelic-specific behavioral changes. As expected, the NMDA antagonists induced hyperlocomo2on and ataxia, while 5-HT2A agonists induced head-twitch responses and amphetamine induced hyperlocomo2on and stereotypy (p<0.001; Figure, 1C, S1 and S2). We found no evidence for a behavior that was similarly and specifically altered by both 5-HT2A agonists and NMDA antagonists. Importantly, this means that psychedelic-specific changes in neuronal ac2vity cannot be explained trivially by changes in motor behavior.
UNIT AC+VITY
In total, 365 units were iden2fied from 169 recording sites in 9 freely behaving animals (Figure). When comparing the ac2vity of all neurons from all structures to baseline, we observed that the spontaneous firing rates of a large propor2on of cells on this global scale were affected by the drugs: 65% were inhibited and only 8% were excited with 5-HT2A agonists, while a more balanced response was found for the NMDA antagonists (40% were inhibited and 40% were excited at the p < 0.01 level; Figuresand). As a comparison, in response to amphetamine, 32% were inhibited and 46% were excited (Figure). We further inves2gated the specific effects of each drug on different cell popula2ons. Ater grouping of the recording sites based on structural and func2onal similarity (Table), 7 groups had enough cells to warrant further analysis of popula2on modula2on (≥4 cells from both 5-HT2A and NMDA experiments; Table); olfactory cortex, orbitofrontal cortex, medial prefrontal cortex, temporal associa2on area, sensorimotor cortex, ventral striatum and integra2ve thalamus. Cell classifica2on based on waveform features was performed on all these structures except thalamus, allowing for a puta2ve division into principal cells (PC), interneurons (IN) and an intermediate group of cells that could not reliably be assigned to either group (uniden2fied cells; X) (see Figure). For the 5-HT2A agonists the domina2ng effect on firing rates was inhibi2on in each of these structures, both in terms of net standardized rates and number of significantly up/down modulated cells (see Figure-G, let panels). For the NMDA antagonists the effect on firing rates was more mixed, with many popula2ons responding with both inhibi2on and excita2on (see Figure, right panel). However, in terms of net standardized rates, we observed that INs were excited and PCs were inhibited in most structures (see Figure, right panel). To look for rate modula2ons specific to the psychedelic state, we grouped all types of psychedelic drugs and compared the firing rates to the non-psychedelic control compound amphetamine. Only neurons in integra2ve thalamus showed a significant difference in this comparison, with a modest reduc2on of -0.3 (z-scored) in the psychedelic state compared to an increase of 1.6 in the amphetamine state (p < 0.05). Further, the rate modula2ons for 5-HT2A and NMDA were not significantly different in integra2ve thalamus (we also controlled that the baselines for all drugs were sta2s2cally similar; p > 0.05). However, on a systems level the NMDA modula2ons throughout all the recorded structures resembled the amphetamine modula2ons more than the 5-HT2A modula2ons. This is illustrated by Figurewhere IN modula2on is plobed against PC modula2on for each structure. NMDA and amphetamine had a similar pabern of simultaneous IN excita2on and PC inhibi2on, while 5-HT2A had simultaneous inhibi2on in both IN and PC. The rela2ve dissimilarity of the clusters was quan2fied by measuring their Bhabacharyya distance): 5HT2A-NMDA was the most dissimilar pair with a distance of 6.6, while NMDA-amphetamine was the most similar pair with a distance of 0.43. The distance for the 5HT2A-amphetamine pair was 0.98. In summary, 5-HT2A agonists caused widely distributed inhibi2on in both principal cells and interneurons. NMDA antagonists caused similar inhibitory effects in principal cells, but opposite, excitatory effects in interneurons. A psychedelic-specific modula2on was found in integra2ve thalamus, The scale bar indicates 10 ms. The signal was low pass filtered at 500 Hz. B. LFP spectrograms from the ventral striatum during administra2on of ketamine (top), LSD (middle) and amphetamine (bobom). A clear increase in high-frequency oscilla2ons around 150 Hz (HFOs) is evident ater injec2on of ketamine or LSD, but not ater amphetamine. The color scale is in units of dBfractal, i.e. decibels normalized to the fractal noise background. The transla2onal movement speed is shown below each spectrogram. C. LFP power spectra averaged over 2me and treatment groups (baseline=grey, 5-HT2A agonists=pink, NMDA antagonist=red). The 2me periods used were -35 to -5 minutes for baseline and 30 to 60 minutes for drug treatment rela2ve to injec2on (also indicated in panel B as horizontal bars above each spectrogram). D. Distribu2on of HFO detec2on rates in different structures, showing that high detec2on rates are much more common in 5-HT2A and NMDA condi2ons. Each value is the average detec2on rate for a condi2on in a structure during one recording session. E. Summary of HFO detec2on rates in different structures, showing a similar pabern of HFO prevalence for 5-HT2A agonists and NMDA antagonists. HFOs were classified as persistent (red; more than 90% detec2ons in more than 33% of sessions), prevalent (orange; more than 50% detec2ons in more than 33% of sessions), occasional (yellow, if not in any other class) or absent (grey; more than 5% detec2ons in less than 5% of sessions). but on a systems level the effects of the NMDA antagonists resembled the effects of amphetamine more than those of the 5-HT2A agonists.
LFP OSCILLA+ONS
To extend on previous results on psychedelics-induced HFOs, we recorded simultaneously from several cor2cal and subcor2cal regions and characterized the extent of HFOs, as well as the interac2ons between HFOs in different regions.
EXTENT AND PREVALENCE OF HFOS
As expected, both 5-HT2A agonists and NMDA antagonists consistently caused strong increases in HFOs around 150 Hz, as exemplified in Figureand 3B. HFOs increased both in terms of their amplitude (Figure) and their prevalence (as measured by the detec2on rate; Figureand Figure). In contrast, amphetamine increased broad gamma ac2vity (30-80 Hz) without causing HFOs (Fig., bobom, Figureand Figure). In data pooled from all structures and individuals there was a clear shit towards higher HFO detec2on rates in 5-HT2A and NMDA condi2ons, but there was also a large spread in the detec2on rate distribu2ons (Figure), indica2ng varia2on across structures. To further illustrate the extent and prevalence of HFOs in different structures we defined 4 prevalence classes based on the distribu2on of detec2on rates (see Figure); for each structure and condi2on, HFOs were classified as persistent (>90% detec2ons in >33% of sessions), prevalent (>50% detec2ons in >33% of sessions), occasional (if not in any other class) or absent (>5% detec2ons in <5% of sessions). According to this classifica2on, 5-HT2A agonists and NMDA antagonists increased HFO prevalence in an iden2cal pabern, while the pabern for the non-psychedelic compound amphetamine was very similar to the baseline pabern (Figure).
HFO DYNAMICS AND INTERACGONS BETWEEN STRUCTURES
Next, we addressed the ques2on if psychedelics-induced HFOs are generated in a single structure and propagate to other regions, or whether they are generated in a more distributed fashion. We did this by examining how HFOs in different structures were related in terms of their amplitude, frequency and phase. Qualita2vely, the HFOs had a spindle-like amplitude modula2on and strong co-modula2on within structures, while co-modula2on between structures seemed to be weaker and more complex (Figure). The autocorrelograms of the instantaneous amplitudes were consistent with a typical recorded synchronously from the olfactory cortex (OC; blue), the ventral striatum (vStr; red), the medial prefrontal cortex (mPFC; yellow) and the orbitofrontal cortex (OFC; purple) during LSD treatment. Strong co-modula2on of the HFO amplitude can be seen both within and between structures, as well as examples of independent modula2on within and between structures. B. Autocorrelograms of the instantaneous amplitude of all channels with clear HFOs during the psychedelic state (5-HT2A or NMDA). The black line is the average. Most autocorrelograms had a single clear peak (FWHM=50±18 ms) consistent with spindle-like amplitude modula2on with a spindle length of about 50 ms. C. Cross-correla2ons of the instantaneous amplitude between all channel pairs with clear HFOs during the psychedelic state (5-HT2A or NMDA). Pairs with both channels in the same anatomical structure ("Within") had higher cross-correla2ons on average than pairs with channels in different structures ("Between"). However, some pairs showed between-structure co-modula2ons that were similar in strength to the within-structure values. D. 2D histogram showing the rela2onship between HFO frequencies in pairs of structures. Each data point comes from an 8 s 2me window. The high count on the diagonal shows that the frequency is very similar in all structures at any given 2me, despite a high degree of frequency modula2on. E. Example traces showing a single HFO spindle recorded synchronously from 7 electrodes in the olfactory bulb (OB; top), the ventral striatum (VS; middle) and the orbitofrontal cortex (OFC; bobom) during treatment with LSD. Electrode posi2ons are shown to the let. Ver2cal dashed lines are aligned to the peaks of the top OB trace (black) to facilitate comparisons of peak 2mes between electrodes. Polar plots to the right show histograms of the phase difference of each electrode rela2ve to the black OB electrode (based on the whole drug treatment period 30-60 minutes ater injec2on). F. Scaber plot showing mean phase differences and 𝜅 values for the phase difference distribu2ons of each electrode pair (n=6237). Most pairs had a non-random phase rela2onship (86% with 𝜅>1). Of those, 95% had a mean phase difference close to 0 ( ; blue dots) and 2% had an inverted phase ( ; red dots). G. Swarm plot of phase differences rela2ve to the olfactory bulb for all electrode pairs (n=1686) grouped on structure. Each black dot is one electrode pair and red crosses indicate the median for the structure. A posi2ve value means that the structure leads the olfactory bulb. Asterisks indicate that medians are significantly different from zero at the p<0.05 (*), p<0.01 (**) and p<0.001 (***) levels (Wilcoxon signed rank). OB=olfactory bulb, OC=olfactory cortex, vStr=ventral striatum, dStr=dorsal striatum, OFC=orbitofrontal cortex, mPFC=medial prefrontal cortex. H. Examples of phase inversion in LFP traces from two nearby electrodes in the anterior olfactory nucleus (AON) and two nearby electrodes in the orbitofrontal cortex (OFC). This indicates the presence of local current dipoles between each electrode pair. Both raw LFP traces ("Wide") and bandpass filtered traces ("Narrow") are shown. Note that the AON pair is not recorded synchronously with the OFC pair in this example. J. Example Granger causality spectra for one bipolar measurement in medial prefrontal cortex (mPFC) and one bipolar measurement in ventral striatum (vStr) during an NMDA antagonist experiment (top) and a 5-HT2A experiment (bobom). The blue spectra show the causality of mPFC on vStr, while the red spectra show the causality in the opposite direc2on. Solid lines show the drug treated period and dashed lines show the corresponding baselines. The dobed ver2cal lines indicate the HFO frequency in the corresponding recording. K. Median Granger causality values (calculated from the spectrum peak in the HFO band) for all co-recorded structures with clear HFOs. Black squares indicate missing data. Asterisks indicate that medians are significantly different from zero at the p<0.05 (*), p<0.01 (**) and p<0.001 (***) levels (Wilcoxon signed rank). Only vStr→mPFC showed a causality significantly above zero. spindle length of about 50 ms (the autocorrela2on func2on had a smooth peak with FWHM=50±18 ms; Figure). In some electrodes, the autocorrelograms revealed periodic amplitude modula2on at the theta frequency (also seen in Figure). Cross-correlograms of the instantaneous amplitudes confirmed that the amplitude modula2on was generally more strongly correlated within structures than between structures (median correla2on coefficient 0.82±0.10 and 0.40±0.15, respec2vely; p<0.001, Wilcoxon rank sum; Figure). However, some pairs of structures showed co-modula2ons that were similar in strength to the within-structure values (see the group above 0.7 in Figure). These strong amplitude correla2ons were found between olfactory cortex, ventral striatum and orbitofrontal cortex, and between medial prefrontal cortex and ventral striatum. The oscilla2on frequency was remarkably similar in different brain structures when comparing simultaneous oscilla2ons at different sites, despite large varia2ons between individuals and between condi2ons (see Figure). Such frequency co-modula2on could occur if a single dominant source of oscilla2ng firing rates entrains all other structures via synap2c transmission. This may be considered a likely scenario given previous data showing the importance of the olfactory bulb in the genera2on of HFOs. However, further analyses of phase rela2onships did not support this hypothesis: Figureshows a single HFO spindle recorded simultaneously from 7 electrodes in the olfactory bulb, the ventral striatum and the orbitofrontal cortex. In this example, the HFOs in the olfactory bulb had almost zero phase difference, while the HFOs in the ventral striatum led the olfactory bulb by about 0.5 radians. Orbitofrontal HFOs had similar phases, ranging between 0 and 0.5 radians in this example. These small -but oten non-zero -phase differences were a general finding: In electrode pairs with detectable HFOs (median amplitude >5 μV, n=6237), most pairs had a non-random phase difference (86% with kappa>1), and 95 % of those pairs had an absolute phase difference smaller than pi/4 (see the blue group in Figure). When we compared the phase of the olfactory bulb to different structures, we saw small devia2ons from zero in all structures (ranging from 0.001 to 0.45 radians, corresponding to temporal delays of <1 ms; see Figure). of The observed near-zero phase lags are not consistent with a single source propaga2ng synap2cally or via volume conduc2on. An alterna2ve explana2on is that the HFOs are generated by a system of several self-sustaining but weakly interac2ng oscillators located in mul2ple structures. Such systems are known to have stable states with near-zero phase lags; similar to a standing wave. To find support for the presence of local HFO generators outside the olfactory bulb, we looked for phase inversions in measurements from adjacent electrodes, since an inverted phase indicates that a local current dipole is present between the electrodes. Such inversions were indeed found in about 2% of intrastructural electrode pairs (defined as an absolute phase difference larger than 3𝜋/4 and ; see the red group in Figure) and in several structures, including the olfactory bulb, the olfactory cortex, the orbitofrontal cortex, the medial prefrontal cortex and the ventral striatum (Figure). To further map the network of influences between structures, we calculated the Granger causality between structure pairs in the frequency domain. The Granger causality spectra oten had clear peaks either in the HFO band or in the gamma band. During baseline, the peak of the Granger causality was mostly in the low gamma band, while it was mostly in the HFO band during the psychedelic state (Figure). The mean Granger causality in the HFO band (HFO frequency ±10 Hz) increased by 97% compared to baseline (p<0.001, Wilcoxon rank sum), while it was not significantly changed in the gamma band. A structure-by-structure analysis revealed that the Granger causality was clearly strongest from ventral striatum to mPFC (Figure). In summary, the most specific neurophysiological correlate of psychedelic drug ac2on was the enhancement of widespread phase-synchronized HFOs around 150 Hz. In the psychedelic state, we have iden2fied a network consis2ng of the olfactory bulb, the olfactory cortex, the ventral striatum, the orbitofrontal cortex and the medial prefrontal cortex that are 2ghtly coupled in the HFO band. Vollenweider & Preller, 2020). Some models explicitly state a direct link between firing rates and func2onal changes, like the thalamocor2cal ga2ng model, in which increased mPFC excitability disinhibits the thalamus and reduces its ability to gate sensory informa2on. Other models have indirectly linked the ac2va2on of certain cell popula2ons to the disintegra2on of canonical network states, as seen for example in the reduced BOLD correla2on between nodes in the default mode network. Exis2ng models have however highlighted apparently contradictory evidence and the field has suffered from an almost complete lack of in vivo electrophysiology data from awake animals that could link the pharmacological effects on single cells with changes in informa2on processing in the brain as a whole. With simultaneous large-scale microwire recordings in mul2ple structures in freely behaving animals, we here show that different classes of psychedelics affect firing rates differently, while they cause similar changes in popula2on dynamics in the form of aberrantly strong HFOs. We also show that the HFOs facilitate func2onal coupling between brain structures, both in terms of phase synchroniza2on and Granger causality. This could have major effects on the exchange and integra2on of informa2on across these neuronal systems. Our results cast doubts on models sugges2ng that specific changes in firing rates are directly linked to the psychedelic state, since the appearance of HFOs is largely independent of popula2on firing rates. This is a surprising finding given the dominant hypothesis that both classic and dissocia2ve psychedelics exert their effect via increased release of glutamate in cor2co-limbic circuits (de Gregorio et of. Intui2vely, increased glutamate release should lead to increased excita2on. However, biological neuronal networks have several homeosta2c mechanisms to regulate the balance between excita2on and inhibi2on, which makes it difficult to predict how changes in glutamate signaling affect the overall behavior of the network. Indeed, the seminal studies repor2ng increased AMPA-dependent excitatory currents in pyramidal cells did not observe simultaneous increases in firing rates of those neuronsand, more generally, it is well known that psychedelics are not proconvulsant. Perhaps it is frui ‡ul to focus less on the excitatory role of glutamate, and more on how it may change the strength and temporal dynamics of effec2ve synap2c coupling. Theore2cal work has shown that the strength and dynamics of synap2c coupling play crucial roles in determining network behavior and, in par2cular, if a network has stable periodic states. Intriguingly, both classic and dissocia2ve psychedelics decrease NMDA-mediated synap2c currents, while they increase AMPA-mediated currents (as men2oned above). However, AMPA channels have a much shorter deac2va2on 2me constant than NMDA channels (AMPA: 2-5 ms, NMDA: 50-100 ms), which will lead to drama2cally faster temporal dynamics of the excitatory postsynap2c currents. This decreased "temporal smoothing" could in turn enable stable oscillatory states with a shorter period than would otherwise be possible. While this might be a possible explana2on for the appearance of local oscillatory states, the long-range HFO synchroniza2on observed in the current study is s2ll perplexing. Intui2vely, it seems unlikely that such fast oscilla2ons can synchronize across long distances considering the sizeable delays caused by the propaga2on of ac2on poten2als and the delayed ac2va2on of chemical synapses. On the other hand, gap junc2ons and ephap2c coupling could influence neighboring neurons almost instantaneously, but have very short range. However, mathema2cal analysis of idealized coupled oscillators has shown that stable synchronous states can exist with only local connec2vity and even with delayed influences. Interes2ngly, such systems oten display a surprising complexity, where mul2ple stable synchronous states can co-exist and have different synchroniza2on frequencies. Taken together, the current study represents a first step towards bridging psychedelics-induced physiological phenomena at the single cell level, via networks, to global brain states. Increasing our mechanis2c understanding of how this class of substances induce psychedelic states will be essen2al to develop improved therapies for several neuropsychiatric condi2ons.