Subacute effects of the psychedelic ayahuasca on the salience and default mode networks
=======================================================================================

* Lorenzo Pasquini
* Fernanda Palhano-Fontes
* Draulio B. Araujo

## Abstract

**Background** Neuroimaging studies have just begun to explore the acute effects of psychedelics on large-scale brain networks’ functional organization. Even less is known on the neural correlates of subacute effects taking place days after the psychedelic experience. This study explores the subacute changes of primary sensory brain networks and networks supporting higher-order affective and self-referential functions 24h after a single session with the psychedelic ayahuasca.

**Methods** We leveraged task-free functional magnetic resonance imaging data one day before and one day after a randomized placebo-controlled trial exploring the effects of ayahuasca in naïve healthy participants (21 placebo/22 ayahuasca). We derived intra- and inter-network functional connectivity of the salience, default mode, visual, and sensorimotor networks, and assessed post-session connectivity changes between the ayahuasca and placebo groups. Connectivity changes were associated with Hallucinogen Rating Scale scores assessed during the acute effects.

**Results** Our findings revealed increased anterior cingulate cortex connectivity within the salience network, decreased posterior cingulate cortex connectivity within the default mode network, and increased connectivity between the salience and default mode networks one day after the session in the ayahuasca group compared to placebo. Connectivity of primary sensory networks did not differ between-groups. Salience network connectivity increases correlated with altered somesthesia scores, decreased default mode network connectivity correlated with altered volition scores, and increased salience-default mode network connectivity correlated with altered affect scores.

**Conclusion** These findings provide preliminary evidence for subacute functional changes induced by the psychedelic ayahuasca on higher-order cognitive brain networks that support interoceptive, affective, and self-referential functions.

Keywords
*   Ayahuasca
*   default mode network
*   functional connectivity
*   interoception
*   salience network

## Introduction

Classic psychedelics belong to a family of substances that lead to altered states of consciousness, with a broad range of effects on sensory, cognitive, autonomic, interoceptive, self-referential, and emotional systems (Nichols, 2016). Assisted sessions with psychedelic substances such as ayahuasca and psilocybin have recently gained prominent use for the treatment of affective disorders (Carhart-Harris, Bolstridge, et al., 2016; Griffiths et al., 2016; Grob et al., 2011; Osório et al., 2015; Palhano-Fontes et al., 2019; Ross et al., 2016; Sanches et al., 2016).

Ayahuasca is an indigenous preparation that contains leaves of *Psychotria viridis* and the bark of *Banisteriopsis caapi*, two plants vastly found in the Amazon basin (Mckenna et al., 1984). The admixture contains the psychedelic tryptamine N,N-dimethyltryptamine (N,N-DMT) and monoamine oxidase inhibitors. The first has a particularly high affinity to serotonin and sigma-1 receptors while the latter renders N,N-DMT orally active by inhibiting the enzyme monoamine oxidase in the gut (Barker, 2018; Fontanilla et al., 2009; Mckenna et al., 1984; Riba, Rodríguez-Fornells, et al., 2001). Recent clinical trials have explored the antidepressant effects of ayahuasca (Osório et al., 2015; Palhano-Fontes et al., 2019; Sanches et al., 2016). In an open-label trial, 17 patients with treatment-resistant depression were submitted to a single session with ayahuasca. Significant antidepressant effects were observed one day after the session and persisted for 21 days (Osório et al., 2015; Sanches et al., 2016). Consistent results were found in a separate randomized placebo-controlled trial, where a significant rapid antidepressant effect of ayahuasca was observed already one day after the session, compared to placebo (Palhano-Fontes et al., 2019). Yet, little is known about the underlying neural mechanisms by which ayahuasca may modulate affect, mood, and internal representations of oneself.

Neuroimaging techniques provide a powerful tool to explore brain changes elicited by psychedelic agents (Dos Santos et al., 2016; Dos Santos et al., 2020). Task-free functional magnetic resonance imaging (tf-fMRI, also referred to as resting-state fMRI) measures synchronous, low frequency (<0.1 Hz) fluctuations in blood oxygen level-dependent (BOLD) activity among distant gray matter regions (Allen et al., 2011; Fox et al., 2005; Smith et al., 2009) and has been used to delineate distinct intrinsic large-scale brain networks supporting primary sensory and higher cognitive and emotional functions (Fox et al., 2005; Grayson and Fair, 2017; Smith et al., 2009; van den Heuvel and Hulshoff Pol, 2010). Importantly, human emotions and social behavior are modulated by a distributed set of brain regions including the anterior cingulate cortex, bilateral anterior insula, and subcortical and limbic structures that regulate affect, interoception, and self-awareness by interacting with the autonomic nervous system (Craig, 2009; Critchley and Harrison, 2013; Lorenzo Pasquini et al., 2019; L. Pasquini et al., 2019; Zhou and Seeley, 2014). These regions coactivate and are structurally connected, cohesively forming the salience network, a large-scale brain system involved in homeostatic behavioral guidance and supporting socioemotional functions by integrating visceral and sensory information (Rankin et al., 2006; Seeley et al., 2007; Sturm, Brown, et al., 2018; Uddin, 2015). A second prominent large-scale brain system is the default mode network, which is primarily composed of the posterior cingulate cortex, precuneus, inferior parietal lobe, parahippocampal gyrus, and medial prefrontal cortex (Buckner et al., 2008). This network has been repeatedly associated with self-referential processes such as mind-wandering, introspection, and memory retrieval under healthy and pathological conditions (Buckner et al., 2008; Raichle et al., 2001).

Recent neuroimaging studies have explored the acute effects of psychedelics on functional organization of human brain networks (Carhart-Harris et al., 2012; Carhart-Harris, Muthukumaraswamy, et al., 2016; Palhano-Fontes et al., 2015; Roseman et al., 2014; Viol et al., 2017, 2019), with reports of increased functional integration at the global brain level (Petri et al., 2014; Tagliazucchi et al., 2016; Viol et al., 2017), and connectivity increases primarily affecting the primary visual cortex and frontal areas (Carhart-Harris, Muthukumaraswamy, et al., 2016; De Araujo et al., 2012). Only one study has reported salience network connectivity decreases under acute conditions (Lebedev et al., 2015), while most studies report decreased connectivity in important hubs of the default mode network, such as the posterior cingulate cortex (Carhart-Harris et al., 2012; Lebedev et al., 2015; Palhano-Fontes et al., 2015).

A recent study assessing subacute neural changes 24 hours after a single dose of ayahuasca reported decreased negative connectivity between the anterior cingulate cortex, a salience network hub, and key default mode network regions such as the posterior cingulate cortex and medial temporal lobes. These connectivity changes correlated with increased self-compassion scores and other measures of enhanced psychological capacities assessed 24 hours after the session, providing a biological basis for the post-acute or “after-glow” stage of psychedelic effects (Majic et al., 2015; Sampedro et al., 2017). Likewise, subacute effects of psilocybin in patients with treatment-resistant depression one day after dosing led to increased connectivity between hubs of the salience and default mode networks such as the subgenual cingulate and posterior cingulate cortices, while connectivity decreases were observed in patients between the prefrontal and parahippocampal cortices, two important default mode network regions (Carhart-Harris et al., 2017). Little is known, however, on inter- and intra-network subacute changes induced by ayahuasca on the salience and default mode networks of naïve healthy participants hours to days after a psychedelic session, and how functional changes may relate to the altered state of consciousness elicited during the psychedelic experience.

To address this critical gap, we leveraged tf-fMRI data assessed before and 24 hours after a randomized placebo-controlled trial exploring the subacute effects of ayahuasca in naïve healthy participants. We hypothesized that 24 hours after ayahuasca administration, we would observe changes in higher-order affective and self-referential networks but not primary sensory networks. Furthermore, we predicted that changes in higher-order affective and self-referential networks would relate to the acute effects of ayahuasca. We applied a seed-based connectivity approach to map the salience, default mode, visual, and sensorimotor networks and assessed subacute inter- and intra-network connectivity changes one day after the session with ayahuasca compared to placebo.

## Methods and Materials

