This mega-analysis integrated 11 resting-state MRI datasets on acute effects of psilocybin, lysergic acid diethylamide, mescaline, N,N-dimethyltryptamine and ayahuasca, and found that psychedelics commonly increased connectivity between higher-order and sensory brain networks. It also found altered links involving the thalamus, caudate, putamen and cerebellum, with only modest and variable reductions within networks.
Psychedelic drugs are re-emerging as promising scientific and clinical tools. However, despite a rapidly expanding literature on their therapeutic value, the neural mechanisms underlying psychedelic effects remain unclear. Resting-state functional magnetic resonance imaging studies of acute psychedelic effects, conducted independently by several research groups, have so far yielded fragmented and sometimes inconsistent findings. Here, to help facilitate greater convergence, we conducted a 'mega-analysis' integrating 11 independent resting-state functional magnetic resonance imaging datasets across five psychedelic drugs (psilocybin, lysergic acid diethylamide, mescaline, N,N-dimethyltryptamine and ayahuasca) from research groups spanning three continents and five countries. By applying a uniform preprocessing pipeline and a Bayesian hierarchical modeling framework, we discovered several common features in the induced alterations to brain function across drugs and sites. Most prominently, we identified a core signature of increased functional connectivity between transmodal (default, frontoparietal and limbic) and unimodal networks (visual and somatomotor), with subnetwork specificity. Furthermore, key subcortical regions (thalamus, caudate and putamen) and the cerebellum exhibited altered coupling with sensorimotor networks. In contrast to several single-site reports, Bayesian modeling revealed weak-to-moderate and selective reductions in within-network functional connectivity, with substantial variability across drugs and networks. Together, these findings extend past work by demonstrating that psychedelics reconfigure large-scale cortical organization while selectively engaging subcortical circuitry. This study provides the most comprehensive synthesis of psychedelic brain action to date, helping resolve inconsistencies and offering a probabilistic map of how psychedelics alter large-scale brain organization. We hereby provide a cornerstone to benchmark and shepherd future psychedelic neuroimaging research.
Papers cited by this study that are also in Blossom
Vollenweider, F. X., Preller, K. H. · Nature Reviews Neuroscience (2020)
Rucker, J., Iliff, J., Nutt, D. J. · Neuropharmacology (2017)
Nutt, D. J., Erritzoe, D., Carhart-Harris, R. L. · Cell (2020)
Andersen, K. A. A., Carhart-Harris, R. L., Nutt, D. J. et al. · Acta Psychiatrica Scandinavica (2020)
Johnson, M. W., Hendricks, P. S., Barrett, F. S. et al. · Pharmacology and Therapeutics (2019)
Girn and colleagues describe a field in which acute psychedelic effects on the brain have been studied with resting-state functional MRI, but the literature has remained fragmented, methodologically variable and sometimes inconsistent. Previous single-site studies had suggested that classic psychedelics may reduce within-network functional connectivity and increase connectivity between large-scale networks, yet the exact network-level pattern varied across datasets, drugs and analysis pipelines. As a result, it was difficult to identify which findings were genuinely reproducible across studies and which might reflect site-specific or methodological differences. The authors set out to address this problem with an international mega-analysis that pooled resting-state functional MRI data across 11 datasets from five psychedelic drugs: psilocybin, LSD, mescaline, DMT and ayahuasca. Their aim was to characterise shared and drug-specific changes in brain circuit function during the acute psychedelic state, using a uniform preprocessing pipeline and Bayesian hierarchical modelling to quantify both effect size and uncertainty. The study is framed as a consortium-based effort to provide a more convergent, probabilistic map of psychedelic brain effects and to establish a benchmark for future neuroimaging research in this area.
The researchers conducted a mega-analysis of resting-state functional MRI datasets collected during the acute effects of classic psychedelics in healthy adults. Eleven independent datasets were ultimately included, with access obtained through contact with the primary investigators of earlier studies. The extracted text indicates that the core inclusion criteria were acute resting-state fMRI under a classic psychedelic and a healthy adult sample. Most datasets were double-blind, randomised, placebo-controlled trials, usually with within-subject crossover designs and counterbalanced order; one dataset used a fixed-order design without placebo. The extracted text notes that 13 datasets were considered in total, but 11 were included in the analysis. All imaging data were processed through the same centralised fMRIPrep-based pipeline, with each site running the standardised preprocessing locally. The preprocessing included coregistration, motion correction, slice-time correction, spatial normalisation to MNI space, and generation of standard confound measures such as framewise displacement, DVARS and global signal. The authors excluded subjects and sessions with poor registration, segmentation problems or excessive motion; they removed six subjects entirely and 41 connectomes in total, leaving 519 connectomes for analysis. The final sample was reported as 267 unique participants and more than 500 scanning sessions. To reduce artefacts, the researchers applied four denoising pipelines: aCompCor with and without global signal regression (GSR), and aggressive ICA-AROMA with and without GSR. The main text focuses on the aCompCor results, while ICA-AROMA results were placed in the supplementary repository. Functional connectivity was calculated as Pearson correlations between parcels. The cortical parcellation used the Schaefer atlas grouped into Yeo 17 networks, and subcortical and cerebellar regions were extracted using standard atlases, including a dedicated thalamic parcellation. Analytically, the study first conducted descriptive comparisons of drug-placebo differences in between-network and within-network connectivity, including parcel-level measures of within-network integration and between-network integration. It then applied Bayesian hierarchical models to estimate drug effects for each network and network pair, with study- and drug-level variation explicitly modelled. The Bayesian approach generated posterior distributions rather than relying on binary significance testing, allowing the authors to quantify uncertainty and consistency across datasets and drugs.
Across the pooled datasets, the researchers analysed 267 participants and more than 500 resting-state scans after motion-related exclusions. The descriptive, non-Bayesian analyses showed a broadly similar pattern across the psychedelic conditions: increased between-network connectivity, especially between transmodal association networks such as the default mode and frontoparietal networks and unimodal or sensorimotor networks such as visual, somatomotor and dorsal attention networks. There were also prominent increases in coupling between sensorimotor networks and subcortical regions, particularly the putamen, caudate, thalamus and cerebellum, although these patterns varied somewhat depending on whether GSR was used. In contrast, within-network connectivity tended to decrease across many networks, with the strongest reductions in visual, somatomotor and some default mode subnetworks. At the drug level, psilocybin and LSD showed highly similar patterns, with widespread increases in between-network connectivity and decreases in sensorimotor connectivity. DMT produced the largest apparent response and resembled an amplified version of the psilocybin/LSD pattern, although the dataset was small. Mescaline showed a broadly similar but somewhat more selective pattern. Ayahuasca was the most idiosyncratic, with weaker or absent evidence for the same broad increase in transmodal-unimodal coupling seen for the other drugs. The extracted text notes that some mean effects were visible in matrices, but that these did not always translate into high-confidence Bayesian effects. The Bayesian hierarchical analysis sharpened the picture by identifying the most consistent cross-drug effects. The clearest and most reproducible positive effects were increases in connectivity between unimodal networks, especially visual and somatomotor systems, and transmodal networks including default mode and frontoparietal subnetworks. Strong posterior support was also found for increased coupling between the caudate and visual or somatomotor networks, and more generally between the dorsal striatum and cortical networks. By contrast, thalamic coupling was less robust: although some descriptive analyses suggested increased thalamocortical connectivity, the Bayesian posteriors were more mixed and often overlapped with zero. For within-network connectivity, the Bayesian results indicated weak-to-moderate reductions rather than widespread network disintegration. Many posterior distributions overlapped with zero, and the most reliable reductions were concentrated in visual and somatomotor subnetworks. Effects in default mode, frontoparietal, dorsal attention and salience networks were smaller and less consistent. The authors also report that GSR tended to push estimates in a more negative direction for within-network connectivity, but this was not uniform across all networks and drugs.
Girn and colleagues argue that their mega-analysis provides the most comprehensive synthesis to date of how classic psychedelics alter human brain circuit function. They interpret the main result as evidence for a reproducible cross-drug signature: psychedelics increase functional coupling between transmodal association circuits and unimodal sensorimotor circuits, and also increase coupling between these cortical systems and key subcortical regions, especially the caudate and putamen. In their view, this supports previous ideas that psychedelics relax or flatten the brain’s intrinsic hierarchical organisation, although the functional connectivity data do not establish directionality. The authors position their findings as extending earlier work while also refining it. They say the study confirms the general observation that psychedelics increase whole-brain integration, but identifies which inter-network links are most robust and shows that these effects are more specific than some prior broad claims suggested. They note that previous studies often reported widespread within-network disintegration, whereas their Bayesian analysis found only limited and selective evidence for such reductions. They also highlight the value of subnetwork and subcortical analyses, which revealed distinctions that were previously obscured at coarser scales. They discuss several drug-specific patterns. LSD and psilocybin were described as nearly identical at the network level, mescaline as broadly similar but more selective, DMT as showing the strongest perturbations in raw terms, and ayahuasca as the most idiosyncratic. However, they caution that interpretation of drug-specific differences is limited because the DMT and ayahuasca datasets were small and because the studies differed in dose, timing, route of administration and pharmacology. The authors acknowledge important limitations. These include variation in scanner field strength, voxel size and repetition time across datasets, residual site-specific noise and confounding, head motion under psychedelic conditions, and heterogeneity in study design. They also note that most datasets were double-blind randomised controlled crossover studies, but one lacked placebo and another used a fixed-order design. They emphasise that blinding can be difficult with psychedelics, though they believe expectancy effects are less likely to drive resting-state BOLD findings than subjective or clinical outcomes. Overall, they argue that the consortium-based Bayesian approach helps reconcile fragmented literature and provides a more reliable foundation for future multisite, standardised psychedelic neuroimaging studies.
By carrying out a systematic mega-analysis pooled across 11 independent rsfMRI datasets, we were able to delineate the acute effects of psychedelics on human brain function. Results, encompassing 267 unique participants and more than 500 individual brain scanning sessions, painted a sharp picture of drug-mediated alterations in FC across brain regions and networks. For the sake of completeness and transparency, we present results for the aCompCor denoising pipeline with and without global signal regression (GSR). Results with the independent component analysis (ICA) with automatic removal of motion artifacts (ICA-AROMA) pipeline (with and without GSR) are presented in the Supplementary GitHub repository atBOLD_psychedelics_consortium.
To characterize the overall nature of psychedelic-induced changes in FC in a synoptic exploration, we first examined mean drug-placebo After a half-a-century-long 'psychedelic research winter', psychedelic drugs have resurfaced as drivers of scientific insight and clinical innovation in mental health. Characterized by their ability to induce wide-ranging changes to conscious experience, the so-called 'classic psychedelics,' which exert their primary effects through the shared mechanism of 5-HT2A receptor agonism, include psilocybin (the active compound in magic mushrooms), lysergic acid diethylamide (LSD), mescaline (the psychoactive alkaloid in certain cacti), N,N-dimethyltryptamine (DMT; found in many plant species and is the central ingredient of the psychedelic brew ayahuasca). These compounds have demonstrated strong therapeutic potential, with a combined total of over a dozen randomized, placebo-controlled clinical trials reporting efficacy in treating depression, end-of-life distress, generalized anxiety disorder, tobacco addiction and alcoholism. Reflecting this rapid expansion of clinical interest, a current search for 'psychedelic therapy' on clinicaltrials.gov returns more than 400 active trials exploring these compounds as a potential treatment for psychiatric and neurological conditions. Ongoing trials span all four of the classic psychedelic classes which, despite the presence of pharmacological differences between them, are assumed implicitly to overlap in key mechanisms. The growing promise of psychedelic-assisted therapy necessitates a deeper understanding of the neurobiological mechanisms underlying their effects, including the possibility of mechanistic overlap across distinct drug agents. Psychedelic neuroscience research has accelerated rapidly over the past decade, spanning from cellular-molecular investigations of neuronal morphology and plasticityto investigations of large-scale functional networks in humans. The latter investigations have relied primarily on (task-free) resting-state functional magnetic resonance imaging (rsfMRI)-a method measuring the spontaneous coordinated activity among brain regions. This approach has begun to shed light on the acute changes in neural activity and connectivity induced by psychedelic compounds and their relation to psychological effects. In addition to providing important mechanistic data to help inform regulatory decisions, progress in this area is promising for the development of precision medicine approaches for psychedelic therapy, including personalized treatment approaches and more reliable early-stage prognostics. Progress also has potential to shed light on brain circuits relevant to mental health, as well as the neurobiological and psychological sequelae of brain circuit perturbation, which can help propel drug-based and neuromodulatory therapeutic innovation. Despite compounding research interest, psychedelic rsfMRI studies in humans have so far been scattered into isolated efforts from separate research laboratories. These independent efforts have examined the changes in functional connectivity (FC; the covariation in two or more brain regions over time) induced by all of the classic psychedelics. Findings from this work suggest that psychedelics decrease FC within and increase FC between most large-scale cortical networks-a finding first observed by Roseman and colleagueswith psilocybin, and which has since been reported for additional drugs and datasets. However, the specific network-level changes underlying these broad connectivity shifts have shown considerable heterogeneity across studies. For example, a direct comparison of three existing datasets (two LSD, one psilocybin) failed to identify a single between-network FC increase consistently observed across all datasets. Another notable inconsistency arises from 'global FC' analyses, which compute the mean whole-brain FC of individual regions to assess their overall integration with the rest of the brain. Strikingly, two independent research groups applying this approach to psychedelic neuroimaging data reported topographically opposite effect: one group found increased global FC in transmodal cortex and decreases in unimodal cortex across psilocybin, DMT and LSD datasets), whereas the other reported a nearly inverse pattern for both LSD and psilocybin. Accordingly, it is difficult to extract consistently confirmable conclusions from the current literature. This ambiguity probably stems, at least in part, from the methodological and analytical variability ('researcher's degreesdifferences across the 11 eligible datasets. Figurepresents these effects, averaged across cortical networks, subcortical regions and cerebellar regions. Region-wise drug-induced changes in cortical FC are shown in Supplementary Fig.. Note that these results are meant to offer a descriptive account of the observed mean differences and are before any formal modeling. Between-network FC changes. Figuredisplays the mean druginduced changes in between-network FC, averaged across cortical networks and key subcortical regions. Both with and without GSR (lower and upper triangle, respectively; Fig.), increases in between-network FC were strongest between transmodal association networks, such as the default network (DN) (DN A ) and frontoparietal network (FPN) (FPN A and FPN B ) and unimodal/heteromodal sensory networks, such as visual (VIS) networks (VIS A and VIS B ), somatomotor (SMN) networks (SMN A and SMN B ) and dorsal attention networks (DAN) (DAN A and DAN B ). Prominent increases in FC with sensorimotor networks were also evident for subcortical regions-most notably, the putamen (PUT) and caudate (CAU) without GSR, and PUT, CAU, the thalamus (ventral thalamic nuclei (vTHA) and dorsal thalamic nuclei (dTHA)) and cerebellum (CEREB) with GSR. Furthermore, several reductions in FC in between-network emerged. With GSR, this was observed predominantly between VIS and SMN networks, with scattered decreases elsewhere. Without GSR, reductions in FC between sensory networks were also observed, as well as between limbic networks and association networks, and and outliers (1.5× interquartile range as whiskers). c,d, Mean drug-placebo difference in the mean FC of each parcel with all parcels within its Yeo-Schaefer network ('Within-network integration') and the mean FC of each parcel with all parcels outside of its network ('Between-network integration'). Red and blue indicate higher and lower FC change, respectively. Kernel density plot display the distribution of within-(blue) and between-(green) network FC changes. Zero is marked with a vertical black line. Results are shown without GSR (c) and with GSR (d). Apparent group-level trends may not correspond to high-confidence Bayesian effects (reported in 'Bayesian posterior inference' below), which reflect both effect size and consistency across subjects and datasets. between the amygdala (AMY) and hippocampus (HIP) and most cortical networks. Within-network FC changes. Figuredisplays within-network FC changes across cortical networks, with results separately shown for pipelines with (bottom) and without (top) GSR. For both pipelines, all networks showed a mean decrease in within-network FC, with the single exception of limbic (LIM) A (LIM A ) with GSR. The strongest reductions were observed in VIS (VIS A and VIS B ), SMN (SMN A and SMN B ) and DN C networks. Subcortical results (Supplementary Fig.) revealed decreased integration within all regions (that is, across the subparcels that comprise each structure), both with and without GSR. Parcel-wise within-and between-network integration. Figure,d display the region-wise mapping of mean parcel-wise within-network (left) and between-network (right) integration for each denoising pipeline. Both with and without GSR, reductions in within-network FC were observed across the cortex and subcortex-with frontal, temporal and left-lateralized lateral parietal regions showing stronger reductions with GSR. With regard to between-network FC, with GSR, cortical and subcortical regions both showed increases in predominantly between-network FC. Strongest increases were observed in lateral parietal and frontal regions, as well as along the cingulate gyrus. Weak decreases were observed in medial temporal and orbitofrontal regions, as well as the temporal sulcus. Without GSR, increases in between-network FC were more circumscribed and observed predominantly in lateral parietal and frontal regions, as well as dorsomedial prefrontal cortex. In the subcortex, increases were predominantly in the CEREB as well as the mostly right-lateralized CAU and PUT. Prominent reductions in between-network FC were observed in medial temporal and orbitofrontal regions. Minimal hemisphere-specific effects in absolute global FC were observed, as shown for cortical networks in Supplementary Fig.and for subcortical regions in Supplementary Fig..
Next we examined the drug-induced changes in mean FC within-and between-network for each drug separately, with and without GSR (Fig.). Region-wise drug-induced changes in cortical FC are shown in Supplementary Fig.. Psilocybin and LSD. Psilocybin (six datasets; n = 106 total) and LSD (four datasets; n = 119 total) revealed FC changes highly similar to each other and to the all-drugs results. With and without GSR both compounds showed widespread increases in between-network FC, particularly among transmodal association networks (FPN A and FPN B , DN A ), as well as LIM A and unimodal/heteromodal sensorimotor networks (VIS, SMN, DAN and salience (SAL)). Increased FC between subcortical regions and sensorimotor networks was also observed and was particularly prominent for LSD. Decreased FC between sensorimotor regions was also observed for both drugs. For LSD, GSR led to further decreases in between-network FC, involving LIM networks as well as temporoparietal network (TPar), HIP, AMY and nucleus accumbens (NAc). N,N-dimethyltryptamine. DMT (one dataset; n = 16) exhibited the largest drug response effect of the drugs and featured a pattern largely resembling an amplified version of the psilocybin, LSD and all-drugs average. Of note, frontoparietal (FPN A , FPN B , FPN C ) and DN (DN A , DN B ) networks, as well as CAU, vTHA and dTHA, showed particularly pronounced increases in FC with unimodal/heteromodal sensorimotor networks (VIS, SMN, DAN, SAL). DMT also induced strong FC reductions within and between VIS (VIS A and VIS B ), SMN (SMN A and SMN B ) and DAN A networks, as well as between globus pallidus (GP), PUT, CAU, vTHA, dTHA and CEREB. LIM networks showed strong reductions in FC with the rest of the brain after GSR. Mescaline. Mescaline (one dataset; n = 31) exhibited a pattern that moderately resembled that of psilocybin, LSD, DMT and all drugs. Without GSR, we observed widespread increases in between-network FC, with the most prominent effects occurring between SAL B , LIM A , FPN A , FPN networks (FPN A , FPN B , FPN C ), DN A , DN B and sensorimotor networks (VIS, SMN and DAN), broadly mirroring effects seen in psilocybin and LSD. Also, similar to psilocybin and LSD, strong increases between each of the NAc, GP and dTHA, and cortical regions were observed. With GSR, FC increases were predominantly attenuated, with pronounced changes persisting between sensorimotor and association networks, as well as involving PUT and CAU. Ayahuasca. Ayahuasca (one dataset; n = 9) exhibited a relatively idiosyncratic pattern of FC changes. In contrast to the other examined drugs, results without GSR revealed reduced FC between most networks and did not reveal a pattern of increased FC between unimodal and transmodal networks. With GSR, prominent reductions in FC were observed between unimodal/heteromodal sensorimotor networks (VIS, SMN, DAN, SAL), as well as between LIM A , DN C , HIP and AMY and sensorimotor networks. Increases were scattered, with prominent decreases between several subcortical regions (PUT, CAU, vTHA, dTHA) and sensorimotor networks. Without GSR, prominent decreases were also observed involving LIM networks, SAL B , FPN C and other scattered network pairs.
After qualitatively assessing the broad topography of psychedelicinduced changes in large-scale network connectivity, we next applied Bayesian hierarchical inference (Methods) to formally quantify the strength and uncertainty of these effects across drugs and studies. Although some mean effects seemed prominent in average FC matrices, not all reached high-confidence levels in the Bayesian analyses. This reflects the model's sensitivity to intersubject variability and noise, highlighting only those effects with robust, consistent support across drugs and datasets. Unlike traditional frequentist approaches aiming at categorical yes/no answers on effects, our Bayesian framework allows for graded probabilistic inference by answering the question: 'How sure are we that these two brain regions/networks change in FC as a result of psychedelic drug effects?' We constructed Bayesian models for each network (intra-network FC) and network pair (between-network FC), deriving full posterior distributions for each drug's estimated effects. Robust cross-psychedelic effects are those that have distributions that (1) are distant from zero, (2) show relatively low dispersion and (3) show overlap/consistency in divergence from zero across drugs. We avoid strict thresholds as these would be, by necessity, arbitrary and in opposition to Bayesian thinking, and instead allow readers to make their own judgments on the basis of the posterior distributions provided. Bayesian posterior distributions for the most prominent Apparent group-level trends may not correspond to high-confidence Bayesian effects (reported in 'Bayesian posterior inference' below), which reflect both effect size and consistency across subjects and datasets. between-network FC changes (based on our noninferential descriptive analyses) are shown in Fig.(additional results available in the Supplementary GitHub Repository atpsychedelics_consortium). Across network pairs, posterior inference revealed a consistent pattern of increased between-network coupling, with effect magnitudes and certainty varying by drug, network pair and preprocessing choice. Distributions were most commonly overlapping for LSD and psilocybin and, to a lesser extent, mescaline. Distributions for these three drugs also demonstrate the least dispersion, explained, probably at least partly, by larger sample sizes in these drug groups. Less dispersion or tighter posterior distributions increases confidence that similarities between LSD and psilocybin may be meaningful and reliable. Ayahuasca and DMT tended to exhibit the least certainty in effects, probably owing to the converse-namely, single studies and relatively small sample sizes. Without GSR, the strongest and most consistently positive posterior shifts were observed for CAU coupling with unimodal networks, including CAU-VIS A , CAU-VIS B , CAU-SMN A and CAU-SMN B . Robust positive shifts were also evident for cross-network coupling between VIS and transmodal subnetworks, most prominently DN A -VIS B , DN B -VIS A , FPN B -VIS A and FPN A -VIS B , as well as for SMN-frontoparietal coupling, particularly FPN B -SMN A . In contrast, THA coupling showed posterior distributions that were consistently negatively skewed across all VIS and SMN subnetworks. With GSR, posterior distributions across these between-network pairings were generally shifted in a more positive direction, although the magnitude of this shift varied across drugs and network pairs. Regarding within-network FC effects, Bayesian parameter distribution plots are shown in Fig.for all 17 (sub)networks and the subcortical regions with the most prominent within-network effects (based on our noninferential descriptive analyses). Across networks, posterior distributions generally indicated weak-to-moderate reductions in within-network FC, with substantial variability in both effect magnitude and uncertainty across drugs and preprocessing pipelines. Although many distributions overlapped with zero, several networks exhibited consistent posterior shifts toward decreased within-network FC. Drug response effects were most reliable for psilocybin and LSD, which showed relatively narrow posterior distributions and consistent
Posterior distribution plots from Bayesian modeling for the largest nonzero effects in drug-mediated changes in between-network FC for each drug. These posterior distributions provide principled, probabilistically grounded estimates of both the uncertainty (as indicated by distribution width) and magnitude (as indicated by the peak location on the x axis) of drug-induced effects for each modeled drug class. The x axis represents change in Pearson's product-moment correlation coefficient (r) between conditions (drug-placebo) and the y axis represents Bayesian probability in unspecified units. directionality across several networks. Mescaline exhibited broadly similar but more variable effects, whereas DMT and ayahuasca showed wider, more dispersed posteriors, reflecting greater uncertainty, probably attributable to smaller sample sizes. Without GSR, posterior distributions were shifted toward negative values across sensorimotor networks (VIS A , VIS B , SMN A , SMN B ) and select hetero/transmodal networks (SAL A , DN C , TPar). On the whole, transmodal networks, including default and frontoparietal subnetworks, showed smaller and less consistent shifts, with posteriors frequently overlapping zero. With GSR, posterior distributions tended to be shifted further in the negative direction. However, this shift was nonuniform across drugs and networks, and several transmodal networks continued to exhibit substantial posterior overlap with zero. Subcortical regions, including the caudate, putamen and thalamus, also exhibited modest negative shifts in within-region coupling.
There has been an upsurge of research activity on the neural underpinnings of the acute psychedelic experience in humans over the last decade. These efforts have produced a variety of relatively scattered findings, spanning several independent research groups across the globe. Owing to a variety of factors, this expanding literature is marked by heterogeneous and sometimes discrepant findings. Although much preclinical and wet laboratory science is challenging to homogenize and formally integrate, the brain imaging community has the key advantage of agreement on data formats, acquisition procedures and reference atlases. This study represents the most comprehensive synthesis to date of how classic psychedelics modulate human brain circuit function. To reach this goal, we created an international consortium-the 'BOLD Psychedelic Consortium'-that aggregated and processedpsychedelic-mediated changes uniformly in rsfMRI functional coupling profiles from 11 datasets spanning four psychedelic drugs, totaling 267 participants and over 500 connectomes. By implementing a fully Bayesian hierarchical modeling framework, we provide probabilistically grounded evidence for both shared and drug-specific alterations in FC within and between large-scale cortical networks and subcortical regions. Our findings reaffirm previous observations of psychedelic-induced enhanced whole-brain integration yet refine them by highlighting the changes in inter-network coupling links that exhibit the greatest robustness (defined here as high posterior confidence and low dispersion) across studies and by revealing effect specificity at the subnetwork and subcortical levels. We further note that, although certain patterns were apparent in mean FC matrices-such as reduced within-network coupling and increased THA-unimodal coupling-they did not always yield high-confidence posteriors. This reflects a core strength of the Bayesian approach: it considers jointly both the magnitude and uncertainty of an effect, discounting patterns with high between-subject variability or inconsistent replication across studies. In doing so, it avoids overinterpretation of apparent group-level trends and foregrounds the most robust and generalizable findings. As our core conclusion, we found that psychedelics most robustly increase functional integration between select pairs of transmodal and unimodal subnetworks, as well as between key subcortical regions (PUT, CAU) and both unimodal and transmodal cortical areas. These conclusions are supported by posterior distributions showing high confidence (narrow dispersion) and consistent divergence from zero across compounds. In addition, our results challenge previous claims of widespread within-network disintegration, as Bayesian hierarchical modeling revealed limited and selective evidence for reductions in within-network FC, with few effects showing consistent divergence from zero. Instead, psychedelics seem to selectively reconfigure large-scale network-network interactions while modulating subcortical-cortical connectivity. Comparable analyses using ICA-AROMA-based denoising pipelines are reported in the Supplementary GitHub repository atedelics_consortium. These showed broadly similar patterns of FC change, although with reduced effect magnitudes and notable sensitivity to the inclusion of GSR-highlighting the trade-offs between denoising stringency and preservation of neural signal. In sum, despite substantial methodological heterogeneity across studies, a core set of reproducible, cross-drug network changes emerged, pointing toward conserved neurobiological mechanisms across psychedelic drugs. As the central finding of the present work, our analysis identified a brain signature of increased functional coupling between transmodal association circuits (FPN A , FPN B , FPN C , DN A , DN B ) and unimodal/heteromodal sensorimotor circuits (VIS A , VIS B , SMN A , SMN B , DAN A , DAN B ). Notably, our Bayesian modeling revealed many high-confidence posterior probability distributions supporting these increases, particularly for LSD, psilocybin and mescaline, with some effects observed in amplified form for DMT. The stark reaction of transmodal systems under psychedelic effects is consistent with findings from high-resolution in vivo positron emission tomography imaging that have shown that the primary mechanistic target of psychedelics, the 5-HT2A receptor, is expressed most densely in transmodal (as well as VIS) cortices, and that these regions are among those showing statistically increased 5-HT2A occupancy after psilocybin administration. Given that the unimodal-transmodal axis reflects a central cortical hierarchy separating sensorimotor from abstract cognitive processing, the observed increase in their coupling suggests a flattening of this intrinsic processing hierarchy (cf. refs. 14,21). Although our FC data cannot determine directionality, these effects may arise from either increased top-down transmodal constraint or increased bottom-up influence from unimodal regions, both of which are plausible and supported by previous work. The former interpretation aligns with known patterns of task-dependent coupling between frontoparietal and sensory networks in goal-directed behavior. The latter is consistent with evidence that 5-HT2A receptor activation disrupts transmodal cortical synchrony, potentially relaxing hierarchical precision weighting, as posited by the 'relaxed beliefs under psychedelics' model of psychedelic effectsand also consistent with the cortico-striato-thalamocortical gating model. The pattern of increased unimodal-transmodal (and transmodal-transmodal) coupling under psychedelics, first shown by Roseman et al., has been reported in several previous studies of LSD, psilocybin and DMT-all of which were included in the present analyses. However, the precise spatial topography of this effect has varied across studies. Here we expanded upon these past studies and identified, with fine-grained subnetwork specificity and probabilistic likelihood, the specific unimodal transmodal network pairs that exhibited the most reliable increases in connectivity across drugs and datasets. In particular, we found that regions within DN and FPN (sub)networks may show the most reliable increases in coupling with unimodal networks, reinforcing the notion that psychedelics may flatten lines of communication in the brain's intrinsic neural processing hierarchy. Bayesian modeling also identified a robust pattern of intensified coupling between unimodal cortical networks and the dorsal striatum, particularly the CAU and PUT. Both the CAU and PUT receive dense convergent input from visual, motor and association cortices and play a central role in action selection, sensorimotor integration and the contextual modulation of perception and behavior. Increased striatal coupling with unimodal cortex may therefore reflect altered weighting of corticostriatal interactions that link sensory input to motor and behavioral output, a process that is plausibly engaged during psychedelic states. These findings loosely align with cortico-striato-thalamocortical accounts of psychedelic action, which implicate changes in basal ganglia signaling in shaping large-scale information flow. However, support for thalamic involvement was comparatively weaker in the present data. Although qualitative mean effects suggested increased coupling between thalamic regions (vTHA and dTHA) and unimodal cortex, these effects did not emerge as reliable in the Bayesian analyses. Instead, thalamocortical posteriors largely overlapped with zero, with a tendency toward weak negative shifts without GSR and moderately positive shifts with GSR. This pattern contrasts with previous reports emphasizing psychedelic-induced thalamic disinhibition and increased afferent flow to cortex, although important methodological differences-for example, analytic framework, parcellation strategy and treatment of global signal-limit direct comparability. Taken together, our findings indicate that striatal involvement represents the clearest and most reproducible subcortical feature of the observed connectivity changes. Regarding within-network functional integration, our Bayesian findings point to a selective and modest set of effects, rather than widespread psychedelic-induced network 'disintegration' as has been reported previously (for example, refs. 13,14,27). Across networks, Bayesian posterior distributions revealed weak-to-moderate reductions in within-network FC, with substantial variability in effect magnitude and uncertainty across drugs, networks and preprocessing pipelines. Of note, many posterior distributions overlapped with zero, indicating that within-network effects were neither uniform nor universally robust. Within-network reductions were most consistent in sensorimotor (VIS and SMN) subnetworks, both with and without GSR. In contrast, hetero-and transmodal networks-including the DN, FPN, DAN and SAL-showed smaller and less reliable effects, with posterior mass frequently centered near zero. Notable exceptions include DN C and SAL A . Subcortical regions exhibited modest reductions in within-region coupling only after GSR. Broadly, the inclusion of GSR systematically shifted within-network estimates in a more negative direction; however, this effect was nonuniform across networks and drugs. An important point to note is that past work has applied primarily an ICA dual regression approach to assess within-network FC changes, whereas we applied an inter-regional FC approach. This may have influenced the observed effects (or lack thereof). Nonetheless, our finding of limited and heterogenous within-network effects with our analysis approach and Bayesian framework suggests that past work may have overestimated the robustness of this effect. Finally, our results also underscore the importance of fine-grained subnetwork analyses, which here revealed meaningful dissociations within networks that were obscured previously. Although a core cross-drug signature was evident, our Bayesian modeling also revealed distinct connectivity profiles across individual psychedelics. We focus here on drug-induced effects that were broadly consistent across GSR and no GSR pipelines. LSD and psilocybin-our largest sample-size datasets-displayed virtually identical network-level alterations, consistent with their comparable pharmacological properties and phenomenological profiles. Mescaline in turn exhibited a broadly similar pattern but with more selective enhancements in transmodal-unimodal integration and subcorticalcortical integration. DMT qualitatively exhibited the strongest network perturbations across all drugs. This trend was reflected in our Bayesian posterior distributions, though to a somewhat lesser extent than expected from raw effect size comparisons, probably due to high interindividual variability and the small sample size of the DMT dataset (n = 15). Ayahuasca, which contains DMT but also monoamine oxidase inhibitors, exhibited the most idiosyncratic connectivity profile, probably due to both its pharmacological complexity and its extremely small dataset (n = 9), which limits interpretability. Overall, although our study provides the most comprehensive evaluation to date of common psychedelic-induced connectivity changes, it remains challenging to draw strong conclusions regarding drug-specific effects. A complex mix of methodological and pharmacological factors probably contribute to the observed differences across drugs, including dosage, route of administration, receptor binding affinities (for example, 5-HT2A, 5-HT1A and others), time since drug administration and motion-related artefacts. Future studies should aim to directly compare different psychedelics under matched experimental conditions with standardized administration protocols and larger sample sizes to more precisely delineate the unique versus shared neural signatures of each compound (cf. ref. 30). Despite the advantages of pooling data across several psychedelic neuroimaging studies, our mega-analytic approach carries important limitations. First, the included datasets varied in several fMRI acquisition parameters. Most notably, scanner field strength included 1.5 T (ayahuasca), 7 T (Maastricht psilocybin) and 3 T (all remaining datasets). Voxel size and repetition time also varied across studies (2-3 mm isotropic voxels; 2-3 s repetition time). Although our uniform preprocessing pipeline was designed to mitigate these differences, such variability inevitably introduces noise that may obscure more subtle effects. Second, although pooling subjects across several studies enhances confidence in the generalizability and robustness of findings, it does not guarantee data quality. Site-specific artifacts, residual physiological noise and unmeasured confounds may still propagate through the analysis, even with harmonized preprocessing and denoising. Head motion, in particular, remains a well-known concern in resting-state fMRI and is especially relevant here, given that participants tend to move more under the influence of psychedelics. To address this, we computed correlations between individual-level changes in framewise displacement (FD) and corresponding changes in FC. As shown in the Supplementary Information, these correlations were weak-to-moderate, varied across drugs and did not resemble the spatial pattern of the observed drug effects. This suggests that our main findings are unlikely to be driven solely or systematically by motion-related artifacts. Third, as detailed in Table, studies varied in dosage, route of administration and latency of scanning relative to administration-all of which may also influence the observed effects. Although our Bayesian hierarchical modeling approach helps mitigate the impact of such heterogeneity, some residual influence is likely to persist. Finally, although our mega-analytic design was aimed at a harmonized analysis across studies, variability in study design remains. Although most datasets employed counterbalanced within-subject double-blind, randomized, controlled trials (DB-RCT) designs, one study lacked a placebo control and another used a fixed-order designboth of which may introduce confounds related to novelty, time or expectancy. Furthermore, although double-blinding aims to reduce bias, the strong psychological effects of psychedelics make true blinding difficult to maintain. Nevertheless, we believe that such expectancy and unblinding effects are less likely to substantially influence resting-state BOLD signal compared to subjective or clinical outcomes. Still, we acknowledge that design heterogeneity may contribute to residual variance in our findings. As the field advances, future collaborative efforts would benefit from prospectively harmonized multisite studies with standardized study designs, acquisition protocols, dosing regimens and participant selection criteria. Our confederated effort is an instrumental step in bringing together the fragmented landscape of psychedelic research, yielding a more integrated whole that will in turn provide more clarity on the path forward. By coalescing rsfMRI data from several research groups across North America, South America and Europe, we identified a core, cross-drug signature of enhanced connectivity between transmodal association and unimodal sensorimotor networks, alongside enhanced selectively enhanced subcortical-cortical coupling and reduced subcortical-subcortical integration. These findings support neurocognitive models of psychedelic-induced hierarchical relaxation, while adding nuances to previous work by revealing fine-grained subnetwork and subcortex-specific effects. Contrary to some past findings, our mega-analytic approach did not find consistent evidence for within-network disintegration in the psychedelic state. Ultimately, cutting across methodological heterogeneity and additional limitations, our analysis revealed a robust cross-drug neural fingerprint of psychedelic states. As psychedelic research enters a new era of clinical and neuroscientific exploration, our findings provide a critical foundation upon which future investigations can build, ensuring that the field progresses toward greater rigor, reliability and translational impact.
Any methods, additional references, Nature Portfolio reporting summaries, source data, extended data, supplementary information, acknowledgements, peer review information; details of author contributions and competing interests; and statements of data and code availability are available at.
This mega-analysis was designed and carried out as a community-wide effort that sought to integrate all available rsfMRI datasets examining the acute effects of classic serotonergic psychedelics in healthy adults. Eleven independently acquired fMRI neuroimaging datasets were ultimately included (Table). Access to these datasets was obtained by contacting the first and senior authors of the primary publications for previous acute psychedelic neuroimaging studies. Our two core inclusion criteria were: (1) rsfMRI acquired during the acute effects of a classic psychedelic and (2) a healthy adult participant sample. To our knowledge, no eligible datasets were excluded for reasons other than feasibility of data sharing. One site-Copenhagen University Hospital (principal investigator P. Fisher)-was unable to contribute due to restrictions related to general data protection regulation. As of our final outreach in August 2024, we were not aware of any additional eligible datasets that were excluded due to lack of contact or nonresponse. Of the 13 datasets, 12 were DB-RCTs that included a matched placebo condition. Most employed within-subject crossover designs with counterbalanced order of drug and placebo administration. One datasetused a fixed-order design without a placebo condition. Tablesummarizes the dosage, timing, and study design of each dataset. Full details on each dataset can be found in the original publications.
All resting-state fMRI data underwent an identical preprocessing protocol using the fMRIprep software v.22.1.1. This standardized, uniform preprocessing protocol was centrally coordinated and executed locally in an automated fashion in each laboratory where the data originated. For each of the BOLD runs found per subject, our analytical protocol performed the following preprocessing steps. First, a reference volume and its skull-stripped version were generated using a custom methodology of fMRIPrep. The BOLD reference was then coregistered to the T1w reference using bbregister (FreeSurfer), which implements boundary-based registration. Coregistration was configured with six degrees of freedom. Head motion parameters with respect to the BOLD reference (transformation matrices, and six corresponding rotation and translation parameters) were estimated before any spatiotemporal filtering using McFlirt (FSL). BOLD runs were slice-time corrected using '3dTshift' from AFNI. The BOLD timeseries were resampled onto their original, native space by applying the transforms to correct for head motion. We then resampled the ensuing voxelwise BOLD timeseries from individual subject space into standard anatomical reference space, generating preprocessed BOLD runs in the Montreal Neurological Institute (MNI) 152NLin2009cAsym 2-mm space. fMRIPrep then calculated the following confounds automatically based on the preprocessed BOLD: FD, spatial root mean square of the data after temporal differencing (DVARS)and the global signal (mean signal within a whole-brain gray matter (GM) mask). FD and DVARS were calculated for each functional run, both using their implementations in Nipype. All fMRIPrep outputs were generated by each site following fMRIPrep's standardized quality control protocols. This included verification of: (1) accurate brain extraction; (2) correct tissue segmentation into GM, white matter (WM) and cerebrospinal fluid (CSF); (3) precise coregistration between functional and anatomical images, with alignment confirmed in several planes; (4) appropriate normalization to MNI-space and (5) complete and accurate brain masking. Corresponding checks were guided by fMRIPrep-generated HTML reports and aligned with best practices from the Human Connectome Project. Each site conducted these checks before sharing preprocessed data with the central analysis team. Subjects showing poor registration, segmentation errors or mask artifacts were excluded. To automatically remove high-motion subjects based on the obtained fMRIPrep QC outputs, we followed past work and excluded subjects which featured >15% of timepoints with FD > 0.5 (refs. 13,14,47,58). A total of six subjects were fully excluded from this study (five psilocybin, one DMT), corresponding to 12 total connectomes. An additional 29 connectomes were removed due to excessive motion at the session level (13 psilocybin, 7 LSD, 1 mescaline, 1 DMT and 7 placebo). Thus, a total of 41 connectomes were removed, resulting in a total of 519 remaining for the analyses.
Following preprocessing, we applied four distinct pipelines to remove physiological, scanner-related and motion-related artifacts from the neuroimaging data, yielding four different sets of denoised brain signal volumes and, therefore, four sets of results. These pipelines included: (1 and 2) the anatomical CompCor (aCompCor) approach of Behzadi et al., with and without GSR, and (3 and 4) the 'aggressive' ICA-AROMA approach of Pruim et al., also with and without GSR. Pipelines 1 and 2 (aCompCor) are reported in the main text, whereas Pipelines 3 and 4 (ICA-AROMA) are included as a more stringent secondary approach in the Supplemental GitHub Repository at. For the aCompCor-based pipelines, with and without GSR, a set of physiological regressors were extracted to allow for latent factor component-based noise correction. Principal components were estimated after high-pass filtering the preprocessed BOLD timeseries using a discrete cosine filter (128-s cut-off). Probabilistic masks were generated in anatomical space for CSF and WM, and combined for an additional CSF + WM mask. Unlike the original implementation of aCompCor, which erodes the nuisance masks in BOLD space by contours of two voxels, we employed an anatomically informed refinement whereby a dilated GM mask, extracted from FreeSurfer's 'aseg' segmentation, was subtracted to avoid partial volume contamination. This GM-removed CSF + WM mask was resampled into BOLD space and binarized at a 0.99 threshold. Components were extracted separately from each of CSF, WM and CSF-WM combined masks, with the top k components retained such that their timeseries cumulatively explained 50% of the variance within each mask. The remaining components were discarded from further analysis. These aCompCor-derived regressors were then used to denoise the BOLD signal with and without GSR. In contrast, the ICA-AROMA-based pipelines, with and without GSR, involved applying ICA as a basis to classify and remove motion-related artifacts automatically. This was performed on MNI-space BOLD volumes that were first smoothed using a 6-mm FWHM Gaussian kernel (smoothing was applied only for ICA-based noise-versus-no-noise classification, not for final analysis). The fMRI timeseries data were decomposed into spatially independent components (ICs), which were then distinguished as signal or noise by a supervised pattern-learning algorithm using three features: (1) temporal correlation with motion parameters, (2) high-frequency spectral content and (3) spatial overlap with CSF or brain edges. ICA-AROMA includes both 'nonaggressive' and 'aggressive' variants; we implemented the aggressive variant, which regresses only the noise ICs, thereby providing a more stringent denoising procedure compared to aCompCor. In contrast, the nonaggressive approach would have retained shared variance between signal and noise components, potentially leaving residual artifacts. By employing aggressive ICA-AROMA, we aimed to maximize noise removal; thus, using this approach as a conservative approach to compare with aCompCor. ICA-AROMA was applied with and without GSR. Following denoising, all FC matrices were assessed using a standardized procedure. Each matrix was plotted using identical color scales and thresholds to ensure consistent evaluation across participants and studies. As part of our protocol, we checked and confirmed the presence of canonical resting-state network structure-specifically that within-network connectivity was visibly higher than between-network connectivity for key neural systems such as the default mode, VIS and SMN networks. Matrices were also assessed in terms of the direction andrelative magnitude of drug-induced changes and were confirmed to be consistent with those reported in the original studies for each dataset.
To facilitate comparability to other studies and aid neuroscientific interpretation, we extracted the denoised BOLD timeseries from cortical and subcortical regions using field-standard publicly available atlases. Cortical timeseries were extracted using the Schaefer et al.local-global parcellation (100 region resolution and were grouped into large-scale subnetworks based on the 17 network assignments of Yeo et al.). Subcortical timeseries were extracted using the Tian et al.subcortical parcellation. The 'S2' parcellation was used for all nonthalamus subcortical regions (24 parcels total) and the 'S3' parcellation was used for the thalamus (14 parcels total), given the special interest in the latter in psychedelic research. CEREB timeseries were extracted using the Buckner et al.cerebellar parcellation (17 parcels in total).
Interregional FC was calculated as the Pearson's product-moment correlation r between all parcels. Between-network FC corresponds to FC between parcels of different networks as defined based on the 17-network assignments of Yeo et al., whereas within-network FC corresponds to FC between parcels of the same network. As visualized in Fig.,d below, we additionally computed parcel-wise measures that we refer to as 'within-network integration' and 'between-network integration', following Siegel et al.. Defined here, 'within-network integration' corresponds to the mean FC of each parcel with all parcels within its Yeo-Schaefer network, whereas 'between-network integration' corresponds to the mean FC of each parcel with all parcels outside of its network. In this sense, it is a variant of previous global FC or 'global brain connectivity' approaches (for example, refs. 20,22).
We carried out a Bayesian mega-analysis to fully quantify the uncertainty of the different sources of variation at play, as a natural choice of method. This step of our analysis workflow enabled the discovery and probabilistic characterization of the common neural features underlying psychedelic drug effects as captured by human brain imaging.
Most previous neuroscience studies on psychedelic drugs attempted to draw sharp boundaries for regional drug effects using classical null hypothesis testing and P value thresholds. To overcome these limitations, the present work adopts Bayesian principles as a formal estimation framework to quantify the degree of change in functional coupling links under the psychedelic state. In contrast to frequentist post hoc estimations, the Bayesian regime provides a direct quantification of uncertainty around model parameters by appropriately handling all considered sources of variation. This enables more careful estimation of even very subtle drug effects. As an increasing trend, Bayesian mega-analysis has been used in medicine in high-stake areas such as to re-evaluate drug effects across several RCTs. Bayesian hierarchical regression provides several key advantages relevant to the present aims: (1) it enables principled treatment of uncertainty by deriving full Bayesian posterior distributions for all measures, (2) it allows for continuous statements with explicitly quantified degrees of confidence (for example, 'How sure are we that these two drugs lead to similar changes in FC?') and (3) it does not require correction for multiple comparisons when models are independently specified, which is a key advantage as we can run one separate Bayesian model per brain feature. More broadly, although frequentist models are designed to support binary decision-making under controlled assumptions (for example, rejecting a null hypothesis), Bayesian models are better suited for probabilistic reasoning about effect magnitudes and uncertainty. Rather than making dichotomous claims of presence or absence, our approach provides a graded, uncertainty-aware estimate of drug-related FC changes across diverse datasets. This distinction in modeling goals is especially relevant in the context of heterogeneous multisite neuroimaging data, where quantifying between-study variation is as important as identifying consistent effects. Interpretation. Our analytical approach aimed to answer a distinct question: 'How certain are we that a particular psychedelic drug in a certain brain region leads to a drug-induced change in functional coupling, and how strong is this effect?' Accordingly, our models estimate the complete shape of the effect uncertainty of regional drug effects in the form of principled Bayesian posterior distributions. That is, our study did not ask 'Is there a strict categorical difference in inter-regional or inter-network FC between drug and placebo states?' Rather, we sought to directly quantify the population uncertainty distributions of functional coupling effects in relation to psychedelic drug influence, rather than focusing exclusively on mean FC differences. This modeling approach thus allows for probabilistic insight into the presence and strength of psychedelic effects across brain circuits, as well as their variability across drugs and studies. A fully specified generative model of the neural responses further opens the possibility to identify graded or overlapping shifts in brain network organization, consistent with the idea that drug-induced mental states may differ in degree rather than in kind. Model specification. The random-effects probability model with parameters that vary by study and by drug, using a subject-wise MRI-derived brain measure as outcome and the experimental condition as input variable, had the following form: Full Bayesian model specification. where α denotes the intercept parameter corresponding to each study j out of 11 total studies (datasets), bundled by the hyperprior hyper α , β denotes the slope parameter explaining the presence of the psychedelic state in connectivity differences corresponding to drug k out of the 5 overall drugs, bundled by the across-drug hyperprior hyper β , drug_condition denotes whether the observed resting-state scan was acquired during the placebo or drug-influenced state and y conn_link captures a region-region functional coupling strength, initially computed by Pearson's correlation across the duration of a given resting-state brain scan in a particular participant, and then rescaled by Fisher z transformation (arctanh()) as outcome for this model estimation. As such, this model was inputted as many datapoints as we had resting-state scans (four-dimensional timeseries) from all 11 studies.
Create a free account to open full-text PDFs.
Kwan, A. C., Olson, D. E., Preller, K. H. et al. · Nature Medicine (2022)
Jaster, A. M., González-Maeso, J. · Molecular Psychiatry (2023)
McCulloch, D. E-W., Knudsen, G. M., Barrett, F. S. et al. · Neuroscience and Biobehavioral Reviews (2022)
Girn, M., Rosas, F. E., Daws, R. E. et al. · Trends in Cognitive Sciences (2023)
Linguiti, S., Vogel, J. W., Sydnor, V. J. et al. · Neuroscience and Biobehavioral Reviews (2023)
Roseman, L., Leech, R., Feilding, A. et al. · Frontiers in Human Neuroscience (2014)
Carhart-Harris, R. L., Muthukumaraswamy, S., Roseman, L. et al. · PNAS (2016)
Timmermann, C., Roseman, L., Haridas, S. et al. · PNAS (2023)
Müller, F., Dolder, P. C., Schmidt, A. et al. · NeuroImage (2018)
Stoliker, D., Novelli, L., Vollenweider, F. X. et al. · Biological Psychiatry (2024)
Madsen, M. K., Stenbaek, D. S., Arvidsson, A. et al. · European Neuropsychopharmacology (2021)
Tagliazucchi, E., Roseman, L., Kaelen, M. et al. · Current Biology (2016)
Girn, M., Roseman, L., Bernhardt, B. et al. · NeuroImage (2022)
Preller, K. H., Burt, J. B., Adkinson, B. et al. · Biological Psychiatry (2020)
Preller, K. H., Burt, J. B., Adkinson, B. et al. · eLife (2018)
Palhano-Fontes, F., Andrade, K. C., Tófoli, L.F. et al. · PLOS ONE (2015)
Madsen, M. K., Fisher, P. M., Burmester, D. et al. · Neuropsychopharmacology (2019)
Barrett, F. S., Zhou, Y., Carbonaro, T. M. et al. · Frontiers in Neuroergonomics (2022)
Muthukumaraswamy, S. D., Carhart-Harris, R. L., Moran, R. J. et al. · Journal of Neuroscience (2013)
Carhart-Harris, R. L., Friston, K. J. · Pharmacological Reviews (2019)
Doss, M. K., Madden, M. B., Gaddis, A. et al. · Brain (2021)
Preller, K. H., Razi, A., Zeidman, P. et al. · PNAS (2019)
Avram, M., Müller, F., Rogg, H. et al. · Biological Psychiatry (2022)
Carhart-Harris, R. L., Erritzoe, D., Williams, T. et al. · PNAS (2012)
Barrett, F. S., Doss, M. K., Sepeda, N. D. et al. · Scientific Reports (2020)
Ley, L., Holze, F., Arikci, D. et al. · Neuropsychopharmacology (2023)