LSDPsilocybinPlacebo

Serotonergic psychedelic drugs LSD and psilocybin reduce the hierarchical differentiation of unimodal and transmodal cortex

This study (2022) applied a non-linear dimensionality reduction technique previously used to map hierarchical connectivity gradients to assess cortical organization in the LSD and psilocybin state from two previously published pharmacological resting-state fMRI datasets. The principal gradient of cortical connectivity was significantly flattened under both drugs relative to their respective placebo conditions. This study is the first to show that macroscale connectivity gradients are sensitive to acute pharmacological manipulation.

Authors

  • Bernhardt, B.
  • Carhart-Harris, R. L.
  • Girn, M.

Published

NeuroImage
individual Study

Abstract

Lysergic acid diethylamide (LSD) and psilocybin are serotonergic psychedelic compounds with potential in the treatment of mental health disorders. Past neuroimaging investigations have revealed that both compounds can elicit significant changes to whole-brain functional organization and dynamics. A recent proposal linked past findings into a unified model and hypothesized reduced whole-brain hierarchical organization as a key mechanism underlying the psychedelic state, but this has yet to be directly tested. We applied a non-linear dimensionality reduction technique previously used to map hierarchical connectivity gradients to assess cortical organization in the LSD and psilocybin state from two previously published pharmacological resting-state fMRI datasets (N = 15 and 9, respectively). Results supported our primary hypothesis: The principal gradient of cortical connectivity, describing a hierarchy from unimodal to transmodal cortex, was significantly flattened under both drugs relative to their respective placebo conditions. Between-condition contrasts revealed that this was driven by a reduction of functional differentiation at both hierarchical extremes - default and frontoparietal networks at the upper end, and somatomotor at the lower. Gradient-based connectivity mapping indicated that this was underpinned by a disruption of modular unimodal connectivity and increased unimodal-transmodal crosstalk. Results involving the second and third gradient, which respectively represent axes of sensory and executive differentiation, also showed significant alterations across both drugs. These findings provide support for a recent mechanistic model of the psychedelic state relevant to therapeutic applications of psychedelics. More fundamentally, we provide the first evidence that macroscale connectivity gradients are sensitive to an acute pharmacological manipulation, supporting a role for psychedelics as scientific tools to perturb cortical functional organization.

Unlocked with Blossom Pro

Research Summary of 'Serotonergic psychedelic drugs LSD and psilocybin reduce the hierarchical differentiation of unimodal and transmodal cortex'

Introduction

Research on serotonergic psychedelics such as LSD and psilocybin has reported substantial changes in whole-brain functional organisation and dynamics, including reduced segregation of large-scale networks, increased global functional connectivity, and greater complexity of regional activity. A recent theoretical account—the RElaxed Beliefs Under Psychedelics (REBUS) model—posits that psychedelics increase the influence of low-level sensory and limbic inputs on high-level representations encoded in transmodal cortex, implying a reduction in hierarchical differentiation between unimodal and transmodal regions. However, whether psychedelics acutely alter macroscale hierarchical organisation has not been directly tested. Girn and colleagues set out to test this idea by applying gradient-mapping techniques to pharmacological resting-state fMRI datasets acquired under LSD and psilocybin. Their primary hypothesis was that both drugs would produce a contraction of the principal cortical connectivity gradient—an axis running from unimodal sensorimotor regions to transmodal association cortex—reflecting reduced hierarchical differentiation and increased unimodal–transmodal crosstalk. Secondary aims included examining drug effects on the second and third gradients, which capture sensory and executive differentiation respectively, and relating gradient changes to self-reported subjective effects such as ego-dissolution and complex imagery. The study therefore investigates whether macroscale connectivity gradients are sensitive to acute serotonergic pharmacological manipulation and whether such changes relate to phenomenology.

Methods

The study analysed previously published within-subject resting-state BOLD fMRI datasets for LSD and psilocybin. For LSD, 20 subjects were scanned on two occasions (LSD and placebo) separated by two weeks, with 75 μg LSD or saline administered intravenously; two non-music resting scans (scans 1 and 3) were used. For psilocybin, 15 subjects underwent two scanning sessions (psilocybin and placebo) with a 2 mL psilocybin infusion in 10 mL saline or saline alone; a single 12-minute scan was used and the post-infusion half of each session was analysed. Subjects were blind to session order; experimenters and analysts were not. Ego-dissolution ratings were collected using intra-scanner visual analogue scales and complex imagery was assessed from the 11-factor altered states of consciousness (ASC) questionnaire completed after each day. Preprocessing and denoising followed an identical pipeline for both datasets. Excessive motion (>15% volumes with FD ≥ 0.5) led to exclusions: four participants in the LSD dataset and six in the psilocybin dataset were removed, and one additional LSD participant left the scanner due to anxiety, yielding final samples of 15 (LSD) and 9 (psilocybin). Standard steps included removal of initial volumes, de-spiking, slice-timing and motion correction, brain extraction, rigid and non-linear registration to MNI space, scrubbing (FD = 0.4), 6 mm spatial smoothing, band-pass filtering (0.01–0.08 Hz), linear and quadratic detrending, regression of six motion parameters and three anatomical nuisance signals; global signal regression was not performed. Quality control indicated no distance-dependent motion confounds after denoising. Gradient-mapping analyses were conducted using the BrainSpace toolbox. Surfaces were downsampled to 10,000 vertices, and for each subject a 10,000 × 10,000 Pearson correlation matrix was computed, z-transformed and thresholded row-wise at 90% sparsity to retain strongest connections. Cosine similarity was computed on the thresholded matrix and diffusion map embedding (α = 0.5) was applied to obtain low-dimensional embedding components (gradients). Individual gradients were aligned to a group template using Procrustes rotation. Group contrasts and behavioural associations were evaluated with surface-based linear models implemented in SurfStat. To probe hierarchical specificity, vertex-wise gradient scores were binned into ten percentile ROIs per subject; combined unimodal (lowest three bins) and transmodal (highest three bins) ROIs were created and bin-wise functional connectivity (Pearson r) with these combined ROIs was compared between drug and placebo. Seedmaps for each bin were also computed and compared at vertex-wise and network-wise levels.