We recruited 50 participants from local media and the Internet. All individuals included in the study were cognitively and physically healthy and naïve to ayahuasca. The following exclusion criteria were adopted: (i) diagnosis of current chronic disease (e.g. heart disease, diabetes) based on anamnesis, physical, and laboratory examination; (ii) pregnancy; (iii) history of neurological/psychiatric diseases; (iv) substance abuse (DSM-IV); (v) restrictions to magnetic resonance imaging assessment. All procedures took place at the Onofre Lopes University Hospital (UFRN), Natal-RN, Brazil. The study was approved by the University Hospital Research Ethics Committee (# 579.479). All subjects provided written informed consent before participation. The study was registered at [http://clinicaltrials.gov](http://clinicaltrials.gov) ([NCT02914769](http://medrxiv.org/lookup/external-ref?link_type=CLINTRIALGOV&access_num=NCT02914769&atom=%2Fmedrxiv%2Fearly%2F2020%2F02%2F21%2F19007542.atom)).

Half of the participants (n = 25) received a single low dose of 1 ml/kg of ayahuasca and the other half (n = 25) received 1 ml/kg of a placebo substance. Alkaloids concentrations were quantified by mass spectroscopy analysis. A single ayahuasca batch was used throughout the study, containing on average (mean ± S.D.): 0.36 ± 0.01 mg/ml of N, N-DMT, 1.86 ± 0.11 mg/ml of harmine, 0.24 ± 0.03 mg/ml of harmaline, and 1.20 ± 0.05 mg/ml of tetrahydroharmine. The batch was prepared and provided by a branch of the Barquinha church, Ji-Paraná, Brazil. The liquid used as placebo was designed to simulate organoleptic properties of ayahuasca, such as a bitter and sour taste, and a brownish color. It contained water, yeast, citric acid, zinc sulfate, and caramel colorant. The presence of zinc sulfate also produced modest gastrointestinal distress allowing to control for the non-specific effects of ayahuasca (Palhano-Fontes et al., 2019). Analogously to two previous studies investigating the subacute effects of psychedelics (Carhart-Harris et al., 2017; Sampedro et al., 2017), baseline tf-fMRI assessments occurred one day before dosing, and a second tf-fMRI assessment occurred 24 hours after dosing (Figure 1A). Blood samples were drawn at baseline and one day after the session.

![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/02/21/19007542/F1.medium.gif)

[Figure 1.](http://medrxiv.org/content/early/2020/02/21/19007542/F1)

Figure 1. Study pipeline.
**(A)** Participants were randomly assigned to the placebo (n = 21; 1 ml/kg of ayahuasca adjusted to contain 0.36 mg/kg of DMT) or ayahuasca group (n = 22; 1 ml/kg of a placebo). Subscales of the Hallucinogenic Rating Scale (HRS) were assessed during the acute effects of ayahuasca/placebo. All participants were assessed with tf-fMRI one day before and one day after dosing. **(B)** Map of the salience, default mode, visual, and sensorimotor networks at baseline overlaid onto the MNI template (warm colors; extent probability threshold of p < 0.01 FWE corrected, voxel-wise threshold of p < 0.001 FWE corrected for multiple comparisons). Individual brain network maps were derived by seeding bilateral region-of-interest (in blue). Left is on the left, bar reflects t-values.

To assess the acute psychedelic experience, we used the Hallucinogenic Rating Scale (HRS), which was initially designed to assess the psychedelic effects of N,N-DMT (Strassman et al., 1994). It is divided into six subscales: (i) intensity, which reflects the strength of the overall experience; (ii) somesthesia, which assesses somatic effects including interoception, visceral, and tactile effects; (iii) affect, which assesses comfort, emotional, and affective changes; (iv) perception, which assesses visual, auditory, gustatory, and olfactory effects; (v) cognition, which assesses changes in thoughts; and (vi) volition, which is related to the subject’s capacity to willfully interact with his/her ‘self’ (Riba, Rodri, et al., 2001; Strassman et al., 1994). We used the validated HRS version translated to Brazilian Portuguese (Mizumoto et al., 2011). Participants completed the HRS at the end of the dosing session (approximately 4 hours after ayahuasca or placebo ingestion).

### Neuroimaging data acquisition

tf-fMRI scanning was performed on a 1.5 Tesla scanner (General Electric, HDxt). Functional MRI data were acquired using an EPI sequence with the following parameters: TR = 2000 ms; TE = 35 ms; flip angle = 60°; FOV = 24 cm; matrix = 64 ⨯ 64; slice thickness = 3 mm; gap = 0.3 mm; number of slices = 35; volumes = 213. T1-weighted images were also acquired using a FSPGR BRAVO sequence with the following parameters: TR = 12.7 ms, TE = 5.3 ms, flip angle = 60°, FOV = 24 cm, matrix = 320 ⨯ 320; slice thickness = 1.0 mm; number of slices = 128.

### Neuroimaging data preprocessing

Of those 50 participants that met our initial inclusion criteria, three participants in the ayahuasca group and four in the placebo group were excluded from further analyses. One subject was excluded due to loss to follow up resulting in the missing of the second scan; one subject due to obvious scanning artifacts as missing of the parietal cortex in the field of view; one participant had to be excluded due to a cystic lesion identified after visually inspecting the anatomical scans after the session; four subjects had to be excluded due to high blood glucose levels indicative of diabetes, which were revealed to the experimentator about a week after completion of the experiment.

The first three fMRI volumes were discarded to allow T1 saturation. Preprocessing steps were performed using CONN ([https://web.conn-toolbox.org/](https://web.conn-toolbox.org/)) (Whitfield-Gabrieli and Nieto-Castanon, 2012) and included: slice timing correction, head motion correction, and spatial smoothing (6-mm FWHM Gaussian kernel) (**Supplementary Figure S1**). Functional images were coregistered to the subject’s anatomical image, normalized into the Montreal Neurological Institute (MNI) template, and resampled to 2 mm3 voxels. Motion artifact was examined using the Artifact Detection Toolbox (ART; [http://www.nitrc.org/projects/artifact_detect/](http://www.nitrc.org/projects/artifact_detect/)). Volumes were considered outliers if the global signal deviation was superior to five standard deviations from the mean signal or if the difference in frame displacement between two consecutive volumes exceeded 0.9 mm. Physiological and other spurious sources of noise were removed using the CompCor method (Behzadi et al., 2007). CompCor is a method used for the reduction of non-physiological noise (head movement, cerebrospinal fluid) in BOLD data. Anatomical data is used to define regions of interest composed primarily of white matter and cerebrospinal fluid. Principal components are derived from these noise-related regions-of-interest and are then included as nuisance parameters within general linear models for BOLD fMRI time-series data denoising. Residual head motion parameters (3 rotation and 3 translation parameters), outlier volumes, white matter, and cerebral spinal fluid signal were regressed out. Finally, a temporal band-pass filter of 0.01 Hz to 0.1 Hz was applied. We evaluated the amount of head motion in our sample by individually estimating the head mean framewise displacement at both sessions (**Supplementary Table S1**). Participants showed head mean framewise displacement way below the widely used threshold of 0.50 mm (Carhart-Harris et al., 2017) and close to the stringent threshold of 0.25 mm recommended by experts in the field (Parkes et al., 2018; Power et al., 2012, 2014; Satterthwaite et al., 2012, 2013).

### Intra-network functional connectivity analyses

For each subject, average BOLD activity time courses were extracted using bilateral regions-of-interest from the Brainnetome Atlas (seed numbers 187 and 188 for the anterior cingulate cortex, seed numbers 175 and 176 for the posterior cingulate cortex, seed numbers 203 and 204 for the occipital cortex, and seed numbers 59 and 60 for the precentral cortex, [http://atlas.brainnetome.org/](http://atlas.brainnetome.org/)) (Fan et al., 2016). Salience, default mode, visual, and sensorimotor network functional connectivity maps were generated through voxel-wise regression analyses using the extracted time-series as regressors (Seeley et al., 2007). Baseline one-sample t-test maps derived from the seed-based approach were visually explored by experts and subsequently spatially correlated with publicly available maps of the visual (Rseed = 0.44), sensorimotor (Rseed = 0.49), default mode (Rseed = 0.69) and salience networks (Rseed = 0.64) (Smith et al., 2009), and interactively compared to maps available on NeuroVault ([https://neurovault.org](https://neurovault.org)) (Gorgolewski et al., 2015). To control for our method of choice, the seed-based approach, we replicated the analysis using group independent component analysis (ICA) implemented through the GIFT toolbox ([http://icatb.sourceforge.net](http://icatb.sourceforge.net)). Briefly, preprocessed tf-fMRI data of all subjects were decomposed into 20 spatially independent components within a group-ICA framework (Calhoun et al., 2001; Pasquini et al., 2015). Data were concatenated and reduced by two-step principal component analysis, followed by independent component estimation with the infomax-algorithm. This procedure results in a set of averaged group components, which are then back reconstructed into the subject’s space. As in the seed-based approach, baseline one-sample t-test maps derived from the ICA-based approach were spatially correlated with publicly available maps of the default mode network (RICA = 0.63) (Buckner et al., 2008) and salience network (RICA = 0.56) (Smith et al., 2009).

### Inter-network functional connectivity analyses

Average time series of resting BOLD activity were derived before and after the session from the one-sample t-test maps of the salience, visual, sensorimotor, and default mode networks identified in the seed-based analyses (**Figure 1B**). Pairwise Pearson’s correlation coefficients were used to estimate the level of functional connectivity between networks at baseline and after the ayahuasca session.

### Statistical analyses

Using individual brain network connectivity maps across both samples at baseline, a voxel-wise one-sample t-test was implemented in SPM12 ([https://www.fil.ion.ucl.ac.uk/spm/software/spm12/](https://www.fil.ion.ucl.ac.uk/spm/software/spm12/)) to assess significant functional connections to our seeds of choice (t = 6.5, extent probability threshold of p < 0.01 FWE corrected, voxel-wise threshold of p < 0.001 FWE corrected for multiple comparisons; **Figure 1B**). These maps were used as masks for subsequent voxel-wise statistical analyses. Change maps of brain network functional connectivity were derived by individually subtracting the functional connectivity maps post-dosing from the ones at baseline. Voxel-wise two-sample t-tests were used to assess differences in brain network functional connectivity between placebo and ayahuasca. Control analyses were performed by adding the difference in head mean framewise displacement between both scans as a covariate of no interest. We explicitly refrained from comparing functional connectivity maps at baseline across the two groups in accordance with the recommendation from the Consolidated Standards of Reporting Trials Group (CONSORT, [http://www.consort-statement.org/](http://www.consort-statement.org/)), which encompasses various initiatives developed by journal editors, epidemiologist, and clinical trialists to alleviate the problems arising from inadequate reporting of randomized control trials. If not reported otherwise, all voxel-wise findings are reported at the joint cluster and extent probability thresholds of p < 0.05 and p < 0.01 uncorrected. Analogously to the intra-network connectivity approach, changes in inter-network connectivity were assessed by subtracting the baseline inter-network functional connectivity value from the post-session value. The change in inter-network functional connectivity was subsequently compared across groups via two-sample t-tests (p < 0.02 uncorrected).

Since the use of liberal statistical thresholds in neuroimaging studies has been recently shown to be prone to inflate false-positives results (Eklund et al., 2016), we additionally assessed group differences in functional connectivity through nonparametric permutation tests as implemented using the FSL package Randomise ([https://fsl.fmrib.ox.ac.uk/fsl/fslwiki/Randomise](https://fsl.fmrib.ox.ac.uk/fsl/fslwiki/Randomise)) (p < 0.05 threshold-free cluster enhancement uncorrected and FWE corrected for multiple comparisons) (Smith and Nichols, 2009; Winkler et al., 2014). Two-sample t-tests were used to assess group differences in demographic variables, head mean framewise displacement, and HRS subscales (p < 0.05, corrected for multiple comparisons if not specified otherwise). Average functional connectivity changes were derived from clusters identified with the voxel-wise two-sample t-tests using an uncorrected joint cluster and extent probability threshold of p < 0.05 for the default mode network and a more stringent joint cluster and extent probability threshold of p < 0.01 for the salience network. Partial correlation coefficients corrected for the difference in head mean framewise displacement between the second and first scan were used to associate intra- and inter-network functional connectivity changes with HRS subscale scores (p<0.05 if not specified otherwise).

### Data and Code availability

The data that support the findings of this study and the code used to generate the findings are available on request from the corresponding author D.B.A. Statistical t-maps shown in this study are publicly available on NeuroVault (Gorgolewski et al., 2015) ([https://identifiers.org/neurovault.image:312874](https://identifiers.org/neurovault.image:312874))([https://identifiers.org/neurovault.image:312875](https://identifiers.org/neurovault.image:312875)). The data are not publicly available due to them containing information that could compromise the privacy of participants in the study.

## Results

Forty-three (43) ayahuasca naïve participants assigned either to ayahuasca (n = 22) or placebo (n = 21) passed our neuroimaging data quality check (**Table 1**). All participants were assessed with tf-fMRI one day before and one day after dosing (**Figure 1A**), and salience, default mode, visual, and sensorimotor networks maps were individually derived by seeding bilateral regions-of-interest (**Figure 1B, Supplementary Table S2**). Participants were comparable in terms of age, gender distribution, and head mean framewise displacement at both sessions (**Table 1**). We found significant between-group differences for all HRS subscales, with significant increases in the ayahuasca compared to the placebo group in levels of altered somesthesia (p < 0.0005), affect (p < 0.005), perception (p < 0.0005), cognition (p < 0.0005), volition (p < 0.009), and intensity (p < 0.0005) (**Table 2**). All HRS subscale comparisons survived multiple correction testing (p < 0.008) with the exception of the volition subscale.

View this table:
[Table 1.](http://medrxiv.org/content/early/2020/02/21/19007542/T1)

Table 1. Sample characteristics.
&Chi-square statistics, chi-square test used to test group differences in gender distribution. Paired t-test used to test differences in head motion between the first and second sessions: At = 1.43, p = 0.17; Bt = 0.41, p = 0.69.

View this table:
[Table 2.](http://medrxiv.org/content/early/2020/02/21/19007542/T2)

Table 2. Hallucinogenic rating scale (HRS).
Subscales assessed during the ayahuasca/placebo sessions. All comparisons survive multiple correction testing (p < 0.008) with the exception of the volition subscale.

### Subacute salience and default mode connectivity changes

We next explored differences in brain network functional connectivity between the ayahuasca and the placebo groups. A voxel-wise two-sample t-test revealed significant increased functional connectivity in the ayahuasca compared to the placebo group within the salience network (t = 1.3; joint cluster and extent probability thresholds of p < 0.05 and p < 0.01) (**Figure 2A** red and yellow maps, **Supplementary Table S3**). Functional connectivity increases with our seed of interest were located bilaterally within the anterior cingulate cortex and superior frontal gyrus, with a tendency towards the left hemisphere (**Figure 2A** first three brain slides and **Supplementary Table S3**). Conversely, we found functional connectivity decreases within the default mode network in the ayahuasca group compared to placebo, predominantly affecting our seed of choice in the posterior cingulate cortex (**Figure 2C** and **Supplementary Table S3**). Average levels of functional connectivity changes derived from the clusters identified in **Figures 2A** and C are schematized as box plots (**Figure 2B** and D). Similar results were derived by adding the difference in head mean framewise displacement as a covariate of no interest in the two-sample t-test models (**Supplementary Figure S2** and **Supplementary Table S4**), suggesting no influence of head motion in functional differences found across groups. Similar functional connectivity differences to our seeds of choice in the salience network and default mode network were also found when using nonparametric permutation-based tests (**Supplementary Figure S5** and **Supplementary Table S5**). Furthermore, consistent with our seed-based findings, ICA revealed patterns of increased anterior cingulate cortex and superior frontal gyrus functional connectivity for the salience network (t = 1.7; joint cluster and extent probability thresholds of p < 0.05, **Supplementary Figure S4 A-B** and **Supplementary Table S6**) and a trend in decreased posterior cingulate cortex functional connectivity of the default mode network in the ayahuasca group (t = 1.3; joint cluster and extent probability thresholds of p < 0.1, **Supplementary Figure S4 C-D** and **Supplementary Table S6**). Inter-network functional connectivity analyses (**Figure 3**) revealed increased functional connectivity between the salience and default mode networks in the ayahuasca group compared to the placebo group (t = 2.5; p < 0.02). Intra-network functional connectivity of primary sensory networks (visual and sensorimotor) and inter-network functional connectivity of primary sensory networks with the default and salience networks did not differ significantly between ayahuasca and placebo, suggesting some level of specificity of the subacute functional changes elicited by ayahuasca on the salience and default mode networks.

![Figure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/02/21/19007542/F2.medium.gif)

[Figure 2.](http://medrxiv.org/content/early/2020/02/21/19007542/F2)

Figure 2. Intra-network functional connectivity changes.
**(A)** Significant salience network functional connectivity increases within the anterior cingulate cortex for the ayahuasca compared to the placebo group. Joint extent and cluster probability thresholds of p < 0.05 (red) and p < 0.01 (yellow). Left is on the left, bar reflects t-values. **(B)** Box plot schematizing group differences in functional connectivity (average functional connectivity levels derived from the cluster identified in panel A with a joint cluster and extent probability thresholds of p < 0.01). **(C)** Significant default mode network functional connectivity decreases within the posterior cingulate cortex for the ayahuasca compared to the placebo group. Joint cluster and extent probability thresholds of p<0.05 (red). Left is on the left, bar reflects t-values. **(D)** Box plot schematizing group differences in functional connectivity (average functional connectivity levels derived from the cluster identified in panel C with a joint cluster and extent probability thresholds of p < 0.05).

![Figure 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/02/21/19007542/F3.medium.gif)

[Figure 3.](http://medrxiv.org/content/early/2020/02/21/19007542/F3)

Figure 3. Inter-network functional connectivity changes.
**(A)** Significant connectivity increases between the salience and default mode networks in the ayahuasca (in red) compared to the placebo group (in cyan). Connectivity between the **(B)** salience network and visual network, **(C)** salience network ad sensorimotor network, **(D)** default mode network and visual network, **(E)** default mode network ad sensorimotor network, **(F)** and visual network and sensorimotor network did not significantly differ between groups. r = Pearson’s correlation coefficient used to assess inter-network functional connectivity. *p < 0.02

### Ayahuasca induced subacute changes in functional connectivity correlate with acute alterations in somesthesia, affect, and volition

Average functional connectivity changes were derived from the previously identified clusters in the salience network and default mode network, and associated with HRS subscores assessed during the acute effects in an explorative way. Partial correlation analyses corrected for the difference in head mean framewise displacement between the second and first scanning session revealed significant positive associations between increases in salience network functional connectivity and somesthesia in the ayahuasca group but not in the placebo group (PRaya = 0.55, p < 0.01; PRpla = 0.01, p = 0.95) (**Figure 4**). Positive trend correlations were found between increased salience network functional connectivity and affect in the ayahuasca group and volition in the placebo group (**Supplementary Table S7**). Subacute changes in default mode network functional connectivity showed a significant negative correlation only with volition in the ayahuasca group and a trend in placebo (PRaya= - 0.59, p < 0.005; Rpla= −0.43, p = 0.06) (**Figure 4, Supplementary Table S7**). Further, subacute changes in functional connectivity between the salience and default mode networks correlated significantly with affect in the ayahuasca but not in the placebo group (PRaya= 0.47, p < 0.03; Rpla= −0.06, p = 0.81) (**Figure 4, Supplementary Table S7**). Trend correlation were found also to cognition in the ayahuasca and to perception in the placebo group (**Supplementary Table S7**). We subsequently estimated multiple linear regression models over all participants, using somesthesia as the dependent variable, and within salience network connectivity change, the difference in head mean framewise displacement between sessions, and a categorical group variable as predictors. This analysis revealed a main effect of group (β = −0.51, t = −3.00, p < 0.005) indicating that the two slopes associating within salience network connectivity change and somesthesia scores significantly differ across the ayahuasca and the placebo groups. An analogous analysis exploring the association between affect and functional connectivity between the salience and default mode networks revealed a similar main effect of group (β = - 0.39, t = −2.64, p < 0.02) indicating significant group differences in the slopes. Using volition as the dependent variable and within default mode network functional connectivity change as predictor revealed that the main effect of group was not significant (β = −0.13, t = −0.88, p = 0.383), indicating that the slopes do not differ across groups.

![Figure 4.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/02/21/19007542/F4.medium.gif)

[Figure 4.](http://medrxiv.org/content/early/2020/02/21/19007542/F4)

Figure 4. Ayahuasca induced subacute changes in functional connectivity correlate with altered somesthesia, volition, and affect.
Boxplot showing group differences in **(A)** somesthesia, **(B)** volition and **(C)** affect assessed with the HRS during the acute effects of ayahuasca (in red) and placebo (in cyan). Scatter plots and partial correlation coefficients corrected for the difference in head mean framewise displacement between the second and first scan associating **(D)** somesthesia, **(E)** volition, or **(F)** affect with averaged levels of functional connectivity changes identified within the salience network, within the default mode network, and between the salience and default mode networks (ayahuasca red circles and line, placebo cyan triangles and line). +p < 0.1, ^p < 0.03, *p < 0.01, **p < 0.005, \***|p < 0.0005.

## Discussion

Neuroimaging techniques provide a unique opportunity to study the neural correlates of altered states of consciousness induced by psychedelic agents (Carhart-Harris et al., 2012; Carhart-Harris, Muthukumaraswamy, et al., 2016; De Araujo et al., 2012; Palhano-Fontes et al., 2015; Riba et al., 2004; Viol et al., 2017; Vollenweider et al., 1997). In this study, ayahuasca naïve participants were randomly assigned to a single session with either ayahuasca or placebo, and functional connectivity of different large-scale brain networks was assessed with tf-fMRI one day before and one day after the dosing session. We found (i) salience network connectivity increases within the anterior cingulate cortex and between the anterior cingulate cortex and the superior frontal gyrus in ayahuasca compared to placebo; (ii) increased connectivity between the salience and default mode networks in the ayahuasca group; (iii) while default mode network connectivity was decreased within the posterior cingulate cortex in the ayahuasca group. Ayahuasca induced subacute functional connectivity increases in the salience network correlated with altered levels of somesthesia, which reflects somatic changes including interoceptive, visceral, and tactile effects induced by psychedelics (Riba, Rodri, et al., 2001). Reduced default mode network connectivity correlated with altered levels of volition, proposed to reflect the subject’s capacity to willfully interact with his/her ‘self’ during the psychedelic experience (Riba, Rodri, et al., 2001). Increased connectivity between the salience and default mode networks correlated with altered levels of affect reflecting emotional responses during the acute psychedelic session (Riba, Rodri, et al., 2001). No subacute inter- and intra-network connectivity differences were detected in primary sensory brain networks, suggesting specificity of the observed changes in the salience and default mode networks. Although little is known on the long-term functional impact of psychedelics on large-scale brain networks, these findings suggest that ayahuasca has sustained effects on neural systems supporting interoceptive, affective, and self-referential processes (Brewer et al., 2013; Craig, 2009; Critchley and Harrison, 2013; Raichle, 2015; Seeley et al., 2007; Uddin, 2015).

### Effects of psychedelic substances on the salience network and default mode network

N,N-DMT, one of the major constituents of ayahuasca, has high affinity to the 5HT2A serotonergic receptor (Barker, 2018), while other constituents such as harmine, harmaline and tetrahydroharmine act as monoamine oxidase inhibitors and serotonin reuptake inhibitors partially by interacting with the serotonin transporter (5-HTT) (Mckenna et al., 1984; Riba, Rodríguez-Fornells, et al., 2001). A recent molecular neuroimaging study involving positron emission tomography acquired in a large sample of healthy individuals revealed *in vivo* multimodal density maps of serotonin transporter and major serotonin receptors in the human brain (Beliveau et al., 2017). 5HT2A receptors displayed a widespread pattern of distribution spanning the visual cortex, temporal, parietal and frontal regions overlapping with major hubs of the salience, default mode, and task-control networks (Buckner et al., 2008; Raichle et al., 2001; Seeley et al., 2007), with density of serotonin transporter being highest within frontotemporal regions overlapping with the insula and anterior cingulate cortex. Previous studies investigating psychedelics have reported decreased within salience network functional coupling during the acute effects of psilocybin (Lebedev et al., 2015). Psilocybin-induced reduced functional connectivity was associated with “ego-dissolution”, a construct frequently reported as altered by psychedelics, characterized by the feeling that the border between one’s self and the external world is dissolving (Goodman, 2002; Griffiths et al., 2011; Lyvers and Meester, 2012; Trichter et al., 2009). Similar studies have reported increased salience network entropy levels during the acute effects of LSD (Lebedev et al., 2016). Ketamine, an N-methyl-D-aspartate receptor antagonist, is a common dissociative anesthetic, which at subanesthetic doses has rapid and sustained antidepressant effects (Ionescu et al., 2018). In both healthy participants and patients with depression, acute ketamine administration has been shown to dampen connectivity of key regions of the salience and default mode networks such as the anterior cingulate and posterior cingulate cortices (Bonhomme et al., 2016; Lehmann et al., 2016; Scheidegger et al., 2012), whereas connectivity between these regions has been shown to increase post-session in patients (Abdallah et al., 2017). A study investigating the structural correlates of long-term ayahuasca use found increased cortical thickness within the anterior cingulate cortex (a region standing out throughout our analyses) of regular ayahuasca users compared to controls, while cortical thickness of default mode network regions such as the precuneus and posterior cingulate cortex was decreased in long-term ayahuasca users (Bouso et al., 2015). The salience network has been proposed to orchestrate dynamic switching between mentalizing states anchored on the default mode network and externally-driven attentional states anchored on the task-control network in order to guide behavior by segregating relevant internal and extrapersonal stimuli (Menon and Uddin, 2010; Uddin, 2015; Zhou and Seeley, 2014). Intriguingly, increased global coupling and functional connectivity between salience and default mode network nodes have been found under acute psilocybin administration (Carhart-Harris et al., 2013) and subacutely with ayahuasca (Sampedro et al., 2017), but future work, ideally combining multimodal molecular and functional neuroimaging, is needed in order to elucidate how mid- and long-term functional interactions of both networks are changed by psychedelic agents.

### Impact of psychedelics on interoception, affect, and self-referential processes

Psychedelic substances modulate blood pressure, body temperature, and heart rate (Holze et al., 2019) by interacting with the sympathetic and parasympathetic autonomic nervous system, and have shown to impact emotional and affective functions by improving clinical symptoms in mood disorders (Carhart-Harris, Bolstridge, et al., 2016; Griffiths et al., 2016; Grob et al., 2011; Palhano-Fontes et al., 2019; Ross et al., 2016; Sanches et al., 2016). By which mechanisms may ayahuasca mediate changes in interoception, affect, and self-referential processes? Hubs of the salience network, such as the insula and the anterior cingulate, have been consistently associated with emotional processing (Craig, 2009; Critchley and Harrison, 2013; Etkin et al., 2011; Ochsner and Gross, 2005; Seeley et al., 2007; Touroutoglou et al., 2012) and dysfunctions in these regions underlie depression and anxiety in various affective disorders (Williams, 2016). Increased connectivity in the salience network has been associated with early heightened emotional contagion in preclinical Alzheimer’s disease (Fredericks et al., 2018; Sturm et al., 2013), while widespread salience network functional and structural degeneration is observed in behavioral variant frontotemporal dementia, a neurodegenerative disease characterized by loss of empathy, socioemotional symptoms, and autonomic dysfunctions (Rankin et al., 2006; Seeley et al., 2012; Sturm, Sible, et al., 2018). Based on these studies, the salience network has been proposed to support socioemotional-autonomic processing through its interoceptive afferents in the anterior insula processing autonomic activity streams regarding the “moment-to-moment” condition of the body (Craig, 2009; Critchley and Harrison, 2013; Uddin, 2015; Zhou and Seeley, 2014). The anterior cingulate cortex subsequently receives integrated anterior insula input and serves to mobilize visceromotor responses to salient socioemotional stimuli in order to guide behavior (Craig, 2009; Critchley and Harrison, 2013; Uddin, 2015; Zhou and Seeley, 2014). Critically, our findings of subacute default mode network functional connectivity decreases with ayahuasca are in contrast with a previous study reporting subacute default mode network connectivity increases with psilocybin in treatment-resistant depression (Carhart-Harris et al., 2017). Previous works have reported both default mode network connectivity increases and decreases in major depression (Berman et al., 2011; Greicius et al., 2007; Sheline et al., 2009; Yan et al., 2019), with heightened connectivity being associated with altered self-referential thoughts such as rumination and negative internal representations (Berman et al., 2011; Sheline et al., 2009). Although in healthy subjects, our findings provide preliminary evidence that ayahuasca may have a sustained effect on dampening functional connectivity of systems involved in self-referential processes associated with the characteristic symptoms of depression (Brewer et al., 2013; Sheline et al., 2009). We propose that future studies combining longitudinal tf-fMRI and autonomic recordings in healthy and neuropsychiatric populations could shed light on the acute and long-term impact of psychedelic agents on brain networks supporting interoceptive, affective, and self-referential processes.

### Limitations

An important limitation of our study consists in the use of liberal voxel-wise statistical thresholds when assessing group differences in subacute functional connectivity across distinct brain networks. Recent research has shown that liberal cluster-extend thresholds are prone to inflate false-positive results, questioning the validity of weakly significant neuroimaging findings (Eklund et al., 2016). The low power attained in our study may have been caused by the small number of participants involved and by the moderate strength of the MRI scanner (1.5 Tesla), urging for future studies involving larger sample sizes and state-of-the-art neuroimaging acquisition. Further, weak effects on functional connectivity were to be expected after a single ayahuasca session in a sample of cognitively healthy participants. Future studies implementing longitudinal sessions and the inclusion of neuropsychiatric populations could help overcome this limitation. We performed, however, several control analyses in order to mitigate methodological and statistical concerns. First, we assessed the reliability of our findings by controlling for the used methodological approach, which revealed similar functional connectivity changes when comparing findings from the seed-based versus ICA approach. In order to address concerns related to the liberal statistical threshold, we applied nonparametric permutation tests replicating the main findings in this study. Second, we assessed the impact of head movement on intra-network functional connectivity differences across groups and on the association between altered states of consciousness and subacute functional connectivity changes through the use of partial correlation analyses. These analyses revealed that head movement in the scanner, a common confounder in tf-fMRI studies, did not significantly impact our findings. However, findings of decreased functional connectivity, particularly within the default mode network need to be interpreted with a grain of salt, since these did not survive more stringent statistical thresholds. Further, the volition subscale of the HRS has been shown to have a low internal consistency (Riba, Rodri, et al., 2001). Given the exploratory nature of the correlational analyses and the associated danger of false-positive inflation, we advise caution in interpreting the association between altered acute levels in the volition scale and subacute default mode network functional connectivity decreases induced by ayahuasca. In particular, both the ayahuasca and the placebo groups exhibited a comparable relationship between changes in default mode network functional connectivity and volition. Although studies suggest that the anterior cingulate cortex, rather than posterior brain regions, has a critical role in modulating generalized effects (Sikora et al., 2016; Wager et al., 2004), our findings could potentially reflect general placebo effects.

## Conclusions

Psychedelic substances have sustained effects on affect and self-referential processes in healthy and clinical populations days to weeks after dosing. The novelty of our study resides in elucidating the subacute effects of the psychedelic ayahuasca on functional organization of the salience and default mode networks, two brain systems distinctly involved in interoceptive, affective, and self-referential functions. While primary sensory networks did not show subacute changes in functional connectivity, increased functional connectivity of the salience network one day after the session with ayahuasca related to altered acute somesthesia levels, decreases in default mode network functional connectivity related to altered levels of volition, while salience network-default mode network connectivity increases related to altered affect levels. Our findings suggest that ayahuasca may have long-lasting effects on mood by modulating those neural circuits supporting interoceptive, affective, and self-referential functions.

## Data Availability

Data is available upon reasonable request to the corresponding author.

## Financial Disclosures

The authors have no conflict of interest to declare.

## Acknowledgments

L.P. was supported by the German Academic Foundation. The study was funded by the Brazilian federal agencies CNPq (grants #466760/2014 & #479466/2013) and CAPES (grants #1677/2012 & #1577/2013).

*   Received September 27, 2019.
*   Revision received February 18, 2020.
*   Accepted February 21, 2020.


*   © 2020, Posted by Cold Spring Harbor Laboratory

This pre-print is available under a Creative Commons License (Attribution-NonCommercial-NoDerivs 4.0 International), CC BY-NC-ND 4.0, as described at [http://creativecommons.org/licenses/by-nc-nd/4.0/](http://creativecommons.org/licenses/by-nc-nd/4.0/)

## References

1.  Abdallah CG, Averill LA, Collins KA, et al. (2017) Ketamine Treatment and Global Brain Connectivity in Major Depression. Neuropsychopharmacology: official publication of the American College of Neuropsychopharmacology 42(6): 1210–1219.
    
    

2.  Allen EA, Erhardt EB, Damaraju E, et al. (2011) A baseline for the multivariate comparison of resting-state networks. Frontiers in systems neuroscience 5: 2.
    
    

3.  Barker SA (2018) N, N-Dimethyltryptamine (DMT), an Endogenous Hallucinogen: Past, Present, and Future Research to Determine Its Role and Function. Frontiers in neuroscience 12: 536.
    
    

4.  Behzadi Y, Restom K, Liau J, et al. (2007) A component based noise correction method (CompCor) for BOLD and perfusion based fMRI. NeuroImage 37(1): 90–101.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neuroimage.2007.04.042&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17560126&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F02%2F21%2F19007542.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000248526700010&link_type=ISI) 

5.  Beliveau V, Ganz M, Feng L, et al. (2017) A High-Resolution In Vivo Atlas of the Human Brain’s Serotonin System. The Journal of Neuroscience 37(1): 120–128.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Njoiam5ldXJvIjtzOjU6InJlc2lkIjtzOjg6IjM3LzEvMTIwIjtzOjQ6ImF0b20iO3M6Mzk6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDIvMjEvMTkwMDc1NDIuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

6.  Berman MG, Peltier S, Nee DE, et al. (2011) Depression, rumination and the default network. Social cognitive and affective neuroscience 6(5): 548–555.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/scan/nsq080&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20855296&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F02%2F21%2F19007542.atom) 

7.  Bonhomme V, Vanhaudenhuyse A, Demertzi A, et al. (2016) Resting-state Network-specific Breakdown of Functional Connectivity during Ketamine Alteration of Consciousness in Volunteers. Anesthesiology 125(5): 873–888.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1097/ALN.0000000000001275&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27496657&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F02%2F21%2F19007542.atom) 

8.  Bouso JC, Palhano-Fontes F, Rodríguez-Fornells A, et al. (2015) Long-term use of psychedelic drugs is associated with differences in brain structure and personality in humans. European Neuropsychopharmacology 25(4): 483–492.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.euroneuro.2015.01.008&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25637267&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F02%2F21%2F19007542.atom) 

9.  Brewer J, Garrison K and Whitfield-Gabrieli S (2013) What about the ‘self’ is processed in the posterior cingulate cortex? Frontiers in human neuroscience 7: 647.
    
    

10. Buckner RL, Andrews-Hanna JR and Schacter DL (2008) The brain’s default network: anatomy, function, and relevance to disease. Annals of the New York Academy of Sciences 1124: 1–38.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1196/annals.1440.011&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18400922&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F02%2F21%2F19007542.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000257139500002&link_type=ISI) 

11. Calhoun VD, Adali T, Pearlson G, et al. (2001) Group ICA of functional MRI data: separability, stationarity, and inference. In: Proc. Int. Conf. on ICA and BSS.
    
    

12. Carhart-Harris RL, Erritzoe D, Williams T, et al. (2012) Neural correlates of the psychedelic state as determined by fMRI studies with psilocybin. Proceedings of the National Academy of Sciences 109(6): 2138–2143.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMDoiMTA5LzYvMjEzOCI7czo0OiJhdG9tIjtzOjM5OiIvbWVkcnhpdi9lYXJseS8yMDIwLzAyLzIxLzE5MDA3NTQyLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

13. Carhart-Harris RL, Leech R, Erritzoe D, et al. (2013) Functional connectivity measures after psilocybin inform a novel hypothesis of early psychosis. Schizophrenia bulletin 39(6): 1343–1351.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/schbul/sbs117&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23044373&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F02%2F21%2F19007542.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000326189700022&link_type=ISI) 

14. Carhart-Harris RL, Muthukumaraswamy S, Roseman L, et al. (2016) Neural correlates of the LSD experience revealed by multimodal neuroimaging. Proceedings of the National Academy of Sciences: 2–7.
    
    

15. Carhart-Harris RL, Bolstridge M, Rucker J, et al. (2016) Psilocybin with psychological support for treatment-resistant depression: an open-label feasibility study. The Lancet Psychiatry 3(7): 619–627.
    
    

16. Carhart-Harris RL, Roseman L, Bolstridge M, et al. (2017) Psilocybin for treatment-resistant depression: fMRI-measured brain mechanisms. Scientific reports 7(1): 13187.
    
    

17. Craig ADB (2009) How do you feel--now? The anterior insula and human awareness. Nature reviews. Neuroscience 10(1): 59–70.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nrn2555&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19096369&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F02%2F21%2F19007542.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000261934500015&link_type=ISI) 

18. Critchley HD and Harrison NA (2013) Visceral influences on brain and behavior. Neuron 77(4): 624–638.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neuron.2013.02.008&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23439117&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F02%2F21%2F19007542.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000315561300006&link_type=ISI) 

19. De Araujo DB, Ribeiro S, Cecchi GA, et al. (2012) Seeing with the eyes shut: Neural basis of enhanced imagery following ayahuasca ingestion. Human brain mapping 33(11): 2550–2560.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/hbm.21381&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21922603&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F02%2F21%2F19007542.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000310264500004&link_type=ISI) 

20. Dos Santos RG, Balthazar FM, Bouso JC, et al. (2016) The current state of research on ayahuasca: A systematic review of human studies assessing psychiatric symptoms, neuropsychological functioning, and neuroimaging. Journal of psychopharmacology 30(12): 1230–1247.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1177/0269881116652578&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27287824&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F02%2F21%2F19007542.atom) 

21. Dos Santos RG and Hallak JEC (2020) Therapeutic use of serotoninergic hallucinogens: A review of the evidence and of the biological and psychological mechanisms. Neuroscience & Biobehavioral Reviews. DOI: 10.1016/j.neubiorev.2019.12.001.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neubiorev.2019.12.001&link_type=DOI) 

22. Eklund A, Nichols TE and Knutsson H (2016) Cluster failure: Why fMRI inferences for spatial extent have inflated false-positive rates. Proceedings of the National Academy of Sciences 113(28): 7900–7905.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMToiMTEzLzI4Lzc5MDAiO3M6NDoiYXRvbSI7czozOToiL21lZHJ4aXYvZWFybHkvMjAyMC8wMi8yMS8xOTAwNzU0Mi5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

23. Etkin A, Egner T and Kalisch R (2011) Emotional processing in anterior cingulate and medial prefrontal cortex. Trends in cognitive sciences 15(2): 85–93.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.tics.2010.11.004&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21167765&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F02%2F21%2F19007542.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000287295500005&link_type=ISI) 

24. Fan L, Li H, Zhuo J, et al. (2016) The Human Brainnetome Atlas: A New Brain Atlas Based on Connectional Architecture. Cerebral cortex 26(8): 3508–3526.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/cercor/bhw157&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27230218&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F02%2F21%2F19007542.atom) 

25. Fontanilla D, Johannessen M, Hajipour AR, et al. (2009) The hallucinogen N,N-dimethyltryptamine (DMT) is an endogenous sigma-1 receptor regulator. Science 323(5916): 934–937.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEyOiIzMjMvNTkxNi85MzQiO3M6NDoiYXRvbSI7czozOToiL21lZHJ4aXYvZWFybHkvMjAyMC8wMi8yMS8xOTAwNzU0Mi5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

26. Fox MD, Snyder AZ, Vincent JL, et al. (2005) The human brain is intrinsically organized into dynamic, anticorrelated functional networks. Proceedings of the National Academy of Sciences 102(27): 9673–9678.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMToiMTAyLzI3Lzk2NzMiO3M6NDoiYXRvbSI7czozOToiL21lZHJ4aXYvZWFybHkvMjAyMC8wMi8yMS8xOTAwNzU0Mi5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

27. Fredericks CA, Sturm VE, Brown JA, et al. (2018) Early affective changes and increased connectivity in preclinical Alzheimer’s disease. Alzheimer’s & Dementia 10: 471–479.
    
    

28. Goodman N (2002) The serotonergic system and mysticism: could LSD and the nondrug-induced mystical experience share common neural mechanisms? Journal of psychoactive drugs 34(3): 263–272.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1080/02791072.2002.10399962&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12422936&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F02%2F21%2F19007542.atom) 

29. Gorgolewski KJ, Varoquaux G, Rivera G, et al. (2015) NeuroVault.org: a web-based repository for collecting and sharing unthresholded statistical maps of the human brain. Frontiers in neuroinformatics 9: 8.
    
    

30. Grayson DS and Fair DA (2017) Development of large-scale functional networks from birth to adulthood: A guide to the neuroimaging literature. NeuroImage 160: 15–31.
    
    

31. Greicius MD, Flores BH, Menon V, et al. (2007) Resting-State Functional Connectivity in Major Depression: Abnormally Increased Contributions from Subgenual Cingulate Cortex and Thalamus. Biological psychiatry 62(5): 429–437.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.biopsych.2006.09.020&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17210143&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F02%2F21%2F19007542.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000249042800009&link_type=ISI) 

32. Griffiths RR, Johnson MW, Richards WA, et al. (2011) Psilocybin occasioned mystical-type experiences: immediate and persisting dose-related effects. Psychopharmacology 218(4): 649–665.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s00213-011-2358-5&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21674151&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F02%2F21%2F19007542.atom) 

33. Griffiths RR, Johnson MW, Carducci MA, et al. (2016) Psilocybin produces substantial and sustained decreases in depression and anxiety in patients with life-threatening cancer: A randomized double-blind trial. Journal of psychopharmacology 30(12): 1181–1197.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1177/0269881116675513&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27909165&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F02%2F21%2F19007542.atom) 

34. Grob CS, Danforth AL, Chopra GS, et al. (2011) Pilot Study of Psilocybin Treatment for Anxiety in Patients With Advanced-Stage Cancer. Archives of general psychiatry 68(1): 71.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/archgenpsychiatry.2010.116&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20819978&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F02%2F21%2F19007542.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000285927800014&link_type=ISI) 

35. Holze F, Vizeli P, Müller F, et al. (2019) Distinct acute effects of LSD, MDMA, and D-amphetamine in healthy subjects. Neuropsychopharmacology: official publication of the American College of Neuropsychopharmacology. DOI: 10.1038/s41386-019-0569-3.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41386-019-0569-3&link_type=DOI) 

