Validation of mature cortical neuron properties.

(A) hiPSCs were differentiated in vitro into glutamatergic cortical neurons over a neuronal progenitor step since the cerebral cortex is a key region for psychedelic effects and psychiatric disorders. Mature 40-day old cortical neurons were glutamatergic (vGLUT2 expression) and expressed cortical layer markers like TBR1, CTIP2 and the neuronal markers NeuN, MAP2 and TAU. cFOS expression is a sign for neuronal activity. GFAP expression indicated a low amount of astrocytes. Scale bar: 50 pm, for GFAP staining 100 pm.

Psilocin-induced increase in BDNF level was 5-HT2A receptor and PKC- and endocytosis mediated and induced activation of m-BDNF/TrkB-associated downstream pathway.

(A-B) Representative image of a neuronal network for pre-treatment condition (Ctrl) and 24 hrs after a 10 min 10 μM short psilocin trigger (Psi). Scale bar: 50 pm, (B) close up: 2 pm. (C) BDNF density was significantly increased 24 hrs after a 10 min 10 μM psilocin trigger (Psi). Four Ctrl cell lines were included in the analysis (Ctrl with N = 613 neurites; Psi with N = 529 neurites). (D) Representative image 24 hrs after the simultaneous treatment with psilocin and ketanserin (Psi + Ket) and single treatment with ketanserin (Ket), scale bar: 50 pm. (E) Ketanserin co-treatment (Psi + Ket) and ketanserin monotreatment (Ket) significantly reduced BDNF density compared to the 24 hrs monotreatment psilocin condition, suggesting a 5-HT2A-mediated process. Ketanserin monotreatment provoked a significant reduction in BDNF density compared to ketanserin co-treatment with psilocin. Two Ctrl cell lines were included in the analysis, each with one biological batch (Ctrl with N = 99 neurites; Psi with N = 102 neurites; Psi + Ket with N = 102 neurites; Ket with N = 102 neurites). (F) Representative image 24 hrs after the simultaneous treatment of 10 μM psilocin with chelerythrine (Psi + Chel) or dynasore (Psi + D), scale bar: 50 pm. (G) Chelerythrine (Psi + Chel, selective PKC inhibitor) or dynasore (Psi + D, inhibition of clathrin-coated vesicle invagination) co-treatment significantly reduced BDNF density compared to the psilocin monotreatment condition (Psi). Two Ctrl cell lines were included in the analysis, each with one biological batch (Ctrl with N = 90 neurites; Psi with N = 90 neurites; Psi + Chel with N = 90; Psi + D with N = 90 neurites). (H) Phosphorylated AKT (pAKT) protein level was increased 72 hrs after psilocin exposure (Psi), reversed by ketanserin co-treatment (Psi + Ket), (I) both effects were not significant. Ctrl cell line 3 was included in the analysis (Ctrl with N = 4 data points, Psi with N = 4 data points, Psi + Ket with N = 4 data points), each with two biological batches. For all analyses Kruskal-Wallis-Test for independent samples was calculated. Post hoc Wilcoxon rank sum test. Bonferroni-correction, adjusted p < .05, mean ± SD. Significance levels against the respective control and for multiple group comparisons are *p < .05.

Psilocin displays fast and enduring changes of the genetic landscape.

