Integrated Metabolomic and Lipidomic Analysis Reveals the Neuroprotective Mechanisms of Bushen Tiansui Formula in an Aβ1-42-Induced Rat Model of Alzheimer's Disease

Bushen Tiansui Formula (BSTSF) is a traditional Chinese medicine prescription. It has been widely applied to treat Alzheimer's disease (AD) in the clinic; however, the mechanisms underlying its effects remain largely unknown. In this study, we used a rat AD model to study the effects of BSTSF on cognitive performance, and UPLC-MS/MS-based metabolomic and lipidomic analysis was further performed to identify significantly altered metabolites in the cerebral cortices of AD rats and determine the effects of BSTSF on the metabolomic and lipidomic profiles in the cerebral cortices of these animals. The results revealed that the levels of 47 metabolites and 30 lipids primarily associated with sphingolipid metabolism, glycerophospholipid metabolism, and linoleic acid metabolism were significantly changed in the cerebral cortices of AD rats. Among the altered lipids, ceramides, phosphatidylethanolamines, lysophosphatidylethanolamines, phosphatidylcholines, lysophosphatidylcholines, phosphatidylserines, sphingomyelins, and phosphatidylglycerols showed robust changes. Moreover, 34 differential endogenous metabolites and 21 lipids, of which the levels were mostly improved in the BSTSF treatment group, were identified as potential therapeutic targets of BSTSF against AD. Our results suggest that lipid metabolism is highly dysregulated in the cerebral cortices of AD rats, and BSTSF may exert its neuroprotective mechanisms by restoring metabolic balance, including that of sphingolipid metabolism, glycerophospholipid metabolism, alanine, aspartate, and glutamate metabolism, and D-glutamine and D-glutamate metabolism. Our data may lead to a deeper understanding of the AD-associated metabolic profile and shed new light on the mechanism underlying the therapeutic effects of BSTSF.


Introduction
Alzheimer's disease (AD) is a common neurodegenerative disorder of the central nervous system characterized by progressive memory loss, cognitive impairment, abnormal behavior, and personality disorders. Dementia, including that related to AD, is the fifth leading cause of death worldwide, and 40-50 million people are thought to be affected by this condition [1], making it a major and increasing global health challenge. However, an effective treatment for AD remains elusive [2].
Traditional Chinese medicines (TCMs) have been used in the treatment of dementia for thousands of years. The Bushen Tiansui Formula (BSTSF, also known as "Naoling decoction") is derived from Sagacious Confucius' Pillow Elixir, a classic Chinese medicinal formula mainly used to treat cognitive decline [3]. This formula comprises six herbs, including Epimedium acuminatum Franch. (Yinyanghuo), Fallopia multiflora (Thunb.) Harald. (Heshouwu), Polygala tenuifolia Willd. (Yuanzhi), Acorus tatarinowii Schott. (Shichangpu), Plastrum testudinis (Guiban), and Ossa draconis (Longgu). We recently showed that BSTSF could improve learning and memory deficits in AD model rats through regulating serum lipid metabolism and the amino acid metabolic pathway [4]. However, the mechanism underlying the therapeutic effects of BSTSF, especially its effects on metabolic stress and impaired lipid metabolism in brain tissue, is poorly understood.
Senile plaques, neurofibrillary tangles, and lipid granule accumulation were the three defining neuropathological features in the cerebral cortex of AD patients identified in the original analysis by Alois Alzheimer [5]. Subsequently, a large number of studies have shown that beta-amyloid (Aβ) plaques accumulate in the cerebral cortex in the early stages of AD [6][7][8]. These observations highlight that the cerebral cortex is the main pathological region in brain tissue of AD. The prefrontal cortex (PFC) is implicated in cognitive processes including working memory, temporal processing, decision-making, flexibility, and goal-oriented behavior [9]. Alterations in prefrontal cortex (PFC) function and abnormalities in its interactions with other brain areas (i.e., the hippocampus) have been related to Alzheimer's disease (AD) [10]. Several potential biomarkers have been proposed for preclinical AD and are mainly related to lipid metabolism. Perturbations of sphingolipid metabolism in brain tissue are consistently associated with endophenotypes across preclinical and prodromal AD, indicating that sphingolipids may be biologically relevant biomarkers for early AD detection [11]. Moreover, phosphatidylcholine breakdown may be mediated by phospholipase A2, leading to significantly elevated levels of glycerophosphocholine in human cerebrospinal fluid [12], while amino acids such as valine, arginine, and histidine are also associated with AD [13,14]. These observations suggest that severe metabolic disorder and dysregulated lipid metabolism may have an important role in the pathogenesis of AD.
In this study, an ultra-high-performance liquid chromatography-mass spectrometry-(UHPLC-MS-) based metabolomic and lipidomic analysis was performed in cerebral cortices of control and AD model rats. Furthermore, the therapeutic effects of BSTSF and the mechanism underlying its ameliorating effects on the pathogenesis of AD were also explored for the first time using a metabolomic strategy. This study may provide a basis for a better understanding of the AD metabolic profile and novel insight into the clinical utility of BSTSF.