36. Ionescu DF, Felicione JM, Gosai A, et al. (2018) Ketamine-Associated Brain Changes: A Review of the Neuroimaging Literature. Harvard review of psychiatry 26(6): 320–339.
    
    

37. Lebedev AV, Lövdén M, Rosenthal G, et al. (2015) Finding the self by losing the self: Neural correlates of ego-dissolution under psilocybin. Human brain mapping 36(8): 3137–3153.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/hbm.22833&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26010878&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F02%2F21%2F19007542.atom) 

38. Lebedev AV, Kaelen M, Lövdén M, et al. (2016) LSD-induced entropic brain activity predicts subsequent personality change. Human brain mapping 37(9): 3203–3213.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/hbm.23234&link_type=DOI) 

39. Lehmann M, Seifritz E, Henning A, et al. (2016) Differential effects of rumination and distraction on ketamine induced modulation of resting state functional connectivity and reactivity of regions within the default-mode network. Social cognitive and affective neuroscience 11(8): 1227–1235.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/scan/nsw034&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27075438&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F02%2F21%2F19007542.atom) 

40. Lyvers M and Meester M (2012) Illicit use of LSD or psilocybin, but not MDMA or nonpsychedelic drugs, is associated with mystical experiences in a dose-dependent manner. Journal of psychoactive drugs 44(5): 410–417.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1080/02791072.2012.736842&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23457892&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F02%2F21%2F19007542.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000312747100008&link_type=ISI) 

