Multi-metric evaluations of acute psychedelic effects on fMRI brain entropy
This brain imaging study (n=28) examined how psilocybin affected several fMRI-based measures of brain entropy in healthy volunteers. Some entropy metrics rose with subjective and objective drug effects, but many did not change and the measures were only weakly related to each other, suggesting they do not capture a single brain process.
Authors
- Gitte Knudsen
- Patrick Fisher
- Dea Stenbæk
Published
Abstract
A prominent theory of psychedelics is that they increase brain entropy. Thirteen studies have evaluated psychedelic effects on fMRI brain entropy, each applying a distinct measure. Here we evaluated these metrics in an independent 28-participant healthy cohort with 121 pre- and post-psilocybin fMRI scans. We assessed relations between brain entropy and objective and subjective psychedelic drug effects using linear mixed-effects models. All metrics were evaluated using two parcellation strategies and 7 denoising pipelines. We observed consistent significant positive associations for Shannon entropy of the spatial eigendistribution of the time by voxel matrix, path-length, instantaneous correlations, brain-state switching, and sample entropy at short time-scales. We consistently did not observe significant effects for 8 of 14 entropy metrics and observe inconsistent positive effects for Lempel-Ziv complexity of the BOLD signal. Brain entropy quantifications showed limited inter-measure correlations. Our observations support a nuanced acute psychedelic effect on brain entropy, empirically demonstrating that these metrics do not reflect a singular construct.
Research Summary of 'Multi-metric evaluations of acute psychedelic effects on fMRI brain entropy'
βBlossom's Take
Introduction
Psychedelics are thought to produce altered states of consciousness partly through serotonin 2A receptor agonism, and clinical trials up to Phase IIb have suggested potential benefit in some psychiatric conditions when used with psychological support. In parallel, fMRI studies have been used to investigate whether acute psychedelic effects can be captured as changes in brain entropy, but earlier work had examined a set of 14 different entropy-related metrics in only a few datasets, often with small samples, different drugs or routes of administration, and no independent re-evaluation of the full set of measures. The authors frame this as an important gap because the Entropic Brain Hypothesis proposes that psychedelic phenomenology reflects increased entropy in functional brain signals, yet it was unclear whether the various published metrics actually measure a single construct or yield consistent findings. McCulloch and colleagues therefore set out to test 14 previously reported fMRI brain entropy metrics in an independent cohort of healthy participants given psilocybin, with concurrent measures of subjective drug intensity, plasma psilocin level, and estimated 5-HT2A receptor occupancy. They hypothesised that entropy would increase after psilocybin, and they also wanted to examine whether observed associations were robust to different parcellation schemes and denoising pipelines, and how strongly the different entropy measures correlated with one another. The study is positioned as an independent replication and comparison of candidate acute psychedelic biomarkers, with direct relevance to the Entropic Brain Hypothesis.
Methods
The study used a single-blind crossover design in 28 healthy volunteers, recruited from a database of people interested in psychedelic research. Participants were screened for neurological, somatic and psychiatric illness, and the protocol was approved by the relevant Danish ethics and medicines authorities and registered at ClinicalTrials.gov. The extracted text states that the data were collected between 2018 and 2021, and that no formal sample size calculation was performed, although the sample was larger than in previous publications in this area. Participants received a single oral dose of psilocybin at 0.2-0.3 mg/kg, with a mean dose of 19.7 mg, or 20 mg ketanserin; the ketanserin scans were outside the scope of this paper and are not presented. After drug administration, resting-state fMRI scans were acquired at approximately 40, 80, 130 and 300 minutes post-dose, producing 121 scan sessions in total. After each scan, participants rated current subjective drug intensity on a 0-10 scale, blood was drawn for plasma psilocin measurement, and 5-HT2A receptor occupancy was estimated from the plasma psilocin level using a previously reported Hill-Langmuir-based relation. The researchers refer to the combined subjective, pharmacokinetic and receptor-occupancy measures as PsiFx. Pre-processing was applied uniformly across metrics using SPM12 and CONN, including slice-timing correction where appropriate, unwarping, realignment, coregistration, segmentation, normalisation, smoothing, linear detrending, aCompCor, motion regression, scrubbing of artefactual volumes, and band-pass filtering from 0.008-0.09 Hz. Because data came from two scanners, some time series were temporally downsampled, and cerebellar regions were removed from atlases when not consistently in view. The main analyses regressed each entropy metric on PPL, SDI and estimated occupancy in separate linear mixed-effects models with a subject-specific random intercept, adjusting for motion, age, sex and scanner. Regional metrics were corrected for family-wise error across regions using maxT permutation correction, while whole-brain metrics used permutation p-values. The authors also repeated analyses across two parcellation strategies and seven denoising pipelines to assess robustness, and examined pairwise correlations between entropy measures. The paper then details how each entropy metric was operationalised, including static connectivity measures, dynamic connectivity measures, dynamic activity measures, and Lempel-Ziv complexity measures.
Results
Participants showed substantial subjective drug intensity and plasma psilocin after psilocybin, as expected. In the primary analysis, using the original parcellation and a standard denoising pipeline, the authors found significant positive associations between PsiFx and five of the 14 entropy metrics, one inconsistent positive association, and null results for eight metrics. The five metrics with consistent positive relations across denoising pipelines were geodesic entropy, normalised global spatial complexity (NGSC), dynamic conditional correlation (DCC) entropy, meta-state complexity, and sample entropy at scale 1. Lempel-Ziv temporal BOLD complexity (LZct) showed inconsistent positive effects across some pipelines, whereas spatial BOLD complexity (LZcs) was not associated with PsiFx. The authors note that the metrics were generally weakly correlated with one another. Among the individual metrics, out-network connectivity distribution entropy, degree distribution entropy, von Neumann entropy, intra-network synchrony distribution entropy, motif-connectivity distribution entropy, LEiDA-state Markov-rate, integration/segregation-state distribution entropy, and spatial BOLD complexity were not significantly associated with PsiFx in the primary analysis. Geodesic entropy showed a significant positive association, with weak-to-moderate correlations for PPL, occupancy and SDI; this replicated across most denoising pipelines except global signal regression. NGSC was significantly and moderately positively associated with PsiFx and remained robust across all denoising pipelines. Meta-state complexity was weakly positively associated with occupancy and SDI, but not clearly with PPL, and remained significant across most pipelines. DCC entropy showed the strongest and broadest signal: it was positively associated with PsiFx in 35 of 36 network-to-network connections, with moderate-to-strong effect sizes, although the pattern was sensitive to temporal filtering and could weaken or reverse with some preprocessing choices. Sample entropy showed a short-scale positive pattern at scale 1 in 7 of 17 networks and a long-scale negative pattern at scale 5 in 14 of 17 networks, but these results were sensitive to denoising and did not appear at intermediate scales. Temporal Lempel-Ziv complexity was associated with occupancy and SDI but not clearly with PPL in the primary model, and its significance varied across pipelines. When the authors re-ran the analyses using a common parcellation, the overall pattern was similar: geodesic entropy, meta-state complexity, DCC entropy, and short-scale sample entropy remained the main positive findings, whereas other metrics remained null. Scanner moderation analyses suggested that some associations varied by scanner, including geodesic entropy, sample entropy at scale 5, and von Neumann entropy, whereas DCC entropy and NGSC appeared more robust. In the correlation analysis between whole-brain entropy measures, some metrics were positively correlated, some negatively correlated, and several conceptually related measures clustered together, such as NGSC with von Neumann entropy, and LZcs with LZct. The extracted text also states that the authors considered the possibility that the original studies reported values that were not always mathematically plausible for some motif-based metrics, and that several prior findings were not reproduced here.
Discussion
The authors interpret the study as showing that acute psilocybin effects on fMRI entropy are real but heterogeneous rather than uniform across all entropy measures. They argue that they independently replicated increased geodesic entropy, increased NGSC, and the short- versus long-scale pattern previously reported for sample entropy, and they also identified a strong positive DCC entropy signal that had not been reported before. Meta-state complexity and temporal Lempel-Ziv complexity showed additional evidence of association across some preprocessing pipelines. At the same time, eight of the 14 metrics were not associated with psilocybin measures, which the authors say suggests that only some published entropy metrics are sensitive to acute psychedelic effects. Relative to earlier research, they note that their findings converge with some previous psilocybin, LSD and ayahuasca reports, but also fail to replicate a number of earlier positive results. They emphasise that this likely reflects the fact that the various metrics are not measuring a single underlying construct called “brain entropy”. Instead, the authors present the metrics as capturing distinct properties of functional brain activity, connectivity, or complexity. They also stress that the Entropic Brain Hypothesis is broad and not specific about which entropy measure should be expected to change, making the present mixed findings difficult to classify as simple support or refutation. The authors discuss several limitations. The study used two scanners with different sequences, requiring some temporal downsampling, although each participant was scanned on only one scanner. There was no placebo condition and no active non-psychedelic comparator, so the specificity of the observed effects for psychedelics cannot be determined. Motion increased with psilocybin intensity, and although motion was adjusted for and scrubbed, the authors could not exclude residual confounding. Scan durations were relatively short for some sessions, which may affect the stability of entropy estimates, and multi-band acquisition may have influenced signal quality. They also did not directly measure physiological variables such as respiration or heart rate, and they note that psychedelic effects on cerebral blood flow and neurovascular coupling could complicate interpretation of BOLD-based entropy measures. For implications, the authors suggest that future work should include active comparators such as stimulants, ketamine, MDMA and cannabinoids to determine whether these entropy changes are specific to serotonergic psychedelics or reflect broader altered brain states. They also recommend that future psychedelic fMRI studies measure subjective intensity at scan time and collect plasma drug levels, and they point to the need for methodological transparency, independent replication and further work on preprocessing sensitivity. They additionally mention a MATLAB-based toolbox developed for calculating the metrics studied here, which they present as a resource for future research.
Conclusion
The authors conclude that psilocybin produces consistent positive effects on a subset of fMRI entropy and complexity metrics, especially NGSC, DCC entropy, geodesic entropy, meta-state complexity and short-scale sample entropy, but not on many other published measures. They state that the limited inter-correlation between metrics and their differential sensitivity to denoising indicate that these measures do not represent a single unitary construct of brain entropy. Overall, they describe their findings as nuanced support for the Entropic Brain Hypothesis and as identifying candidate biomarkers for future clinical psychedelic studies, while emphasising the need for independent replication and transparent methodology.
View full paper sections
A R T I C L E I N P R E S S ARTICLE IN PRESS
In the current study, we sought to evaluate the acute effects of 0.2-0.3 mg/kg psilocybin administration on these 14 brain entropy metrics in an independent Dataet of 28 healthy individuals. Based on the EBH, we hypothesised that brain entropy metrics would be increased following psilocybin administration. Participants completed a 5 or 10min resting-state fMRI scan a single time before, and multiple times following psilocybin administration (121 total scan sessions). All scans for each participant were performed on one of two scanners. Each scan was accompanied by a self-report measure of subjective drug intensity (SDI) and a blood sample to quantify plasma psilocin level (PPL), from which brain serotonin 2A (5-HT2A) receptor occupancy (Occ 2A ), was estimated based on its relation to PPL established in a previous study from our lab. We evaluated the relations between each entropy metric and SDI, PPL and Occ 2A using linear mixed-effect models with a subject-specific random intercept and correction for age, sex, scanner, and motion (see Methods Supplementary Material for more details). We report uncorrected permutation p-values (p perm ) for each metric producing a single whole-brain value. For entropy metrics producing regional values, we report family-wise error rate corrected permutation p-values (p FWER ) with correction across all regions within each metric using maxT correction. We report as significant those metrics which were significantly (i.e., p < 0.05) associated with SDI, PPL and Occ 2A , collectively referred to as "PsiFx". Standardised effect sizes are reported (Pearson's rho). This evaluation was repeated across two parcellation strategies and seven preprocessing pipelines to explore the robustness of effects to pre-processing decisions. Finally, we explored the intercorrelation between brain entropy metrics to characterise their associations.
RESULTS
Participants showed substantial SDI and PPL following drug administration as anticipated (Supplementary Figure). See Figurefor a summary of entropy metrics. Primary analyses used the spatial parcellation reported in the original publication (with cerebellum removed) and a generic denoising pipeline (12 motion parameters, aCompCor, interpolation (scrubbing) of artefactual volumes, and bandpass filtering in the range 0.008-0.09Hz). To explore moderating effects of denoising pipelines on brain entropy metrics, we considered six variations 1) including global-signal regression (GSR), 2) removal of the low-pass filter (0.09 Hz), 3) applying a narrow bandpass filter (0.03-0.07 Hz), 4) regressing 24 motion parameters, 5) omitting scrubbing, and 6) a stricter scrubbing threshold. For a summary of results across denoising pipelines see Supplementary Data S8 and for a full breakdown across all pipelines see Supplementary Data S9. Here we show a significant and consistent relation between objective (i.e., PPL) and subjective (i.e., SDI) measures of psilocybin effects and five of the 14 entropy metrics evaluated, inconsistent but positive effects for one, and consistently null effects for eight metrics. Those demonstrating a consistent relation with psilocybin effects across denoising pipelines were distributed across entropy classes: two features of static connectivity, two of dynamic connectivity and one of dynamic activity. No single entropy metric stands out as the clearest candidate for future research. Our observation that the metrics are generally weakly correlated with each other emphasises that they represent distinct phenomena, i.e., "brain entropy" is not a singular construct.
OUT-NETWORK CONNECTIVITY DISTRIBUTION ENTROPY
Shannon entropy of regional distribution of out-network correlation coefficients was not significantly associated with psilocybin effects (PsiFx) in any of the 181 non-cerebellar brain regions after controlling for multiple comparisons (p FWER > 0.07 for all regions for at least one of PsiFx, Supplementary Data S2). Across additional denoising pipelines, only one region had a moderate significant positive relation and one had a weak significant negative relation, each for only a single pipeline (Supplementary Data S8, S9).
DEGREE DISTRIBUTION ENTROPY
Shannon entropy of global degree distribution at a correlation coefficient threshold corresponding to a mean degree of 27 was not associated with any of PsiFx (t(90) = -0.03, p perm = 0.99, rho = 0.00, 95% CI [-0.0054, 0.0052]; Figure). We also did not observe significant effects for thresholds producing a mean degree between 1 and 48 (Supplementary Data S3). Across additional denoising pipelines, there were sparse weak-to-moderate significant positive associations in three of six pipelines, with little overlap in which degrees were significant (Supplementary Data S8, S9).
GEODESIC DISTRIBUTION ENTROPY
The entropy of path-length (geodesic) distribution was significantly positively associated with PsiFx at the a priori described threshold producing mean degree 27 (t(90) = 2.04, p perm = 0.04, rho = 0.23, 95% CI [0.0001, 0.0069]; Figure) . The associations were weak to moderate (Pearson's rho = 0.39, 0.27, and 0.23 for PPL, Occ 2A and SDI, respectively). Significant weak to moderate positive associations with PsiFx were also observed across thresholds producing mean degrees from 22 to 38 (Supplementary Data S3). This weak-to-moderate positive effect
S ARTICLE IN PRESS
was consistent across five of six additional denoising pipelines, the exception being GSR (Supplementary Data S8, S9).
VON NEUMANN ENTROPY
The Von Neumann entropy of correlation matrices was not significantly associated with PsiFx in the primary analysis (t(90) = 0.14, pperm = 0.90, rho = 0.02, 95% CI [-0.0007, 0.0008]; Figureand Supplementary Data S3)Supplementary Data. However, weak-to-moderate significant positive effects were observed in three of six denoising pipelines for some PsiFx, with the strongest relation following GSR (Pearson's rho = 0.43, 0.52 and 0.36 for PPL, Occ 2A and SDI, respectively) (Supplementary Data S8, S9).
NORMALIZED GLOBAL SPATIAL COMPLEXITY
Normalized Global Spatial Complexity (NGSC) was significantly moderately positively associated with PsiFx (t(90) = 4.30, pperm < 0.0003, rho = 0.47, 95% CI [0.0006, 0.0016]; Figureand Supplementary Data S3; Pearson's rho = 0.47, 0.44, and 0.52 for PPL, Occ2A and SDI)Supplementary DataThis relation was consistent across all denoising pipelines (Supplementary Data S8, S9).
INTRA-NETWORK SYNCHRONY DISTRIBUTION ENTROPY
Intra-network synchrony distribution entropy was not significantly associated with PsiFx in any of nine networks in the primary analysis (All p FWER > 0.98, Supplementary Data S2, Supplementary Figure). Following the removal of the lowpass filter, the visual and somatomotor networks showed a significant moderate negative relation with all PsiFx (Supplementary Data S8, S9).
MOTIF-CONNECTIVITY DISTRIBUTION ENTROPY
The four-ROI motif-connectivity state distribution was not significantly associated with PsiFx at any window length from 15 to 150s, except a single weak association at window-length 100s (t(90) = 1.96, pperm = 0.05, rho = 0.24, 95% CI [-0.000, 0.035]; Pearson's rho = 0.30, 0.25, and 0.24 for PPL, Occ2A, and SDI) surrounded by nonsignificant findings (Supplementary Figureand Supplementary Data S3, S9). Across only four denoising pipelines, only few individual window lengths were significantly associated with PsiFx, some of which were positive and some negative (Supplementary Data S8, S9).
LEIDA-STATE MARKOV-RATE ENTROPY
LEiDA-state Markov-rate was not significantly associated with PsiFx (t(90) = 0.17, p perm = 0.94, rho = 0.02, 95% CI [-0.0029, 0.0034]; Figure, Supplementary Data S3, S9)Supplementary Data. A weak positive relation with SDI and Occ 2A was observed only following GSR (Pearson's rho = 0.26 and 0.21, respectively) (Supplementary Data S8, S9).
META-STATE COMPLEXITY
Meta-state complexity was positively weakly associated with Occ 2A (t(90) = 2.14, p perm = 0.03, rho = 0.22, 95% CI [0.00, 0.13]) and SDI (t(90) = 3.08, pperm = 0.003, rho = 0.33, 95% CI [0.002, 0.011]), but not PPL (t(90) = 1.84, pperm = 0.08, rho = 0.20, 95% CI [-0.0002, 0.0050], Figure, Supplementary Data S3). Supplementary DataThe weak positive associations remained significant across five of the six additional denoising pipelines (Supplementary Data S8, S9).
INTEGRATION/SEGREGATION-STATE DISTRIBUTION ENTROPY
Integration sub-state entropy was not significantly associated with PsiFx in either the primary or additional pipelines (t(90) = 0.83, p perm = 0.40, rho = 0.09, 95% CI [-0.0014, 0.0033]; Figure, Supplementary Data S3, S8, S9) Supplementary Data.
DYNAMIC CONDITIONAL CORRELATION DISTRIBUTION ENTROPY
Dynamic conditional correlation (DCC) entropy was significantly positively associated with PsiFx in 35 of 36 network-network connections (18/36 p FWER < 0.0001, i.e., observed data superseded all permutations, 29/36 p FWER < 0.001, 35/36 p FWER < 0.05; Figure; Supplementary Data S2). Associations were moderate to strong (t(90) range: 2.30 to 9.35, Pearson's rho range: 0.35 to 0.78, Supplementary Data S2). The one association with at least one non-significant relation was for edges within the motor cortex (t(90) = 6.83, p FWER = 0.13, rho = 0.30). The positive associations replicated in four of six denoising pipelines. DCC entropy effects were particularly sensitive to temporal filtering: employing a narrow bandpass filter removed the associations and removing the low-pass filter changed to negative the association for 8 of 36 network-network connections. (Supplementary Data S8, S9).Supplementary DataSupplementary DataSupplementary Data
MULTI-SCALE SAMPLE ENTROPY
At scale 1, (i.e., no time-series compression), sample entropy was significantly positively associated with PsiFx (p FWER < 0.05) in 7 of 17 networks (i.e., Central Visual, Dorsal Attention A, Control A, B and C, Default-Mode A and C). This finding replicated only in two of six additional denoising pipelines. When the low-pass filter was removed, negative associations were observed (Supplementary Data S8, S9). At scales 2, 3, and 4, no associations were significantly associated with PsiFx (p FWER > 0.05). At scale 5, sample entropy was significantly negatively associated with PsiFx in 14 of 17 networks; Control A and C and Default-Mode A (p FWER < 0.001), Somatomotor A, Dorsal Attention A and B, Salience-Ventral-Attention B, Limbic B, Control B, Temporal-Parietal, Default-Mode B and C (p FWER < 0.05, Figure, Supplementary Data S2). Scale 1 associations were weak to moderate (t(90) range: 2.39 to 4.57, Pearson's rho range: 0.26 to 0.47), as were Scale 5 associations (t(90) range: -4.92 to -2.53, Pearson's rho range: -0.49 to -0.27). The scale 5 finding only replicated for one additional denoising pipeline, i.e., removal of the lowpass filter (Supplementary Data S8, S9).
SPATIAL AND TEMPORAL DYNAMIC BOLD COMPLEXITY
Temporal BOLD complexity (LZct) was not associated with PPL (t(90) = 1.53, p perm = 0.14, rho = 0.17, 95% CI [-0.0001, 0.0006]; Figure) but was associated with Occ2A (t(90) = 2.20, p perm = 0.03, rho = 0.23, 95% CI [0.001, 0.016]) and SDI (t(90) = 2.67, p perm = 0.009, rho = 0.30, 95% CI [0.0002, 0.0013]). Associations were weak (Pearson's rho = 0.23, 0.30 and 0.17 for Occ2A, SDI, and PPL respectively). Spatial BOLD complexity (LZcs) was not significantly associated with PsiFx (t(90) = -0.21, p perm = 0.82, rho = -0.02, 95% CI [-0.013, 0.011]; Figure, Supplementary Data S3). LZct was significantly moderately associated with all PsiFx in two of six denoising pipelines, i.e., no scrubbing and strict scrubbing. LZcs was not associated with PsiFx for any pipeline (Supplementary Data S8, S9).Supplementary DataSupplementary Data
EFFECT OF PARCELLATION
Main analyses were performed using the same atlas as the original publication. To evaluate PsiFx on brain entropy metrics using a common parcellation, all analyses using the original denoising pipeline were re-run using an atlas combining the Schaefer-100 (7 Yeo networks) and Tian-16 subcortical atlases, see Supplementary Data S5 for detailed results. Geodesic entropy showed a weak positive association with all PsiFx at mean degrees 31 to 38. Meta-state complexity was weakly positively associated with Occ, but not PPL nor SDI. LEiDA-state Markov-rate was weakly negatively associated with all PsiFx and DCC distribution was weak-to-strongly associated with all PsiFx across most network edges. Sample Entropy was weak-to-moderately associated with PsiFx at scale 1 but no significant associations were observed at longer scales, although the trend of increased entropy at short scales and decreased entropy at long-scales was maintained. No parcellation was applied for NGSC (voxel-wise) or Motifconnectivity Distribution Entropy (uses predefined seeds). All other metrics were not significantly associated with PsiFx.
MODERATING EFFECT OF SCANNER
To evaluate scanner effects on observed associations from the main analysis pipeline, we fit linear mixed models estimating the moderating effect of scanner, see Supplementary Data S6 for detailed Results. We observed a significant moderating effect of scanner on the relation between PPL and entropy for geodesic entropy over the range of thresholds at which we observe significant associations with PsiFx. The nature of this interaction was that for scanner A the effect was closer to zero than for scanner B. We also show a significant moderating effect of scanner for Sample Entropy scale 5 across many ROIs that are significantly associated with PsiFx. However, the observed effect is the same direction for scanner A and B, only numerically stronger for scanner A. We observed a significant moderating effect of scanner for Von Neumann entropy such that the associations with PsiFx were opposite for the two scanners. We did not observe a significant moderating effect of scanner for DCC entropy, nor NGSC, supporting the robustness of these metrics. For entropy metrics that were not associated with PsiFx in the main model, we observed moderating effects of scanner for I/S-state distribution, some window-lengths of motifconnectivity distribution, some ROIs of out-network connectivity, degree distribution, and some sample entropy scale 4 ROIs.
CORRELATION BETWEEN WHOLE-BRAIN ENTROPY QUANTIFICATIONS
We estimated the correlation between whole-brain entropy metrics to explore their association with one another across all included scans from the main analysis pipeline. Some metrics were positively and some negatively
DISCUSSION
Overview Recent studies have reported acute psychedelic effects on functional brain entropy, but to date none of these metrics have been independently reevaluated. In this study we evaluated 14 previously reported entropy metrics in an independentsample of 28 healthy participants scanned with BOLD fMRI before and several times after psilocybin with concomitant measurements of subjective drug intensity and plasma psilocin level. We observed statistically significant psilocybin effects that echoed previous reports for three brain entropy metrics: geodesic entropy, wherein we replicate increased entropy at previously reported thresholds; sample entropy, wherein we replicate a previously observed increase in entropy at short scales and decrease at long scales; and NGSC where we replicate previously observed increases. We observed a strong positive relation between PsiFx and brain entropy measured by Dynamic Conditional Correlation analyses that has not been previously reported. Meta-state complexity showed significant relations with PsiFx across most preprocessing pipelines, and temporal Lempel-Ziv complexity showed evidence for associations with PsiFx across some pipelines. For 8 of 14 brain entropy metrics previously reported, we did not observe a significant association with psilocybin measures, and we see limited correlation between entropy metrics. These mixed findings underscore the importance of corroborating outcomes in independent Dataets. Although we observe some evidence supporting the entropic brain hypothesis, these variable findings underscore the broadness of this theory and the need to more clearly establish which brain entropy metrics of functional brain imaging signals are acutely affected by psychedelics.
NORMALISED GLOBAL SPATIAL COMPLEXITY
We found that NGSC, the entropy of (spatial) eigenvalues of z-scored voxel-wise data, was significantly positively associated with PsiFx across denoising pipelines. This is consistent with the original report from Siegel et. al.,who report increases following oral psilocybin and intravenous LSD using the publicly available Dataet. This indicates that psilocybin increases the uniformity of variance distribution across spatial principal components. Not noted previously, NGSC and Von Neumann entropy are nearly mathematically equivalent, Von Neumann entropy uses a parcellated, region-wise instead of voxel-wise, correlation matrix. Unsurprisingly, they are strongly positively correlated in our sample (Pearson's rho = 0.80). Despite this, Von Neumann entropy was not associated with PsiFx, except after GSR, implying that GSR affects parcellation more than voxel-wise connectivity matrices. Although reducing spatial dimensionality may be necessary for comparison of brain regions across participants, it may inadvertently average away complexity observable at the voxel level. Future studies exploring these trade-offs could better highlight when one method or the other is preferred. Our results provide evidence for NGSC as a sensitive biomarker of psychedelic drug effects, robust to denoising strategy. However, we are not aware of any fMRI studies evaluating the effect of other psychoactive drugs on NGSC so cannot infer specificity.
GEODESIC ENTROPY
We report a significant positive association between the Shannon entropy of the distribution of path lengths across the whole brain as previously reportedand all three psilocybin metrics evaluated: PPL, Occ 2A , and SDI. We observed statistically significant associations at a range of correlation coefficient thresholds that produce graphs with a mean degree from 22 to 38, the previous study reported significant differences at thresholds producing mean degrees from 24 to 35. Characteristic path length is a description of the number of edges that must be traversed to get from any one brain region to another, a putative measure of capacity for information flow. Our results suggest that one of the effects of psilocybin on the brain can be described as a broadening of the histogram of path lengths across region-to-region connections. Notably, this does not imply that the average path length is shorter or longer, only that there is a wider distribution of these across the whole brain. Our convergent results are encouraging considering that the previously reported Dataet used a different drug (ayahuasca, which contains MAOIs as well as the psychedelic N,N-dimethyltryptamine) and different imaging parameters, suggesting robustness of the metric. This association with PsiFx was relatively robust to pre-processing strategies, yet was sensitive to GSR, suggesting that acutely increased global connectivity often reported following psychedelic drug administration may directly affect geodesic entropy. We are not aware of other studies evaluating geodesic entropy, so comparison to other drugs or psychiatric conditions are not yet possible and should be evaluated in future studies.
DYNAMIC CONDITIONAL CORRELATION DISTRIBUTION ENTROPY
We observed a statistically significant positive relation between PsiFx and DCC entropy for all within or between network relations except within the motor network. DCC entropy is a measure of the width of the distribution of instantaneous connectivity values for any region-region edge across each scan. The previous study found no change in DCC distribution at one-week and one-month post administration of psilocybinand importantly did not evaluate acute effects. We observed moderate-to-strong correlations with PsiFx (i.e., Pearson's rho with PPL > 0.7 for three edges, all including DMN (DMN-DMN, DMN-Frontoparietal, DMN-Medial Frontal), and Pearson's rho > 0.5 for 28/36 network edges). The strength of these associations is remarkable, perhaps as large as any previously reported fMRI effect of psychedelic action, suggesting that DCC distribution may be a strong candidate neural correlate of acute psychedelic effects and among the strongest correlations observed in pharmaco-fMRI. Our results suggest that psilocybin increases the variability of connectivity between regions across time across almost all region-region pairs, which are summarised into networks. Furthermore, this association was robust to most preprocessing strategies. Notably, the association became weakly negative in some network edges when the low-pass filter was removed, suggesting that PsiFx effects on DCC entropy manifest mostly in relatively slow frequencies and are obscured by higher frequencies. As above, we are not aware of other pharmaco-fMRI studies evaluating DCC distribution. Notwithstanding, the sheer magnitude of the observed associations suggests DCC distribution may be a sensitive marker for acute psychedelic effects on the brain and so we encourage independent replication.
MULTI-SCALE SAMPLE ENTROPY
We observed a significant positive relation between PsiFx and scale-1 sample entropy (i.e., temporal resolution = 2-seconds) in seven out of 17 brain networks. Conversely, we observed a significant negative relation between scale 5 sample entropy (temporal resolution = 10-seconds) in 14 out of 17 networks. However, these associations at both scales were affected by the denoising pipeline, only maintaining a positive association for two and one additional pipelines respectively and even negative associations at the short scale after removing the low-pass filter. Multi-scale sample entropy measures the irregularity of a signal over its entire length. Increased sample entropy in
S ARTICLE IN PRESS
most networks at scale-1 and decreased sample entropy at scale-5 align in both cases with the original observation following LSD. fMRI-measured multi-scale sample entropy has been shown to be increased in the default-mode, visual, motor and lateral-prefrontal networks following caffeineand decreased at scale 1 during sleep, although certain parameters used in their calculations were different to those employed here. As such, it is possible that the effects that we, and Lebedev and colleagues, observed may reflect differences in wakefulness and may thus be nonspecific to psychedelic effects. Positive symptoms of schizophrenia have been positively associated with sample entropy at scales 1 and 2, and negatively associated in certain brain regions at scales 3, 4, and 5. This is consistent with our observations and is also phenomenologically consistent, as the high-dose psychedelic state has some overlap with some positive symptoms of schizophrenia e.g., verisimilitude, alterations in visual perception (though psychedelics do not normally produce 'true' hallucinations, i.e., sensory appearances indistinguishable from reality, as are present in schizophrenia), and sense of self. We are aware of the problematic history of psychedelic 'psychotomimetic' research and urge caution in overinterpretation of this apparent convergence. Our convergent results with Lebedev and colleagues are intriguing considering that the original paper reported effects following intravenous LSD administration whereas we administered psilocybin orally. However, the diverging effects of denoising pipelines suggests that this metric is susceptible to retained signal frequencies and global signal. Beyond denoising, regional SampEn has also been shown to be positively associated with systolic blood pressure in HCP data across subjects. Although a between-subject association does not necessarily translate to within-subject drug effects, psilocybin acutely raises blood pressureso residual blood-pressure-related variance is a plausible contributor to the SampEn changes we observed. Taken together, increases in sample entropy at short temporal resolutions may be a candidate biomarker for psychedelic effects, if they cannot be explained by, e.g., wakefulness state or physiological confounds.
LEMPEL-ZIV COMPLEXITY
Intriguingly, Lempel-Ziv complexity of two measures (meta-state complexity and temporal BOLD complexity (LZct)) were significantly positively associated across several denoising pipelines. Meta-state complexity was significantly weak-moderately associated with at least one PsiFx in all but one pre-processing pipeline. Furthermore, recent work has shown a negative relation between meta-state complexity and Global Control Energy (GCE), representing the amount of input required to shift from one brain-state to the next. GCE has been shown to decrease during DMT infusion and associated with LZct of simultaneously collected EEG signals. This suggests a potential bridge between fMRI meta-state complexity and EEG LZct through the lens of Network Control Theory. These convergent findings support meta-state complexity as an insightful metric of psychedelic effects on brain activity. However, the original study of meta-state complexity does not report a statistical analysis of intravenous LSD nor intravenous psilocybin effects, but the data are publicly available and do not support a significant effect of either drug. The original study of BOLD complexity reports an increase in spatial BOLD complexity following LSD, but not psilocybin and does not report any findings pertaining to temporal BOLD complexity. One MEG and four EEG studies have reported increased LZc following psychedelic administration 54-58 , providing convergence for its utility as a marker of psychedelic effects, however one reports increases in LZct in the absence of subjective drug effects, indicating a potential epiphenomenon. Given this intriguing convergence across modalities, future studies should evaluate whether there is a relation between LZct of the EEG and BOLD timeseries. Across previous studies analysing regional timeseries, there is inconsistency in the quantification of LZc in the temporal (LZct) and spatial (LZcs) domain. Our borderline statistically significant associations were observed for LZct only. It is our perspective that LZct is more sensible and should be used in future studies as it preserves region-specific temporal information, whereas LZcs is sensitive to arbitrary region order.
NULL FINDINGS
We consistently did not observe a significant association with PsiFx for 8 of the 14 brain entropy metrics considered here (Figure). Of these, the original studies reported either increased entropy following psychedelic drug administration, no effect, or did not formally evaluate the effect of psychedelic drug administration(Figure, Supplementary Data S1). Our observed entropy estimates for hippocampal-ACC motif entropy are markedly different from those previously reported. We are concerned that the originally reported values are not mathematically possible, see the Supplementary Text and Supplementary Figurefor a detailed consideration. Although our null findings with respect to these metrics does not establish that they have no relation to acute psychedelic drug effects, they imply a smaller relation that limits their utility as biomarkers of acute psychedelic effects. The discrepancy between our observations and those reported previously underscores the need to replicate or corroborate findings in independent cohorts to validate initial reports. Our inability to replicate previous findings may be due to greater statistical power and different statistical models i.e., linear regression with PsiFx. Incongruence may also be attributed to differences in data collection. All previous studies reporting acute effects on brain entropy following a range of routes of administration and across psychedelic drug classes though the acute effect appear similar, we cannot rule out differences due to drug or route of administration. However, if entropic brain effects are not consistent across drugs this would indicate that these metrics are not useful neural correlates of the psychedelic experience. We encourage all future fMRI studies evaluating psychedelic effects on brain function to measure subjective drug intensity at time of scanning and to collect plasma samples for quantification of plasma drug levels as described in the psychedelic fMRI consensus paper.
INTER-CORRELATION BETWEEN ENTROPY METRICS
Despite the large set of brain entropy metrics that have been reported previously, no studies have considered whether these measures are inter-correlated. We observed positive associations between geodesic and degree distribution entropies, which are based on the same graph-theory representation of connectivity, between LZcs and LZct which are conceptually very much related, and between the mathematically similar NGSC and Von Neumann entropy as discussed above. Notably, we observed four pairs of brain entropy metrics that were significantly negatively correlated. This highlights the importance of specificity in describing "brain entropy." Many of these metrics represent distinctly different constructs, their individual meaning and collective representation of psychedelic effects is muddied by superficially considering them all metrics of "brain entropy." Future studies should be cognisant of this variable relation in considering whether findings are consistent or convergent across studies.
IMPLICATIONS FOR THE ENTROPIC BRAIN HYPOTHESIS
We evaluate this set of metrics in part to address uncertainties around The Entropic Brain Hypothesis, a prominent model that suggests 1) there is a relation between the magnitude of entropy of a brain activity signal and the subjective quality of that conscious state, 2) subjective psychedelic effects stem from a drug-induced increase in brain function signal entropy. (See Figureof Carhart-Harris et al., 2014). Herein we have observed that only certain types of fMRI brain-entropy are associated with psychedelic drug effects. Notably, the EBH does not articulate a specific entropy quantification, our mixed findings highlight a relevant lack of specificity within the hypothesis; it is difficult to clearly say whether we find support for or against it. Our findings do reproduce and identify significant effects, e.g., geodesic entropy, sample entropy, DCC entropy, and NGSC. It is also notable that some metrics advanced as supporting this hypothesis, e.g., LZ complexity, are not measures of entropy, "in the information theoretic sense" as described in the original hypothesis paper, but rather complexity. The limited correlation between measures highlights that the previous studies provide evidence for distinct hypotheses and not necessarily support for a singular model. Extending beyond entropy, future work should consider whether candidate entropy metrics align with other models of psychedelic brain-action such as the REBUSand the thalamic gating hypothesis.
ALTERNATIVE NEUROIMAGING TECHNIQUES
Although we focus here on fMRI quantifications of entropy, it is worth noting that psychedelic effects on brain entropy, specifically LZc, have been applied to five original MEG and EEG Dataets, over eight papers. The entropy measures applied in these studies leverage the high temporal sampling rate that is not clearly applicable to temporally slower fMRI and were not evaluated here. Further, the methods capture different aspects of physiological response to psychedelics. Future work evaluating psychedelic effects on brain entropy using multimodal neuroimaging and evaluating relations between alternative quantifications of brain-entropy will contribute meaningfully to the field.
LIMITATIONS
Brain imaging data were acquired on two different MRI scanners with different sequences (e.g., different TRs) requiring temporal downsampling of some data to match the other. However, each participant was scanned on only one scanner, enabling us to map within-subject changes onto PsiFx independent of scanner differences. Our study did not include a placebo condition, but we did acquire a pre-drug scan with which we estimated brain entropy metrics in the absence of psilocybin effects. Due to the large within-subject variability in fMRI outcomes in participants scanned several days apart, pre-drug vs post-drug scans performed on the same day may be superior to placebo scans performed weeks apart for evaluating drug effects because it limits this within-subject variance component. However, we did not administer negative control substances that would allow investigation of specificity. Consequently, we cannot determine which metrics are specific to psychedelics. Future research aiming to identify psychedelic-specific entropy biomarkers should include active comparators from drug classes including stimulants (amphetamine, caffeine), NMDA agonists (ketamine, nitrous oxide), MDMA, and cannabinoids (delta-9-THC). Such comparative studies would clarify whether observed entropy changes reflect serotonergic 5-HT2A agonism specifically or broader alterations in brain state complexity common to various psychoactive substances. PPL and SDI were associated with increased motion in the scanner, see Supplementary Figure; although we included an estimate of motion as a covariate in our models, employed scrubbing, motion correction and denoising strategies, and show that motion was not positively associated with any whole-brain entropy metrics, we cannot rule out that motion confounds our reported effects. It has been reported in many groups that head motion is increased following psychedelic drug administration so this is not a limitation unique to our data. Most fMRI scans were 10 minutes long, though some were only five. This may not be long enough to derive stable estimates of brain entropy metrics, e.g., previous studies have recommended >13 minutes for single-echo fMRI. Our data include scan sessions of three different volume lengths due to differences in TR and scan duration; notably, all scans for each subject were the same length. Entropy quantifications that preserve the temporal ordering of the time-series, such as sample entropy and the Lempel-Ziv complexity measures may be sensitive to time-series length. Future studies resolving optimal scan parameters for estimating stable brain entropy estimates would benefit this emerging field. Approximately half of the scans analysed herein utilised a multi-band acceleration protocol that may negatively affect signal-to-noise, though these effects may be less pronounced for task-free imaging as performed here. We corrected for estimated physiological noise using aCompCor but did not measure physiological effects such as by using eye-tracking or statistically model physiological effects such as changes in respiration, heart rate, or vasoconstriction which are affected by psilocybin and may have confounded our findings. Psilocybin has been shown to increases cerebral blood flow, and LSD has been shown to decrease it. Furthermore, the psychedelic DOI has been shown to alter neurovascular coupling i.e., the relation between neural activity and hemodynamic response, in mice. These haemodynamic effects of psychedelics may alter the BOLD signal in a way that confounds the interpretation of the underlying neural effects. Future research should aim to quantify psychedelic effects on neurovascular coupling in humans and the effect of this on signal complexity. Our statistical models assume a close temporal relation between brain entropy and PsiFx (measured adjacent to scans), thus, if changes in brain entropy occur after PsiFx, they would not be well captured.
DENOISING
A prevailing challenge in fMRI research is how best to handle the enormous flexibility in denoising strategies. Here we explored this space by evaluating the moderating effect of different denoising strategies on brain entropy metrics. Most of our results were robust to atlas or parcellation choice. The associations with PsiFx of some metrics (DCC entropy, NGSC, meta-state complexity and geodesic entropy) were robust to most denoising strategies, whereas other metrics were sensitive to denoising strategy (e.g., Sample entropy). For those metrics which remained significantly associated with PsiFx across pipelines, the strength of some associations varied across pipelines. Notably, removing the low-pass filter and applying a relatively narrow bandpass filter both substantively affected the statistical relations between PsiFx and brain entropy metrics. This is consistent with previous reports that preprocessing decisions can influence observed effects on fMRI outcome measures, which underscores the need for future studies in large, normative Dataets to evaluate brain entropy metric characteristics across this denoising multiverse. This is all the more relevant to advance their predictive or prognostic utility in clinical cohorts. For the purposes of this manuscript, we show that our main denoising pipeline was very similar to all previously applied in this space (See Supplementary Data S7) and thus believe that our results are directly comparable with previously reported findings. We also present here the Copenhagen Brain Entropy toolbox (CopBETolsen/CopBET), a Matlab-based toolbox containing functions to calculate the entropy metrics reported here, supporting the reproducibility of this and future studies.In conclusion, we observed consistent positive effects of psilocybin administration on five considered complexity/entropy metrics: NGSC, DCC entropy, geodesic entropy, meta-state complexity, and sample entropy (scale 1). DCC entropy showed remarkably strong effects, and NGSC demonstrated robust associations across all denoising pipelines. LZct showed inconsistent effects and eight showed very limited or no associations. Limited inter-correlations between metrics and differential denoising sensitivity indicate that these measures capture distinct features of complexity and entropy of BOLD signals, rather than a unitary construct. Our findings provide nuanced support for the Entropic Brain Hypothesis, indicating that not all entropy metrics may reliably capture psychedelic effects. These results establish candidate biomarkers for clinical psychedelic studies and underscore the critical importance of independent replication and methodological transparency in psychedelic neuroimaging research.
METHODS
Twenty-eight healthy volunteers participated in the study (10 female, mean age ± SD : 33 ± 8) and were recruited from a database of individuals interested in participating in a study involving psychedelics. Sex of participants was determined by the sex described in their medical records. A detailed description of the study design can be found in the Supplementary Text and has been reported previously. The study protocol was approved by the ethics committee of the capital region of Copenhagen (H-16026898) and the Danish Medicines Agency (EudraCT no.: 2016-004000-61). The study was registered at clincaltrials.gov (NCT03289949). Data presented here were collected between 2018 and 2021. A subset of the functional brain imaging data presented here has been included in different studies reported previously. Details of recruitment, procedures during the psilocybin session, ethical approvals, MRI acquisition and quality control, are described in the Supplementary Text. Analyses were preregistered on the 3rd of August 2022 (). Some analyses that met our inclusion criteria (i.e., fMRI studies investigating entropy changes pertaining to psychedelics) were identified after preregistration and were added. No statistical methods were used to pre-determine sample sizes but our sample sizes are larger than all previous publications (see Figure).
DATA COLLECTION
After obtaining written informed consent and screening for neurological, somatic and psychiatric illness, participants completed a single-blind, cross-over study design wherein participants received a single 0.2-0.3 mg/kg dose of psilocybin (mean ± SD dose: 19.7 ± 3.6 mg, administered in units of 3 mg capsules) or 20 mg of ketanserin. Data from ketanserin scans are outside the scope of the current evaluation and not presented here. After drug administration, participants completed MRI scan sessions including resting-state fMRI (see Supplementary Text for details) approximately 40, 80, 130, and 300 minutes after administration. Following each scan, participants were asked, "On a scale from 0 to 10 how intense is your experience right now" to measure SDI and a venous blooddraw used to quantify PPL (see Supplementary Text for details). After each resting-state fMRI scan, participants were asked if they had fallen asleep (no participants reported doing so). Occ 2A , i.e., occupancy of psilocybin at the 5-HT2A receptor is closely related to PPL and SDI 2 . Here we applied the previously reported parameter estimates relating PPL to occupancy based on the Hill-Langmuir equation: where Occ max refers to the maximum measurable occupancy, C p refers to the measured concentration of the ligand in plasma (i.e., PPL), and EC 50 refers to the concentration in plasma at which occupancy is equal to 50% of Occ max (fixed parameters used to compute Occ 2A : EC 50 = 1.95 µg/L and Occ max = 76.6%).
PRE-PROCESSING
Pre-processing and denoising was uniform across all entropy metrics despite differences in the pipelines of the original publications. Our pipeline included slice-timing correction (where applicable), unwarping, realignment, co-registration of structural scans to functional data, segmentation, normalisation, and smoothing in SPM12. Two MR-scanners were used to acquire the data, and some functional data were temporally downsampled so that the sampling frequency was consistent across scan sessions. Denoising in CONN v.19b 76 included linear detrending, aCompCor, 12-motion (three translations, three rotations and their first derivatives) and artefact-flagged volume regression (z > 4 SDs or motion > 2 mm using ART), band-pass filtering (0.008-0.09 Hz) and parcellation. Cerebellar ROIs were removed from included atlases as they were not consistently within the field of view. See Supplementary Text for more details.
ENTROPY QUANTIFICATIONS
All below entropy quantifications are replications of previously reported methods. Where a metric has specific hyperparameters, we ensured that they were applied as in the previous report via communications with the original authors.
ENTROPY OF STATIC CONNECTIVITY
Four studies evaluated the entropy of static connectivity given by the matrix of Pearson correlation coefficients, R, computed from 𝑁-regional time-series data, 𝑁 being the number of ROIs in the atlas used by the study.
OUT-NETWORK CONNECTIVITY DISTRIBUTION ENTROPY
Following a graph-theory framework, ROIs from the 200-region Craddock-atlaswere partitioned into communities using the Louvain modularity algorithm applied to the average connectivity matrix across scan sessions. The "Out-network Connectivity", referred to as "diversity coefficient" in the original publication and Brain Connectivity Toolbox, of an ROI was calculated for each scan session as the Shannon entropy of the distribution of connectivity estimates between a given ROI and the set of ROIs assigned to a different network.
DEGREE DISTRIBUTION ENTROPY
Degree refers to the number of non-zero elements in any given row of a thresholded matrix. ROI-specific degrees are computed based on R, the Pearson correlation matrix between ROIs, with 𝑁=105 using the Harvard-Oxford-105 atlas. Both this analysis and geodesic distribution entropy use the absolute correlation values. The thresholding for this analysis occurred in two steps. In the first step, any correlation for which the corresponding pvalue was above 0.05 was set to 0. In the second step the goal is to reach a pre-specified mean degree across rows. In order to achieve this, a threshold below which all absolute values are set to 0 is gradually increased until the mean number of non-zero elements is at the desired level. Here we applied a scan-specific threshold that produced a mean degree of 27 because this was the threshold that produced the largest effect in the original publication. This means that each scan may have a different absolute threshold value, but identical mean degree. The final entropy quantification is simply the Shannon entropy of the distribution of degrees across ROIs. We also calculated entropy for mean degrees of 1 up to the point at which for any given scan session an increase in absolute threshold did not produce an increase in mean degree, i.e., 48. This also applies to the geodesic distribution entropy described below.
GEODESIC ENTROPY
Again using absolute correlation values, the matrix was thresholded using only the mean-degree criteria and not the p-value threshold. The matrix was then binarised, setting all non-zero elements to 1. The shortest path length was then computed as the fewest edges one must traverse to go from one node to another. The Shannon entropy of the distribution of path lengths from each node to all other nodes was then calculated. Geodesic entropy was evaluated for correlation coefficient thresholds up to a mean degree of 53. where 𝜆 are the eigenvalues of the scaled correlation matrix ρ=R/N. The von Neumann entropy may also be defined as S(ρ) = -tr(ρlogρ), where log represents the matrix logarithm.
NORMALISED GLOBAL SPATIAL COMPLEXITY
NGSC was calculated via the eigenvalues of the connectivity matrix computed by a principal component analysis (PCA) of the z-scored time by voxel data matrix. Only voxels within a cortical gray matter mask were used. The eigenvalues were then rescaled by their sum to represent a probability distribution, and the Shannon entropy calculated and rescaled by the maximal value, log(m), where m is the number of non-zero eigenvalues equal to the number of degrees of freedom of the PCA. Notably, disregarding the rescaling by log(m), NGSC is exactly equal to Von Neumann entropy except that the latter used an atlas to reduce the spatial dimensionality of data.
ENTROPY OF DYNAMIC CONNECTIVITY
Intra-network Synchrony Distribution Nine brain networks were defined according to a previous study 82 : auditory, dorsal attention, default mode, left and right frontoparietal, motor, salience, visual 1 and visual 2. For a given network, for a given time point, the variance across voxels within the network was evaluated. The Shannon entropy was then calculated on the histogram of the variance estimates over time.
MOTIF-CONNECTIVITY DISTRIBUTION
Dynamic functional brain connectivity was evaluated in four regions (10mm diameter spheres) located at bilateral hippocampi, MNI coordinates: right: (26, -21, -16), left: (-34, -22, -16), and anterior cingulate cortices, right: (4, 35, 18), left: (-2, 23 ,28) using a non-overlapping sliding window approach with varying window lengths (15-150s). In each window, the partial correlation coefficient and corresponding p-value was calculated for every region pair, controlling for the remaining regions and the motion framewise displacement time-series. These time-series were standardised before windowing. The 4 x 4 partial correlation matrix was binarised for every window, according to a corrected significance threshold p=0.0083 (i.e., 0.05/6, where 6 is the number of region pairs). A probability distribution of the frequency of each of the 64 possible graph structures was established and the Shannon entropy was calculated.
LEIDA-STATE MARKOV-RATE
Notably, this entropy metric was not applied to evaluate psychedelic effects in the original paper. Rather, the authors provided a computational framework wherein parameters were learned by optimising this entropy measure. For each scan session, Leading Eigenvector Dynamics Analysis (LEiDA) 83 was applied to the time-series of 90 AAL atlas regions 84 . The phase series was computed using the Hilbert transform and, for each time point, a phase coherence matrix was estimated based on the cosine of the difference between pairwise instantaneous phases. The phase coherence matrices were decomposed using the eigenvalue decomposition and the first eigenvector was retained for every time point. The set of eigenvectors was clustered using K-means with K = 3 states. Subsequently, the transition probability matrix was computed for each scan session. The entropy rate of the transition matrix, where, 𝑝 is the leading eigenvector of 𝑃. The final entropy measure is given as 𝑆 = ∑ 𝑆 𝑖 𝑙𝑜 ⁄ 𝑔 2 (𝐾) 𝐾 𝑗=1
DYNAMIC CONDITIONAL CORRELATION DISTRIBUTION
Regional time-series were evaluated for each of the regions described in the Shen 268 region atlas. Windowless framewise correlation coefficients were calculated for all edges using the Dynamic Conditional Correlation (DCC) toolbox. Subsequently, the probability distribution over each ROI-to-ROI DCC time-series was established, and the Shannon entropy was calculated. Each ROI was assigned to one of eight networks: default mode, fronto-parietal, medial-frontal, motor, subcortical-cerebellar, visual association, visual 1, and visual 2. Each ROI-to-ROI pair was assigned to its respective network-to-network association (e.g., motor-to-motor, default mode-to-motor) and the mean entropy of each network-to-network association was calculated. Although the original publication applies bin-width correction, they do not report an effect of bin width and we report findings using MATLAB's histcounts function, which automatically calculates bin-width. Thus, we did not implement bin-width correction.
META-STATE COMPLEXITY
Regional time-series were evaluated for each of the regions described in the Lausanne 463 region atlas. BOLD time-series across all scan sessions were clustered using K-means into K = 4 states using the Pearson correlation distance metric. The clustering procedure was repeated 200 times with random initialisations and the best repeat in terms of K-means loss was extracted. The four states were grouped into two meta-states because the clustering procedure typically produces sign-symmetric states. Each volume was assigned to meta-state 0 or 1 and the Lempel-Ziv complexity (LZ76 exhaustive algorithm) of this binary sequence was calculated. Integration/Segregation-state Distribution Regional time-series were evaluated for each region described in the Schaefer 200 region atlas, augmented with 32 subcortical regions from the Tian atlas. A sliding-window correlation analysis was performed using a window defined by convolving a rectangular window of size 44 seconds with a temporal Gaussian kernel (FWHM = 3s). The correlation matrix was established for each window (stride of 1), and the Louvain modularity algorithmwas applied to estimate the module degree z-score and participation coefficient for each region. The Louvain modularity algorithm was repeated 100 times to ensure an optimal assignment. K-means clustering with K = 2 states was applied to a cartographic profile, i.e., a two-dimensional unnormalised histogram of these measures, using the correlation distance and 500 replications. The Shannon entropy was computed on the probability distribution of state occurrences.
ENTROPY OF REGIONAL DYNAMICS
Multi-scale Sample Entropy Networks were defined using the Yeo 17-network atlas. Sample entropy is defined as the negative logarithm of the conditional probability that if two vectors with length 𝑚 (set to 2) are dissimilar below a threshold distance 𝑟 (set as 0.3), then vector pairs with length 𝑚 + 1 will also have distance below the threshold. Scales 1-5 were evaluated for each network, meaning that each time-series was split into non-overlapping windows of length (scale) 𝑠 volumes and the means of each window were concatenated to form a condensed time-series upon which sample entropy was calculated.
BOLD COMPLEXITY
Regional time-series were evaluated for each of the regions described in the Schaefer 1000 region atlas. BOLD time-series for each ROI were first Hilbert-transformed. The amplitude of the Hilbert series was then binarised around the mean amplitude for that region, i.e., assigned as 1 if greater than the mean and 0 if less. These binarised time-series were combined into an 𝑇 × 𝑁 matrix, where 𝑁 = 1000 is the number of regions and 𝑇 and is the number of time points.This matrix was collapsed into a single vector to compute 1) the Lempel-Ziv complexity over time (LZct, LZ78 algorithm) wherein regional time-series were concatenated or 2) Lempel-Ziv complexity over space (LZcs) wherein time-adjacent region-series were concatenated. LZct represents a calculation of the temporal entropy of each ROI, whereas LZcs represents a calculation of the spatial entropy at each timepoint. The original publicationreported only LZcs, but LZct is also described inwhom Varley and colleagues reference as the source of their methods.
STATISTICAL MODEL
Effects of psilocybin on brain entropy metrics were estimated using a linear mixed effects model with relevant R packages, i.e., predictmeans (v1.0.6), lme4 (v1.1.30), nlme (v3.1.157), lmerTest (v3.1.3) and LMMstar (v0.7.6). We regressed each metric against each of the three measures (PPL, SDI, or Occ 2A ) separately with a subject-specific random intercept and adjusting for motion, age, sex, and scanner. A test statistic for the association between metric and measure was obtained using the Wald statistic. To ensure adequate control of the family-wise error rate (FWER) across regions within each of the 14 metrics, (e.g., 17 networks for one time scale of multi-scale sample entropy), we calculate p FWER adjusted using the maxT test methodin a permutation framework similar to, employing 10000 permutations. As such, if observed data superseded all permutations, the p-value is reported as p < 0.0001. "Motion" reflects the framewise displacement computed using the Artifact Detection Toolbox (ART) (see Supplementary Text) and "scanner" controls for MR scanner, of which there were two. We do not adjust p-values across metrics, nor across SDI, PPL and Occ 2A ; unadjusted p-values are reported for non-regional metrics as p perm . We defined findings as statistically significant if they were associated with all three psilocybin effects, SDI, PPL and Occ 2A (collectively summarised "PsiFx") at p perm < 0.05 for non-regional metrics or p FWER < 0.05 for regional metrics. Effect sizes are reported as Pearson's correlation coefficient between the partial residuals of the entropy metrics (adjusted for covariates using the mixed-model described above) and each of PsiFx. The strength of Pearson's correlation coefficients for significant associations are described as weak (≤0.3), moderate (>0.3 and ≤0.6), or strong (>0.6) as previously defined.
MODERATING EFFECT OF SCANNER
Our data were collected on one of two MRI scanners. In order to investigate whether scanner choice had an impact on the estimated relation between PPL and entropy moderating effects of scanner were explored in separate models that included the scanner-x-PPL interaction as an additional covariate, i.e., entropy ~ PPL + FD + Age + Sex + scanner + PPL:scanner.
EFFECT OF PARCELLATION
To explore parcellation effects on outcomes, all entropy metrics were evaluated using the Schaefer 100 region atlas with 16 subcortical regions from the Tian atlas. For metrics using network definitions, the Yeo 7-network atlas was applied as a common atlas.
CORRELATION BETWEEN METRICS
Simple Pearson correlation coefficients were calculated between each whole-brain entropy metric pair using the common atlas described above. Of the two graph theory metrics requiring thresholding, the threshold producing a mean degree of 27 was used. For the motif-connectivity distribution, 15 and 100 second windows were selected to represent fast and slow dynamics, respectively. P-values were adjusted using Bonferroni correction. All scans remaining after pre-processing were used in these analyses.
EFFECT OF DENOISING PIPELINES
To explore the effect of denoising decisions on the associations between PsiFx and brain entropy metrics, analyses were repeated for six additional pre-processing pipelines. Each pre-processing pipeline was run on the data parcellated as described in the section "Effect of parcellation" i.e., 116 ROIs assigned to seven networks. Each pipeline changed one variable from the main pipeline. These were as follows: 1) adding global signal regression, 2) removing the low-pass 0.09 Hz filter (i.e., not removing high-frequency signal), 3) expanding the 12-motion regressors to include squares of the derivatives (i.e., Volterra expansion), 4) not regressing out flagged volumes 5) regressing flagged volumes with a stricter threshold (z >3 or motion>0.5mm), () applying a narrower bandpass filter (0.03-0.07 Hz).
Full Text PDF
Full Paper PDF
Create a free account to open full-text PDFs.
Study Details
- Study Typeindividual
- Populationhumans
- Characteristicsopen labelbrain measures
- Journal
- Compound
- Topics
- Authors
- APA Citation
Related Papers
References (40)
Papers cited by this study that are also in Blossom
Madsen, M. K., Fisher, P. M., Burmester, D. et al. · Neuropsychopharmacology (2019)
Nichols, D. E. · Pharmacological Reviews (2016)
Anderson, B. T., Danforth, A. L., Daroff, R. et al. · EClinicalMedicine (2020)
Bogenschutz, M. P., Ross, S., Bhatt, S. R. et al. · JAMA Psychiatry (2022)
Carhart-Harris, R. L., Giribaldi, B., Watts, R. et al. · New England Journal of Medicine (2021)
Carhart-Harris, R. L., Bolstridge, &. M., Day, C. M. J. et al. · Psychopharmacology (2017)
Goodwin, G. M., Aaronson, S. T., Alvarez, O. et al. · New England Journal of Medicine (2022)
Holze, F., Gasser, P., Müller, F. et al. · Biological Psychiatry (2023)
Johnson, M. W., Garcia-Romeu, A., Griffiths, R. R. · The American Journal of Drug and Alcohol Abuse (2016)
Griffiths, R. R., Johnson, M. W., Richards, W. A. et al. · Psychopharmacology (2011)
Show all 40 referencesShow fewer
McCulloch, D. E-W., Madsen, M. K., Jensen, P. S. et al. · Frontiers in Pharmacology (2022)
Schmid, Y., Liechti, M. E. · Psychopharmacology (2017)
McCulloch, D. E-W., Knudsen, G. M., Barrett, F. S. et al. · Neuroscience and Biobehavioral Reviews (2022)
Doss, M. K., Madden, M. B., Gaddis, A. et al. · Brain (2021)
Carhart-Harris, R. L. · Neuropharmacology (2018)
Carhart-Harris, R. L., Leech, R., Shanahan, M. et al. · Frontiers in Human Neuroscience (2014)
Kringelbach, M. L., Cruzat, J., Cabral, J. et al. · PNAS (2020)
Lebedev, A. V., L€ Ovd En, M., Rosenthal, G. et al. · Human Brain Mapping (2015)
Lebedev, A. V., Kaelen, M., L€ Ovd En, M. et al. · Human Brain Mapping (2016)
Luppi, A. I., Carhart-Harris, R. L., Roseman, L. et al. · NeuroImage (2021)
Singleton, S. P., Luppi, A. I., Carhart-Harris, R. L. et al. · Nature Communications (2022)
Varley, T. F., Carhart-Harris, R., Roseman, L. et al. · NeuroImage (2020)
Viol, A., Palhano-Fontes, F., Onias, H. et al. · Scientific Reports (2017)
Barrett, F. S., Doss, M. K., Sepeda, N. D. et al. · Scientific Reports (2020)
Doss, M. K., May, D. G., Johnson, M. W. et al. · Scientific Reports (2020)
Carhart-Harris, R. L., Muthukumaraswamy, S., Roseman, L. et al. · PNAS (2016)
Nichols, D. E., Walter, H. · Pharmacopsychiatry (2020)
Singleton, S. P., Timmermann, C., Luppi, A. I. et al. · Communications Biology (2025)
Schartner, M., Carhart-Harris, R. L., Barrett, A. B. et al. · Scientific Reports (2017)
Timmermann, C., Roseman, L., Haridas, S. et al. · PNAS (2023)
Timmermann, C., Roseman, L., Schartner, M. et al. · Scientific Reports (2019)
Pallavicini, C., Cavanna, F., Zamberlan, F. et al. · Journal of Psychopharmacology (2021)
Murray, C., Frohlich, J, Haggarty, C. J., Tare, I. et al. · Neuropsychopharmacology (2024)
Carhart-Harris, R. L., Friston, K. J. · Pharmacological Reviews (2019)
Avram, M., Rogg, H., Korda, A. et al. · Frontiers in Psychiatry (2021)
Holze, F., Becker, A. M., Kolaczynska, K. E. et al. · Clinical Pharmacology and Therapeutics (2022)
Larsen, K., Lindberg, U., Ozenne, B. et al. · Journal of Cerebral Blood Flow and Metabolism (2025)
McCulloch, D. E-W., Larsen, K., Johansen, A. et al. · MedRvix (2025)
Madsen, M. K., Stenbaek, D. S., Arvidsson, A. et al. · European Neuropsychopharmacology (2021)
Olsen, A. S., Ozenne, B., Madsen, M. K. et al. · NeuroImage (2022)