Exploring Psilocybe spp. mycelium and fruiting body chemistry for potential therapeutic compounds

Introduction

Psilocybin and psilocin are the primary hallucinogenic alkaloids found in psychedelic mushrooms (Wieczorek et al., 2015) from the basidiomycete genera Psilocybe, Inocybe, Panaeolus, Gymnopilus, Copelandia, Hyboloma, Pluteus, Conocybe, and Panaeolina [also referred to as hallucinogenic, entheogenic, “magic”, medicinal, neurotropic, psychoactive, sacred, and saint mushrooms (Guzmán, 2008)]. Psilocin exerts a neurologic effect by binding to serotonin subtype 2A receptors (5-HT_2A) and producing neuropsychological effects, including oceanic boundlessness, anxious ego dissolution, visionary restructuring, auditory alterations, and reduction of vigilance; collectively, these comprise the “psychedelic experience” (Hasler et al., 2004). On oral administration in humans, psilocybin is rapidly converted to psilocin via intestinal alkaline phosphatase, where nearly 100% of the bioactive molecule is found in circulation 30 min to 6 h after ingestion (Dinis-Oliveira, 2017). Psilocin is able to pass through the blood–brain barrier (Dinis-Oliveira, 2017) and bind to 5-HT_2A receptors in the central nervous system with nanomolar affinity (Rickli et al., 2016), leading to downstream hallucinogenic effects (López-Giménez and González-Maeso, 2017).

Psilocybin is synthesized using tryptophan derived either directly from amino acids or from the shikimate acid pathway as its primary substrate, with the latter being the most common in both fruiting bodies and mycelia (Mahmood, 2013). In brief, L-tryptophan is decarboxylated to tryptamine by the psilocybin decarboxylase (PsiD) enzyme where tryptamine is then hydroxylated at position 4 to yield 4-hydroxy-tryptamine by psilocybin hydroxylase (PsiH). Subsequently, 4-hydroxy-tryptamine may be phosphorylated by psilocybin kinase to form norbaeocystin, which is then twice-methylated at the amine group by the psilocybin methyltransferase enzyme to form psilocybin (PsiM) (Fricke et al., 2017; Torrens-Spence et al., 2018) (Figure 1). Further methylation of the amine group by PsiM yields aeruginascin, a trimethyl ammonium analog of psilocybin found originally in Inocybe aeruginascens (Jensen et al., 2006). The above-mentioned indole alkaloids may exhibit biological activity and have been detected in samples of Psilocybe cubensis.

FIGURE 1

Figure 1 Biosynthesis pathway of psilocybin in Psilocybe cubensis. (Adapted from Fricke et al., 2017). PsiD, psilocybin decarboxylase; PsiH, psilocybin hydroxylase; PsiM, psilocybin methyl transferase; PsiK, psilocybin kinase; SAM, S-adenosyl-L-methionine; SAH, S-adenosyl homocysteine.

Each stage of the fungal life cycle creates a unique fingerprint of fungal metabolites which may have their own unique therapeutic potential. For example, Hericium erinaceus (lion’s mane) mycelia preferentially expresses p450 enzymes involved in the synthesis of diterpenoid “erinacines” (Li et al., 2018), whereas fruiting bodies exhibit higher concentrations of polyketide “hericenones”, but a low erinacine content (Ma et al., 2010). Erinacines have been cited as being anti-inflammatory and reactive oxygen species scavengers, whereas hericenones are promoters of nerve growth factor (NGF) secretion and neurite outgrowth (Phan et al., 2014). Mycelium has been targeted for the development of pharmaceuticals with culture broths as a means of increasing fungal biomass for the purpose of isolating bioactive compounds (Toth et al., 1983; Ikekawa et al., 1985; Kawagishi et al., 1996; Okuyama et al., 2004; Park et al., 2004; Mygind et al., 2005; Kawagishi et al., 2006; Jung et al., 2008; Singh et al., 2011; John et al., 2013; Phan et al., 2015). Some examples of these compounds are an antibiotic compound from Pseudoplectania nigrella (Mygind et al., 2005), immunomodulatory compounds from Aspergillus nidulans (Singh et al., 2011), and an antitumor agent isolated from Flammulina velutipes (Ikekawa et al., 1985).

The multivariate analysis of metabolomic data performed on Ganoderma lucidum (reishi or lingzhi) revealed the differential accumulation of metabolites, with 48, 25, 364, and 25 metabolites being unique to the fermentation broth, mycelium, fruiting body, and spore powder, respectively (Xie et al., 2020). Several cultivation methods, substrates, and media may be utilized to produce the fungal material, with expected differences found in the metabolomic profiles. In this study, we analyze the difference between liquid culture and grain spawn substrate for mycelial colonization. Although the genome of P. cubensis has recently been drafted, with 32 contigs mapping 97.6% of the genome by benchmarking universal single-copy orthologs (BUSCO) completeness score (McKernan et al., 2021), to date, no metabolomic mapping of Psilocybe spp. growth stages has been published. Psilocybe spp. mycelia have markedly lower concentrations of psilocybin and psilocin (Gartz and Moller, 1989; Gotvaldová et al., 2021); however, unique compounds were detected in a sample of mycelia when compared with fruiting bodies of the same species using ultra-high-pressure liquid chromatography coupled with mass spectrometry (UPLC-MS) data (Bradshaw Alexander et al., 2022), and, therefore, a more comprehensive analysis of targeted and untargeted fungal metabolite profiles from both fruiting body and mycelial datasets is needed.