Results

Diffusion map embedding revealed multiple orthogonal gradients per subject; the report focuses on the first three gradients that explained the largest variance. The principal gradient replicated prior work by spanning unimodal somatomotor regions at one extreme to transmodal default-network and superior frontal regions at the other. Variance explained by this gradient was: LSD placebo mean 12.9%, LSD mean 12.6%; psilocybin placebo mean 13.1%, psilocybin mean 11.9%. Quantitatively, the principal gradient contracted under both drugs: LSD showed a significant reduction in the difference between maximum and minimum principal gradient values (t28 = -2.5, p = 0.01) and psilocybin showed a larger contraction (t16 = -3.9, p = 0.001). Vertex-wise and network-wise contrasts indicated that both unimodal-proximal and transmodal-proximal regions became less differentiated along this axis. Specific topographical changes varied somewhat by drug but displayed broad convergence. Under LSD, vertex-wise increases in principal gradient scores appeared in somatomotor cortex and medial/lateral visual regions, while decreases were observed in precuneus, premotor cortex, superior and inferior frontal gyri, superior parietal lobule, and Wernicke's area (FWE p < 0.05, cdt < 0.01). Psilocybin produced vertex-wise increases in somatomotor, auditory, and lateral visual cortex, and decreases in superior frontal gyrus and inferior parietal lobule (FWE p < 0.05, cdt < 0.01). Network-wise, LSD increased gradient scores in visual and somatomotor networks and decreased scores in default, frontoparietal control, and dorsal attention networks (p < 0.05 FDR-corrected). Psilocybin showed a network-wise increase in the somatomotor network and decrease in the default network (p < 0.05 FDR-corrected). Control analyses found no significant correlation between principal gradient differences and motion. Gradient-based connectivity mapping supported the hypothesis of increased unimodal–transmodal crosstalk. Lower percentile bins (unimodal side) exhibited significantly reduced connectivity with the combined unimodal ROI and less negative connectivity with the combined transmodal ROI, while higher percentile bins (transmodal side) showed increased connectivity with the unimodal ROI but not with the transmodal ROI (p < 0.05, FDR-corrected or uncorrected). Seedmap analyses indicated that lower bins decreased FC with somatosensory and visual networks and increased FC with frontoparietal control regions, notably lateral prefrontal and lateral parietal cortex (FWE p < 0.05, cdt < 0.05 vertex-wise; network-wise p < 0.05 FDR or uncorrected). Higher bins tended toward increased FC with visual, dorsal attention, and salience networks. The second gradient separated visual cortex from somatomotor/auditory cortex and explained: LSD placebo mean 8.8%, LSD mean 8.6%; psilocybin placebo mean 9.8%, psilocybin mean 10.3%. The third gradient separated executive/task-positive regions from other cortex and explained: LSD placebo mean 7.7%, LSD mean 6.0%; psilocybin placebo mean 8.4%, psilocybin mean 7.5%. LSD produced vertex-wise increases in lateral/ventral visual cortex and middle/superior frontal gyri and decreases in somatomotor/auditory cortex for the second gradient (FWE p < 0.05, cdt < 0.01), with network-wise increases in frontoparietal control and visual networks and decreases in the somatomotor network (FDR p < 0.05). Psilocybin showed no significant second-gradient effects, which the authors attribute to limited statistical power. For the third gradient, both drugs showed patterns consistent with a contraction reflecting reduced differentiation of executive regions, though reported effects were stronger for LSD. Finally, behavioural associations identified one significant relationship: in the LSD dataset, greater principal gradient values in a left medial prefrontal/anterior cingulate cluster extending into superior frontal gyrus were positively associated with increases in ego-dissolution (FWE p < 0.05, cdt < 0.01). No significant associations were observed for psilocybin or for complex imagery in either dataset.

Discussion