Material and Methods
2.1. Preparation of BSTSF. The six herbs (H. epimedii, P. multiflorum, tortoise plastron, O. draconis, P. tenuifolia, and R. acori graminei) comprising BSTSF were mixed at a ratio of 3 : 3 : 4 : 4 : 2 : 2. To ensure the quality of the herbal medicine, all herbs were obtained from the TCM pharmacy of the Second Xiangya Hospital of Central South University (CSU, Changsha, China) and were authoritatively identified by Professor Suiyu Hu. For details of the preparation and quality control of lyophilized powder, refer to our previous publication [15].

Animals and Experimental Design.
Adult male Sprague-Dawley rats, weighing 180-220 g, used in our experiments, were purchased from the Laboratory Animal Centre of Central South University (Changsha, China). All animal experiments were performed following guidelines from the Committee on the Ethics of Animal Experiments of Central South University. Rats were housed under standard animal house conditions and randomly allocated into one of three groups: sham, AD, and BSTSF. For the AD and BSTSF groups, rats were injected intracerebroventricularly with oligomeric Aβ 1-42 to generate a validated AD model, as we previously described [16]. The sham rats were injected bilaterally with the vehicle into the lateral ventricles. According to our previous study, the BSTSF group was intragastrically administrated with 27 g/kg BSTSF once daily from 1 to 28 days, whereas the sham and AD groups were intragastrically administrated with an equal volume of distilled water. By referring to the calculation formula from the Experiment Methodology of Pharmacology, the conversion factor between rats (200 g) and humans (70 kg) is 0.018; therefore, the calculated gavage dose of BSTSF for rats is about 9 g/kg/d. Our previous study explored the efficacy of three doses (9 g/kg/d, 27 g/kg/d, and 54 g/kg/d) of this prescription, and the result demonstrated that BSTSF owns optimal efficacy when it is administered at three times the regular dose; therefore, a dose of 27 g/kg/d was chosen for the experiments in the current study [15].

Morris Water Maze
Test. The Morris water maze (MWM) test was used to assess the hippocampusdependent learning and memory ability of the rats, as described in our previous study with minor modifications [17]. In brief, a spatial acquisition test was carried out from the 24th to 28th days after Aβ 1-42 infusion, and animals were subjected to a five-day memory acquisition experiment to assess their spatial learning ability. Subsequent spatial probe experiments were conducted on day 29, to determine rat spatial memory retention ability. We applied SuperMaze video tracking and analysis systems to analyze experimental parameters (XR-XM101, Shanghai Softmaze Information Technology Co. Ltd., China).

Sample Preparation.
After the MWM test, the prefrontal cortices were harvested and immediately washed with precooled physiological saline and stored at −80°C for further analysis. Samples (50 mg) were accurately weighed, and the metabolites extracted using 400 μL of a methanol : water (4 : 1, v/v) solution. The mixture was allowed to settle at −20°C and treated using a high-throughput Wonbio-96c tissue crusher (Shanghai Wanbo Biotechnology Co., Ltd., China) at 60 Hz for 6 min, vortexed for 30 s, and sonicated at 40 kHz for 10 min at −20°C. This step was performed three times. The samples were placed at −20°C for 30 min to precipitate the proteins. After centrifugation at 13,000 × g at 4°C for 15 min, the supernatants were carefully transferred  ; other data were analyzed by one-way ANOVA followed by least significant difference tests). * p < 0:05, * * p < 0:01 vs. sham group; # p < 0:05, ## p < 0:05 vs. AD group.