In this article, we employed ultra-high-pressure liquid chromatography coupled with ultraviolet and visible light detector (UHPLC-UV/VIS; TSQ Altis™) and ultra-high-pressure liquid chromatography coupled with high-resolution mass spectrometry (UHPLC-HRMS; Orbitrap Exploris™) data coupled with multivariate analysis to compare the targeted and untargeted (metabolomic) chemical profiles of Psilocybe spp. mycelia grown in liquid culture, mycelia grown on rye grain substrate, and fruiting bodies harvested from bulk substrate. Evidence of the preferential accumulation of indole alkaloids and also the resulting correlations among other fungal metabolites in mycelial samples vs. fruiting bodies is essential in not only the development of Psilocybe fruiting bodies in mycotherapy, but also that of Psilocybe mycelia as a lower psychedelic potential alternative.

Materials and methods

Sample collection

All fungal strains originated as externally sourced spore prints. The spores were germinated and the dikaryotic mycelium was isolated to either single-sector isolates (P-AZ-1, 5051-B.1, P-CY-1, BP-STAR, and BP-A-STAR) or multi-strain isolates (BP-A-MUL, BP-C-MUL, AV-MUL, and LK-MUL). The strains P-AZ-1, P-CY-1, BP-STAR, BP-A-STAR, BP-A-MUL, and BP-C-MUL were identified through the amplification and Sanger sequencing of the ITS1–5.8S–ITS2 region. DNA extraction was performed at Numinus Bioscience using the Invitrogen™ PureLink™ Genomic Plant DNA Purification Kit, and amplification and sequencing was performed at the Slot Laboratory at Ohio State University, Columbus, OH, USA. The strains AV-MUL, LK-MUL, and 5051-B.1 were identified by confirming the morphological characteristics of the fruiting bodies. All fungal strains were preserved as mycelium cultures stored at 4°C in our parental cell line library.

The mycelial cultures were grown on solid yeast malt extract agar (YMEA) or potato dextrose agar (PDA) in 10-cm Petri dishes, and the mycelia showed no preference toward either substrate. The strains were incubated at 26°C (P. cubensis) or 21°C (Psilocybe allenii and Psilocybe cyanescens) until fully colonized. The cultures were visually inspected for sectoring and, in some cases, continuously subcultured until a single-sector isolate was achieved, while in other cases, cultures were not subcultured and were kept as multi-strain isolates. The mycelia and agar slices were transferred to sterilized rye kernel grain spawn jars and left to fully colonize at 26°C. Once fully colonized, the grain spawn jars were used to inoculate grain spawn bags and left to fully colonize at 26°C. These fully colonized spawn bags were then used to inoculate the bulk substrate bags at a ratio of 1 : 3. On full colonization of the bulk substrate bags, fruiting and harvesting occurred. The fruiting bodies were considered mature and ready for harvest once the veil was separated from the cap margin (Supplementary Figure 1).

The liquid cultures were grown by transferring mycelia + agar slices to sterilized malt extract broth in ventilated culture flasks and incubated at 26°C (P. cubensis) or 21°C (P. allenii and P. cyanescens) for approximately 21 days. The bulk substrate, consisting of 600 g of coconut coir, 280 g of vermiculite, and 124 g of gypsum at field capacity for moisture, once fully colonized, was subjected to fruiting conditions. The fruiting conditions consisted of a 12-hour photoperiod at a Kelvin color temperature of approximately 6,400 and a decrease in temperature to 21°C. The bag of bulk substrate, which had a 0.5-micron filter patch to allow for gas exchange, acted as a miniature greenhouse, maintaining moisture and humidity within the bag to allow for sporocarp formation. The sporocarps were harvested approximately 14 days after being subjected to fruiting conditions.

The mushrooms and grain mycelia samples were dehydrated and ground, whereas liquid mycelial cultures were collected from the cell culture flasks and freeze-dried using Toption and Harvest Right Vacuum Freeze-Dryers. For both Orbitrap Exploris and TSQ Altis, the samples were weighed out to 25 mg± 5 mg, dissolved in 5 mL of methanol, and mixed at 1,800 rpm in a Benchmark Scientific BenchMixer™ multi-tube vortexer, followed by a 5-min sonication. The samples were then centrifuged at 3,100 x g for 15 min at −8°C. The supernatant was collected, and the material was once again extracted following the same procedure, after which the supernatants were pooled.

Screening of targeted compounds

For the targeted analysis, samples were analyzed using a TSQ Altis Triple Quadrupole Mass Spectrometer coupled with UHPLC equipped with heated electron spray ionization (H-ESI) for data acquisition. Data points were collected for 25 min for each sample. The scan was carried out by way of selected reaction monitoring (SRM) with positive polarity, collecting 13.3 points per peak. For each compound, the following parameters were included in the SRM: retention (RT), precursor [mass-to-charge ratio (m/z)], product(s) (m/z), collision energy (V), and RF lens (V). The resolution for both Q1 and Q3 was 0.7 FWHM (full width at half maximum). A spiked sample was prepared along with every sample run as part of the quality control process for analytical testing. Each sample run included three replicates for a sample (A, B, and C) and also a procedural blank consisting of methanol and internal standards at the same dilution, which was subjected to the same extraction procedure, and a solvent blank consisting of methanol only. The level of detection (LOD) was 0.00002 mg/g and the limit of quantification (LOQ) varied across the standardized compounds.

The fungal samples were given a unique ID (Supplementary Table 1) along with a certificate of analysis for standardized compounds, including: alpha-glycerylphosphorylcholine (AGPC), arginine (ARG), baeocystin (BAEO), carnitine (CARN), choline (CHOL), ergothioneine (ERGO), glutamate (GLU), glutamine (GLN), histidine (HIS), methoxy-tryptamine (MTRY), N-acteylglucosamine-anhydro (NAGA), nicotinamide (NICO), nicotinic acid (NICA), norbaeocystin (NORB), pantothenic acid (PANT), phenylethylamine (PEA), psilocin (PSC), psilocybin (PSB), serotonin (SERO), trimethylglycine (TMG), trimethyllycine (TML), tryptophan (TRP), aeruginascin (AERU), norpsilocin (NORP), and 4-hydroxy-N,N,N-trimethyltryptamine (4HTT). The samples were run as three technical replicates (n = 3) and averaged. All the targeted data are presented in mg/g.

