Open-label trial
Depression severity
Of the 19 patients recruited, 3 were excluded due to excessive fMRI head-motion. The remaining 16 patients (mean age=42.75, SD=10.15, 4 female) were the final analytical sample (Figure 2a). Baseline BDI scores indicated severe depression (mean BDI=34.81, SD=7.38). As previously reported, rapid and sustained reductions in depression severity were observed post-treatment 31. Relative to baseline, significant BDI reductions were observed at 1 week (mean difference, -21.0 points; 95% CI, -27.30 to -14.71, P<.001) and this was still evident 6 months (mean difference, -14.19 points; 95% CI, -21.29 to -7.09, P<.001).
Decreased brain modularity following psilocybin therapy
Confirming our primary hypothesis prediction, brain network modularity was significantly reduced (Figure 3a) one day after psilocybin therapy (mean difference, -0.29; 95% CI 0.07 to 0.50, P=.012), indicative of an increased integration of brain networks.
Decreased modularity predicts long-term clinical outcomes
We hypothesised that decreased brain network modularity would relate to sustained improvements in depression severity following psilocybin therapy. To test this, we calculated Pearson correlations between the post-treatment brain modularity and the BDI scores from the 3 post-treatment timepoints (1 week, 3 months, 6 months). After false discovery rate (FDR) correction for multiple-comparisons, a strong significant correlation was observed at 6 months (Figure 3b - Pearson, r=0.64; P=.023). Although consistent with this, relationships at 3 months (r=0.46; P=.114) or 1 week (r=0.29; P=.284) did not survive correction. Furthermore, the pre vs post-treatment change in modularity significantly correlated with the change in BDI score at 6 months, relative to baseline (Figure 3c - Pearson, r=0.54; P=.033). These results indicate that decreased brain modularity relates to long-term improvements in depression symptom severity.
Decreased DMN & increased DMN-frontoparietal FC post-treatment
Consistent with previous work and our a priori assumptions, psilocybin therapy was related to significantly (FDR-corrected) decreased DMN network recruitment (Figure 3d - mean difference, -0.54; 95% CI, -0.92 to -0.15, P=.009), and increased integration between the DMN and multiple frontoparietal networks (DMN - EN, mean difference, 0.53; 95% CI, 0.15 to 0.90, P=.01; DMN - SN, mean difference, 0.55; 95% CI, 0.14 to 0.95, P=.01). A post-hoc exploratory analysis of network recruitment and integration indicated a general increase in DMN integration with other higher-order networks (Figure 3e).
Double-blind randomised controlled trial
Of the 59 MDD patients recruited, 29 were randomly allocated to the escitalopram-arm. Of those, 21 patients (mean age [SD], 40.9 [10.1], 6 [29%] female) were included in this imaging sample. 30 patients were randomly allocated to the psilocybin-arm. Of those, 22 patients (mean age [SD], 44.5 [11.0], 8 [36%] female) were included (Figure 2b).
Psilocybin therapy has greater efficacy than escitalopram for treating depression
Decreased depressive symptom severity was significantly greater under psilocybin than escitalopram, indicating superior efficacy of psilocybin therapy vs. escitalopram (Figure 4). This was confirmed within this neuroimaging sample by a significant arm x timepoint ANOVA interaction for the BDI scores (F, 4.47; P=0.005). FDR-corrected pairwise comparisons relative to baseline were significantly different at 2 weeks (mean difference, -8.73; 95% CI = -13.55 to -3.91, P=0.002), 4 weeks (mean difference, -7.79, 95% CI = -13.62 to -1.95, P=0.013) and at 6 weeks (mean difference, -8.78, 95% CI = -15.58 to -1.97, P=0.013), all favouring the psilocybin-arm (see 32 for full sample).
Increased brain network integration is specific to psilocybin therapy
Confirming our primary hypothesis (Figure 5a-b) and replicating the findings of the open-label trial, brain network modularity significantly reduced following psilocybin therapy (mean difference, -0.39; 95% CI = -0.75 to -0.02, P=0.039). Individuals’ decreases in brain network modularity significantly correlated with greater depression recovery at the 6-week primary endpoint (Pearson, r=0.42, P=.025, one-tailed).
Importantly, this replication was specific to the psilocybin-arm; in the escitalopram group (Figure 5d-e), modularity did not change from baseline to week 6 (mean difference, 0.01; -5% CI -0.35 to 0.33, P=0.945) and there was no significant relationship with changes in BDI scores (Pearson, r=0.08; P=0.361, one-tailed).
Depression recovery correlates with increased cognitive network flexibility.
Next, we examined the dynamic flexibility of the brain’s canonical networks. This finer-grained metric summarises how often brain regions change their community allegiance during the course of an fMRI scan. Post-treatment change in network flexibility were correlated with the changes in BDI score. Specifically, increased EN dynamic flexibility related to greater depression recovery at the 6-week primary endpoint for the psilocybin-arm (Pearson, r=-0.76, P=0.001). Significant relationships predominantly involved the EN, SN and dorsal attention networks (Figure 5c). No significant correlations between BDI and dynamic flexibility were observed in the escitalopram-arm (Figure 5f).