(A) Enrichment in differentially expressed genes associated with synapse formation, neuronal plasticity and axonogenesis GO terms 1 day after psilocin administration, and the effect 3 days later. (B) Psilocin induced within 24 hours a first wave of upregulation of selected GO genes based on DESeq2 normalized counts. Log2fold changes of each differentially expressed genes between two conditions are shown as dots. Zscores indicate if more genes in the respective GO term are upregulated or downregulated also indicated by height and color of the bars (C) TPM-normalized mean (red line) of psilocin-induced temporal expression pattern of genes belonging to the indicated GO. Shaded area is indicating SD. (D) Chord plots showing significant genes appearing in at least 4 GO terms (of 6 selected GO terms) 1 day after psilocin administration and appearing in at least 3 GO terms 3 days after psilocin administration. (E). Heat map (z-scaled normalized counts) showing expression of immediate early genes and AMPA/ NMDA receptor genes throughout psilocin administration. (F) Differentially expressed genes, that are up-/ downregulated upon psilocin treatment, showed a reversed effect upon co-treatment with ketanserin. Mean TPM-normalized expression values are shown (red line), shaded area indicates S.D. BP, biological process; CC, Cellular Compartment; KEGG, Kyoto Encyclopedia of Genes and Genomes. TPM, transcripts per kilobase million. Significance levels against the respective control ns: p-adj. > .05, *: p-adj. <= .05, **: p-adj. <= .01, ***: p-adj. <= .001.

Psilocin induced neurite branching.

(A-C) Cells were transduced with AAV CamKIIa p-hCHR2(134a)-mCherry. (A) Representative image for pre-treatment condition (Ctrl) and 48 hrs after a 10 min short psilocin trigger (48 hrs) for mCherry staining showed an increase in neurite intersections for the latter. (B) Singular analysis: 48 hrs after a 10 min 10 μM psilocin trigger the number of primary (25 μm distance from soma) neurite intersections significantly increased compared to the untreated control condition (Ctrl). The number of intersections at 50 μm distance from soma significantly increased 24 hrs and 48 hrs after a 10 min 10 μM psilocin trigger. (C) Significant changes for the calculated total neurite length and for the total number of neurite intersections significantly increased 24 hrs and 48 hrs after a 10 min 10 μM psilocin trigger. Two Control cell lines with two biological batches (Ctrl with N = 43 - 48 neurites; 24 hrs with N = 47 - 48 neurites; 48 hrs with N = 31 - 35 neurites) were included. (D) Sholl analyses summary for results shown in figure (B). For all analyses Kruskal-Wallis-Test for independent samples was calculated. Post hoc Wilcoxon rank sum test. Bonferroni-correction, adjusted p < .05, mean ± SD. Significance levels against the respective control are *p < .05.

Psilocin-induced increase in synaptic strength and synaptogenesis.

(A) Total number of evoked action potentials (eAPs) significantly increased at day 7 after 24hrs permanent psilocin administration (Psi 24hrs; day 7) and was increased 7 days after a 10 minutes psilocin trigger (Psi 10min; day 7). AP amplitudes stayed constant. Representative traces for total number of APs. One Ctrl cell line with 2 biological batches was included in the analysis. Ctrl 3 with N = 28 cells, Psi 10min; day 7 with N = 29 cells, Psi 24hrs; day 7 with N = 34 cells. (B) Increase in sEPSCs amplitude and frequency 7 days after 10 minutes and 24hrs permanent psilocin administration. One Ctrl cell line with 2 biological batches was included in the analysis, Ctrl 3 with N = 19 cells, Psi 10min; day 7 with N = 20 cells, Psi 24hrs; day 7 with N = 21 cells. (C) Increase in mEPSCs amplitude 6 days after permanent 10 μM psilocin administration, Ctrl 3 with N = 16 cells, Psi 24hrs; day 7 with N = 14 cells. Increase in mEPSCs amplitude 10 days after 24hrs permanent psilocin administration, Ctrl 3 with N = 15 cells, Psi 24hrs; day 11 with N = 20 cells. (D) Significant increase in sEPSCs amplitude 6 days after 96hrs permanent 10 μM psilocin administration (Psi 96hrs; day 10). Ctrl 3 with N = 18 cells, Psi 96hrs; day 10 with N = 13 cells. Representative traces for control and psilocin condition. For all experiments Mann-Whitney-U-test for independent samples was calculated, mean ± SEM. Significance levels against the respective control are *p < .05. (E-F) Representative dendritic PSD-95 and synapsin staining for the untreated condition (Ctrl), 4 days after permanent 10 μM psilocin administration (4 days) and 10 days after 4 days permanent 10 μM psilocin administration (10 days). (E) Scale bar: 50 μm (F), closeup: 10 pm. Trend in synapsin (G) density and (H) intensity increase 4 days after permanent psilocin administration. (I) PSD-95 particle per neurite length and (J) intensity per area and (K) synapsin/ PSD- 95 co-localization were significantly increased 4 and 10 days after 4 days permanent 10 μM psilocin treatment compared to an untreated control condition. (G-K) Ctrl with N = 180 neurites, 4 days with N = 180 neurites, 10 days with N = 180 neurites. Two control cell lines each with two biological batches were included. For all analyses Kruskal-Wallis-Test for independent samples was calculated. Bonferroni- correction, adjusted p < .05, mean ± SD. Significance levels against the respective control are *p < .05.

