Ganoderma lucidum promotes sleep through a gut microbiota-dependent and serotonin-involved pathway in mice

Ganoderma lucidum is a medicinal mushroom used in traditional Chinese medicine with putative tranquilizing effects. However, the component of G. lucidum that promotes sleep has not been clearly identified. Here, the effect and mechanism of the acidic part of the alcohol extract of G. lucidum mycelia (GLAA) on sleep were studied in mice. Administration of 25, 50 and 100 mg/kg GLAA for 28 days promoted sleep in pentobarbital-treated mice by shortening sleep latency and prolonging sleeping time. GLAA administration increased the levels of the sleep-promoting neurotransmitter 5-hydroxytryptamine and the Tph2, Iptr3 and Gng13 transcripts in the sleep-regulating serotonergic synapse pathway in the hypothalamus during this process. Moreover, GLAA administration reduced lipopolysaccharide and raised peptidoglycan levels in serum. GLAA-enriched gut bacteria and metabolites, including Bifidobacterium, Bifidobacterium animalis, indole-3-carboxylic acid and acetylphosphate were negatively correlated with sleep latency and positively correlated with sleeping time and the hypothalamus 5-hydroxytryptamine concentration. Both the GLAA sleep promotion effect and the altered faecal metabolites correlated with sleep behaviours disappeared after gut microbiota depletion with antibiotics. Our results showed that GLAA promotes sleep through a gut microbiota-dependent and serotonin-associated pathway in mice.

GLAA potentiates pentobarbital-induced sleep. Using a mouse model of pentobarbital-induced hypnosis, we observed that gavage administration of GLAA for 4 weeks led to significant decreases in sleep latency and increases in sleeping time compared with the solvent control containing 0.05% sodium carboxymethyl acid (0.05% CMC-Na). The average sleep latency of the control group was approximately 300 s, while those of groups administered GLAA at levels ranging from 25 to 100 mg/kg were shortened to approximately 250 s (Fig. 1c). The sleeping time of mice gavaged with 25 mg/kg, 50 mg/kg and 100 mg/kg GLAA increased significantly from 1570.4 s (at 0 mg/kg) to 2816.5 s, 3000 s and 3523.5 s, respectively, showing a certain dose dependency (Fig. 1d). These results suggested that gavage administration of GLAA could significantly promote the sleep quality of mice.
GLAA increases the 5-HT content in the hypothalamus. To determine which neurotransmitters were responsible for the sleep changes in mice, 5-HT, GABA, norepinephrine and dopamine were detected in www.nature.com/scientificreports/ the hypothalamus of mice with a commercial ELISA kit. Only the 5-HT level was significantly increased after administration of 50 mg/kg and 100 mg/kg GLAA (Fig. 1e), indicating that 5-HT may act as a target neurotransmitter for GLAA to promote sleep health in mice.