Screening of untargeted compounds

For the untargeted analysis, the samples were analyzed using an Orbitrap Exploris 120 equipped with H-ESI with a full-scan method. The data points were collected for 22 min for each sample. The resolution of the mass spectrometer was set to 15,000, with a scan range of 80 m/z–700 m/z. The instrument was set to positive polarity. Each sample run included three replicates for a sample (A, B, and C), in addition to a procedural blank and a solvent blank.

The raw data files were returned from Xcalibur™ (Thermo Fisher Scientific), which were then converted to an mzML file format using a ProteoWizard MS converter and loaded into MZmine2.53. The mass detection of raw data was done through a local minima search with the noise level set to 4.0E + 06, and peaks with retention times (RTs) between 0 min and 12 min were included. Throughout data preprocessing, an m/z tolerance of 0.005 was set. An automated data analysis pipeline (ADAP) chromatogram builder was used to generate the chromatograms. The minimum threshold was set to 3.0E + 06, and the minimum peak intensity was set to 4.0E + 06. The chromatogram deconvolution consisted of visualizing the individual spectral scans to capture unique peaks, changing the chromatographic threshold to 20%, the minimum RT range to 0.05 min, the minimum relative height to 15%, the minimum absolute height to 5.0E + 06, and the minimum ratio of the peak top edge to 1. An isotopic peak grouper was used with a maximum charge of 1 and a monotonic shape was applied to the peak. A join aligner allowed buckets to be made to create a two-dimensional aligned chromatogram, and gap-filling was conducted using the peak finder function with a 5.0% intensity tolerance and 0.05 RT tolerance. The joined and gap-filled chromatogram was then exported to CSV with the peak ID, RT, and m/z value. Data pruning consisted of removing the features that did not reach the minimum threshold of 4.0E + 06 in samples or that were abundantly present in the procedural blanks. Peak annotation was performed in Xcalibur using the compound discovery application, where percentage identity scores were created based on MS2 spectra when compared with internal and external compound databases.

Statistical analysis

Principal component analysis (PCA) is an unsupervised clustering method that reduces the dimensionality of multivariate data while preserving most of their variance. All samples were collected from a unique fruiting bag, grain spawn bag, or mycelial culture, and were assigned a unique internal lot number and barcode. Studying the extracted fungal material using the H-ESI full scan revealed a large cloud of mass features, the data structure of which was parsed from the nose using PCA. On generating the PCA plots, all scree plots (data not shown) were observed to confirm that inflection occurs after PC1 and PC2. The analysis of the mycelia, grain spawn mycelia, and fruiting bodies, and also the comparisons between strains of P. cubensis fruiting bodies, used an unsupervised clustering method along with ellipses, and all data matrices were Pareto scaled before analysis.

A univariate analysis of indole alkaloid and minor compound content was used to determine the relative difference in mycochemicals among the groups and further evaluate any emerging correlations in the dataset. ANOVA and Tukey’s honestly significant difference (HSD) were used to generate F-statistics, fold differences, and to identify significant differences among the groups based on targeted or untargeted compounds. For all data groups, a Shapiro–Wilk test was used to ensure normality assumptions were met, and all data groups contained ≥ five replicates.

Results

Targeted analysis of Psilocybe mycelium, grain mycelium, and fruiting bodies

The univariate analysis of indole alkaloids and other fungal metabolites revealed that fruiting bodies preferentially accumulate aeruginascin, alpha-glycerophosphocholine, carnitine, choline, ergothioneine, glutamate, glutamine, histidine, norbaeocystin, pantothenic acid, phenylethylamine, psilocin, psilocybin, trimethylglycine, trimethyllysine, and tryptophan (Figure 2) when compared with mycelium and grain mycelium. The average psilocybin (± SE) content in fruiting bodies was 9.913 mg/g ± 0.389 mg/g, 0.041 mg/g ± 0.014 mg/g for mycelia, and 0.047 mg/g ± 0.023 mg/g for grain mycelium samples. In contrast, no significant differences were found in the concentrations of 4-hyrdoxy-N,N,N-trimethyltryptamine, arginine, methoxy-tryptamine, N-acetyl glucosamine anhydro, nicotinamide, nicotinic acid, and norpsilocin between fruiting bodies and mycelia (Supplementary Table 2).

FIGURE 2

Figure 2 Box plot representing the accumulation of targeted compounds analyzed using TSQ Altis as determined by a standard curve. The data were collected from 14 fruiting bodies, seven mycelia and six grain mycelia samples. * (p < 0.05) in ANOVA and Tukey’s HSD. 4HTT, 4-hydroxytryptamine; AERU, aeruginascin; AGPC, α-glycerylphosphorylcholine; ARG, arginine; BAEO, baeocystin; CARN, carnosine; CHOL, choline; ERGO, ergotheioneine; GLN, glutamine; GLU, glutamic acid; HIS, histidine; MTRY, methoxytryptamine; NAGA, N-acetylglucosamine-anhydro; NICA, nicotinamide; NICO, nicotinic acid; NORB, norbaeocystin; NORP, norpsilocin; PANT, panthotenic acid; PEA, phenylethylamine; PSB, psilocybin; PSC, psilocin; TMG, trimethylglycine; TML, trimethyllysine; TRP, tryptophan.