41. Majić T, Schmidt TT and Gallinat J (2015) Peak experiences and the afterglow phenomenon: When and how do therapeutic effects of hallucinogens depend on psychedelic experiences? Journal of psychopharmacology 29(3): 241–253.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1177/0269881114568040&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25670401&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F02%2F21%2F19007542.atom) 

42. Mckenna DJ, Towers GHN and Abbott F (1984) Monoamine-Oxidase Inhibitors in South-American Hallucinogenic Plants - Tryptamine and Beta-Carboline Constituents of Ayahuasca. Journal of ethnopharmacology 10(2): 195–223.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/0378-8741(84)90003-5&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=6587171&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F02%2F21%2F19007542.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1984SQ64300003&link_type=ISI) 

43. Menon V and Uddin LQ (2010) Saliency, switching, attention and control: a network model of insula function. Brain structure & function 214(5-6): 655–667.
    
    

44. Mizumoto S, da Silveira DX, Barbosa PCR, et al. (2011) Hallucinogen Rating Scale (HRS) – Versão brasileira: tradução e adaptação transcultural Hallucinogen Rating Scale (HRS) – A Brazilian version: translation and cross-cultural adaptation. Rev Psiq Clín. 38(6): 231–237.
    
    

45. Nichols DE (2016) Psychedelics. Pharmacological reviews 68(2): 264–355.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoicGhhcm1yZXYiO3M6NToicmVzaWQiO3M6ODoiNjgvMi8yNjQiO3M6NDoiYXRvbSI7czozOToiL21lZHJ4aXYvZWFybHkvMjAyMC8wMi8yMS8xOTAwNzU0Mi5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

