Schisandrin alleviates the cognitive impairment in rats with Alzheimer’s disease by altering the gut microbiota composition to modulate the levels of endogenous metabolites in the plasma, brain, and feces

Abstract

Schisandrin is one of the main active compounds isolated from the fruit of Schisandrae chinensis Fructus, which is scientifically proven to have beneficial effects on Alzheimer’s disease (AD) treatment at the cellular and whole organism level. However, the oral availability of schisandrin is very low, thus implying that the underlying mechanism of therapeutic effect on AD treatment is yet to be clarified fully. Therefore, we speculated that the therapeutic effect of schisandrin on AD is mainly by regulating the imbalance of the gut microbiota (GM). In this study, behavioral experiments and H&E staining were used to confirm the pharmacological effects of schisandrin on rats with AD. 16S rDNA gene sequencing and feces, plasma, and brain metabolomics techniques were utilized to investigate the therapeutic effects and the underlying mechanisms of schisandrin on cognitive impairment in rats with AD. The results indicated that schisandrin improved cognitive impairment and hippocampal cell loss in rats. The UPLC-QTOF/MS-based metabolomics studies of the feces, plasma, and brain revealed that 44, 96, and 40 potential biomarkers, respectively, were involved in the treatment mechanism of schisandrin. Schisandrin improved the metabolic imbalance in rats with AD, and the metabolic changes mainly affected the primary bile acid biosynthesis, sphingolipid metabolism, glycerophospholipid metabolism, and unsaturated fatty acid biosynthesis. Schisandrin can improve the GM structure disorder and increase the abundance of beneficial bacteria in the gut of rats with AD. The predictive metagenomics analysis indicated that the altered GM was mainly involved in lipid metabolism, steroid hormone biosynthesis, arachidonic acid metabolism, biosynthesis of unsaturated fatty acids, and bacterial invasion of epithelial cells. Spearman’s correlation analysis showed a significant correlation between affected bacteria and metabolites in various metabolic pathways. Overall, the data underline that schisandrin improves the cognitive impairment in rats with AD by affecting the composition of the GM community, thus suggesting the potential therapeutic effect of schisandrin on the brain–gut axis in rats with AD at the metabolic level.

1 Introduction

Alzheimer’s disease (AD), being a neurodegenerative disorder, brought a heavy mental burden and economic pressure to the people worldwide since its onset (Lee et al., 2018). The most compelling hypothesis about the pathophysiology of AD is the accumulation of extracellular amyloid-β (Aβ) plaques and the formation of neurofibrillary hyperphosphorylated tau (p-tau) tangles, eventually leading toward the loss of neurons and synapses (Weiner et al., 2017).

The gut microbiota (GM) is considered as an important companion to help balance the vital life functions of the host through participation in the processes of maintaining the host’s defense and inflammatory responses, digestion, and nutrient bioavailability (Cryan et al., 2020; Pereira et al., 2020). Previous studies suggest that the microbiota and its induced immunity are involved in the pathogenesis and progression of neurodegenerative diseases (Bell et al., 2019; Kundu et al., 2019). Also, there is a clear correlation between the changes in GM and the development of AD. For example, in patients with AD, the abundance of microbiota has changed (Li et al., 2019), and further studies have revealed that the GM may participate in Aβ aggregation (Marizzoni et al., 2020). Therefore, GM can be used as a potential new target for AD treatment, and the pathogenesis of AD and its possible treatment methods can be discovered by studying the GM, which will be of great significance.

Schisandrin is one of the main active ingredients of the Schisandrae chinensis Fructus, popular for its therapeutic effect on AD at the molecular, cellular, and whole organism level (Wei et al., 2018; Zhao et al., 2019). Previous studies conducted on the efficacy of schisandrin had a main focus on the brain, and further research has not been conducted on the overall biological activity of schisandrin and its specific treatment mechanisms. The oral bioavailability of schisandrin is only 7.69%, and the blood–brain barrier absorption is 0.1, thereby hindering its absorption by the blood through oral administration. Herein, we speculated that schisandrin may show its biological activity in the intestines, primarily by affecting the genetically modified environment rather than absorption through blood. Therefore, the systematic study about the effects of schisandrin on the brain–gut axis in rats with AD can further clarify the underlying therapeutic mechanisms of treatment and would be of great significance in evaluating the effect of GM on metabolism.

Metabolomics can systematically analyze small molecules in cells, tissues, organs, or biological systems which helps out in understanding the mechanism of changes being occurred (Alseekh and Fernie 2018). In this regard, the 16S rDNA gene can accurately quantify all the bacterial species in the GM, as it is the most commonly used molecular marker for the systematic classification of prokaryotic microorganisms (Lu et al., 2020). The expression of the core functions of GM can alter with the change of the gut microenvironment, which would have variable effects on many tissues of the host. Metabolomics analysis helps in comparing the metabolic changes in the GM with the metabolic changes in the host (Evans et al., 2014). Considering the complex mode of action of schisandrin, this study mainly uses metabolomics combined with 16S rDNA gene sequencing technology to systematically explore its protective effect on rats with AD for the first time at the level of brain–gut axis. Furthermore, database prediction and correlation statistics methods were also employed in the research process, which provides research directions for future in-depth research.

2 Materials and methods

2.1 Chemical reagents

Methanol, formic acid, and acetonitrile were obtained from Merck AG (Germany), and the water used throughout the experiment was sourced from Watsons. Among other reagents, D-(+) Galactose was purchased from Beijing Soleibao Technology Co., Ltd. (China), Aβ_25-35 was purchased from Shanghai Yuanye Biological Technology Co., Ltd. (China), and schisandrin (purity >98%) was purchased from the National Institute for Food and Drug Control. The chemical structure for schisandrin is provided as a diagram in Supplementary Figure S1.

2.2 Animals

Regarding the animals utilized in this study, male SD rats (8 weeks old, 200–250 g) were purchased from Changsheng Biotechnology Co., Ltd. (China), from them 18 rats were housed in the cages (one cage for every three rats) which were placed in a 23 ± 2°C and 50 ± 5% humidity environment (food and water were freely available). The animal study was carried out in accordance with the guidelines for Animal Experimentation of China Medical University. The Animal Ethics Committee had approved the institutional protocol (CMU2019257). The Morris water maze (MWM) test was carried out in an air-regulated and soundproofed experimental room, and the rats got used to the experimental room at least 1 h beforehand.