Girn and colleagues interpret their results as evidence that serotonergic psychedelics acutely attenuate macroscale cortical hierarchy, consistent with predictions derived from the REBUS model. The principal gradient contracted under both LSD and psilocybin, indicating reduced differentiation between unimodal and transmodal cortex. Gradient-based connectivity analyses suggest this contraction reflects structured increases in unimodal–transmodal crosstalk: unimodal regions became less internally cohesive and more connected to transmodal hubs, while transmodal regions preserved intra-transmodal connectivity and increased connections to unimodal regions. The authors propose that this pattern corresponds to fewer intervening processing steps between sensorimotor and association cortex, which could underlie the experiential blurring of internal and external boundaries reported during psychedelic states. The discussion places these findings in the context of prior literature showing altered static and dynamic functional connectivity and increased signal complexity with psychedelics, arguing that the gradient approach offers a low-dimensional organisational framework that can reconcile heterogeneous network-level findings. Differences between LSD and psilocybin were noted for the second gradient—LSD showed reduced visual differentiation not evident with psilocybin—which aligns with prior reports of stronger visual-cortex whole-brain connectivity changes with LSD. The positive association between principal gradient changes in left mPFC/ACC and ego-dissolution under LSD is highlighted as a potentially meaningful link between hierarchical cortical organisation and a clinically-relevant subjective effect, though the authors explicitly caution that replication is needed given the small samples and the absence of the effect in the psilocybin data. Key limitations acknowledged by the investigators include relatively small final sample sizes, which constrain statistical power and generalisability; difficulties in collecting larger psychedelic datasets due to logistical, funding, and ethical hurdles are noted. The authors also emphasise that, while LSD and psilocybin have complex pharmacology, the findings are consistent with past work implicating 5-HT2A receptor agonism in psychedelic effects, suggesting a possible serotonergic modulation of cortical functional differentiation. They conclude that cortical gradients are sensitive to acute pharmacological manipulation and recommend replication and further research to determine whether acute reductions in cortical hierarchy relate to therapeutic outcomes observed in psychotherapeutically mediated psychedelic treatments. The authors present their results as a foundational contribution to understanding how serotonergic psychedelics reorganise macroscale brain organisation and its relation to conscious experience.

View full paper sections

INTRODUCTION

The past decade has seen a resurgence of scientific interest in serotonergic psychedelic compounds such as lysergic acid diethylamide (LSD), psilocybin, and dimethyltryptamine (DMT)/ayahuasca, primarily motivated by suggestive findings from preliminary clinical trials for depression, end-of-life-anxiety, alcoholism, and tobacco addiction. Mirroring their complex subjective effects, functional neuroimaging investigations have shown significant alterations to whole-brain functional organization and dynamics following psychedelic administration. More specifically, findings suggest that psychedelic administration shifts the brain towards greater global functional integration, as reflected by the reduced functional segregation of large-scale brain networks and increased global functional connectivity. In addition, the complexity of brain dynamics has been shown to increase under psychedelics, as indexed by increases in regional and population-level entropy/complexity, as well as increases in the repertoire or dynamic range of functional connectivity (FC) statesNotably, the recently proposed RElaxed Beliefs Under Psychedelics (REBUS) modelunifies past psychological and neural findings with psychedelics into a common theoretical framework based on hierarchical predictive coding and the Free Energy Principle. A primary hypothesis of this model is that psychedelics elicit their characteristic subjective effects by increasing the sensitivity of high-level representations (e.g., beliefs or assumptions) encoded within transmodal cortex to low-level inputs from sensory or limbic sources). An acute reduction in the hierarchical differentiation of transmodal versus unimodal cortex, consistent with increased crosstalk between these typically segregated functional zones, would provide support for this hypothesized effect. However, the impact of psychedelics on neural hierarchical organization has yet to be directly tested. Gradient-mapping techniques have emerged in recent years as valuable tools to characterize brain functional organization. In contrast to parcellation approaches which decompose the brain into discrete regions/networks, these approaches model the brain as the superposition of eigenmodes describing continuous axes of differentiation. Such techniques have consistently revealed a principal component describing a gradient of FC spanning from unimodal sensorimotor regions to transmodal association regions centered on the default network. This gradient is consistent with tracttracing work identifying hierarchical cortical organization in human and non-human primates, corresponds to cortical myeloarchitectonic and transcriptomic transitions, and represents a functional hierarchy from low-level sensorimotor processing to abstract cognition. Gradient-mapping has now been used to characterize neural hierarchy in diverse contexts, including autism, neonatal development, schizophrenia, and lifespan development. The present study investigates alterations of macroscale cortical hierarchy in the psychedelic state by applying gradient-mapping analyses to pharmacological fMRI datasets collected with two serotonergic psychedelic compounds: LSD and psilocybin. We apply cortex-wide and network-specific analyses to characterize psychedelic-dependent changes. In addition, we leverage the organizational scheme provided by the hierarchical gradient to probe hierarchyspecific changes in whole-brain functional connectivity. We predicted that both LSD and psilocybin will be associated with a contraction of the principal gradient, indicating reduced hierarchical differentiation. Further, we predicted that this hierarchical contraction would be specifically consistent with increased crosstalk between unimodal and transmodal cortices. In addition, we investigated psychedelic-induced changes to the second and third gradient, which have been shown, respectively, to differentiate visual from somatomotor/auditory cortex, and executive/task-positive cortex from the rest of cortex. We also predicted gradient dedifferentiation along these axes, given past work which highlighted alterations to sensory and task-positive connectivity in the psychedelic state. Finally, we examined whether psychedelic-dependent changes in cortical gradients are associated with self-reported subjective effects. Collectively, this study provides the first comprehensive assessment of whether serotonergic drugs LSD and psilocybin pharmacologically alter the primary axes of cortical functional organization. In addition, it constitutes the first direct assessment of whether cortical gradients are sensitive to acute pharmacological manipulation.

LSD.