The preparations of P. cubensis mycelium and grain mycelium were compared for their targeted compound content. Targeted analysis revealed the fruiting body chemistry to be chemically distinct from that of mycelium and grain mycelium, as indicated in a PCA score plot (Figure 3A). The loading plot (Figure 3B) reveled psilocybin influencing the clustering of fruiting bodies while α-GPC and trimethylglicine the mycelia. However, based on targeted data, no significant difference was found between grain mycelia and mycelia. Principal components 1 and 2 captured 75.3% of the total variation in the dataset. Hierarchical cluster analysis of the samples, using an unweighted pair group method with arithmetic mean (UPGMA) and a Euclidean distance function, revealed similar patterns in terms of the relatedness of samples in targeted compounds data and also of the relatedness of the compounds (Supplementary Figure 2).

FIGURE 3

Figure 3 Score (A) and loading (B) plots of Psilocybe cubensis fruiting bodies (scaled and unsupervised). The mycelium (red) consisted of seven Psilocybe mycelium samples of Psilocybe cyanescens, Psilocybe allenii, and P. cubensis. The grain mycelium (green) consisted of six samples of P. cubensis and P. cyanescens. The fruiting bodies (black) consisted of 18 samples from P. cubensis. The samples were extracted and quantified using TSQ Altis, and the data were averaged over three technical replicates (n = 3) for each data point.

Psilocybin is present in mycelia and grain mycelia, but the isolation of indole alkaloids from fungal materials will be most effective by the extraction of fruiting bodies. There was no significant difference in psilocin, norpsilocin, and aeruginascin content among the life stages. An ANOVA of the individual compounds revealed that levels of alpha-glycerylphosphorylcholine (α-GPC), N-acetylglucosamine, and trimethylglycine accumulation were found to be higher in mycelia and grain mycelia than in fruiting bodies, according to Pareto-scaled data (p < 0.05). It is known that α-GPC is a class of phosphocholines that delivers choline to the brain and is a precursor to acetylcholine, and it is also a neurotransmitter and modulator involved in arousal, attention, memory, and motivation (Sigala et al., 1992). In addition, α-GPC has been implicated in enhanced physical motor performance (Bellar et al., 2015; Marcus et al., 2017) and enhanced cognition (Lee et al., 2017), while N-acetylglucosamine is the monomeric unit of the polymer chitin which is abundant in medicinal mushroom preparations (Chen et al., 2010). Finally, trimethylglycine (betaine) is essential to the methylation of methionine in the homocysteine methyltransferase cycle (Brütting et al., 2021) and also to cellular osmoregulation (Lever and Slow, 2010). Considering the low psilocybin content of mycelia and the accumulation of beneficial fungal metabolites in them, Psilocybe mycelia may be good candidates for mycotherapeutic development.

Untargeted analysis of Psilocybe mycelium, grain mycelium, and fruiting bodies

The preprocessing of Orbitrap Exploris 120 full-scan data yielded an aligned data matrix consisting of 978 mass features across samples of all types. Each feature is denoted by an RT and m/z (Figure 4). Only 31.2% of variability in the dataset was captured by principal components 1 and 2. However, the enhanced separation of individual samples was observed within and among the groups. The fruiting body chemistry is different from that of mycelium and grain mycelium (Figure 4A). A greater separation of mycelium and grain mycelium confidence was observed; however, these differences were insubstantial. Hierarchical cluster analysis using an unweighted pair group with arithmetic mean (UPGMA) of samples and a Euclidean distance function revealed similar patterns in terms of the relatedness of the samples in targeted compounds data and also of the relatedness of metabolites (Supplementary Figure 3). The clustering revealed that fruiting bodies, mycelia, and grain mycelia exhibit a unique suite of fungal metabolites.

FIGURE 4

Figure 4 Score (A) and loading (B) plots of Psilocybe fruiting bodies (non-scaled and unsupervised). The mycelium (red) consisted of seven Psilocybe mycelium samples of Psilocybe cyanescens, Psilocybe allenii, and Psilocybe cubensis. The grain mycelium (green) consisted of six samples of P. cubensis and P. cyanescens. The fruiting bodies (black) consisted of 18 samples from P. cubensis. The principal components were determined based on the loaded 978 mass features returned by Orbitrap Exploris H+ESI. The mass features were pruned with method and procedural blanks and then they were Pareto scaled.

The separation within fruiting bodies appeared to be influenced by the mass features RT5.14_285.10045 and RT1.16_360.15077, which are precursor peaks for psilocybin and trehalose, respectively, as determined by MS2 spectra. These mass features appeared to polarize the data in the PC1 loadings plot when the mycelia and fruiting bodies were compared. The further isolation and quantification of fungal sugars such as α-glucans and β-glucans (Pérez-Vendrell et al., 1995), and of monosaccharides and disaccharides may be beneficial in future studies. These results may further expand on the potential inverse relationship between Psilocybe psilocybin and carbohydrate content.

An analysis of Morchella (morel) based on UPLC-ESI-MS/MS data revealed that mycelium and fruiting bodies also displayed different metabolomes with upregulated and downregulated metabolic pathways (Yang et al., 2021). A comprehensive study of Isaria cicadae (cordyceps) revealed that growth media and environmental conditions influence mycelial and fruiting body metabolomes, such that fruiting bodies may accumulate higher antioxidant compounds due to their exposure to oxygen (He et al., 2019). This pattern may also hold true for Psilocybe species and should be considered in the further development of research on psychedelic mushrooms.

Discussion