46. Ochsner KN and Gross JJ (2005) The cognitive control of emotion. Trends in cognitive sciences 9(5): 242–249.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.tics.2005.03.010&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15866151&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F02%2F21%2F19007542.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000229466500010&link_type=ISI) 

47.  Osório F de L, Sanches RF, Macedo LR, et al. (2015) Antidepressant effects of a single dose of ayahuasca in patients with recurrent depression: a preliminary report. Brazilian Journal of Psychiatry 37(1): 13–20.
    
    

48. Palhano-Fontes F, Andrade KC, Tofoli LF, et al. (2015) The psychedelic state induced by Ayahuasca modulates the activity and connectivity of the Default Mode Network. PloS one 10(2).
    
    

49. Palhano-Fontes F, Barreto D, Onias H, et al. (2019) Rapid antidepressant effects of the psychedelic ayahuasca in treatment-resistant depression: a randomized placebo-controlled trial. Psychological Medicine 49(4): 655–663.
    
    

50. Parkes L, Fulcher B, Yücel M, et al. (2018) An evaluation of the efficacy, reliability, and sensitivity of motion correction strategies for resting-state functional MRI. NeuroImage 171: 415–436.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/J.NEUROIMAGE.2017.12.073&link_type=DOI) 

51. Pasquini L, Scherr M, Tahmasian M, et al. (2015) Link between hippocampus’ raised local and eased global intrinsic connectivity in AD. Alzheimer’s & Dementia 11(5): 475–484.
    
    