2.3 Animal experiment design

Aβ_25-35 peptide solution was incubated in sterile water at 37°C for 7 days before injection. A total of 18 rats were randomly divided (Excel’s RAND function) into three groups (blank group, model group, and schisandrin administration group (Sch group), n = 6). From day 1 after the completion of adaptive feeding, the rats from the model and Sch groups were administered with D-galactose (120 mg kg⁻¹) through intraperitoneal injection and were given AlCl₃ (80 mg kg⁻¹) by oral gavage, whereas the blank group rats were given an equal volume of normal saline at the same time, and this whole administration was administered once a day until the MWM test was conducted. From the sixth week onward, rats from the Sch group were treated with schisandrin at the rate of 10 mg kg⁻¹, calculated according to the conversion algorithm of human and rat body surface area. The rats from the model and blank groups received an equal volume of normal saline. This administration continued until the ninth week. On the last day of the eighth week, all the rats were kept fasted for 24 h before the surgery. On day 1 of the ninth week, all the rats were anaesthetized with 10% chloral hydrate and were fixed on the brain stereotaxic instruments. According to the stereotactic map of the brain, 3 μL of Aβ_25-35 (1 mg ml⁻¹) was injected into each bilateral hippocampal CA1 area of rats from the model and Sch groups, whereas, the rats in the blank group were injected with an equal volume of saline (the anterior fontanelle is the zero starting point. The rats were drilled 3.6 mm behind the anterior fontanelle and 2.6 mm on the left and right sides of the median raphe by the stereotaxic apparatus. With a depth of 2.9 mm and a slow injection of 3 μL Aβ_25-35, the rat scalp tissues were sutured with sterile surgical sutures). The rats from the Sch group were continuously given schisandrin for 1 week after the cerebroventricular injection test.

2.4 Morris water maze test

After the drug treatment, the rats were trained on the MWM, in which the main body of the water maze device was a circular pool (diameter, 150 cm; height, 60 cm), and it was divided into four quadrants artificially. The escape platform was located at the first quadrant (1 cm below the water surface). From day 1–5, the rats under observation were randomly placed into the water maze with their backs facing toward the pool wall. The recording time was set to 60 s, during which the rats climb to the platform and stays there for 3 s, after which the time recording stopped and started counting as the escape latency (if the platform was not found for more than 60 s, the rats were guided to see the venue and allowed to stay on the platform for 10 s, the escape latency was counted as 60 s). The rats from all the groups were trained daily in four quadrants. The platform was removed from the first quadrant on day 4 to conduct the space exploration experiment. The number of times was recorded for each rat crossing the target platform position, and that number was considered as final statistics which the rat has made crossings in the first quadrant in 60 s. The experimental arrangement plan is provided in Supplementary Figure S2.

2.5 Sample collection and processing

2.5.1 Sample collection

The samples were collected after the MWM tests were over. Fresh feces samples were collected in sterile EP tubes directly from the anus, and were stored at −80°C for metabolomics analysis and 16S rDNA gene sequencing. On the next day, the rats were anesthetized by intraperitoneal injection of 10% chloral hydrate (0.25 mg 100 g⁻¹), and plasma samples were collected from the abdominal aorta into EP tubes coated with sodium heparin. The centrifuged plasma supernatant was stored at −80°C (centrifuged at 1,600 × g for 10 min at 4°C) for metabolomics analysis. Thereafter, the rats were sacrificed by decapitation, and their brains were immediately removed. One half of each brain was soaked in 4% paraformaldehyde solution, and the other half was stored at −80°C for Ultra Performance Liquid Chromatography/Quadrupole time-of-flight mass spectrometer (UPLC-QTOF/MS)-based analysis.

2.5.2 Sample processing

Each 200 µl aliquot of the plasma sample was added with 600 µL of cooled methanol. The mixtures were then vortexed (60 s) and placed at −20°C (20 min) to allow the proteins to precipitate. Thereafter, the mixtures were centrifuged (6,200 × g, 10 min, 4°C), and all the supernatants transferred to EP tubes were allowed to evaporate until dryness with the stream of nitrogen (35°C). Thereafter, 200 µl of 80% methanol solution was added to the dried residue. After the mixture was completely dissolved, it was centrifuged (6,200 ×g, 10 min, 4°C), and the supernatant was taken for analysis.

Regarding the feces sample, each 50 mg feces sample was mixed with 1,500 µl cooled methanol, and the mixture was homogenized (60 HZ, 30 s) and vortexed (60 s) to allow the proteins to precipitate. Then, the mixture was placed at −20°C (20 min) and centrifuged twice (6,200 ×g, 10 min, 4°C) and all the supernatants were filtered through a 0.22 µm microporous millipore membrane. Thereafter, it was centrifuged (6,200 ×g, 10 min, 4°C), and the supernatant was taken for analysis.

For the brain samples, each 50 mg brain samples was mixed with 1,000 µL of 50% cooled methanol solution, and the mixture was homogenized (35 HZ, 60 s, 4°C) and subsequently vortexed (60 s). Then, the sample was placed at −20 °C (20 min) to allow the proteins to precipitate. The mixtures were then centrifuged (6,200 ×g, 10 min, 4°C), and all the supernatants transferred to EP tubes were allowed to evaporate until dryness with the stream of nitrogen (35°C). After this, 200 µL of 80% methanol solution was added to the dried residue, and the mixture was completely dissolved and centrifuged (6,200 ×g, 10 min, 4°C), and the supernatant was taken for analysis.

When all the samples were processed separately, quality control (QC) samples with a volume of 120 µl were configured at the same time. The QC sample was a mixture of all samples.

2.6 H&E staining