Although the metabolites from the indole alkaloid synthesis pathway have an established pharmacology, their bioactivity is not limited to this class of compounds, as other selected fungal metabolites, amino acids, lipids, and sugars may positively or negatively modulate indole alkaloid involvement in various pathologies (Brown and Gordon, 2003; Evidente et al., 2014). The crude extract preparations of Psilocybe mushrooms may exhibit synergistic effects when compared with pure psilocybin due to the binding of the orthosteric site on 5-HT_2A receptors (Glatfelter et al., 2022), or the potential affinity for off-target receptors, which may be additive, such as inhibiting the degradation or reuptake of endogenous serotonin. Allosterism has been described for serotonin receptors (Fasciani et al., 2020), and this site may be targeted by metabolites that exist in fruiting body or mycelial extracts. This may be similar to the relationship between Δ⁹-tetrahydrocannabinol (THC) and cannabidiol (CBD) from Cannabis, where CBD is a negative allosteric modulator of THC at the cannabinoid receptor 1 (Laprairie et al., 2015). The synergism between the mushroom's active compound and other fungal metabolites may enhance the therapeutical effect described with pure psilocybin. Another example is the anti-depressant synergistic effects of active compounds from Hypericum perforatum (St. John’s Wort) (Simmen et al., 2001). Hypericin and other active compounds showed inhibitory activity at corticotrophin-releasing factor receptors as well as gamma-aminobutyric acid type A GABA_A, serotonin, and opioid receptors.

Active sites of tryptamine derivatives in Psilocybe spp. match those of serotonin receptors (5-HT_2A) (Almaula et al., 1996). The binding of 5-HT_2A G-protein-coupled receptors (GPCRs) and the resulting hallucinogenic effects elicited by psilocybin/psilocin are determined by distinct downstream signaling events, biased agonism (Karaki et al., 2014; Pottie et al., 2020), the binding of off-target receptors, and the formation of heterocomplexes between the 5-HT_2A and D_2L and mGlu2/3 receptors (Albizu et al., 2011; Inserra et al., 2021). The benefits of psilocybin are dependent on a delicate cascade of signaling spurred by the activation of GPCRs and further recruitment of G subunits and β-arrestins (López-Giménez and González-Maeso, 2018). The effects of Psilocybe metabolite allosteric interactions at the level of the receptor and potential synergism by binding of dopamine and glutamate receptors have yet to be investigated. The cumulative effects of crude extracts substantiate why traditional medicine practitioners believe these preparations have better therapeutic activity than the compounds isolated from them (Koltai and Namdar, 2020).

The levels of indole alkaloids psilocybin, baeocystin, norbaeocystin, and tryptophan were significantly higher in fruiting bodies than in the mycelium. Overall, mycelium and grain mycelium could not be distinguished from each other in the targeted and untargeted analysis; however, fruiting body chemistry was found to be significantly different in the univariate analysis. The retention times and mass-to-charge ratios were used to identify dominant peaks in MS2 spectra; however, many mass features remain unknown in this analysis. As Psilocybe preparations become more prevalent in the recreational and psychedelic-assisted therapy space, a well-assembled fungal mass spectra library to identify potential bioactive or harmful compounds will be a useful tool for researchers and regulators. Mycelia and grain mycelia cultures are easy to expand, and fungal metabolites in them remain in high-enough proportions, with the level of psilocybin accumulation remaining low. In this study we provide the analysis necessary to propose Psilocybe mycelia as a non-intoxicating mycotherapeutic with chemistries significantly different than those of fruiting bodies.

Data availability statement

The data presented in the study are deposited in the Center for Open Science repository and can be accessed at https://osf.io/jctyk.

Author contributions

AW: Conceptualization, Data curation, Formal Analysis, Investigation, Writing – original draft. MB: Data curation, Investigation, Methodology, Writing – original draft. SN: Data curation, Investigation, Methodology, Writing – original draft. JA: Conceptualization, Project administration, Supervision, Writing - review & editing. SS: Conceptualization, Funding acquisition, Project administration, Resources, Supervision, Writing – review & editing.

Funding

The author(s) declare financial support was received for the research, authorship, and/or publication of this article. All contributing authors were employed by Numinus Bioscience. This study was privately funded by Numinus Bioscience. Numinus Bioscience had provided all funding to carry out research activities and materials for the study but had no role in in the study design, collection, analysis, interpretation of data, the writing of this article or the decision to submit it for publication. All authors declare no other competing interests.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/ffunb.2023.1295223/full#supplementary-material

Supplementary Figure 1 | Fully colonized fruiting block of P. cubensis pre-harvest.

Supplementary Figure 2 | Heat map of fruiting body (n = 18), mycelial (n = 7), and grain mycelial samples (n = 6) and known targeted compounds. The Euclidean distance function was used to generate dendrograms, and columns and rows were clustered using an unweighted pair group with arithmetic mean (UPGMA) in R.

Supplementary Figure 3 | Heat Map of Fruiting Body (n=18), Mycelial (n=7) and Grain Mycelial samples (n=6) and untargeted compounds. The Euclidean distance function was used to generate dendrograms, and columns and rows were clustered using an unweighted pair group with arithmetic mean (UPGMA) in R.

Abbreviations

5-HT_2A, serotonin receptor subtype 2A; α-GPC, L-alpha-glycerylphosphorylcholine; ADAP, automated data analysis pipeline; CBD, cannabidiol; FWHM, full width at half maximum; MS2, tandem mass spectrometry; PDA, potato dextrose agar; PsiD, psilocybin decarboxylase; PsiH, psilocybin hydroxylase; PsiM, psilocybin methyl transferase; PsiK, psilocybin kinase; m/z, mass-to-charge ratio; H-ESI, heated electron spray ionization; PC, principal component; PCA, principal component analysis; RT, retention time; SRM, selected reaction monitoring; THC Δ⁹-, tetrahydrocannabinol; UHPLC-UV/VIS, ultra-high-pressure liquid chromatography coupled with ultraviolet and visible light detector; UHPLC-HRMS, ultra-high-pressure liquid chromatography coupled with high-resolution mass spectrometry; YMEA, yeast malt extract agar.

References

Albizu L., Holloway T., Gonzalez-Maeso J., Sealfon S. C. (2011). Functional crosstalk and heteromerization of serotonin 5-HT_2A and dopamine D2 receptors. Neuropharmacology 61 (4), 770–777. doi: 10.1016/j.neuropharm.2011.05.023