52. Pasquini L., Nana AL, Toller G, et al. (2019) Salience network atrophy links neuron type-specific degeneration to loss of empathy in frontotemporal dementia. bioRxiv. biorxiv.org. Available at: [https://www.biorxiv.org/content/10.1101/691212v2.abstract](https://www.biorxiv.org/content/10.1101/691212v2.abstract).
    
    

53.  Pasquini Lorenzo, Toller G, Staffaroni A, et al. (2019) State and trait characteristics of anterior insula time-varying functional connectivity. NeuroImage 208: 116425.
    
    

54. Petri G, Expert P, Turkheimer F, et al. (2014) Homological scaffolds of brain functional networks. Journal of the Royal Society 11(101).
    
    

55. Power JD, Barnes KA, Snyder AZ, et al. (2012) Spurious but systematic correlations in functional connectivity MRI networks arise from subject motion. NeuroImage 59(3): 2142–2154.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neuroimage.2011.10.018&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22019881&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F02%2F21%2F19007542.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000299494000017&link_type=ISI) 

56. Power JD, Mitra A, Laumann TO, et al. (2014) Methods to detect, characterize, and remove motion artifact in resting state fMRI. NeuroImage 84: 320–341.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neuroimage.2013.08.048&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23994314&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F02%2F21%2F19007542.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000328868600030&link_type=ISI) 