Neuroimaging data from an already published datasetwas used for the present analyses. The participants and the acquisition protocol are described in greater detail in the original publication and in the Supplementary Information. In a within-subjects design, resting-state BOLD fMRI data were acquired in 20 subjects for each of LSD and placebo conditions, with conditions spaced 2 weeks apart. The order was counterbalanced across subjects; subjects were blind to this order, but the researchers and those analyzing the data were not. For the placebo condition, subjects were given 10mL of saline and for the LSD condition they were given 75μg of LSD in 10-mL saline via intravenous injection. Three scans were collected on each scan day, 2 without music and 1 without music. The two nomusic scans were used in the present analyses (scans 1 and 3). Ego-dissolution scores were collected via intra-scanner visual analogue scale ratings, whereas complex imagery scores were derived from the 11-factor altered states of consciousness (ASC) questionnairewhich was completed at the end of each scan day.

PSILOCYBIN.

Neuroimaging data from an already published datasetwas used for the present analyses. The participants and the acquisition protocol are described in greater detail in the original publication and in the Supplementary Information. In a within-subjects design, resting-state BOLD fMRI data were acquired in 15 subjects for each of psilocybin and placebo conditions, with conditions spaced 2 weeks apart. The order was counterbalanced across subjects; subjects were blind to this order, but the researchers and those analyzing the data were not. For the placebo condition, subjects were given 10mL of saline and for the psilocybin condition they were given 2mL of psilocybin in 10mL saline, via intravenous injection. Each scanning day consisted of one 12-minute scan with infusion beginning at 6 minutes following the start of the scan. The post-infusion half of the scan for each condition was used in the present analyses. Ego-dissolution scores were collected via intra-scanner visual analogue scale ratings, whereas complex imagery scores were derived from the 11-factor altered states of consciousness (ASC) questionnairewhich was completed at the end of each scan day.

ETHICS STATEMENT

Data collection for both LSD and psilocybin datasets were approved by the National Research Ethics Service committee London-West London and was conducted in accordance with the revised declaration of, the International Committee on Harmonization Good Clinical Practice guidelines, and National Health Service Research Governance Framework. Imperial College London sponsored the research, which was conducted under a Home Office license for research with schedule 1 drugs.

NEUROIMAGING DATA PREPROCESSING AND DENOISING

Both datasets underwent an identical preprocessing protocol, as described in detail elsewhere. Subjects with excessive head motion were discarded from analyses (>15% volumes with FD >= 0.5). This resulted in 4 exclusions in the LSD dataset, and 6 exclusions in the psilocybin dataset. The excessive head motion was found in scans conducted during the drug conditions. One additional subject in the LSD dataset exited the scanner due to intra-scanner anxiety, leaving a final sample of 15 for the LSD dataset and 9 for the psilocybin dataset. The following pre-processing and denoising steps were performed on the BOLD resting-state fMRI data for each dataset: removal of the first three volumes, de-spiking (3dDespike, AFNI), slice time correction (3dTshift, AFNI), motion correction (3dvolreg, AFNI) by registering each volume to the volume most similar to all others, brain extraction (BET, FSL); rigid body registration to anatomical scans, non-linear registration to a 2mm MNI brain (Symmetric Normalization (SyN), ANTS), scrubbing (FD = 0.4), spatial smoothing (FWHM) of 6mm, bandpass filtering between [0.01 to 0.08] Hz, linear and quadratic de-trending (3dDetrend, AFNI), regression of 6 motion parameters, and regression of 3 anatomical nuisance regressors (ventricles, draining veins, and local white matter). Global signal regression was not performed. Quality control tests confirmed the lack of distance-dependent motion confounds in the denoised data.

STATISTICAL ANALYSIS