PubMed Abstract | CrossRef Full Text | Google Scholar

Almaula N., Ebersole B. J., Zhang D., Weinstein H., Sealfon S. C. (1996). Mapping the binding site pocket of the serotonin 5-hydroxytryptamine2A receptor: Ser3. 36 (159) provides a second interaction site for the protonated amine of serotonin but not of lysergic acid diethylamide or bufotenin. J. Biol. Chem. 271 (25), 14672–14675. doi: 10.1074/jbc.271.25.14672

PubMed Abstract | CrossRef Full Text | Google Scholar

Bellar D., LeBlanc N. R., Campbell B. (2015). The effect of 6 days of alpha glycerylphosphorylcholine on isometric strength. J. Int. Soc. Sports Nutr. 12 (1), 42. doi: 10.1186/s12970-015-0103-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Bradshaw Alexander J., Backman T. A., Ramirez-Cruz V., Forrister D. L., Winter J. M., Guzman-Davalos L., et al. (2022). DNA authentication and chemical analysis of psilocybe mushrooms reveal widespread misdeterminations in fungaria and inconsistencies in metabolites. Appl. Environ. Microbiol. 88 (24), e01498–e01422. doi: 10.1128/aem.01498-22

PubMed Abstract | CrossRef Full Text | Google Scholar

Brown G. D., Gordon S. (2003). Fungal β-glucans and mammalian immunity. Immunity 19 (3), 311–315. doi: 10.1016/S1074-7613(03)00233-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Brütting C., Hildebrand P., Brandsch C., Stangl G. I. (2021). Ability of dietary factors to affect homocysteine levels in mice: a review. Nutr. Metab. 18 (1), 1–9. doi: 10.1186/s12986-021-00594-9

CrossRef Full Text | Google Scholar

Chen J.-K., Shen C.-R., Liu C.-L. (2010). N-acetylglucosamine: production and applications. Mar. Drugs 8, 2493–2516. doi: 10.3390/md8092493

PubMed Abstract | CrossRef Full Text | Google Scholar

Dinis-Oliveira R. J. (2017). Metabolism of psilocybin and psilocin: clinical and forensic toxicological relevance. Drug Metab. Rev. 49 (1), 84–91. doi: 10.1080/03602532.2016.1278228

PubMed Abstract | CrossRef Full Text | Google Scholar

Evidente A., Kornienko A., Cimmino A., Andolfi A., Lefranc F., Mathieu V., et al. (2014). Fungal metabolites with anticancer activity. Natural Product Rep. 31 (5), 617–627. doi: 10.1039/C3NP70078J

CrossRef Full Text | Google Scholar

Fasciani I., Petragnano F., Aloisi G., Marampon F., Carli M., Scarselli M., et al. (2020). Allosteric modulators of G protein-coupled dopamine and serotonin receptors: a new class of atypical antipsychotics. Pharmaceuticals 13 (11), 388. doi: 10.3390/ph13110388

PubMed Abstract | CrossRef Full Text | Google Scholar

Fricke J., Blei F., Hoffmeister D. (2017). Enzymatic synthesis of psilocybin. Angewandte Chemie Int. Edition 56 (40), 12352–12355. doi: 10.1002/anie.201705489

CrossRef Full Text | Google Scholar

Gartz J., Moller G. K. (1989). Analysis and cultivation of fruit bodies and mycelia of Psilocybe bohemica. Biochemie und Physiologie der Pflanzen 184 (3-4), 337–341. doi: 10.1016/S0015-3796(89)80023-X

CrossRef Full Text | Google Scholar

Glatfelter G. C., Pottie E., Partilla J. S., Sherwood A. M., Kaylo K., Pham D. N. K., et al. (2022). Structure–activity relationships for psilocybin, baeocystin, aeruginascin, and related analogues to produce pharmacological effects in mice. ACS Pharmacol. Trans. Sci. 5 (11), 1181–1196. doi: 10.1021/acsptsci.2c00177

CrossRef Full Text | Google Scholar

Gotvaldová K., Hájková K., Borovička J., Jurok R., Cihlářová P., Kuchař M. (2021). Stability of psilocybin and its four analogs in the biomass of the psychotropic mushroom Psilocybe cubensis. Drug testing Anal. 13 (2), 439–446. doi: 10.1002/dta.2950

CrossRef Full Text | Google Scholar

Guzmán G. (2008). Hallucinogenic mushrooms in Mexico: an overview. Economic Bot. 62 (3), 404–412. doi: 10.1007/s12231-008-9033-8

CrossRef Full Text | Google Scholar

Hasler F., Grimberg U., Benz M. A., Huber T., Vollenweider F. X. (2004). Acute psychological and physiological effects of psilocybin in healthy humans: a double-blind, placebo-controlled dose-effect study. Psychopharmacol. (Berl) 172 (2), 145–156. doi: 10.1007/s00213-003-1640-6

CrossRef Full Text | Google Scholar

He Y., Zhang W., Peng F., Lu F., Zhou H., Bao G., et al. (2019). Metabolomic variation in wild and cultured cordyceps and mycelia of Isaria cicadae. Biomed. Chromatogr. 33 (4), e4478. doi: 10.1002/bmc.4478

PubMed Abstract | CrossRef Full Text | Google Scholar

Ikekawa T., Maruyama H., Miyano T., Okura A., Sawasaki Y., Naito K., et al. (1985). PROFLAMIN, A NEW ANTITUMOR AGENT PREPARATION, PHYSICOCHEMICAL PROPERTIES AND ANTITUMOR ACTIVITY. Japanese J. Cancer Res. GANN 76 (2), 142–148.

Google Scholar

Inserra A., De Gregorio D., Gobbi G. (2021). Psychedelics in psychiatry: neuroplastic, immunomodulatory, and neurotransmitter mechanisms. Pharmacol. Rev. 73 (1), 202–277. doi: 10.1124/pharmrev.120.000056