Data Processing.
After UPLC-MS/MS analysis, the raw data were imported into Progenesis QI 2.3 (Nonlinear Dynamics, Waters) for peak detection and alignment. The preprocessing results generated a data matrix that consisted of the retention time (RT), mass-to-charge ratio (m/z) values, and peak intensity. After filtering, half of the minimum metabolite values were imputed for specific samples in which the metabolite levels fell below the lower limit of quantitation and each metabolic feature was normalized by sum. The internal standard was used for data QC (reproducibility), and metabolic features with a QC relative standard deviation ðRSDÞ > 30% were discarded. Following normalization and imputation, statistical analysis was performed on log-transformed data to identify significant differences in metabolite levels between comparable groups. The mass spectra of these metabolic features were identified by using was applied to obtain an overview of the metabolic data, general clustering, trends, and outliers. All the metabolite variables were scaled to unit variances before PCA. Orthogonal partial least squares discriminant analysis (OPLS-DA) was used to determine global metabolic changes between comparable groups. All the metabolite variables were scaled using Pareto scaling before the OPLS-DA. Model validity was evaluated from model parameters R2 and Q2, which provide information for the interpretability and predictability, respectively, of the model and avoid the risk of overfitting. Variable importance in the projection (VIP) was calculated in the OPLS-DA model. p values were estimated with paired Student's t-tests for single-dimensional statistical analysis. Significance among groups was assumed with VIP scores > 1 and p values < 0.05.