Gradient-mapping. Cortical gradients were computed using the BrainSpace toolbox (;) as implemented in MATLAB. Surfaces were first downsampled from fsaverage 5 space (20,484 vertices) to 10,000 vertices for computational efficiency. For each subject, a 10,000 x 10,000 connectivity matrix was then computed by calculating the pairwise Pearson's correlation between all vertices. As has been done previously (e.g.,, this matrix was ztransformed and thresholded row-wise at 90% sparsity in order to retain only the strongest connections. Cosine similarity was then computed on the thresholded z-matrix in order to generate a similarity matrix which captures the similarity in whole-brain connectivity patterns between vertices. This similarity matrix is required as input to the diffusion map embedding algorithm we employed here. The use of cosine similarity as the similarity metric of choice is consistent with past work using this approach. Diffusion map embedding, a non-linear manifold learning technique from the family of graph Laplacians, was applied to similarity matrices in order to identify gradient components at the individual subject level. The technique estimates a low-dimensional set of embedding components (gradients) from a high-dimensional similarity matrix, where each embedding can intuitively be viewed as a dimension of FC pattern similarity covariance. In the embedding space, vertices that are strongly connected (as weighted by FC pattern similarity) by many connections or a few very strong connections are closer together, whereas vertices with little or no connections are farther apart. Euclidean distance between two points in the embedding space is equivalent to the diffusion distance between probability distributions centered at those points (hence the name of the algorithm), each of which are equivalent to 'difference in gradient value' as referred to in the main text. The algorithm is controlled by a single parameter α, which controls the influence of density of sampling points on the manifold (α = 0, maximal influence; α = 1, no influence). Diffusion map embedding is specifically characterized by α = 0.5, which allows the influence of both global and local relationships between data points in the estimation of the embedding space. Following past work, to enable comparisons across subjects, Procrustes rotation was performed to align individual-subject embedding (gradient) components to an all-subjects group average embedding component template. This rotation ensures that gradient axes are matched across subjects. Group contrasts and behavioural associations were conducted using surface-based linear models, as implemented in the SurfStat toolbox (;. Gradient-based connectivity mapping. In order to further explicate how observed gradient differences related to interregional FC, we examined unimodal versus transmodal functional connectivity changes as stratified by gradient scores. More specifically, we wanted to determine whether unimodal regions (as defined based on the hierarchical gradient) preferentially exhibited increased FC with transmodal regions, as opposed to non-specific global increases. To do this, we first separated vertex-wise gradient scores into 10 percentile binsfor each subject. These 10 bins were then used as ROIs in FC analyses. We were interested in whether each bin preferentially increased its connectivity with unimodal versus transmodal cortex after psychedelic administration. A combined unimodal cortex ROI was created as the combination of the first three percentile bins and a combined transmodal cortex ROI was defined as the last three percentile bins. Unimodal-specific FC and transmodal-specific FC was computed for each of the 10 bins, by computing the Pearson's correlation (r) between each bin and the combined unimodal ROI and the combined transmodal ROI. T-tests were applied to evaluate drug vs. placebo differences at each bin for each of unimodal-specific and transmodalspecific FC. To complement this analysis and examine the spatial distribution of FC associated with each gradient percentile bin, we additionally computed seedmaps for each bin and compared across drug and placebo conditions at both vertex-wise and network-wise levels. This enabled a more detailed look at unimodal-specific versus transmodal-specific changes in whole-brain FC.

DATA AND CODE AVAILABILITY

All analyses were conducted in MATLAB using custom code and included functions within the BrainSpace and SurfStat toolboxes described above. Data is freely available at.

RESULTS

We applied gradient-mapping techniques to characterize differences in macroscale cortical gradients in each of LSD and psilocybin states relative to their respective placebo states. The diffusion-map embedding algorithm used here is, compared to other gradient-mapping approaches, robust to noise, computationally inexpensive, and governed by a single parameter controlling the influence of the sampling density on the manifold (see Methods for details). Applying this algorithm with standard settings revealed 98 mutually orthogonal gradient components per subject in each of the LSD dataset and psilocybin dataset. We presently include discussion of the first three gradients revealed by this approach, given that these explain the greatest variance and have been highlighted in past work.

PRINCIPAL GRADIENT

The principal gradient of cortical connectivity revealed in the present datasets replicates past findingsof a putatively hierarchical axis of FC similarity variance spanning from unimodal regions centered in somatomotor cortex on one end to transmodal regions centered on the default network and superior frontal gyrus on the other (Figure). This gradient explained the greatest amount of variance in each of LSD placebo (mean 12.9% variance explained), LSD (mean 12.6% variance explained), psilocybin placebo (mean 13.1% variance explained), and psilocybin (mean 11.9% variance explained) conditions. his ve support for our hypothesis of reduced hierarchical differentiation in the psychedelic state. To quantitatively assess this, we calculated the difference between each subject's maximum and minimum principal gradient value and compared these differences across conditions . Results confirmed the presence of a significant contraction in both LSD (t 28 = -2.5, p=0.01) and psilocybin (t 16 = -3.9, p=0.001) conditions, reflective of reduced functional differentiation along this hierarchical unimodal-transmodal axis. In order to examine overall trends in psychedelic-dependent changes in hierarchical differentiation, we then visualized the principal gradient for each dataset as a scatter plot, color coded for increases (red) and decreases (blue) in each drug condition relative to placebo (Figure). This revealed that both unimodal-proximal regions and transmodal-proximal regions become less differentiated along this axis, suggesting a contraction on both sides of the gradient. Next, we quantitatively assessed between-condition differences in gradient score values at a vertex-(Figureand 1D Mid-Bottom/Bottom) and network-wise (Figureand 1C Bottom) level. Unthresholded maps at the vertex level are included to display the overall consistency in topology between LSD and psilocybin states, indicating that these distinct serotonergic drugs affect the principal gradient in a similar manner . LSD vertexwise contrasts revealed LSD-dependent principal gradient increases in somatomotor cortex, as well as in both medial and lateral visual regions. LSD-dependent principal gradient decreases were found in the precuneus, premotor cortex, superior and inferior frontal gyrus, superior parietal lobule, and Wernicke's area ( FWE p<0.05, cdt<0.01). Psilocybin vertex-wise contrasts revealed psilocybin-dependent principal gradient increases in somatomotor and auditory cortex, as well in visual cortex predominantly on the lateral surface. Psilocybin-dependent gradient decreases were found in the superior frontal gyrus and inferior parietal lobule ( FWE p<0.05, cdt<0.01). Network-wise differences were assessed according to thenetwork parcellation scheme. LSD results indicated a significant increase in gradient score within visual and somatomotor networks, and a significant decrease within default, frontoparietal control, and dorsal attention networks (p<0.05 FDR corrected). Psilocybin results indicated a significant increase in gradient score within the somatomotor network and decrease in the default network (p<0.05 FDR corrected). All between-condition differences were consistent with the gradient scores of significant regions and networks approaching zero, indicating reductions in differentiation along this axis. These results therefore offer further quantitative support for the qualitative trend seen in Figureand 1C Mid-Bottom. Namely, that the LSD state is characterized by a pulling-together of both unimodal (visual and somatomotor) and transmodal (default, frontoparietal control, and dorsal attention) networks in gradient space -reflective of relatively symmetrical reduction in differentiation along this hierarchical axis of cortical connectivity. Control analyses indicate that principal gradient score differences are not significantly correlated with motion (Supplementary Figure).