57. Raichle ME (2015) The brain’s default mode network. Annual review of neuroscience 38: 433–447.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1146/annurev-neuro-071013-014030&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25938726&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F02%2F21%2F19007542.atom) 

58. Raichle ME, MacLeod AM, Snyder AZ, et al. (2001) A default mode of brain function. Proceedings of the National Academy of Sciences 98(2): 676–682.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czo4OiI5OC8yLzY3NiI7czo0OiJhdG9tIjtzOjM5OiIvbWVkcnhpdi9lYXJseS8yMDIwLzAyLzIxLzE5MDA3NTQyLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

59. Rankin KP, Gorno-Tempini ML, Allison SC, et al. (2006) Structural anatomy of empathy in neurodegenerative disease. Brain 129(11): 2945–2956.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/brain/awl254&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17008334&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F02%2F21%2F19007542.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000241783800015&link_type=ISI) 

60. Riba J, Rodrı A, Strassman RJ, et al. (2001) Psychometric Assessment of the Hallucinogen Rating Scale. Drug Alcohol Depend 62: 215–223.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0376-8716(00)00175-7&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11295326&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F02%2F21%2F19007542.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000168188700008&link_type=ISI) 

61. Riba J, Rodríguez-Fornells A, Urbano G, et al. (2001) Subjective effects and tolerability of the South American psychoactive beverage Ayahuasca in healthy volunteers. Psychopharmacology 154(1): 85–95.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s002130000606&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11292011&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F02%2F21%2F19007542.atom) 

62. Riba J, Anderer P, Jané F, et al. (2004) Effects of the South American psychoactive beverage ayahuasca on regional brain electrical activity in humans: a functional neuroimaging study using low-resolution electromagnetic tomography. Neuropsychobiology 50(1): 89–101.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1159/000077946&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15179026&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F02%2F21%2F19007542.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000221802600015&link_type=ISI) 

63. Roseman L, Leech R, Feilding A, et al. (2014) The effects of psilocybin and MDMA on between-network resting state functional connectivity in healthy volunteers. Frontiers in human neuroscience 8: 204.
    
    

64. Ross S, Bossis A, Guss J, et al. (2016) Rapid and sustained symptom reduction following psilocybin treatment for anxiety and depression in patients with life-threatening cancer: a randomized controlled trial. Journal of psychopharmacology 30(12): 1165–1180.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1177/0269881116675512&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27909164&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F02%2F21%2F19007542.atom) 

65. Sampedro F, de la Fuente Revenga M, Valle M, et al. (2017) Assessing the psychedelic ‘ after-glow ‘ in ayahuasca users: post-acute neurometabolic and functional connectivity changes are associated with enhanced mindfulness capacities. Int. J. Neuropsychopharmacol. 20(9): 698–711.
    
    

66. Sanches RF, de Lima Osório F, Dos Santos RG, et al. (2016) Antidepressant Effects of a Single Dose of Ayahuasca in Patients With Recurrent Depression: A SPECT Study. Journal of clinical psychopharmacology 36(1): 77–81.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1097/JCP.0000000000000436&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26650973&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F02%2F21%2F19007542.atom) 