For the H&E staining process, the rats’ brain was extracted from 4% paraformaldehyde solution and embedded in paraffin. Coronal sections of 4 μm thick were prepared for staining, and following steps were performed: section dewaxing, hematoxylin stain (5 min), water rinsing (30 min) for blue return, eosin staining for 2 min, dehydrating, and sealing. After that, each rat was randomly selected in four parts to take a photograph (the magnification of visualized damaged brain area was 20 times). The neurons were quantified in five randomly selected hippocampal CA1 micrographs from each segment, and the number of CA1 cells was counted. Only the cells with apparent nuclei and nucleoli were measured, while the other cells without prominent nuclei and nucleoli were excluded from the quantitative range.

2.7 UPLC-QTOF/MS-based metabolomics

Regarding the metabolomics analysis, metabolic profiling was analyzed on the UPLC-QTOF/MS system (Waters Corporation, Milford, MA, United States) in the presence of positive (ES+) and negative ionization sources (ES-) mode (Zhou et al., 2019). Water ACQUITY UPLC H-Class/Xevo used for UPLC analysis. For this purpose, 3 µL samples were injected into the Acquity HSS (High class) T3 column (2.1 mm × 100 mm, 1.8 μm, Waters). During this procedure, the UPLC analysis conditions were as follows: column temperature, 35°C; flow rate, 0.4 ml min⁻¹; and sample manager temperature, 4°C. For the metabolomics analysis of feces, the prevailing conditions were as: gradient system (ES+): methanol (A) and a mixture of 0.1% formic acid in water (B) (0–2 min, A, 5%; 2–3 min, A, 5–20%; 3–16, A, 20–80%; 16–17 min, A, 80–95%; 17–20 min, A, 95%; 20–20.1 min, A, 95–5%; 20.1–23 min, A, 5%); The gradient system (ES-): acetonitrile (A) and water (B) (0–2 min, A, 5%; 2–3 min, A, 5–20%; 3–16 min, A, 20–80%; 16–17 min, A, 80–95%; 17–20 min, A, 95%; 20–20.1 min, A, 95%–5%; and 20.1–23 min, A, 5%).

The parametric conditions for the brain metabolomics analysis were as follows: the gradient system (ES+): acetonitrile (A) and a mixture of 0.1% formic acid in water (B) (0–2 min, A, 2%; 2–6 min, A, 2–40%; 6–18 min, A, 40–98%; 18–21 min, A, 98%; 21–21.1 min, A, 98–2%; 21.1–23 min, A, 2%); and the gradient system (ES-): acetonitrile (A) and water (B) (0–2 min, A, 2%; 2–6 min, A, 2–60%; 6–14 min, A, 60–98%; 14–17 min, A, 98%; 17–17.1 min, A, 98–2%; 17.1–20 min, A, 2%).The conditions in the plasma metabolomics analysis were as follows: the gradient system (ES+): acetonitrile (A) and a mixture of 0.1% formic acid in water (B) (0–2 min, A, 5%; 2–18 min, A, 5–95%; 18–21 min, A, 95%; 21–21.1 min, A, 95–5%; 21.1–23 min, A, 5%); and the gradient system (ES-): acetonitrile (A) and water (B) (0–2 min, A, 5%; 2–18 min, A, 5–95%; 18–21 min, A, 95%; 18–20 min, A, 5%; 21–21.1 min, A, 95–5%; 21.1–23 min, A, 5%).

Similarly, the MS parameters were as follows: nebulizing gas, high-purity nitrogen (N₂); collision gas, ultra-high pure helium (He); desolvent gas flow, 800 (L⁻¹Hr); source temperature, 120°C; desolvent temperature, 500°C; conical air, 50 L⁻¹ h; collision energy, 20–45 v; capillary voltage: plasma sample, 2.0 kV; brain sample, 2.5 kV; feces sample, 1.5 kV; and scanning range, m/z 100–1,200 Da. Leucine keorphin (200 pg μL⁻¹) has a mass of 556.2771 in ES+ and 554.2615 in ES-. At the start of sequencing, we ran 10 QC samples and inserted a QC sample between every six sequences to further ensure the stability of the instrument.

2.8.16S rDNA gene sequencing analysis

The total DNA of feces samples was extracted using E.Z.N.A. ^® Stool DNA Kit. The V3–V4 region of the 16S rDNA gene was amplified by PCR. The bacterial primers were 341F (5′-CCTACGGGNGGCWGCAG-3′) and 805R (5′-GACTACHVGGGTATTCTAATCC-3′). Thereafter, the PCR products were purified (AMPure XT beads, Beckman Coulter Genomics, Danvers, MA, United States) and quantified (Qubit, Invitrogen, United States). Agilent 2100 Bioanalyzer (Agilent, United States) and Library Quantification Kit for Illumina (Kapa Biosciences, Woburn, MA, United States) were employed to evaluate the size and quantity of the amplicon library, respectively. The NovaSeq PE250 platform was used for sequencing. Quality filtering was performed according to fqtrim (v0.94). Vsearch software (v2.3.4) was used to filter the chimeric sequences.

2.9 Data analysis