GRADIENT-BASED CONNECTIVITY MAPPING

In order to further explicate the observed reduction in hierarchical organization in the psychedelic state, we sought to determine whether the observed changes in gradient values were specifically consistent with increased cross-talk between unimodal and transmodal cortices, as hypothesized by the abovementioned REBUS model. This interest was further motivated by the fact that the principal gradient changes found in both datasets each represent a movement towards zero. As such, the results are consistent with the changes simply corresponding to less of a loading on this gradient in general and not necessarily to increases/decreases in unimodal-or transmodal-specific connectivity. To evaluate whether the gradient changes are indicative of structured and hierarchically specific changes in FC, we applied a gradient-based connectivity mapping approach. We constructed ROIs based on percentile bins along the principal gradient and evaluated drug-induced unimodal versus transmodal changes in FC. In particular, bin ROIs were correlated with each of a combined unimodal cortex ROI (three lowest bins) and combined transmodal cortex ROI (three highest bins; see Methods for more details). Between-condition comparisons of unimodal versus transmodal-specific FC supported our hypotheses of increased crosstalk between unimodal and transmodal cortices. For both LSD and placebo, multiple lower percentile bins (corresponding to the unimodal side of the gradient) showed significantly reduced FC with unimodal cortex and less negative FC with transmodal cortex (p<0.05, FDR-corrected or uncorrected; Figure). In addition, higher percentile bins (corresponding to the transmodal side of the gradient) displayed significantly increased FC with unimodal cortex but not with transmodal cortex. This suggests that the observed gradient changes have their basis in structured, hierarchically specific changes in FC that are consistent with increased unimodal-transmodal crosstalk. To further probe this effect, we also examined the spatial distribution of hierarchically specific changes in FC. In particular, we computed whole-brain seedmaps based on each gradient bin ROI and compared across conditions at both the vertex-wise and network-wise level (Figure). Results were largely consistent across drugs and indicate that lower percentile bins exhibit significantly reduced FC with somatosensory and visual networks and increased connectivity with the frontoparietal control network, particularly in lateral prefrontal and lateral parietal cortex ( FWE p<0.05, cdt<0.05 vertex-wise, p<0.05 FDR or uncorrected network-wise). This effect was largely sustained in the median percentile bin, with greater differences involving the default and limbic networks for LSD. With respect to higher percentile bins, with the exception of the default network for the eighth bin of psilocybin, results indicated more of a trend toward increased FC, particularly involving visual, dorsal attention, and salience networks ( FWE p<0.05, cdt<0.05 vertex-wise, p<0.05 FDR or uncorrected network-wise).

SECOND AND THIRD GRADIENT

Having found support for our hypothesis of a contraction in macroscale functional hierarchy in the LSD state that is consistent with greater unimodal-transmodal crosstalk, we additionally examined the second and third gradient of macroscale functional organization. The second gradient represents an axis of FC similarity variance that separates the visual cortex on one end from the somatosensory/auditory on the other (Figureand). This gradient explained the second most variance in each of LSD placebo (mean 8.8% variance explained), LSD (8.6 variance explained), psilocybin placebo (9.8% variance explained), and psilocybin (mean 10.3% variance explained) conditions. The third gradient represents an axis of FC similarity variance that separates executive control/task-positive regions from the rest of cortex (Figureand). This gradient explained the second most variance in each of LSD placebo (mean 7.7% variance explained), LSD (mean 6% variance explained), psilocybin placebo (mean 8.4% variance explained), and psilocybin (mean 7.5% variance explained) conditions. LSD second gradient vertex-wise contrasts revealed LSD-dependent increases in lateral and ventral visual cortex and the middle and superior frontal gyrus bilaterally, while decreases were found in somatomotor and auditory cortex ( FWE p<0.05, cdt<0.01). Network-wise contrasts revealed significant increases in frontoparietal control and visual networks and decreases in the somatomotor network (FDR p<0.05). Psilocybin vertex-and network-wise contrasts revealed no significant differences, likely due to a lack of power in this dataset. The significant qualitative difference between LSD and the trend-level effects in psilocybin is notable, however, and suggests differential effects on sensory differentiation across these drugs. Visual inspection of).

GRADIENT MANIFOLDS

To visualize the relationship between the three gradients examined in this study and how they differ across drug and placebo conditions, we created gradient manifold scatter plots, color coded for each of the Yeo et al. 2011 7 networks (Figure). Qualitative examination of the LSD manifold plots (Figureand 5B) reveals an overall less diffuse embedding space distribution in the LSD state, marked by notable contractions on the unimodal aspect of the principal gradient and the visual aspect of the second gradient. Qualitative examination of the psilocybin manifold plots (Figureand 5D) similarly shows a significant contraction on the unimodal side of the principal gradient, but less contraction of the visual network along the second gradient relative to LSD. Psilocybin also displays a greater contraction of transmodal (default and frontoparietal control network) nodes in the principal gradient, and, interestingly, a greater contraction of the third gradient primarily involving regions within salience and dorsal attention networks.