GLAA reduces the LPS content and increases the PG content in serum.
LPS and PG are two major components of the bacterial cell wall and two indicators of body inflammation 14 . LPS exists in the cell wall of only gram-negative bacteria, and the PG content in the cell wall of gram-positive bacteria is much higher than that in the cell wall of gram-negative bacteria 17,18 . In our study, administration of GLAA significantly reduced the LPS content but raised the PG content in serum (Fig. 1f), indicating that GLAA may change the composition of the gut microbiota.
GLAA changes the composition of the gut microbiota. Next-generation sequencing of 16S rRNA V3-V4 regions was conducted to investigate the faecal microbiota composition. As shown by the Chao1 index, the microbiota of the control and 100 mg/kg GLAA groups exhibited the same richness (Fig. 2a); as shown by the Simpson index, there was also no significant difference in microbiota diversity between the control group and the 100 mg/kg GLAA group (Fig. 2a). However, the principal coordinate analysis (PCoA) results showed that the microbiota profiles of the two groups were clearly separated from each other (Fig. 2b). PERMANOVA further confirmed that the compositions were different between the control and 100 mg/kg GLAA groups (p = 0.001).
To analyse the effect of 100 mg/kg GLAA on the intestinal microbiota, discriminative taxa were identified by LEfSe analysis at multiple levels, and the most different results (linear discriminant analysis (LDA) score > 3.5) were retained to investigate the alteration in the microbiota structure (Fig. 2c). Compared with the control group, the 100 mg/kg group exhibited enrichment of the family Bifidobacteriaceae and its affiliated genus Bifidobacterium and species Bifidobacterium animalis, the family Lactobacillaceae and its affiliated genus Lactobacillus and species Lactobacillus reuteri, the family Porphyromonadaceae and its affiliated genus Odoribacter, the family Erysipelotrichales and its affiliated genus Turicibacter, and the genus Rikenella. In addition, administration of 100 mg/kg GLAA depleted the family Ruminococcaceae and its affiliated genus Oscillibacter and the genera Anaerotignum and Roseburia. Finally, 100 mg/kg GLAA significantly enriched Bifidobacterium, Rikenella, Odoribacter and Turicibacter at the genus level, as well as Bifidobacterium animalis at the species level after p value correction (p adj < 0.05, Fig. 2d).
The correlation of discriminative bacterial genera with the sleep behavioural index and concentration of 5-HT in the hypothalamus were determined using Spearman's rank correlation analysis (Fig. 2e). Bifidobacterium, Odoribacter, Rikenella and Turicibacter at the genus level, as well as Bifidobacterium animalis at the species level, were negatively correlated with sleep latency. Bifidobacterium and Odoribacter at the genus level, as well as Bifidobacterium animalis at the species level, were positively correlated with sleeping time and the concentration of 5-HT in the hypothalamus. These results indicated that the gut microbiota may be involved in the process by which GLAA promotes sleep in mice.
GLAA alters gut metabolites. To explore the relationship between alterations in gut metabolites and sleep, faecal samples were submitted to metabolome analysis using LC-MS/MS. In total, 13,023 metabolites were identified. 1037 metabolites with score > 0.5 and a level 2 identification were retained (supplemental material 1) 19 . Principal component analysis (PCA) showed that the control group and 100 mg/kg group clustered separately from each other in both the original identification metabolome profiles and the retained metabolome profiles, suggesting that the metabolome profiles of these two groups (Fig. 3a, Supplementary Fig. S2a) are quite different. Similar results were obtained from the orthogonal partial least squares discriminant analysis (OPLS-DA) (Fig. 3b, Supplementary Fig. S2b).
Spearman's rank correlation analysis showed that most metabolites were correlated with the sleep behavioural indexes and concentration of 5-HT in the hypothalamus (Fig. 3c). d-ribose and daidzein showed negative correlation with sleep-promotion index and hypothalamic 5-HT concentration. On the contrary, galactosylsphingosine, 2-hydroxy-7-hydroxymethylchromene-2-carboxylate, 2-hydroxy-2-ethylsuccinic acid and indole-3-carboxylic acid were positively correlated with sleep-promotion index and hypothalamic 5-HT concentration. Hydroxypyruvic acid, malonic acid, avermectin B1b aglycone, steviol, d-tagatose, γ-l-glutamyl-l-2-aminobutyrate, His-Phe-Met, 5-androstenetriol, d-(+)-raffinose, Leu-Asp-Ser, 1-Nitro-7-glutathionyl-8-hydroxy-7,8-dihydronaphthalene and Leu-Arg-Ser-Lys were correlated with at least one sleep behavioural index, which suggested that they were positively correlated with sleep-promotion effect. Conversely, Cys-Arg-Phe, acetylphosphate, Leu-Arg-Ser-Lys and betaine exhibited negatively correlated with sleep-promotion effect. www.nature.com/scientificreports/ Antibiotic-induced depletion of the gut microbiota inhibits the sleep-promoting effect and the increase in 5-HT caused by GLAA. To explore whether the sleep promotion effect of GLAA was dependent on the gut microbiota, we depleted the gut microbiota in mice using antibiotics according to a previous report 20 . After 12 days of depletion of the gut microbiota using an antibiotic concoction, GLAA was gavaged from the 13th day to the 40th day, together with the antibiotic concoction. Notably, neither antibiotic treatment nor GLAA gavage had a significant effect on the body weight of mice (Fig. 4b).
Then, a mouse model of the pentobarbital hypnosis test was established. All the control mice fell asleep within 600 s, while all mice administered 100 mg/kg GLAA fell asleep within 300 s. However, mice pre-treated with antibiotic concoction exhibited only 40% sleep incidence during the experimental period. Even after administration of 100 mg/kg GLAA, the incidence of sleep in mice pre-treated with antibiotic concoction was only 50% (Fig. 4a). Furthermore, pre-treatment with antibiotic concoction inhibited the increase in 5-HT induced by 100 mg/kg GLAA, as observed in mice (Fig. 4a). Therefore, depletion of the gut microbiota prevented the sleep promotion effect of GLAA, indicating that gut microbiota might play a pivotal role in GLAA-mediated promotion of sleep.
Our results showed that antibiotic treatment depleted more than 85% of the operational taxonomic units (OTUs) of gut bacteria, with a significant decrease in the Chao1 and Simpson indexes (Fig. 4c); moreover, the significant decrease in the LPS content in serum further indicated a decrease in microbial quantity in the gut (Fig. 4c). The remaining bacteria mainly included Proteobacteria (99.35%) at the phylum level ( Fig. 4d), Enterobacteriaceae (99.06%) at the family level ( Supplementary Fig. S3a), and Klebsiella (81.4%) and Proteus (13.8%) at the genus level ( Supplementary Fig. S3b). In contrast, in the control and 100 mg/kg groups without antibiotics, the microbiota presented a normal distribution of Firmicutes (63.6% and 56.96%), Bacteroidetes (23.5% and 32.38%), Proteobacteria (9.6% and 6.47%), Deferribacteres (2.9% and 2.3%), and Actinobacteria (0.1% and 1.64%), with other phyla accounting for less than 0.2% (Fig. 4d). The above-mentioned microbial markers altered by GLAA were all depleted after antibiotic treatment, and this depletion was not reversed by GLAA in the antibiotic-gavaged mice (Fig. 4e).
Among the top 30 faecal metabolites altered by 100 mg/kg GLAA, the alterations in 13 metabolites were totally or partially prevented by depletion of the gut microbiota (Fig. 4f). Depletion of the gut microbiota prevented the GLAA-induced elevation of acetylphosphate, 2-hydroxy-7-hydroxymethylchromene-2-carboxylate, Pro-Val-Val, Leu-Arg-Ser-Lys and indole-3-carboxylic acid. It also prevented the GLAA-induced depletion of hydroxypyruvic acid, Asp-Ile-Glu and daidzein. Notably, except for esprocarb, the remaining 12 of the abovementioned 13 metabolites were strongly correlated with sleep latency, sleeping time and/or the concentration of 5-HT in the hypothalamus.
GLAA regulates the serotonergic synapse pathway. Genome-wide transcriptional profiling was performed to analyse the changes in gene transcription in the hypothalamus of mice. In total, 830,782,772 valid sequencing reads were generated from 21 samples (5-6 samples per group), and 55,450 expressed genes were identified after mapping to the mouse genome. The PCA results showed that the transcriptional profiles were clearly separated from each other (Fig. 5a).
All the significant genes mapped by KEGG Mapper revealed that administration of 100 mg/kg GLAA regulated 10 pathways, namely, nitrogen metabolism; d-arginine and d-ornithine metabolism; protein digestion and absorption; prion diseases; NOD-like receptor signalling pathway; glycine, serine and threonine metabolism; Principal coordinate analysis (PCoA) plot with unweighted UniFrac distances based on operational taxonomic units (OTUs) of the gut microbiota. The image is generated by R package "ade4". (c) The relative abundance of gut bacterial taxa differed between the control group and the 100 mg/kg group and was calculated by LEfSe analysis (LDA score > 3.5). The image is generated on line at the Huttenhower Lab (http:// hutte nhower. sph. harva rd. edu/ galaxy). The taxa enriched in the control group are shown in red, and the taxa in the 100 mg/kg group are shown in green. (d) The discriminative bacteria altered by 100 mg/kg GLAA (LDA score > 3.5 and p adj < 0.05). (e) Correlations of the discriminative bacteria, sleep behaviour and the concentration of 5-HT in the hypothalamus were determined using Spearman's rank correlation analysis. The image is generated by R package "corrplot". The results with p values < 0.05 are shown in the heat map. The colour key and circle size indicate the strength of the correlation (r value). Red indicates positive correlations; blue indicates negative correlations. All the data were calculated from 10 samples per group. The data are shown as the means ± SEMs; # p adj < 0.05, ## p adj < 0.01. Except specifically mentioned, the images are generated by GraphPad Prism 8. G + Gram-positive, GLAA the acidic part of the alcohol extract of G. lucidum mycelia, 5-HT 5-hydroxytryptamine. www.nature.com/scientificreports/ cytokine-cytokine receptor interaction; alcoholism; glutamatergic synapse; and serotonergic synapse (Fig. 5d). These pathways were mostly involved in sleep (glutamatergic synapse and serotonergic synapse) and immunity (prion diseases, NOD-like receptor signalling pathway and cytokine-cytokine receptor interaction) regulation 23,24 . Then, we specifically examined the serotonergic synapse, the output of which is 5-HT. After gavage administration of 100 mg/kg GLAA, the transcription levels of Tph2, Itpr2 and Gng13 in the serotonergic synapse pathway increased significantly, and depletion of the gut microbiota also blocked these increases. Tryptophan is converted by the enzyme tryptophan hydroxylase (TPH) to 5-hydroxytryptophan (5-HTP), and 5-HTP is further converted to 5-HT in presynaptic cells 25 . 5-HT is secreted into postsynaptic cells and binds to 5-HT receptors (mainly 5-HT1R, 5-HT2R and 5-HT5R) based on the KEGG pathway, thus activating downstream signalling pathway-related genes (including Itpr2 and Gng13) and promoting sleep (Fig. 5e,f). www.nature.com/scientificreports/

Discussion
Sleep problems have become a hot topic of global concern. G. lucidum has been used as a tranquilizing agent for the treatment of restlessness, insomnia, and palpitation in China for hundreds of years 26 . Chu et al. reported that 80 mg/kg and 120 mg/kg water extracts of G. lucidum prolonged the sleeping time of mice by 50% and 60%, respectively, while a 40 mg/kg dose did not significantly increased the sleeping time 13 . In this study, for the first time, the acidic part of the ethanol extract of G. lucidum was proven to promote sleeping time by 79.35% (25 mg/ kg), 91.03% (50 mg/kg) and 124.37% (100 mg/kg). The sleep promotion effect of GLAA in this study seemed to be better than that of the water extract of G. lucidum according to published data, although this conclusion was limited by comparisons between different studies. The ethanol extract of G. lucidum was reported to contain approximately 130-140 triterpene acids, 20 sterols, fewer alkaloids and other compounds 27,28 . The preparation of GLAA used in our study is also the usual method to enrich triterpene acids that are one of the main medicinal ingredients of the ethanol extract of G. lucidum, playing important roles in the anti-inflammatory, antioxidant and antitumour activities of this mushroom [29][30][31] . The compositions of GLAA are also very complicated. It was reported that the main component of this part was triterpene acids, although the normal content of triterpene acids in different sources of G. lucidum ranges from 6 to 16% 29 . However, it is not known whether the sleep promoting effect of GLAA is mediated by triterpene acids.
Various neurotransmitters, including GABA, 5-HT, norepinephrine and dopamine, affect different brain nuclei to regulate the switch between wakefulness and sleep 32 . In this study, GLAA was found to promote sleep through a 5-HT associated pathway. The monoamine neurotransmitter 5-HT, also known as serotonin, is important for sleep regulation, and animals deficient in 5-HT synthesis have been shown to stay awake longer and sleep less 32 . Moreover, several researchers have found that increasing the 5-HT content in the brain would prolongs sleeping time in animals, which is consistent with our results. In addition to neurotransmitters, other factors could also affect sleep, such as changes in brain structure, ion channels 33 and growth hormone 21,22 , some of which were also found in our study.
Pentobarbital is a well-known sedative-hypnotic and anaesthetic whose hypnotic effect can be enhanced by calcium channel blockers, such as diltiazem and some plant extracts 34 . Ca 2+ influx and Ca 2+ channels are widely involved in sleep/wake regulation 35 . Consequently, sleep latency and sleeping time in the pentobarbital-induced sleep test are commonly used as indicators for assessing the sedative and hypnotic effects of dietary supplements (most of them are plant extracts) and drugs such as diazepam (a longer-acting benzodiazepine) and 5-HTP 32 . Several studies have reported that serotonergic synapse played negative roles in calcium currents/channels 35 . Christopher et al. found that serotonin inhibited calcium transients in mesopontine cholinergic neurons of the laterodorsal tegmental nucleus 36 . Douglas et al. showed that serotonin inhibited N-and P/Q-type calcium channels in caudal raphe neurons of neonatal rat brain slices 37 . Therefore, the synergistic effect of GLAA and pentobarbital in the present study may be related to 5-HT-mediated calcium blockade.
Many studies have shown that sleep under normal physiological conditions may be influenced by the gut microbiota 38 . PG and LPS are two bacterial cell wall components that are enriched in gram-positive bacteria and gram-negative bacteria, respectively 38 . Our results showed that both PG and LPS levels were altered in serum after gavage administration of GLAA, indicating that GLAA-induced sleep may be related to the gut microbiota. In addition, GLAA may also affect serum LPS/PG levels by regulating gut barrier integrity and preventing translocation into the circulation according to the results of other G. lucidum specimens 39 . Using a mouse model in which the gut microbiota was depleted with antibiotics, we further confirmed that the gut microbiota was necessary for GLAA-mediated promotion of sleep in mice. Interestingly, Bifidobacterium and Odoribacter, which were enriched in the GLAA group, were highly correlated with sleep behaviour and 5-HT in our study. Bifidobacterium longum was reported to upregulate serotonin transporter expression in intestinal epithelial cells in vitro 40 . Bubie et al. found that variation in Odoribacter abundance was the likely regulator of sleep architecture, and Odoribacter abundance was correlated with sleep phenotypes in mice 41 . Therefore, these results provide a promising prospect for the study of the role and mechanism of intestinal microorganisms in the promotion of sleep.
Sleep restriction and irregular sleep patterns have been associated with metabolic impairments 42 . Our results showed that more than half of faecal metabolites were correlated with sleep latency, sleeping time and 5-HT. Some of them are associated with sleep. d-ribose was found to improve in energy, sleep, mental clarity, pain intensity and well-being of Fibromyalgia and chronic fatigue syndrome 43 . Steviol could increase the uptake of myocardial and cerebral glucose that play roles in sleep behaviour disorder 44,45 . Indole-3-carboxylic acid is synthesized from tryptophan that can also convert into 5-hydroxytryptophan, the precursor of 5-HT 46 . Acetyl-phosphate can induce Ca 2+ -Ca 2+ exchange in sarcoplasmic reticulum vesicles 47 . Some metabolites are related with function of brain and neurological diseases. Daidzein can affect various neurobiological regulatory mechanisms such as behaviour, cognition, growth, development and reproduction 48 . Betaine (n-trimethylglycine), a common osmolyte, has function in the modulation of hippocampal neurophysiology and neuroprotection 49 . Interestingly, two herbicides dinitramine and esprocarb were found in these metabolites, although they were not correlated with sleep behaviour or indicators. And we suspected that they might come from the bedding material. Additionally, due to the nontargeted metabolome analysis with LC-MS/MS, integrated database and more included compounds in the database of our study, a large number of metabolites have been explored, which provides us with more ideas about sleep, but these metabolites and results need further verification and research.
Recently, bidirectional interactions between the central nervous system and the gastrointestinal tract have attracted increasing attention. Preclinical evidence supports a role of the gut microbiota in behavioural responses associated with pain, emotion, autism, anxiety, depression, and sleep 50 . The microbiota and the brain communicate with each other via various routes, including the immune system, endocrine and neurocrine pathways, the vagus nerve and the enteric nervous system, involving microbial metabolites such as short-chain fatty acids, www.nature.com/scientificreports/ branched-chain amino acids, and peptidoglycans 51 . According to our results, 5-HT and the serotonergic synapse pathway might be among the routes connecting the gut microbiota and sleep signals in the brain.
In summary, in our study, the effect of the acidic part of the alcohol extract of G. lucidum mycelia on sleep was studied to show that GLAA significantly shortened sleep latency and prolonged sleeping time in pentobarbitaltreated mice. This effect was mainly related to an increase in serotonergic synapse signalling in the hypothalamus and was dependent on the gut microbiota of mice.

Materials and methods
Preparation of GLAA. Ganoderma lucidum CGMCC5.616 was provided by Zhejiang Wuyangtang Pharmaceutical Incorporated Company, Ltd. GLAA was prepared as described with minor modifications 52 . Briefly, 100 g of the G. lucidum fruiting body was crushed, passed through a 20-mesh sieve and mixed with 2 L of 95% ethanol (v/v) before ultrasonic extraction. The solution was ultrasonicated twice for 30 min each time at 40 kHz at room temperature. After filtration, the ethanolic extract was dried by rotary evaporation. The dry product was redissolved in 500 mL of water, degreased with 500 mL of petroleum ether, extracted with 500 mL of ethyl acetate three times, and evaporated to dryness. The ethyl acetate extract was redissolved in 500 mL of ethyl acetate and adjusted to pH 9-10 using 1 mol/L NaOH. Then, the alkalized water layer was isolated and acidified to pH 2-3 with 1 mol/L HCl. The acidified water layer was extracted with ethyl acetate three times to obtain crude GLAA. The GLAA was dried and stored at 4 °C. Experimental research on G. lucidum complied with local and Chinese guidelines and legislation.

Animals. Male specific-pathogen-free (SPF) Institute of Cancer Research (ICR) mice were purchased from
Hangzhou Medical College and kept in SPF animal centre of Hangzhou Medical College under controlled light (12 h light-dark cycle), temperature (22-26 °C) and humidity (40-60%) conditions with free access to food and water. Mice that weighed 20-24 g (approximately 6 weeks) were randomly distributed into groups containing ten animals each (five mice per cage) using the standard = RAND() function in Microsoft Excel. The mice were allowed to acclimate to the laboratory conditions for 7 days before pharmacological tests. All experimental procedures were approved by the Animal Experimentation Ethics Committee of Hangzhou Medical College. The experiments were performed in accordance with the criteria of the "Guide for the Care and Use of Laboratory Animals" (NIH publication 86-23, revised 1985) and the ARRIVE guidelines (https:// arriv eguid elines. org).
Pentobarbital-induced hypnosis test. The pentobarbital hypnosis test was performed as described previously 13 . GLAA was resolved in 0.05% CMC-Na. For the GLAA assay, forty mice were assigned to 4 groups and orally administered GLAA (25, 50 or 100 mg/kg) or 0.05% CMC-Na as a solvent control from 8:00 a.m. to 9:00 a.m. every day for 4 weeks. For the antibiotic depletion experiment, forty mice were assigned to 4 groups: pre-treated saline (control and 100 mg/kg groups) or antibiotics (control + AB and 100 mg/kg + AB groups) for 40 days. On the 13th to the 40th day, the 100 mg/kg and 100 mg/kg + AB groups were gavaged with 100 mg/kg GLAA, while the other two groups (control and control + AB) were administered 0.05% CMC-Na.
On the 28th day of the GLAA assay and the 40th day of the antibiotic depletion experiment, 1 h after the last administration of GLAA, 50 mg/kg pentobarbital sodium was injected subcutaneously. The animals were considered asleep if they stayed motionless when positioned on their back and when they lost their righting reflex. The time interval between injection of pentobarbital and start of sleep was considered sleep latency, while prolongation of sleep was described as sleeping time. Faecal samples were collected on the 28th day after administration of GLAA. Mice were sacrificed under carbon dioxide asphyxia 24 h after the pentobarbital-induced hypnosis test. Blood and tissue (brain, spleen, liver, kidney, thymus and colon) were collected. All samples were stored at − 80 °C until use.
Antibiotic-induced depletion of the gut microbiota. Antibiotics treatment was performed according to Refs. 20,53 and divided into two stages. In the first stage, mice were gavaged with amphotericin B (Solarbio Life Science, Beijing, China) at a dose of 1 mg/kg body weight every 12 h for three days. The second stage was from the 4th day to the 40th day, when the water was supplemented with 1 g/L ampicillin (Sangon Biotech, Shanghai, China). Furthermore, an antibiotic combination consisting of 50 mg/kg vancomycin (Sangon Biotech, Shanghai, China), 100 mg/kg neomycin (Sangon Biotech, Shanghai, China), 100 mg/kg metronidazole (Sangon Biotech, Shanghai, China), and 1 mg/kg amphotericin B was administered by gavage every 12 h in the second stage. antibiotics. (f) Antibiotics partially depleted gut metabolites altered by GLAA. The data are shown as the means ± SEMs; *p < 0.05, **p < 0.01, ***p < 0.001, 'ns' indicates that the difference is not significant. All the data are calculated from 8 to 10 samples per group, and the images are generated by GraphPad Prism 8, except specifically mentioned. GLAA the acidic part of the alcohol extract of G. lucidum mycelia, AB antibiotics, 5-HT 5-hydroxytryptamine, LPS lipopolysaccharide. Raw sequence reads were assigned to each sample according to their unique barcode pairs. Overlapping paired-end reads were merged to form tags using FLASH (version 1.2.11) 55 . The following quality control criteria were used: (1) an exact match to at least one end of barcodes and primers; (2) no undetermined bases in the tags; and (3) no more than three mismatches in the overlap region. Operational taxonomic units (OTUs) were clustered with a 97% similarity cut-off using Usearch 11 56 . Alpha diversity was determined using the R packages "rich" and "diversity". Principal coordinate analysis (PCoA), a standard multivariate statistical technique, was performed with the R package "ade4" to explain differences among microbial communities. The longest sequence in each OTU was chosen as the representative sequence for identification using the Ribosomal Database Project (RDP) database v.11.3 57 . The Chao1 index, Shannon index and principal coordinates analysis (PCoA) of unweighted UniFrac beta diversity were calculated by QIIME software (version 1.8.0) to investigate the richness and diversity of each gut microbiome 58 . Linear discriminant analysis effect size (LEfSe) analysis was performed on the Galaxy web platform 59 . The effect size of each differential taxon was estimated by the linear discriminant analysis (LDA) score, and the taxa with LDA scores greater than 3.5 were defined as discriminative taxa.
Metabolic analysis. LC-MS/MS was used for metabolic analysis of 0.1 g of frozen faeces. A mixture containing 0.1 g of frozen faeces, 300 µL of ice-cold water and 3-4 cold 3 mm steel balls (Jingxin, Shanghai, China) was whirled for 5 min. Then, 500 µL of methanol containing 1 ppm dl-2-chlorophenylalanine was added into the mixture without steel balls, and the mixture was further centrifuged at 12,000 rpm and 4 °C for 10 min. Then, 600 µL of supernatant was sucked into another centrifuge tube and concentrated. Then, 100 µL of 5% methanol water was added to the dried product, mixed and centrifuged at 12,000 rpm at 4 °C for 10 min. The original data file obtained by LC-MS analysis was first converted to mzML format by ProteoWizard software 60 . Peak extraction, alignment and retention time correction were performed by the XCMS program 61 . The "SVR" method was used to correct the peak area. Peaks with deletion rates > 50% were filtered from each group of samples. Then, metabolic identification information was obtained by searching the laboratory's in-house database and integrating the public database (HMDB and KEGG) and then combining it with metDNA 62 . Unbiased principal component analysis (PCA), orthogonal partial least squares discriminant analysis (OPLS-DA) and variable importance in projection (VIP) were calculated using SIMCA software (version 14.1) (Sartorius Stedim Biotech, Umeå, Sweden).

Comparative transcriptomic analysis.
Total RNA from the hypothalamus was extracted using TRIzol reagent (Invitrogen, Carlsbad, CA, USA) following the manufacturer's protocol. The total RNA quantity and purity were analysed with a 2100 Bioanalyzer and the RNA 6000 Nano LabChip Kit (Agilent, Palo Alto, CA, USA) with RIN number > 7.0. Poly(A) RNA was purified from total RNA (5 µg) using poly-T oligo-attached magnetic beads using two rounds of purification. Following purification, the mRNA was fragmented into small pieces using divalent cations under elevated temperature. Then, the cleaved RNA fragments were reverse-transcribed to create the final cDNA library in accordance with the protocol for the mRNA-Seq sample preparation kit (Illumina, San Diego, CA, USA); the average insert size for the paired-end libraries was 300 bp (± 50 bp). Next, we performed paired-end sequencing on an Illumina HiSeq 4000 (Illumina, San Diego, CA, USA) following the vendor's recommended protocol. Then, the sample reads were aligned to the UCSC (http:// genome. ucsc. edu/) mouse reference genome using the HISAT package, which initially removes a portion of the reads based on the quality information accompanying each read and then maps the reads to the reference genome. The mapped reads were assembled using StringTie 63 . Then, all transcriptomes from samples were merged to reconstruct a comprehensive transcriptome using Perl scripts. After the final transcriptome was generated, StringTie and edgeR were used to estimate the expression levels of all transcripts. StringTie was used to determine the expression level of mRNAs by calculating fragments per kilobase of exon model per million mapped reads (FPKM) values 64 . The differentially expressed mRNAs and genes were selected with log 2 (fold change) > 1 or log 2 (fold change) < − 1 and with statistical significance (P value < 0.05) by the R package. www.nature.com/scientificreports/ Statistical analysis. The obtained data were analysed using SPSS 22.0 software (IBM Corporation, Armonk, NY, USA). The data were first tested for normal distribution using the Kolmogorov-Smirnov test. The Mann-Whitney U test (for nonnormal distributions) or Student's t test (for normal distributions) was performed to test for differences between two groups. Spearman's rank correlation analysis was conducted using the corrplot package of R. The images are generated by GraphPad Prism 8 (GraphPad Software, La Jolla, CA, USA), R and SIMCA software.

Data availability
All data generated or analysed during this study are included in this published article and its additional information files. The data sets generated during the current study are available in the GenBank Sequence Read Archive Repository under accession number PRJNA731631. www.nature.com/scientificreports/