MarkerLynx within MassLynx software version 4.1 (Waters Corp. Milford, United States) was employed to process all the collected UPLC-QTOF/MS data. Progenesis QI V 2.3 software (Waters, Milford, MA, United States) was used for peak detection and alignment, and EZinfo 3.0 software was used for the multivariate data analysis. To match and distinguish Rt-m/z pairs of different metabolites being selected in each group, an online database named HMDB (Human Metabolome Database, http://www.hmdb.ca/) was consulted and high-precision ion fragments were utilized. Enrichment analysis of the metabolomics pathways was performed through online websites, such as https://www.metaboanalyst.ca/and https://www.kegg.jp/, and the HMDB number of the identified differential metabolites were entered into the Metaboanalyst online data analysis software. Fisher’s exact test was selected as the enrichment method, and the rat Rattus norvegicus (rat) (KEGG) library was selected as the database for pathway analysis.

The data collected from the MWM and H&E staining were expressed as the mean ± standard deviation (SD), and the statistical significance was measured by GraphPad Prism 8.0.1. Spearman’s correlation analysis of the GM and endogenous metabolites was performed by the OmicStudio tools at https://www.omicstudio.cn/tool.

3 Results

3.1 Morris water maze test

The escape latency of each group gradually shortened with the increase in training time. From the third day, the rats in the model group took longer than for those in the blank group to find an escape platform (Figure 1A), indicating that the AD model was successfully established. Similarly, the rats from the Sch group showed a significant effect on the fourth day of training as well. In the space exploration test, residence time for the rats from the model group in the target quadrant (Figure 1B), and the number of target platform crossings (Figure 1C) were significantly reduced as compared to the blank and Sch groups. These results indicated that the rats in the blank and Sch groups had spent more time and stayed longer in the first quadrant, indicating better spatial memorability.

FIGURE 1

3.2 H&E staining

The CA1 pyramidal cell layer of the Sch group was observed to be thicker than that of the model group (Figures 2A–C). The number of neurons in the CA1 of the hippocampus was significantly more in the Sch group than the model group, indicating that schisandrin had a protective effect on the loss of hippocampal neurons in the diseased rats (Figure 2D, p < 0.0001).

FIGURE 2

3.7 Metabolomics research

3.7.1 UPLC-QTOF/MS results

Under the elution program, the base peak images (BPI) showed ideal separation results. The BPIs of the plasma, feces, and brain samples are shown in ES+ and ES- modes, respectively, in Figures 3A-F.

FIGURE 3

3.7.2 Identification of potential differential metabolites

The results from the PLS-DA (Partial least square-discriminant analysis) (Figures 4A-F) showed that the model and blank groups have a clear cluster classification, thus suggesting that there was a certain difference in metabolism between the model and blank group rats. According to OPLS-DA (Orthogonal partial least square-discriminant analysis), marked the compounds with VIP (Variable importance in the projection) > 1 on the S-plot, for which the OPLS-DA figures and S-plot are provided in the supplementary materials (Supplementary Figure S3). Furthermore, 44, 96, and 40 potential biomarkers (Mass error <5, CV (Cross-validated residuals) ≤ 30 (in QC), VIP >1, P (p value) < 0.05) were found in the feces, plasma, and brain samples, respectively, and are listed in Tables 1–3. The metabolites with significant changes in the metabolic level were analyzed by heatmap clustering, which showed that the metabolite levels in the Sch group have different degrees of callback as compared to the model group (Figures 5A-C).

FIGURE 4

FIGURE 5

3.7.3 Analysis of metabolic pathways in the plasma, brain, and feces

The results for the pathway analysis of the plasma and brain are expressed as bubble plots. A total of five metabolic pathways were enriched in the plasma (p ≤ 0.05, Figure 6A), namely, linoleic acid metabolism, sphingolipid metabolism, primary bile acid biosynthesis, glycerophospholipid metabolism, and alpha-linolenic acid metabolism. Similarly, five metabolic pathways were enriched in the brain (p ≤ 0.05, Figure 6B), namely, unsaturated fatty acid biosynthesis, glycerophospholipid metabolism, linoleic acid metabolism, alpha-linolenic acid metabolism, and sphingolipid metabolism. According to the results from pathway analysis of feces (Figure 6C), the metabolites were found mainly involved in the primary bile acid biosynthesis, sphingolipid metabolism, glycerophospholipid metabolism, unsaturated fatty acid biosynthesis, fatty acid biosynthesis, tryptophan metabolism, glycosylphosphatidylinositol (GPI)-anchor set biosynthesis, fructose, mannose metabolism, drug metabolism-cytochrome P450, arachidonic acid metabolism, fatty acid degradation, steroid biosynthesis, and steroid hormone biosynthesis.

FIGURE 6

3.3 16S rDNA gene sequencing analysis

3.3.1 Alpha diversity

Good’s coverage index refers to the microbial coverage, which is an index reflecting whether the sequencing results represent the actual microbial diversity in the sample. In the present study, the slope of the rarefaction curve tended to be 0, indicating that the sequencing was deeply saturated in our samples (Figure 7A).

FIGURE 7

3.3.2 Beta diversity analysis

The PCoA (Principal coordinate analysis) score chart (Figure 7B) revealed that the microbial community structure composition for the model, blank, and Sch groups had apparent clusters, thereby indicating significant differences in the GM of the three groups. The distance between the blank and Sch groups was close, indicating that the changes in the GM community for the Sch group tended to the level of the blank group, which suggests that schisandrin can promote the GM structure imbalance in rats with AD to approach the GM structure of rats in the blank group.

3.4 Species analysis

The results from this analysis revealed that the relative abundance of firmicutes phylum in the Sch and blank groups increased, whereas, the relative abundance of the bacteroides phylum decreased (Figures 8A,B). The results from the top 50 relatively abundant bacterial community indicated that multiple microorganisms showed significant differences at the genus level (p ≤ 0.05) (Figures 8C,D). The content of bacteroides, escherichia-Shigella, ruminococcaceae_unclassified, oscillibacter, and intestinimonas in the model group tyzzerella, alloprevotella, firmicutes_unclassified, and quinella increased as compared to the blank group. However, the content of lactobacillus and muribaculaceae_unclassified decreased. However, the lactobacillus content, ruminococcaceae_UCG-005, muribaculaceae_unclassified, clostridium, and christensenellaceae_unclassified increased in the Sch group as compared to the model group, while the content of escherichia-Shigella, enterorhabdus, bacteroides, ruminococcaceae_unclassified, and alloprevotella decreased in the same group. To further study the changes being occurring in the GM at the genus level, the top 50 bacteria with the highest relative abundance in all the three groups were selected for Spearman’s correlation analysis (Figures 9A-C) (p ≤ 0.05, rho ±0.5). The correlation results showed that the mutual symbiosis relationships among the bacteria in all the three groups has changed significantly. AD severely interfered with the symbiosis of GM in the diseased rats. Schisandrin showed improvement in the imbalance of symbiosis, for example, in the AD model group, eisenbergiella, enterorhabdus, firmicutes_unclassified, marvinbryantia, quinella, candidatus_saccharimonas, and the other GMs have established new symbiotic relationships. The newly established symbiotic relationship of these GMs have disappeared right after the schisandrin treatment. Some GMs, which were also treated with schisandrin, have re-established the same correlation as the blank group compared to the model group. For example, clostridiales_unclassified and oscillibacter and eubacterium_xylanophilum_group, erysipelatoclostridium and eubacterium_nodatum_group have re-established a positive correlation. On the other hand, phascolarctobacterium and romboutsia have established a negative correlation, indicating that schisandrin had strengthen the symbiotic relationships between them. The changes occurred in the correlation between GM indicate that AD might have interfered with the symbiotic relationship between the GM. Schisandrin had improved the imbalance of this relationship to a certain extent, thereby suggesting that this balanced relationship is closely related to the development mechanism of schisandrin to play a protective effect on the rats with AD.

FIGURE 8

FIGURE 9

3.5 Predictive metagenomics analysis of the GM

Based on PICRUSt2, we observed the potential functional changes in the GM of rats from the KEGG database (Figure 10A). A number of metabolic pathways and functions of the GM in the model and Sch groups were significantly enriched, including lipid metabolism, steroid hormone biosynthesis, arachidonic acid metabolism, biosynthesis of unsaturated fatty acids, bacterial invasion of epithelial cells, and steroid biosynthesis.

FIGURE 10

3.6 Relationship between the GM and cognitive impairment of rats with AD

To reveal the pathological significance of the changes being occurring in the GM of rats, top 50 relatively abundant GMs were selected to analyze the correlation between their abundance and behavior (Figure 10B). The correlation results revealed that among the tract bacteria, escherichia-Shigella, rumenococcus_1, lactobacillus, lactobacillus_RF39_unclassified, enterobacter, rumbutzia, lacetospiraceae_NK4A136_group, tyzzerella, lactococcus, enterobacter, quininella, desulfovibrio, and bacteroides found to have a significant correlation with the behavior, lactobacillus and bacteroides were the two most common abundant bacteria among them, thereby suggesting their critical role in the mechanism of schisandrin treatment.

3.8 Potential relationships between the host metabolites and GM

To understand the mechanism of metabolic disorder in the GM caused by the schisandrin in improving the symptoms of diseased rats, the potential differential metabolites enriched in the brain, feces, and plasma sample pathways were selected. Spearman’s correlation coefficient between the metabolites and the top 50 abundant GMs was calculated (Figure 11). The results indicate that most of the GM had a significant correlation with the metabolites in the body. In the plasma samples, citric acid, phenylpyruvic acid, tetradecanoyl-CoA, and 3-Oxo-OPC6-CoA were found positively correlated with christensenellaceae_R-7_group, whereas, christensenellaceae_unclassified, lachnospiraceae_unclassified, oscillibacter, intestinimonas, ruminiclostribacter, intestinimonas, ruminiclostriovic acid, tetradecanoyl-Phenylco Oxo-OPC6-CoA, corticosterone, PC (14:0/22:4 (7Z, 10Z, 13Z, 16Z)), arachidonate, 8Z, 11Z, and 14Z-eicosatrienoyl-CoA were found negatively correlated. In the brain, arachidonate, docosahexaenoic acid, oleic acid, linoleic acid, glycocholic acid, galactosylceramide (d18:1/18:0) were associated with the bacteroides, christensenellaceae_unclassified, phascolarctobacterium, parabacteroides, lachnospiraceae_unclassified, and oscillibacterium, whereas, enterbacillibacteria, edus _ruminantium_group, and peptococcaceae_unclassified were found negatively correlated. In the feces samples, alloprevotella and ruminiclostridium_5 were negatively correlated with myristic acid and palmitic acid, whereas, hydroxylidocaine was positively correlated with LysoPC (P-18:0), dihydrocortisol, and glycocholic acid.

FIGURE 11

4 Discussion

AD is a neurodegenerative disease with insidious onset and reportedly has an excellent therapeutic effect on AD. Schisandrin can inhibit neuroinflammation, Aβ deposition, and helps improving the cognitive impairment in AD. In this study, we have tried to comprehensively explore the pathogenesis of AD and the protective effect of schisandrin on rats with AD through the metabolomics analysis of brain, plasma, and fecal nontargeted combined with 16S rDNA gene sequencing technology systematically from the brain–gut axis (Figure 12). In this regard, an AD rat model was developed by injecting D-gal intraperitoneally, orally administering AlCl₃, and bilaterally injecting Aβ_25-35 intraventricularly to induce the AD. AlCl₃ can induce the rats to produce progressive neurodegeneration similar to that observed in AD (Ulusoy et al., 2015), D-gal is an aging agent that causes cognitive impairment (Cao et al., 2019), and Aβ deposition is the main pathology of AD (Gupta, Sil et al., 2018). Collectively, these three methodologies have a synergistic effect on the establishment of AD rat models (Chiroma et al., 2018; McDonald et al., 2021).

FIGURE 12

Identified metabolites with large fluctuations in the metabolic levels were mainly enriched in unsaturated fatty acid biosynthesis, sphingolipid metabolism, primary bile acid biosynthesis, glycerophospholipid metabolism, and tryptophan metabolism, which were also significantly correlated with the identified GM. Several studies have shown that multiple GMs were associated with the occurrence and development of the aforementioned metabolic pathways. The prediction of GMs pathway further confirmed that the aforementioned metabolic pathways were directly involved in the anti-AD mechanism of schisandrin. Therefore, we focused on the stated metabolic pathways to discuss the mode of action of schisandrin.

The results from the MWM test and H&E staining suggest that schisandrin is directly involved in the improvement of cognitive impairment caused by AD, and could also reduce the loss of hippocampal neurons in rats. The Firmicutes/Bacteroidetes (F/B) ratio of the model group was also decreased, indicating an imbalance in the GM. At the genus level, some microbes showed significant differences among all the three study groups. Muribaculaceae family might have played an essential role in degrading the mucosal layer (Evans et al., 2014; Rooks et al., 2014). Similarly, desulfovibrio can cause a significant reduction in the short-chain fatty acids (SCFAs) levels (Sawin et al., 2015), because the enterobacteriaceae and desulfovibrio bacteria are the main lipopolysaccharide (LPS) endotoxin-producing bacteria (Zhao 2013). Schisandrin inhibits desulfovibrio bacteria and helps in protecting the integrity of intestinal endothelium. Based on the inhibitory effect of schisandrin on desulfovibrio, it can be deduced that it has profound effects on the treatment of AD. Similarly, christensenellaceae is one of the five taxa considered as a sign of healthy intestines. In our case, we have observed that the christensenellaceae content in the Sch group was increased significantly, thereby indicating that schisandrin can restore the intestinal environment to a healthy state by multiplying the content of christensenellaceae in the gut. Furthermore, ruminococcaceae is responsible for the fermentation of dietary fiber and other plant components, and schisandrin has caused a significant increase in the ruminococcaceae content, thus indicating that it can metabolize to produce more SCFAs. Escherichia–shigella is positively correlated with the expression of cytokine and is a typical intestinal pathogen (Chatterjee et al., 2018). We have observed a significant decrease in the escherichia-Shigella content after treating with schisandrin. Taken together, the GM function in diseased rats treated with schisandrin is found to be involved in the bacterial chemotaxis, arachidonic acid metabolism, biosynthesis of unsaturated fatty acids, and bacterial invasion of epithelial cells, which further confirms the improvement of gut inflammation caused by the GM invasion.

It is well known that the supplementation of docosahexaenoic acid (DHA) and eicosapentaenoic acid (EPA) can improve the cognitive impairment in AD patients and these can also reduce the neuroinflammation. Elevated palmitic acid levels are related to the cognitive impairment in humans. Oleic acid improves amyloidosis in AD cells and has neurotrophic effects. Linoleic acid can reverse the inflammatory response of microglia and help in controlling the inflammatory cell response in the brain by targeting the myeloid cells in the gut. Serum dihomogamma-linolenic acid (DGLA), with its anti-inflammatory effects, is associated with the cognitive functions and mild cognitive impairment at much reduced levels in patients with coronary artery disease (Ishihara et al., 2020). Furthermore, stearidonic acid can increase the concentration of EPA in the plasma cells, and arachidonic acid on the other hand is directly or indirectly converted to eicosanoid, thereby leading to cause cell impairment and mediates inflammation. In the present study, schisandrin treated rats had higher levels of DHA, oleic acid, and linoleic acid in the brain, whereas, the palmitic acid was recorded at lower levels. Similarly, in the plasma, Sch group had higher levels of DGLA, stearidonic acid, and lower levels of arachidonate. Regarding the feces, palmitic acid content was reduced and the DHA content was increased. Collectively, our results suggest that schisandrin can improve AD’s inflammatory state, can maintain the metabolic balance, and can impart nutritional effect on the brain.

According to previous studies, the levels of phosphatidylcholine (PC), phosphatidylethanolamine (PE), and phosphatidylinositol (PI) in the neural membranes of different regions of AD patients were reduced as compared to the age-matched control human brains (Söderberg et al., 1991). In this study, the glycerophospholipid metabolic pathways in the brain, plasma, and feces were found disordered, and through the regulatory effect of schisandrin, the metabolic level gradually returned to normal.

Sphingolipid (SL) metabolism is related to the level of amyloid β in the cerebrospinal fluid of AD (Fonteh et al., 2015). Also, the GM of phylum bacteroides can produce SL. Through the prediction analysis of the GM function, it was observed that GM participates in the process of sphingolipid metabolism, and after production sphingolipids enter the host metabolic pathway to modulate the ceramide levels (Johnson et al., 2020). The interaction of sphingolipids and bacteria in the gut stimulates the activity of iNKT cells through CD1b, thereby exerting an immunomodulatory effect (An et al., 2014). In the nervous system, sphingosine-1-phosphate (S1P) can influence the blood–brain barrier when its level decreases in the blood (Vutukuri et al., 2018). In this study, the AD model group showed significant changes in the SL metabolism pathway, such as the decrease in S1P content in the plasma and the decrease in Galactosylceramide (d18:1/18:0) content in the brain. The metabolic level was reversed after schisandrin treatment, suggesting that it can protect diseased rats by changing the SL metabolism level and by improving the influence of microorganisms on SL metabolism.

Cholic acid and alpha-muricholic acid can inhibit the proliferation of Escherichia coli, streptococcus, and other harmful bacteria (Watanabe et al., 2017). Previous findings have elaborated that a decrease in the blood cholic acid is related to the reduction of the overall glucose metabolism in the brain, and the lower cholic acid content and higher chenodeoxycholic acid is associated with more severe brain atrophy (Nho et al., 2019). Similar to previous studies, the cholic acid content in the feces of diseased rats was reduced in our case, and it was adjusted after schisandrin treatment. The harmful bacteria, romboutsia and enterorhabdus, showed a significant negative correlation with cholic acid in feces, which further strengthens that excessive cholic acid might be involved in reducing some harmful bacteria. Furthermore, in the plasma samples, the cholic acid content in the blank and Sch groups was increased, thereby indicating that schisandrin can play a protective role by regulating the bile acid biosynthesis pathway in rats.

Tryptophan participates in the kynurenine pathway, and its metabolites can affect the central nervous system (CNS) through the metabolic pathways. In the model group, the content of the aforementioned metabolites was abnormal, and it was improved after schisandrin treatment was applied. Tryptophan is the only amino acid with an indole structure. Indoleacrylic acid transformed by the GM can promote mucosal homeostasis, and protect mice from bacterial infections, thereby reduce the gut inflammation (Wlodarska et al., 2017). Indoleacrylic acid has anti-inflammatory and antioxidant effect in human peripheral blood mononuclear cells activated by LPS by reducing interleukin (IL)-6 and IL-1β and activating the NRF2-ARE pathway (Joshi and Johnson 2012). In our case, the indoleacrylic acid content in the feces of rats was decreased in the model group. After schisandrin treatment, the indoleacrylic acid content was increased, thus suggesting that schisandrin can protect the gut–blood barrier through anti-inflammatory effect in protecting rats with AD.

5 Conclusion

Although the therapeutic effect of schisandrin had been confirmed in earlier studies, the underlying mechanism to improve the symptoms in rats diagnosed with AD was still unclear. In this context, our study provides convincing results about the mode of action of schisandrin. The metabolomics studies of the plasma, feces, and brain showed that schisandrin had reversed the imbalance of the metabolic level of the brain–gut axis in the diseased rats, and furthermore, it was also involved in the multiple metabolic pathways in AD. Moreover, 16S rDNA gene sequencing results revealed that schisandrin improved the GM imbalance and gut–blood barrier damage caused by gut inflammation in rats. Finally, Spearman’s correlation analysis showed that the changes in intestinal microbes were significantly related to the overall metabolic level in the body. The overall data underline that although schisandrin has low oral bioavailability, it can still improve the metabolic imbalance of rats with AD by strengthening the gut microenvironment.

References

• 1

  AlseekhS.FernieA. R. (2018). Metabolomics 20 years on: What have we learned and what hurdles remain?Plant J.94 (6), 933–942. 10.1111/tpj.13950

  • CrossRef
  • Google Scholar

• 2

  AnD.OhS. F.OlszakT.NevesJ. F.AvciF. Y.Erturk-HasdemirD.et al (2014). Sphingolipids from a symbiotic microbe regulate homeostasis of host intestinal natural killer T cells. Cell156 (1-2), 123–133. 10.1016/j.cell.2013.11.042

  • CrossRef
  • Google Scholar

• 3

  BellJ. S.SpencerJ. I.YatesR. L.YeeS. A.JacobsB. M.DeLucaG. C. (2019). Invited Review: from nose to gut - the role of the microbiome in neurological disease. Neuropathol. Appl. Neurobiol.45 (3), 195–215. 10.1111/nan.12520

  • CrossRef
  • Google Scholar

• 4

  CaoP.ZhangJ.HuangY.FangY.LyuJ.ShenY. (2019). The age-related changes and differences in energy metabolism and glutamate-glutamine recycling in the d-gal-induced and naturally occurring senescent astrocytes in vitro. Exp. Gerontol.118, 9–18. 10.1016/j.exger.2018.12.018

  • CrossRef
  • Google Scholar

• 5

  ChatterjeeR.ShreenivasM. M.SunilR.ChakravorttyD. (2018). Enteropathogens: tuning their gene expression for hassle-free survival. Front. Microbiol.9, 3303. 10.3389/fmicb.2018.03303

  • CrossRef
  • Google Scholar

• 6

  ChiromaS. M.Mohd MoklasM. A.Mat TaibC. N.BaharuldinM. T. H.AmonZ. (2018). d-galactose and aluminium chloride induced rat model with cognitive impairments. Biomed. Pharmacother.103, 1602–1608. 10.1016/j.biopha.2018.04.152

  • CrossRef
  • Google Scholar

• 7

  CryanJ. F.O'RiordanK. J.SandhuK.PetersonV.DinanT. G. (2020). The gut microbiome in neurological disorders. Lancet. Neurol.19 (2), 179–194. 10.1016/S1474-4422(19)30356-4

  • CrossRef
  • Google Scholar

• 8

  EvansC. C.LePardK. J.KwakJ. W.StancukasM. C.LaskowskiS.DoughertyJ.et al (2014). Exercise prevents weight gain and alters the gut microbiota in a mouse model of high fat diet-induced obesity. PLoS One9 (3), 92193. 10.1371/journal.pone.0092193

  • CrossRef
  • Google Scholar

• 9

  FontehA. N.OrmsethC.ChiangJ.CipollaM.ArakakiX.HarringtonM. G. (2015). Sphingolipid metabolism correlates with cerebrospinal fluid beta amyloid levels in Alzheimer's disease. PLoS One10 (5), e0125597. 10.1371/journal.pone.0125597

  • CrossRef
  • Google Scholar

• 10

  GuptaP.SilS.GhoshR.GhoshA.GhoshT. (2018). Intracerebroventricular aβ-induced neuroinflammation alters peripheral immune responses in rats. J. Mol. Neurosci.66 (4), 572–586. 10.1007/s12031-018-1189-9

  • CrossRef
  • Google Scholar

• 11

  IshiharaK.IzawaK. P.KitamuraM.ShimogaiT.KanejimaY.MorisawaT.et al (2020). Serum concentration of dihomo-γ-linolenic acid is associated with cognitive function and mild cognitive impairment in coronary artery disease patients. Prostagl. Leukot. Essent. Fat. Acids158, 102038. 10.1016/j.plefa.2019.102038

  • CrossRef
  • Google Scholar

• 12

  JohnsonE. L.HeaverS. L.WatersJ. L.KimB. I.BretinA.GoodmanA. L.et al (2020). Sphingolipids produced by gut bacteria enter host metabolic pathways impacting ceramide levels. Nat. Commun.11 (1), 2471. 10.1038/s41467-020-16274-w

  • CrossRef
  • Google Scholar

• 13

  JoshiG.JohnsonJ. A. (2012). The nrf2-ARE pathway: A valuable therapeutic target for the treatment of neurodegenerative diseases. Recent Pat. CNS Drug Discov.7 (3), 218–229. 10.2174/157488912803252023

  • CrossRef
  • Google Scholar

• 14

  KunduP.LeeH. U.Garcia-PerezI.TayE. X. Y.KimH.FaylonL. E.et al (2019). Neurogenesis and prolongevity signaling in young germ-free mice transplanted with the gut microbiota of old mice. Sci. Transl. Med.11 (518), eaau4760. 10.1126/scitranslmed.aau4760

  • CrossRef
  • Google Scholar

• 15

  LeeS. R.ParkM. Y.YangH.LeeG. S.AnB. S.ParkB. K.et al (2018). 5α-dihydrotestosterone reduces renal Cyp24a1 expression via suppression of progesterone receptor. J. Mol. Endocrinol.60 (2), 159–170. 10.1530/JME-17-0187

  • CrossRef
  • Google Scholar

• 16

  LiB.HeY.MaJ.HuangP.DuJ.CaoL.et al (2019). Mild cognitive impairment has similar alterations as Alzheimer's disease in gut microbiota. Alzheimers Dement.15 (10), 1357–1366. 10.1016/j.jalz.2019.07.002

  • CrossRef
  • Google Scholar

• 17

  LuW. W.FuT. X.WangQ.ChenY. L.LiT. Y.WuG. L. (2020). The effect of total glucoside of paeony on gut microbiota in NOD mice with sjögren's syndrome based on high-throughput sequencing of 16SrRNA gene. Chin. Med.15, 61. 10.1186/s13020-020-00342-w

  • CrossRef
  • Google Scholar

• 18

  MarizzoniM.CattaneoA.MirabelliP.FestariC.LopizzoN.NicolosiV.et al (2020). Short-chain fatty acids and lipopolysaccharide as mediators between gut dysbiosis and amyloid pathology in alzheimer's disease. J. Alzheimers Dis.78 (2), 683–697. 10.3233/JAD-200306

  • CrossRef
  • Google Scholar

• 19

  McDonaldJ. B.DhakalS.MacreadieI. (2021). A toxic synergy between aluminium and amyloid beta in yeast. Int. J. Mol. Sci.22 (4), 1835. 10.3390/ijms22041835

  • CrossRef
  • Google Scholar

• 20

  NhoK.Kueider-PaisleyA.MahmoudianDehkordiS.ArnoldM.RisacherS. L.LouieG.et al (2019). Altered bile acid profile in mild cognitive impairment and Alzheimer's disease: relationship to neuroimaging and CSF biomarkers. Alzheimers Dement.15 (2), 232–244. 10.1016/j.jalz.2018.08.012

  • CrossRef
  • Google Scholar

• 21

  PereiraF. C.WasmundK.CobankovicI.JehmlichN.HerboldC. W.LeeK. S.et al (2020). Rational design of a microbial consortium of mucosal sugar utilizers reduces clostridiodes difficile colonization. Nat. Commun.11 (1), 5104. 10.1038/s41467-020-18928-1

  • CrossRef
  • Google Scholar

• 22

  RooksM. G.VeigaP.Wardwell-ScottL. H.TickleT.SegataN.MichaudM.et al (2014). Gut microbiome composition and function in experimental colitis during active disease and treatment-induced remission. Isme J.8 (7), 1403–1417. 10.1038/ismej.2014.3

  • CrossRef
  • Google Scholar

• 23

  SawinE. A.De WolfeT. J.AktasB.StroupB. M.MuraliS. G.SteeleJ. L.et al (2015). Glycomacropeptide is a prebiotic that reduces Desulfovibrio bacteria, increases cecal short-chain fatty acids, and is anti-inflammatory in mice. Am. J. Physiol. Gastrointest. Liver Physiol.309 (7), G590–G601. 10.1152/ajpgi.00211.2015

  • CrossRef
  • Google Scholar

• 24

  SöderbergM.EdlundC.KristenssonK.DallnerG. (1991). Fatty acid composition of brain phospholipids in aging and in Alzheimer's disease. Lipids26 (6), 421–425. 10.1007/BF02536067

  • CrossRef
  • Google Scholar

• 25

  UlusoyH. B.SonmezM. F.KilicE.CaliskanK.KaracaB.KaraM.et al (2015). Intraperitoneal administration of low dose aluminium in the rat: how good is it to produce a model for alzheimer disease. Arch. Ital. Biol.153 (4), 266–278. 10.12871/00039829201543

  • CrossRef
  • Google Scholar

• 26

  VutukuriR.BrunkhorstR.KestnerR. I.HansenL.BouzasN. F.PfeilschifterJ.et al (2018). Alteration of sphingolipid metabolism as a putative mechanism underlying LPS-induced BBB disruption. J. Neurochem.144 (2), 172–185. 10.1111/jnc.14236

  • CrossRef
  • Google Scholar

• 27

  WatanabeM.FukiyaS.YokotaA. (2017). Comprehensive evaluation of the bactericidal activities of free bile acids in the large intestine of humans and rodents. J. Lipid Res.58 (6), 1143–1152. 10.1194/jlr.M075143

  • CrossRef
  • Google Scholar

• 28

  WeiB. B.LiuM. Y.ChenZ. X.WeiM. J. (2018). Schisandrin ameliorates cognitive impairment and attenuates Aβ deposition in APP/PS1 transgenic mice: Involvement of adjusting neurotransmitters and their metabolite changes in the brain. Acta Pharmacol. Sin.39 (4), 616–625. 10.1038/aps.2017.135

  • CrossRef
  • Google Scholar

• 29

  WeinerM. W.VeitchD. P.AisenP. S.BeckettL. A.CairnsN. J.GreenR. C.et al (2017). Recent publications from the alzheimer's disease neuroimaging initiative: reviewing progress toward improved AD clinical trials. Alzheimers Dement.13 (4), e1–e85. 10.1016/j.jalz.2016.11.007

  • CrossRef
  • Google Scholar

• 30

  WlodarskaM.LuoC.KoldeR.d'HennezelE.AnnandJ. W.HeimC. E.et al (2017). Indoleacrylic acid produced by commensal peptostreptococcus species suppresses inflammation. Cell Host Microbe22 (1), 25–37. e26. 10.1016/j.chom.2017.06.007

  • CrossRef
  • Google Scholar

• 31

  ZhaoL. (2013). The gut microbiota and obesity: from correlation to causality. Nat. Rev. Microbiol.11 (9), 639–647. 10.1038/nrmicro3089

  • CrossRef
  • Google Scholar

• 32

  ZhaoZ. Y.ZhangY. Q.ZhangY. H.WeiX. Z.WangH.ZhangM.et al (2019). The protective underlying mechanisms of Schisandrin on SH-SY5Y cell model of Alzheimer's disease. J. Toxicol. Environ. Health. A82 (19), 1019–1026. 10.1080/15287394.2019.1684007

  • CrossRef
  • Google Scholar

• 33

  ZhouH.TaiJ.XuH.LuX.MengD. (2019). Xanthoceraside could ameliorate alzheimer's disease symptoms of rats by affecting the gut microbiota composition and modulating the endogenous metabolite levels. Front. Pharmacol.10, 1035. 10.3389/fphar.2019.01035

  • CrossRef
  • Google Scholar