RELATIONSHIP TO SUBJECTIVE MEASURES

Finally, to ascertain the subjective relevance of changes in macroscale gradients during the psychedelic state, we assessed relationships between drug-placebo gradient differences scores for each drug and each of the three gradients, and drug-placebo differences in two self-report measures which index core components of the psychedelic experience: ego dissolution and complex imagery. Ego dissolution corresponds to the experience that one's sense of self and what is usually perceived as 'not self' (e.g. the physical environment and other individuals) becomes blurred or dissolved completely. Complex imagery corresponds to closed-eye visionary experiences of scenes, landscapes, and images, and is contrasted from simple imagery such as geometric patterns and shapes. One significant relationship was found: a positive association between the principal gradient and ego dissolution in the LSD dataset ( FWE p<0.05, cdt<0.01; Figure). This finding indicates that greater hierarchical differentiation of a large cluster within of the left medial prefrontal cortex/anterior cingulate cortex extending laterally into the superior frontal gyrus was significantly associated with increases in ego dissolution scores. No significant results were found with psilocybin or complex imagery.

DISCUSSION

To investigate psychedelic-induced changes in cortical functional organization and test the hypothesis that serotonergic psychedelics attenuate brain hierarchical organization, we characterized macroscale cortical gradients after LSD and psilocybin administration. Our results for both LSD and psilocybin datasets replicated past findingsof a principal gradient spanning a hierarchical axis from unimodal to transmodal cortex. Between-condition contrasts supported our primary hypothesis: relative to placebo conditions, this gradient exhibited a significant contraction in both LSD and psilocybin states, reflective of less differentiated hierarchical organization. This was a relatively symmetrical contraction in which somatomotor regions on the unimodal aspect and default and frontoparietal regions on the transmodal aspect both became less differentiated along this hierarchical axis. Gradient-based connectivity mapping indicated that this contraction has its basis in hierarchically specific changes in FC that are consistent with increased crosstalk between unimodal and transmodal cortices. Overall, there was strong convergence in the topographical changes across the LSD and psilocybin datasets, providing important support for the robustness of our findings given the small sample sizes. Results with the second and third gradient indicate alterations to axes of cortical organization related to sensory differentiation and executive region differentiation, respectively. The topography of changes for these two gradients showed strong overlap across datasets, but with greater significant changes with LSD, likely owing to power differences. One divergence was that LSD displayed reduced differentiation of visual and somatomotor cortex along the second gradient, which was not present in psilocybin. These findings extend previous work suggesting greater whole-brain FC of visual cortices in the LSD-state, a finding that is less consistent with psilocybin. Finally, we observed significant a significant association between LSD-related changes in the principal gradient and self-reported ego-dissolution -a core component of the psychedelic experience. Collectively, these results provide evidence for significant alterations of macroscale cortical gradients following psychedelic administration, marked by attenuation of hierarchical differentiation between unimodal and transmodal cortices as well as dedifferentiation along sensory and executive axes of cortical organization. Past neuroimaging investigations with LSD and psilocybin, including past analyses of the present datasets, have revealed that they elicit a complex set of changes to both static and dynamic FC, as well as to the entropy/complexity of regional timeseries. This work has often attempted to manage this complexity by describing results either in terms of a focus on specific large-scale network interactions (e.g., involving the default network) or in terms of general trends (e.g., towards global integration; Carhart-Harris,. This approach, while useful and necessary given the infancy of the field, can obscure important nuances in the structure of the data. For example, although decreased within-network and increased between-network FC is often listed as a consistent finding with serotonergic psychedelics, the effects with respect to specific networks/network pairs has limited overlap across drugs and datasets. Although this is to be somewhat expected given small sample sizes and analytical discrepancies in this area of research, it suggests the need for approaches which can help organize findings and facilitate comparisons. The cortical gradient approach used here provides a novel perspective on the neural effects of serotonergic psychedelics by collapsing the complex mosaic of increases and decreases in FC into a lowdimensional set of macroscale axes which suggest specific alterations in cortical information processing. This approach can facilitate both functional interpretation and the identification of effects specific to this class of drugs, which, evidently, decreased network segregation is not. As an example, healthy aging is also characterized by decreased within-network and increased between-network FC, yet gradient analyses comparing young and older adults have revealed qualitatively different effects than those observed here. The greatest LSD-dependent increases in the unimodal-transmodal hierarchical gradient were in regions at the lowest end of the cortical hierarchy -i.e., the somatomotor cortex, while the greatest decreases occurred within regions at the highest end -i.e., regions pertaining to the default and frontoparietal control networks such as the ventral posterior cingulate/retrosplenial cortex, the precuneus, and the superior frontal gyrus. This finding provides evidence for a relatively symmetrical contraction of hierarchical organization under LSD. Our gradient-based connectivity mapping approach further indicated that this contraction was underpinned by hierarchically specific changes in FC. In particular, our results revealed that regions on the unimodal side of the brain's hierarchy become less connected amongst themselves and more connected to transmodal regions, while regions on the transmodal side of the hierarchy maintain their connections amongst themselves and become more connected to unimodal regions. This suggests the movement of the brain into a more integrative mode of functioning, wherein unimodal regions are less modular and less differentiated from the globally distributed network of transmodal hubs. Moreover, the reduction in hierarchical differentiation and increased unimodal-transmodal crosstalk observed here is consistent with a reduction in the number of intervening processing steps between low-level sensorimotor cortex and high-level association cortex. This can be understood in the context of resting-state FC investigations of cortical organization which have provided evidence that cortical signals propagate through the macroscale processing hierarchy moving from modular sensory processing, to multi-modal integration, to higher-order processing within a distributed network of transmodal hubs. In addition, both meta-analytic and taskbased analyses have provided evidence that lower levels of the hierarchy pertain to behaviours that are coupled to immediate sensory input, while higher hierarchical zones pertain to perceptually-decoupled, abstract cognitive processes. Thus, in the psychedelic state, the observed levelling of hierarchy may correspond to decreased functional differentiation between sensory and abstract cognitive processing -consistent with the predictions of the REBUS model. This notion of dedifferentiation between concrete and abstract processing is consistent with past reports suggesting that psychedelics can elicit a blurring of the internal-external/subject-object distinction alongside an increased influence of internal mentation on perceptual processing (K. C. R.. Of particular interest is the phenomenon of 'ego-dissolution', which is a central component of the psychedelic experience and is associated with so-called 'mystical-type experiences' which are core predictors of the therapeutic efficacy of psychedelic treatment approaches. It corresponds to an experience of blurred or dissolved boundaries between self and world/other and a reduction in the perceived centrality of one's personal problems and concerns. We observed a significant positive association between LSD principal gradient scores and ego-dissolution within left medial prefrontal cortex/anterior cingulate cortex (mPFC/ACC) -a region that has been consistently implicated in self-related processing. This suggests a potential relationship between mPFC/ACC hierarchical differentiation and ego-dissolution. It is important to note, however, that replication in a larger sample is needed to ascertain the reliability of this effect given the small sample size and lack of replication in the psilocybin dataset. Nonetheless, this finding suggests the potential utility of characterizing changes in unimodal-transmodal differentiation for revealing novel associations between brain function and clinically-relevant dimensions of the psychedelic experience. In addition to changes in cortical hierarchy, we examined changes in the second and third gradient of cortical connectivity. The second gradient represents the differentiation of visual versus somatomotor/auditory processing. We found multiple changes in this gradient with LSD which were not present with psilocybin. Most notable, there was a significant reduction in visual network differentiation along this axis which was trending in the opposite direction with psilocybin. This is consistent with past work that has shown significant increases in V1 wholebrain FC in the LSD state, whereas similar changes are less consistent with psilocybin. The third gradient represents an axis which differentiates the so-called task positive (M. D./multiple demand/executive networkfrom the rest of the cortex. Results with both drugs broadly indicate a contraction of this gradient consistent with reduced differentiation of executive regions from the rest of cortex. Although speculative, this might relate to certain cognitive impairments present while under the influence of these drugs. Finally, to our knowledge, this is the first investigation which has shown significant pharmacologically-induced alterations of macroscale cortical gradients. In particular, we have demonstrated that cortical gradients can be acutely altered by two distinct serotonergic manipulations and, further, that these changes may be paralleled by specific changes in conscious experience. Although LSD and psilocybin exhibit a complex pharmacology, past work has robustly linked their characteristic 'psychedelic' effects and alterations to whole-brain functional dynamics to partial agonism at the 5-HT2A receptor. These findings therefore shed novel light on a potential role for the 5-HT2A receptor in modulating the functional differentiation of unimodal and transmodal systems, as well as in modulating macroscale axes of cortical connectivity more generally. More research is needed to investigate serotonergic influence on global brain dynamics. An important limitation to the present results is the relatively small sample sizes. For this reason, we encourage cautious interpretation of our findings. However, the convergence across drugs found here provides strong evidence for the robustness of the observed effects. Despite consistent evidence suggesting their safety in controlled research settings, collecting large datasets with psychedelics is currently difficult due to hurdles pertaining to funding and ethics board approval. We hope our findings serve as a foundational guide for future research in this nascent field and as motivation for replication in larger samples. The present study extends past findings on the neural underpinning of the psychedelic state by revealing that the whole-brain effects of LSD and psilocybin can be represented as a contraction the brain's macroscale functional hierarchy -directly in line with a recently proposed unified model of psychedelic brain action. We also provided evidence of reductions in sensory and executive region differentiation following psychedelic administration. Future work is needed to ascertain whether reductions in cortical hierarchy in the acute psychedelic experience directly relate to findings of rapid and sustained symptom reductions observed via psychotherapeutically-mediated experiences with these drugs. The findings of this study lend further weight to the view that psychedelics can be powerful tools for investigating brain organization and dynamics.

CONFLICTS OF INTEREST

MG reports receiving consulting fees from Psygen Labs Inc. and EntheoTech Bioscience. RLC-H reports receiving consulting fees from COMPASS Pathways, Entheon Biomedical, Mydecine, Synthesis Institute, Tryp Therapeutics, and Usona Institute. BB, JS, RNS, and LR report no conflicts of interest.

Study Details

Your Library