Primary antibodies used for immunocytochemistry (ICC).

Secondary antibodies used for ICC.

Fluorescent probes.

Primary antibodies used for Western blotting (WB).

IR-dye conjugated secondary antibodies used for WB.

Secondary antibodies were diluted 1:15,000.

Primers for RT-PCR.

Validation of iPSC and neural progenitor properties and HTR2A gene expression.

(A). Timeline: hiPSCs were differentiated in vitro into glutamatergic cortical neurons over a neuronal progenitor step. Maturation of cortical neurons took about 6 weeks. Scheme created in BioRender, Ladewig, J. (2025) https://BioRender.com/v84c084 (B) iPSCs expressed pluripotent marker SOX2 and OCT4. Neuronal progenitors expressed SOX2, PAX6 and NESTIN and were negative for FOXA2 as a midbrain marker. Scale bar: 50 pm. (C) RT-PCR analysis revealed expression of neuronal subtype markers, neurotransmitter receptors, neuronal activity markers, neuronal markers and synapse- associated genes. (D) RT-PCR of the HTR2A gene (also expressed in commercialized samples of fetal and adult brain) confirmed expression of 5-HT2A-R. Scale bar: 50 pm.

Influence of different treatment conditions on BDNF level.

(A) Representative image of a neuronal network for pre-treatment condition (DMSO) and 24 hrs after a 10 min 10 nM, 100 nM, 1 μM or 10 μM trigger showed the strongest increase in BDNF density for the 10 μM treatment condition. Scale bar: 50 pm. (B) BDNF density was exclusively significantly increased 24 hrs after 10 min 10 μM Psilocin trigger (10 pM). Ctrl cell line 1 was included in the analysis, with one biological batch (Ctrl with N = 76 neurites; 24 hrs after a 10 min 10 nM (N = 75 neurites), respectively 100 nM (N = 75 neurites), respectively 1 μM (N = 75 neurites), respectively 10 μM (N = 75 neurites). (C) Representative image of a neuronal network 24 hrs after an “artificial” (10 min 10 pM) and a “physiological” (6 hrs 100nM) treatment condition. Scale bar: 50 pm. (D) BDNF density was significantly increased 24 hrs after a 10 min 10 μM psilocin trigger (10 min 10 pM) and 24 hrs after a 6 hrs 100 nM psilocin stimulation (6 hrs 100 nM). BDNF density was also significantly increased for the” artificial” compared to the “physiological“ stimulation. Two Ctrl cell lines were included in the analysis, each with one biological batch (Ctrl with N = 99 neurites; 10 min 10 μM with N = 102 neurites; 6 hrs 100 nM with N = 102 neurites). (E) Representative image of a neuronal network of the Ctrl condition and 24 hrs after a 10 min 10 μM psilocin trigger showed an increase in axonal TAU+ and dendritic MAP2+ BDNF density. Scale bar: 50 pm. (F) BDNF density was significantly increased for axonal TAU+ and dendritic MAP2+ BDNF density for the psilocin condition. One cell line (Ctrl 3) was included in the analysis, with two biological batches (Ctrl with N = 90 neurites; Psi with N = 90 neurites). (G) BDNF density was significantly increased 24 and still 48 hrs (24 hrs; 48 hrs) after a 10 min 10 μM psilocin trigger. Three Ctrl cell lines were included in the analysis (Ctrl with N = 520 neurites; 24 hrs with N = 478 neurites; 48 hrs with N = 170 neurites). (H) Endogenous phosphorylated (activated) TrKB (pTrkB) receptor, pro- BDNF and m-BDNF protein level increased 24 hrs after psilocin (Psi) administration. (I) Representative WB: phosphorylated AKT (pAKT) protein level, indicating an activation of pro survival AKT pathway, was increased 24 hrs after psilocin (Psi) exposure. For all analyses Kruskal-Wallis-Test for independent samples was calculated. Post hoc Wilcoxon rank sum test. Bonferroni-correction, adjusted p < .05, mean ± SD. Significance levels against the respective control are *p < .05.

Psilocin displays fast and enduring changes on GO terms related to learning, memory and cognition.

(A) Heat map (z-scaled normalized counts) showing expression of genes that are grouped in the selected GO term related to “cognition”. Genes that are also present in the GO term related to “learning and memory” are stated as “present”, indicating great overlay of the genes of the two GO terms. (B-C) Psilocin induced within 24 hours a first “shaking” wave of selected GO genes based on DESeq2 normalized counts. After 3 days those GO genes normalized. (B) Log2fold changes of each gene between two conditions are shown as dots. Zscores indicate if more genes in the respective GO term are upregulated or downregulated indicated by height and color of the bars. (C) TPM-normalized mean (red line) of psilocin-induced temporal expression pattern of genes belonging to the indicated GO. Shaded area is indicating SD. TPM, transcripts per kilobase million. Significance levels against the respective control ns: p-adj. > .05, *: p-adj. <= .05, **: p-adj. <= .01, ***: p-adj. <= .001.

Psilocin-induced increase in synaptic strength.

(A) Increase in mEPSCs amplitude at day 7 after 24hrs permanent psilocin administration. Ctrl 2 with N = 16 cells, Psi 24hrs; day 7 with N = 14 cells. (B) Significant increase in sEPSCs amplitude 6 days after 96 hrs permanent 10 μM psilocin administration (Psi 96hrs; day 10). Ctrl 2 with N = 8 cells, Psi 96hrs; day 10 with N = 9 cells. For all experiments Mann-Whitney-U-test for independent samples was calculated, mean ± SEM. Significance levels against the respective control are *p < .05.

Simplified schematic of signal transduction pathways proposed to mediate psychedelic-induced neuroplasticity in the cortex.

Serotonergic psychedelics activate the 5-HT2A receptor (5-HT2A-R), triggering G-protein-coupled PLC and PLA2 pathways, which, similar to m-BDNF-TrkB activation, increase intracellular Ca2+ and PKC activity. BDNF-TrkB signaling further activates mTOR and NMDAR-dependent plasticity via the PI3K/AKT pathway. Psychedelic-induced glutamate release promotes mTOR activation and BDNF release through AMPARs, enhancing glutamatergic signaling and IEG (e.g., c-FOS) expression. 5- HT2A-R also interacts with PSD-95, while PI3K, ERK1/2, and CaMK activate CREB, stimulating BDNF transcription. Image created using BioRender.com, adapted from24 and used in57 . Abbreviations: 5-HT, serotonin; AKT, protein kinase B; AMPAR, a-amino-3-hydroxy-5-methyl-4- isoxazolepropionic acid receptor; BDNF, brain-derived neurotrophic factor; CaMK, Ca2+/calmodulin- dependent kinase; CREB, cAMP response element-binding protein; ERK1/2, extracellular signal- regulated kinase 1/2; IEG, immediate early gene; MAPK, mitogen-activated protein kinase; mTOR, mammalian target of rapamycin; NMDAR, N-methyl-D-aspartate receptor; PI3K, phosphoinositide 3- kinase; PKC, protein kinase C; PLA2, phospholipase A2; PLC, phospholipase C; PSD-95, postsynaptic density protein 95; TrkB, tropomyosin receptor kinase B.