PubMed Abstract | CrossRef Full Text | Google Scholar

Jensen N., Gartz J., Laatsch H. (2006). Aeruginascin, a trimethylammonium analogue of psilocybin from the hallucinogenic mushroom Inocybe aeruginascens. Planta Med. 72 (07), 665–666. doi: 10.1055/s-2006-931576

PubMed Abstract | CrossRef Full Text | Google Scholar

John P. A., Wong K. H., Naidu M., Sabaratnam V., David P. (2013). Combination effects of curcumin and aqueous extract of lignosus rhinocerotis mycelium on neurite outgrowth stimulation activity in PC-12 cells. Natural Product Commun. 8 (6), 1934578X1300800608. doi: 10.1177/1934578X1300800608

CrossRef Full Text | Google Scholar

Jung J. Y., Lee I. K., Seok S. J., Lee H. J., Kim Y. H., Yun B. S. (2008). Antioxidant polyphenols from the mycelial culture of the medicinal fungi Inonotus xeranticus and Phellinus linteus. J. Appl. Microbiol. 104 (6), 1824–1832. doi: 10.1111/j.1365-2672.2008.03737.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Karaki S., Becaml C., Muratt S., Mannoury la Cour C., Millan M. J., Prézeau L., et al. (2014). Quantitative phosphoproteomics unravels biased phosphorylation of serotonin 2A receptor at Ser280 by hallucinogenic versus nonhallucinogenic agonists. Mol. Cell. Proteomics 13 (5), 1273–1285. doi: 10.1074/mcp.M113.036558

PubMed Abstract | CrossRef Full Text | Google Scholar

Kawagishi H., Shimada A., Hosokawa S., Mori H., Sakamoto H., Ishiguro Y. (1996). Erinacines E, F, and G, stimulators of nerve growth factor (NGF)-synthesis, from the mycelia of Hericium erinaceum. Tetrahedron Lett. 37 (41), 7399–7402. doi: 10.1016/0040-4039(96)01687-5

CrossRef Full Text | Google Scholar

Kawagishi H., Masui A., Tokuyama S., Nakamura T. (2006). Erinacines J and K from the mycelia of Hericium erinaceum. Tetrahedron 62 (36), 8463–8466. doi: 10.1016/j.tet.2006.06.091

CrossRef Full Text | Google Scholar

Koltai H., Namdar D. (2020). Cannabis phytomolecule 'Entourage': from domestication to medical use. Trends Plant Sci. 25 (10), 976–984. doi: 10.1016/j.tplants.2020.04.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Laprairie R. B., Bagher A. M., Kelly M. E. M., Denovan-Wright E. M. (2015). Cannabidiol is a negative allosteric modulator of the cannabinoid CB1 receptor. Br. J. Pharmacol. 172 (20), 4790–4805. doi: 10.1111/bph.13250

PubMed Abstract | CrossRef Full Text | Google Scholar

Lee S. H., Choi Y. B., Kim J. H., Kho A. R., Sohn M., Song H. K. (2017). Late treatment with choline alfoscerate (l-alpha glycerylphosphorylcholine, α-GPC) increases hippocampal neurogenesis and provides protection against seizure-induced neuronal death and cognitive impairment. Brain Res. 1654, 66–76. doi: 10.1016/j.brainres.2016.10.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Lever M., Slow S. (2010). The clinical significance of betaine, an osmolyte with a key role in methyl group metabolism. Clin. Biochem. 43 (9), 732–744. doi: 10.1016/j.clinbiochem.2010.03.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Li I. C., Lee Y. L, Tzeng T. T., Chen W. P, Chen Y. P., Shiao Y. J. (2018). Neurohealth properties of hericium erinaceus mycelia enriched with erinacines. Behav. Neurol. 2018, 5802634. doi: 10.1155/2018/5802634

PubMed Abstract | CrossRef Full Text | Google Scholar

López-Giménez J. F., González-Maeso J. (2017). Hallucinogens and serotonin 5-HT_2A receptor-mediated signaling pathways. Behav. Neurobiol. Psychedelic Drugs p, 45–73. doi: 10.1007/7854_2017_478

CrossRef Full Text | Google Scholar

López-Giménez J. F., González-Maeso J. (2018). Hallucinogens receptor-mediated and serotonin signaling 5-HT pathways. Behav. Neurobiol. Psychedelic Drugs 36, 45. doi: 10.1007/7854_2017_478

CrossRef Full Text | Google Scholar

Ma B.-J., Shen W. J., Yu H. Y., Ruan Y., Wu T. T., Zhao X. (2010). Hericenones and erinacines: stimulators of nerve growth factor (NGF) biosynthesis in Hericium erinaceus. Mycology 1 (2), 92–98. doi: 10.1080/21501201003735556

CrossRef Full Text | Google Scholar

Mahmood Z. A. (2013). Bioactive alkaloids from fungi: psilocybin. Pakistan J. Pharmacol. 23 (3), 349–357.

Google Scholar

Marcus L., Soileau J., Judge L. W., Bellar D. (2017). Evaluation of the effects of two doses of alpha glycerylphosphorylcholine on physical and psychomotor performance. J. Int. Soc. Sports Nutr. 14 (1), 39. doi: 10.1186/s12970-017-0196-5

PubMed Abstract | CrossRef Full Text | Google Scholar

McKernan K., Kane L. T., Crawford S., Chin C. S., Trippe A., McLaughlin S. (2021). A draft reference assembly of the Psilocybe cubensis genome. F1000Research 10, 281. doi: 10.12688/f1000research.51613.2

PubMed Abstract | CrossRef Full Text | Google Scholar