Metabolic Pathway Analysis.
Significantly altered metabolite data were imported into MetaboAnalyst 3.5 (https://www.metaboanalyst.ca) to investigate the therapeutic mechanisms related to BSTSF treatment. The impact value threshold calculated from pathway topology analysis was set to 0.10, and a raw p value < 0.05 was regarded as significant.  Thermo Scientific) was used for peak identification, lipid identification (secondary identification), peak extraction, peak alignment, and quantitative processing. The main parameters were as follows: precursor tolerance, 5 ppm; product tolerance, 5 ppm; and product ion threshold, 5%. Lipid molecules with RSD > 30% were deleted.
2.6. Statistical Analysis. All data are presented as mean ± standard error of the mean (SEM) and were analyzed using SPSS 22.0 software (IBM Corp., Armonk, NY, USA). Student's t-test was carried out for comparisons between two groups, whereas ANOVA was conducted for comparisons of repeated measures. p < 0:05 was defined as indicating a statistically significant difference.

BSTSF Ameliorates Learning and Memory Deficits of
Aβ1-42-Induced AD Rats. First, we examined the effects of BSTSF on learning and memory ability in AD model rats using the MWM test. As shown in Figures 1(a) and 1(b), escape latency time gradually decreased in all groups over time; however, rats in the BSTSF group had significantly lower escape latency than those in the AD group during the last three training days. These results indicate that BSTSF can significantly alleviate the impaired learning ability induced by Aβ  . In the probe test, AD rats crossed the platform fewer times and spent less time in the target quadrant, suggesting that their memory capacity was significantly decreased, while BSTSF treatment significantly reversed these defects in AD rats (Figures 1(c) and 1(d)). Together, these results suggested that BSTSF treatment can ameliorate Aβ 1-42 -induced spatial learning and memory impairment.

Metabolite Identification and Multivariate Statistical
Analysis. In this study, we identified 228 metabolites in a positive-ion mode and 287 in a negative-ion mode. An unsupervised PCA recognition model was generated for the whole dataset to evaluate the clustering trend of the samples with multidimensional data. The clustered QC samples confirmed the repeatability and stability of the instrument and the reliability of the data in the current research (Figure 2(a)). The separation between the AD and sham groups could be clearly observed in the PCA 3D patterns in both negative-and Fold change was calculated as relative intensity obtained from group sham/group AD, and a value less than 1 indicates a decrease in the metabolites of group AD. The levels of potential biomarkers were labeled with "↓" (downregulated) and "↑" (upregulated) ( * p < 0:05). a Change trend compared with the sham group. b Change trend compared with the AD group.
7 Oxidative Medicine and Cellular Longevity positive-ion modes, which demonstrated that the AD model can be successfully induced by Aβ1-42 and there was a severe metabolic disorder in the AD model rats. Furthermore, we also noted a trend towards separation between the BSTSF and AD groups, indicating that metabolism was significantly altered after 28 days of BSTSF administration. A supervised OPLS-DA pattern recognition method was applied to identify the overall metabolic differences between two groups. Fold change was calculated as average relative quantitation obtained from group BSTSF/group AD, and a value less than 1 indicates a decrease in the metabolites of group BSTSF. The levels of potential biomarkers were labeled with "↓" (downregulated) and "↑" (upregulated) ( * p < 0:05). a Change trend compared with the sham group. b Change trend compared with the AD group.

Potential Metabolite Biomarkers in AD and Effect of BSTSF on AD.
Metabolites with VIP scores > 1:0 and p values < 0.05 were defined as potential biomarkers. As shown in Table 1, forty-seven (24 upregulated and 23 downregulated) differential endogenous metabolites were identified between the AD and sham groups. After BSTSF administration, the levels of five metabolites changed significantly, including those of PE(15:0/14:1(9Z)), Cer(d18:0/22:0), fasciculic acid B, citbismine A, and 4-nitrophenol. 14 downregulated metabolites in the AD group compared with the sham group were increased after treatment with BSTSF, and 16 upregulated metabolites were decreased after administration with BSTSF. As shown in Table 2, thirty-four (12 upregulated and 22 downregulated) differential endogenous metabolites were identified between the BSTSF and AD groups. The differential metabolite dataset was imported into R (version 3.4.1) to generate a heatmap. Figures 3(a) and 3(b) show the 47 differential metabolites between the sham and AD groups and the 34 differential metabolites between the AD and BSTSF groups, as well as the relative changes in the concentration of the metabolites in the different groups.

Analysis of Metabolic Pathways.
Metabolic pathway analysis was conducted with MetaboAnalyst 3.5 to further explore the pathogenesis of AD and the possible mechanism by which BSTSF treatment ameliorates this disease. The 47 differential endogenous metabolites between the sham and AD groups were mainly involved in sphingolipid metabolism, glycerophospholipid metabolism, linoleic acid metabolism, and alpha-linolenic acid metabolism ( Table 3). The 34 differential endogenous metabolites between the BSTSF and AD groups were primarily associated with alanine, aspartate, and glutamate metabolism, D-glutamine and D-glutamate metabolism, arginine biosynthesis, and glycosylphosphatidylinositol (GPI) anchor biosynthesis (Table 4). Figure 4 illustrates the differential metabolites classified through the HMDB database. The different colors in each pie chart represent different  Notes: hits represent the matched number of metabolites in one pathway. p represents the original p value calculated from the enrichment analysis. FDR p represents the p value adjusted using the false discovery rate. 10 Oxidative Medicine and Cellular Longevity HMDB classifications, and the area represents the relative proportion of metabolites in the classification. The figure further shows that 60.61% of the differential metabolites between the AD and sham groups were lipid metabolites, indicative of significant changes in the levels of lipid metabolites in the cerebral cortices of AD rats. After BSTSF intervention, 35.71% of the differential metabolites were classified as significantly altered lipids.

HMDB compound classification (sham vs AD)
Lipids and lipid-like molecules: 60.61% Phenylpropanoids and polyketides: 9.09% Organic oxygen compounds: 6.06%  Figure 4: HMDB database classification: (a) differential metabolites between groups sham and AD in HMDB database classification and (b) differential metabolites between groups BSTSF and AD in HMDB database classification.

12
Oxidative Medicine and Cellular Longevity

Lipid Identification and Multivariate Statistical
Analysis. In this study, we identified 1191 lipid species from 31 lipid subclasses. The PCA 3D score plot chart ( Figure 5(a)) shows a trend towards separation between the three groups. The clustered QC samples confirmed the repeatability and stability of the instrument and the reliability of the data. Moreover, OPLS-DA was performed to identify the pattern of separation. As shown in Figure 5(b) (positive-ion mode) and Figure 5(c) (negative-ion mode), there was a clear trend towards separation between the sham and AD groups and between the AD and BSTSF groups. This suggested that lipid metabolism was dysregulated in the cerebral cortices of AD rats, and this dysregulation could be ameliorated by BSTSF administration.

Identification of the Differential Lipid Metabolites in
Cerebral Cortex Samples between the Sham and AD Groups and the BSTSF and AD Groups. VIP scores > 1 and p values < 0.05 were used to determine the significantly altered lipid metabolites between the sham and AD groups and the BSTSF and AD groups. As shown in Table 5, the concentrations of 30 lipids changed markedly after being induced by Aβ1-42 while those of 21 lipids were altered after BSTSF administration. Interestingly, ceramide levels, including those of Cer(d18:1/18:0), Cer(d18:2/20:2), Cer(d36:0), and Cer(d20:1/18:0), and PC(18:0/20:4) and PE(16:0p/22:6) were substantially higher in the AD group than in the sham group, whereas the levels of all these metabolites were significantly reduced in BSTSF-treated rats. PS(20:3/22:6) was downregulated in the AD group compared with the sham group; however, its level was markedly upregulated with BSTSF treatment (Figure 6).

Discussion
In the present study, we found that the cerebral cortex of the AD rats has a distinct metabolomic profile, including the dysregulated sphingolipid metabolism, glycerophospholipid metabolism, linoleic acid metabolism, and alpha-linolenic acid metabolism. The lipidomic analysis indicated that sphingolipids and glycerophospholipids, such as Cer (ceramide), PE (phosphatidylethanolamine), LPE (lysophosphatidylethanolamine), and PC (phosphatidylcholine), are dysregulated in the brain of AD rats. Moreover, the results indicated that BSTSF treatment could restore some dysregulation metabolites and abnormal lipid metabolism in the cerebral cortex of AD rats, especially involving lipid and amino acid metabolism (Figure 7). Lipid homeostasis plays important roles in the central nervous system, including the maintenance of cell membrane structure, signal transduction, and being components of lipid rafts [18][19][20]. Our study identified sphingolipids as a class of lipid metabolites that are closely related to the pathology of AD. As the key intermediates in sphingolipid metabolism, we found that the levels of Cer(d18:0/22:0, d18:1/18:0, d18:2/20:2, d36:0, d20:1/18:0) were markedly upregulated in the AD group when compared with the sham-treated group. Some studies demonstrated that elevated Cer levels in reactive astrocytes were associated with neuroinflammation [21] and can promote the overproduction and aggregation of Aβ through their effects on lipid rafts [22]. Interestingly, the Cer levels were significantly reduced after BSTSF administration and may act as a putative therapeutic target of BSTSF in AD. In addition, sphingomyelins (SMs) were a subclass of sphingolipids, and it also enriched in the central nervous system. They are important constituents of lipid rafts [23] and play a critical role in neuronal cell signaling [24].

13
Oxidative Medicine and Cellular Longevity 14 Oxidative Medicine and Cellular Longevity Importantly, SMs have been associated with amyloidogenic processing of the amyloid precursor protein (APP) [25].
Here, we showed that the levels of SM(d16:1/18:0) were downregulated in AD brains when compared with those from the sham-treated group. Glycerophospholipids were the second major class of lipid metabolites found to be related to AD pathology and BSTSF therapeutic mechanisms in the present study. Concentrations of glycerophospholipids have been associated with the severity of amyloid and neurofibrillary pathology in AD [26]. In this study, the levels of PE  Figure 6: Changes in the relative signal intensities of differential lipids in the rat cerebral cortex from different groups.  Figure 7: Lipid metabolism pathway map depicts the pathomechanism in AD rats and the therapeutic mechanism of BSTSF for AD. Notes: the red and blue letters indicate that the lipids were upregulated or downregulated, respectively, in group AD compared with group sham. The red upward arrow and blue downward arrow indicate that the lipids were upregulated or downregulated, respectively, by BSTSF.
hippocampus of triple transgenic AD mice [27]. Similar to our findings, the level of PC was increased in cerebrospinal fluid from AD patients and was associated with aberrant cerebrospinal fluid Aβ1−42 values, which may be indicative of loss of membrane function and neurodegeneration in the early stages of cognitive dysfunction [14]. Our data confirmed that the concentrations of PC(18:1(11Z)/20:4(5Z,8Z,11Z,14Z)) and PC(18:0/20:4) were increased in the cerebral cortex of the AD rat, which may facilitate the aggregation of Aβ. Treatment with BSTSF may reduce Aβ aggregation through modulation of related PCs. Myelin is enriched in PS in brain tissue, and the docosahexaenoic acid (DHA) content of neuronal PS is of functional importance. Lower PS concentration can reduce the sensitivity of postsynaptic membranes to neurotransmitters such as acetylcholine [28]. PE was also reported to reduce reactive oxygen species levels and protect neuron membranes from oxidative damage [29]. One study found that PS supplementation significantly improved memory, information processing, and the ability to perform daily activities in elderly people [30]. Our data showed that the concentration of PS(20:3/22:6) was markedly upregulated in BSTSFtreated AD rats, which might ameliorate AD pathology by improving sensitivity to neurotransmitters and reducing oxidative damage. The metabolites of sphingolipids and glycerophospholipids both belonged to phospholipid metabolism. So there was a serious disturbance of phospholipid metabolism in brain tissue of AD. Related research confirmed that ApoE proteins are critical determinants of brain phospholipid homeostasis. The phospholipid dysregulation contributes to ApoE4-associated cognitive deficits in AD pathogenesis [31]. APOE4 has been identified as the most prevalent genetic risk factor for AD [32], and it could exacerbate the intraneuronal accumulation of Aβ [33] and plaque deposition in the brain parenchyma [34]. Therefore, BSTSF might treat AD through regulating the APOE4 and phospholipid metabolism.
Additionally, we observed that the level of oleamide was upregulated in AD rats while that of tetracosahexaenoic acid was downregulated after BSTSF administration, and both of these lipids have roles in fatty acid metabolism. Oleamide has been found to accumulate in the cerebrospinal fluid of sleep-deprived cats and mice [35] and is an important regulatory lipid in the brain and central nervous system. Oleamide regulates the sleep-wake cycle, memory, locomotion, and pain perception and exhibits anti-inflammatory, anxiolytic, and neuroprotective properties [36]. Oleamide levels were shown to be increased in the plasma of AD patients [37]. A recent study indicated that DHA is both a product and a precursor of tetracosahexaenoic acid [38], while other investigations have shown that dietary DHA administration can exert protective effects against Aβ production, plaque deposition, and cerebral amyloid angiopathy in an aged mouse model of AD, as well as increasing cerebral blood volume [39,40]. This indicates that BSTSF might regulate DHA and increase tetracosahexaenoic acid treating AD.
In the current study, the level of S-adenosylmethionine (SAM) was significantly upregulated in the cerebral cortex of the AD rat, indicative of aberrant amino acid metabolism in AD. S-Adenosylmethionine is generated through the activity of methionine adenosyltransferases, and more than 50 endogenous substances in the body, including phosphatidylcholine, need SAM as the methyl donor. Homocysteine is formed by SAM demethylation. A related study showed that an elevated total homocysteine baseline serum level during six years was associated with an increased risk of dementia and AD, and the putative mechanism was associated with loss of total brain tissue volume and brain atrophy [41,42]. Therefore, increased SAM concentrations in the AD brain may lead to higher levels of homocysteine, resulting in the loss of total brain tissue volume and brain atrophy. We found that BSTSF ameliorated AD pathology mainly through the alanine, aspartate, and glutamate metabolic pathway and the Dglutamine and D-glutamate metabolic pathway that is involved in regulating amino acid metabolism. N-Acetyl-L-aspartic acid (NAA) and L-aspartic acid are involved in the alanine, aspartate, and glutamate metabolic pathway. NAA is a biomarker for neuronal damage in the human brain during neurodegeneration [43], and reduced cortical levels of NAA were shown to be correlated with clinical scales of dementia severity [44]. L-Aspartic acid also functions as a neurotransmitter. D-Glutamine forms part of the D-glutamate metabolic pathway and is converted from glutamate by astrocytes in the brain [45]. Glutamate, as an excitatory neurotransmitter, is released decarboxylated by neurons to yield the inhibitory neurotransmitter gammaaminobutyric acid (GABA) [46]. Numerous studies have shown that the level of GABA is closely associated with AD [47,48]. D-Glutamine and GABA could be putative therapeutic targets of BSTSF in AD.
Through integrated metabolic and lipidomic analyses, we revealed the metabolic dysregulation in the cerebral cortex of the AD rat, as well as the metabolism-related therapeutic effects of BSTSF on AD. However, this study had some limitations. First, we only obtained a relative quantification of metabolite levels based on untargeted metabolomic analysis. Absolute quantitation of critical metabolites performed by targeted metabolomics is needed in future research. Second, numerous metabolites are regulated by BSTSF; however, the associated molecular mechanisms are not well understood and will be the basis of our future research.

Conclusions
The metabolic disturbance in the cerebral cortex of the AD rat is primarily associated with sphingolipid metabolism, glycerophospholipid metabolism, linoleic acid metabolism, and fatty acid metabolism. Moreover, BSTSF ameliorates the severity of AD by regulating phospholipid metabolism, maintaining fatty acid metabolism, and balancing amino acid metabolism. Our research highlighted some important mechanisms involved in the pathogenesis of AD and revealed the metabolism-related therapeutic effects of BSTSF on AD. 16 Oxidative Medicine and Cellular Longevity