67. Satterthwaite TD, Wolf DH, Loughead J, et al. (2012) Impact of in-scanner head motion on multiple measures of functional connectivity: relevance for studies of neurodevelopment in youth. NeuroImage 60(1): 623–632.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neuroimage.2011.12.063&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22233733&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F02%2F21%2F19007542.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000301218700064&link_type=ISI) 

68. Satterthwaite TD, Elliott MA, Gerraty RT, et al. (2013) An improved framework for confound regression and filtering for control of motion artifact in the preprocessing of resting-state functional connectivity data. NeuroImage 64: 240–256.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neuroimage.2012.08.052&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22926292&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F02%2F21%2F19007542.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000312504200023&link_type=ISI) 

69. Scheidegger M, Walter M, Lehmann M, et al. (2012) Ketamine decreases resting state functional network connectivity in healthy subjects: implications for antidepressant drug action. PloS one 7(9): e44799.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0044799&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23049758&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F02%2F21%2F19007542.atom) 

70. Seeley WW, Menon V, Schatzberg AF, et al. (2007) Dissociable intrinsic connectivity networks for salience processing and executive control. The Journal of Neuroscience 27(9): 2349–2356.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Njoiam5ldXJvIjtzOjU6InJlc2lkIjtzOjk6IjI3LzkvMjM0OSI7czo0OiJhdG9tIjtzOjM5OiIvbWVkcnhpdi9lYXJseS8yMDIwLzAyLzIxLzE5MDA3NTQyLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

71. Seeley WW, Zhou J and Kim E-J (2012) Frontotemporal dementia: what can the behavioral variant teach us about human brain organization? The Neuroscientist 18(4): 373–385.
    
    

72. Sheline YI, Barch DM, Price JL, et al. (2009) The default mode network and self-referential processes in depression. Proceedings of the National Academy of Sciences 106(6): 1942–1947.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMDoiMTA2LzYvMTk0MiI7czo0OiJhdG9tIjtzOjM5OiIvbWVkcnhpdi9lYXJseS8yMDIwLzAyLzIxLzE5MDA3NTQyLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

73. Sikora M, Heffernan J, Avery ET, et al. (2016) Salience Network Functional Connectivity Predicts Placebo Effects in Major Depression. Biological psychiatry. Cognitive neuroscience and neuroimaging 1(1): 68–76.
    
    

74. Smith SM and Nichols TE (2009) Threshold-free cluster enhancement: addressing problems of smoothing, threshold dependence and localisation in cluster inference. NeuroImage 44(1): 83–98.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neuroimage.2008.03.061&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18501637&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F02%2F21%2F19007542.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000262300900010&link_type=ISI) 

75. Smith SM, Fox PT, Miller KL, et al. (2009) Correspondence of the brain’s functional architecture during activation and rest. Proceedings of the National Academy of Sciences 106(31): 13040–13045.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMjoiMTA2LzMxLzEzMDQwIjtzOjQ6ImF0b20iO3M6Mzk6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDIvMjEvMTkwMDc1NDIuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

76. Strassman RJ, Qualls CR, Uhlenhuth EH, et al. (1994) Dose-response study of N,N-dimethyltryptamine in humans: II. Subjective effects and preliminary results of a new rating scale. Archives of general psychiatry 51(2): 98–108.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/archpsyc.1994.03950020022002&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=8297217&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F02%2F21%2F19007542.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1994MW11800002&link_type=ISI) 

77. Sturm VE, Yokoyama JS, Seeley WW, et al. (2013) Heightened emotional contagion in mild cognitive impairment and Alzheimer’s disease is associated with temporal lobe degeneration. Proceedings of the National Academy of Sciences 110(24): 9944–9949.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMToiMTEwLzI0Lzk5NDQiO3M6NDoiYXRvbSI7czozOToiL21lZHJ4aXYvZWFybHkvMjAyMC8wMi8yMS8xOTAwNzU0Mi5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

78. Sturm VE, Brown JA, Hua AY, et al. (2018) Network Architecture Underlying Basal Autonomic Outflow: Evidence from Frontotemporal Dementia. The Journal of Neuroscience 38(42): 8943–8955.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Njoiam5ldXJvIjtzOjU6InJlc2lkIjtzOjEwOiIzOC80Mi84OTQzIjtzOjQ6ImF0b20iO3M6Mzk6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDIvMjEvMTkwMDc1NDIuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

79. Sturm VE, Sible IJ, Datta S, et al. (2018) Resting parasympathetic dysfunction predicts prosocial helping deficits in behavioral variant frontotemporal dementia. Cortex 109: 141–155.
    
    

80. Tagliazucchi E, Roseman L, Kaelen M, et al. (2016) Increased Global Functional Connectivity Correlates with LSD-Induced Ego Dissolution. Current biology: CB 26(8): 1043–1050.
    
    

81. Touroutoglou A, Hollenbeck M, Dickerson BC, et al. (2012) Dissociable large-scale networks anchored in the right anterior insula subserve affective experience and attention. NeuroImage 60(4): 1947–1958.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neuroimage.2012.02.012&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22361166&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F02%2F21%2F19007542.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000302926600003&link_type=ISI) 

82. Trichter S, Klimo J and Krippner S (2009) Changes in Spirituality Among Ayahuasca Ceremony Novice Participants. Journal of psychoactive drugs 41(2): 121–134.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19705674&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F02%2F21%2F19007542.atom) 

83. Uddin LQ (2015) Salience processing and insular cortical function and dysfunction. Nature reviews. Neuroscience 16(1): 55–61.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nrn3857&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25406711&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F02%2F21%2F19007542.atom) 

84. van den Heuvel MP and Hulshoff Pol HE (2010) Exploring the brain network: A review on resting-state fMRI functional connectivity. European neuropsychopharmacology 20(8): 519–534.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.euroneuro.2010.03.008&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20471808&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F02%2F21%2F19007542.atom) 

85. Viol A, Palhano-Fontes F, Onias H, et al. (2017) Shannon entropy of brain functional complex networks under the influence of the psychedelic Ayahuasca. Scientific reports 7(1): 7388.
    
    

86. Viol A, Palhano-Fontes F, Onias H, et al. (2019) Characterizing complex networks using Entropy-degree diagrams: unveiling changes in functional brain connectivity induced by Ayahuasca. Entropy. Available at: [https://www.mdpi.com/1099-4300/21/2/128](https://www.mdpi.com/1099-4300/21/2/128).
    
    

87. Vollenweider FX, Leenders KL, Scharfetter C, et al. (1997) Positron emission tomography and fluorodeoxyglucose studies of metabolic hyperfrontality and psychopathology in the psilocybin model of psychosis. Neuropsychopharmacology 16(5): 357–372.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0893-133X(96)00246-1&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=9109107&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F02%2F21%2F19007542.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1997WT84300006&link_type=ISI) 

88. Wager TD, Rilling JK, Smith EE, et al. (2004) Placebo-induced changes in FMRI in the anticipation and experience of pain. Science 303(5661): 1162–1167.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEzOiIzMDMvNTY2MS8xMTYyIjtzOjQ6ImF0b20iO3M6Mzk6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDIvMjEvMTkwMDc1NDIuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

89. Whitfield-Gabrieli S and Nieto-Castanon A (2012) Conn: a functional connectivity toolbox for correlated and anticorrelated brain networks. Brain connectivity 2(3): 125–141.
    
    

90. Williams LM (2016) Precision psychiatry: a neural circuit taxonomy for depression and anxiety. The Lancet Psychiatry 3(5): 472–480.
    
    

91. Winkler AM, Ridgway GR, Webster MA, et al. (2014) Permutation inference for the general linear model. NeuroImage 92: 381–397.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neuroimage.2014.01.060&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24530839&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F02%2F21%2F19007542.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000335713000035&link_type=ISI) 

92. Yan C-G, Chen X, Li L, et al. (2019) Reduced default mode network functional connectivity in patients with recurrent major depressive disorder. Proceedings of the National Academy of Sciences 116(18): 9078–9083.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMToiMTE2LzE4LzkwNzgiO3M6NDoiYXRvbSI7czozOToiL21lZHJ4aXYvZWFybHkvMjAyMC8wMi8yMS8xOTAwNzU0Mi5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

93. Zhou J and Seeley WW (2014) Network dysfunction in Alzheimer’s disease and frontotemporal dementia: implications for psychiatry. Biological psychiatry 75(7): 565–573.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.biopsych.2014.01.020&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24629669&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F02%2F21%2F19007542.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000333081900008&link_type=ISI)