Mygind P. H., Fischer R. L., Schnorr K. M., Hansen M. T., Sönksen C. P., Ludvigen S. (2005). Plectasin is a peptide antibiotic with therapeutic potential from a saprophytic fungus. Nature 437 (7061), 975–980. doi: 10.1038/nature04051

PubMed Abstract | CrossRef Full Text | Google Scholar

Okuyama S., Sawasaki E., Yokogoshi H. (2004). Conductor compounds of phenylpentane in mycoleptodonoides aitchisonii mycelium enhance the release of dopamine from rat brain striatum slices. Nutr. Neurosci. 7 (2), 107–111. doi: 10.1080/10284150410001710429

PubMed Abstract | CrossRef Full Text | Google Scholar

Park I.-H., Chung S. K., Lee K. B., Yoo Y. C., Kim S. K., Kim G. S, et al. (2004). An antioxidant hispidin from the mycelial cultures ofPhellinus linteus. Arch. Pharmacal Res. 27 (6), 615–618. doi: 10.1007/BF02980159

CrossRef Full Text | Google Scholar

Pérez-Vendrell A. M., Gausch J., Francesch M., Molina-Cano J. L., Brufau J. (1995). Determination of β-(1–3),(1–4)-D-glucans in barley by reversed-phase high-performance liquid chromatography. J. Chromatogr. A 718 (2), 291–297. doi: 10.1016/0021-9673(95)00694-X

PubMed Abstract | CrossRef Full Text | Google Scholar

Phan C.-W., Lee G. S., Hong S. L., Wong Y. T., Brkljača R., Urban S., et al. (2014). Hericium erinaceus (Bull.: Fr) Pers. cultivated under tropical conditions: isolation of hericenones and demonstration of NGF-mediated neurite outgrowth in PC12 cells via MEK/ERK and PI3K-Akt signaling pathways. Food Funct. 5 (12), 3160–3169. doi: 10.1039/c4fo00452c

PubMed Abstract | CrossRef Full Text | Google Scholar

Phan C.-W., David P., Naidu M., Wong K. H., Sabaratnam V. (2015). Therapeutic potential of culinary-medicinal mushrooms for the management of neurodegenerative diseases: diversity, metabolite, and mechanism. Crit. Rev. Biotechnol. 35 (3), 355–368. doi: 10.3109/07388551.2014.887649

PubMed Abstract | CrossRef Full Text | Google Scholar

Pottie E., Dedecker P., Stove C. P. (2020). Identification of psychedelic new psychoactive substances (NPS) showing biased agonism at the 5-HT_2AR through simultaneous use of β-arrestin 2 and miniGαq bioassays. Biochem. Pharmacol. 182, 114251. doi: 10.1016/j.bcp.2020.114251

PubMed Abstract | CrossRef Full Text | Google Scholar

Rickli A., Moning O. D., Hoener M. C., Liechti M. E. (2016). Receptor interaction profiles of novel psychoactive tryptamines compared with classic hallucinogens. Eur. Neuropsychopharmacol. 26 (8), 1327–1337. doi: 10.1016/j.euroneuro.2016.05.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Sigala S., Imperato A., Rizzoelli P., Casolini P., Missale C., Spano P. (1992). L-α-glycerylphorylcholine antagonizes scopolamine-induced amnesia and enhances hippocampal cholinergic transmission in the rat. Eur. J. Pharmacol. 211 (3), 351–358. doi: 10.1016/0014-2999(92)90392-H

PubMed Abstract | CrossRef Full Text | Google Scholar

Simmen U., Higelin J., Berger-Büter K., Schnaffnew W., Lundstrom K. (2001). Neurochemical studies with St. John's wort in vitro. Pharmacopsychiatry 34 (Sup. 1), 137–142. doi: 10.1055/s-2001-15475

PubMed Abstract | CrossRef Full Text | Google Scholar

Singh R. S., Bhari R., Rana V., Tiwary A. K. (2011). Immunomodulatory and therapeutic potential of a mycelial lectin from aspergillus nidulans. Appl. Biochem. Biotechnol. 165 (2), 624–638. doi: 10.1007/s12010-011-9281-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Torrens-Spence M. P., Liu C. T., Pluskal T., Chung Y. K., Weng J. K. (2018). Monoamine biosynthesis via a noncanonical calcium-activatable aromatic amino acid decarboxylase in psilocybin mushroom. ACS Chem. Biol. 13 (12), 3343–3353. doi: 10.1021/acschembio.8b00821

PubMed Abstract | CrossRef Full Text | Google Scholar

Toth O. J., Luu B., Ourisson G. (1983). Les acides ganoderiques tàz : triterpenes cytotoxiques de Ganoderma lucidum (Polyporacée). Tetrahedron Lett. 24 (10), 1081–1084. doi: 10.1016/S0040-4039(00)81610-X

CrossRef Full Text | Google Scholar

Wieczorek P. P., Witkowska D., Jasicka-Misiak I., Poliwoda A., Oterman M., Zielińska K. (2015). “Chapter 5 - bioactive alkaloids of hallucinogenic mushrooms,” in Studies in Natural Products Chemistry. Ed. Atta ur R. (Poland: Elsevier), 133–168.

Google Scholar

Xie C., Yan S., Zhang Z., Gong W., Zhu Z., Zhou Y. (2020). Mapping the metabolic signatures of fermentation broth, mycelium, fruiting body and spores powder from Ganoderma lucidum by untargeted metabolomics. LWT 129, 109494. doi: 10.1016/j.lwt.2020.109494

CrossRef Full Text | Google Scholar

Yang Y., Yang J., Wang H., Jin Y., Liu J., Jia R., et al. (2021). Analysis of primary metabolites of Morchella fruit bodies and mycelium based on widely targeted metabolomics. Arch. Microbiol. 204 (1), 98. doi: 10.21203/rs.3.rs-1015665/v1

PubMed Abstract | CrossRef Full Text | Google Scholar