Original Research

Exp. Biol. Med., 29 August 2025

Sec. Genomics, Proteomics and Bioinformatics

Volume 250 - 2025 | https://doi.org/10.3389/ebm.2025.10452

Integrated analysis of metabolome and microbiome in a mouse model of sodium valproate-induced autism

    SZ

    Shuzhen Zhao 1

    XZ

    Xinyan Zhang 2

    YM

    Yanqiu Miao 3

    XG

    Xueya Gao 4

    QW

    Qiuhua Wan 5

    WQ

    Wei Qiu 6

    HS

    Haixia Si 7

    YH

    Yingjie Han 8

    XD

    Xiao Du 1

    YF

    Yuanyuan Feng 1

    LL

    Lianhua Liu 9*

    YC

    Yuqing Chen 10*

  • 1. Child Rehabilitation Center, Jining Maternal and Child Health Family Planning Service Center, Jining, China

  • 2. Department of Nursing, Jining Maternal and Child Health Family Planning Service Center, Jining, China

  • 3. Department of Child Health, Zoucheng People’s Hospital, Zoucheng, China

  • 4. Department of Pediatrics, Jining Maternal and Child Health Family Planning Service Center, Jining, China

  • 5. Clinical Laboratory, Jining Maternal and Child Health Family Planning Service Center, Jining, China

  • 6. Shandong Yiyang Health Group Baodian Mining Hospital, Zoucheng, China

  • 7. Department of Obstetrics and Gynecology, Jining Maternal and Child Health Family Planning Service Center, Jining, China

  • 8. Department of Obstetrics, Jining No. 1 People’s Hospital, Jining, China

  • 9. Department of Medical Administration, Jining Maternal and Child Health Family Planning Service Center, Jining, China

  • 10. Department of Health, Jining Maternal and Child Health Family Planning Service Center, Jining, China

Article metrics

98

Views

17

Downloads

Abstract

Sodium valproate (SV) has been shown to induce autism in animal models. In this study, the SV method was used to establish a mouse model of autism, and anxiety-like behaviours and learning memory performance were evaluated by behavioural tests. The effects of SV on metabolic profiles and gut microbiota were assessed by integrating gas chromatography-mass spectrometry and 16S ribosomal RNA gene sequencing. Correlations between metabolites and gut microbiota were determined using Spearman correlation coefficient. Behavioral tests, including the three-chambered social assay, repetitive behaviors, open field test, elevated plus-maze test, and novel object recognition test, demonstrated that SV treatment exacerbated anxiety-like behaviors and impeded spatial learning and memory in mice. SV disrupted metabolic pathways in hippocampus, cortex, intestine, and serum, affecting primarily valine, leucine and isoleucine biosynthesis, glycerophospholipid metabolism and glutathione metabolism and so on. SV also altered gut microbiota at the genus level, decreasing the abundances of Dubosiella, Faecalibaculum, Clostridia_UCG-014, Bifidobacterium, and Alloprevotella, while increase the abundances of Lactobacillus, Alistipes, and Lachnospiraceae in intestine. The results of correlation analysis showed that in hippocampus, Bifidobacterium was positively correlated with serine and glycine, while Alistipes was negatively correlated with them. These findings suggested that SV may contribute to the development of autism progression by altering the gut microbiota abundances and metabolite profiles. This may provide new direction for the management of autism.

Impact statement

Sodium valproate (SV) has been shown to induce autism in animal models. In this study, we employed a gas chromatography-mass spectrometry (GC-MS)-based metabolomics approach, complemented by 16S ribosomal RNA (rRNA) sequencing, to elucidate potential associations between gut microbiota components and metabolic pathways following exposure to SV. SV disrupted metabolic pathways in hippocampus, cortex, intestine, and serum, affecting primarily valine, leucine and isoleucine biosynthesis, glycerophospholipid metabolism and glutathione metabolism and so on. SV also altered gut microbiota at the genus level, decreasing the abundances of Dubosiella, Faecalibaculum, Clostridia_UCG-014, Bifidobacterium, and Alloprevotella, while increase the abundances of Lactobacillus, Alistipes, and Lachnospiraceae in intestine. The results of correlation analysis showed that in hippocampus, Bifidobacterium was positively correlated with serine and glycine, while Alistipes was negatively correlated with them. These findings suggest that SV may induce neurotoxicity and promote autism progression by altering the gut microbiota abundances and brain metabolite profiles. Our study provides a more comprehensive understanding of the toxic effects induced by SV in a mouse model of autism.

Introduction

Autism spectrum disorder (ASD) is a neurodevelopmental disorder among children [1]. It is statistically that approximately 20 out of every 10,000 children suffer from ASD worldwide, especially for those as young as 1–3 years old [2]. The clinical symptoms of ASD generally include social communication difficulties, repetitive behaviors or even self-injurious behaviors [3]. In addition, patients with ASD may also suffer from epilepsy and mental retardation [4]. Currently, although the detailed causes of ASD are still unknown, multiple factors such as social influences, environmental insults and genetic aberrations can lead to the development of ASD [5, 6]. In some clinical reports, Purkinje cell integrity loss in the cerebellum, hyperserotonemia and oxidative stress damage are considered as the key pathological findings [79]. Additionally, some researches demonstrated that early life environmental factors such as chemical and drug exposure in the mother and prenatal viral infections may also increase the risk of ASD in the offspring [10, 11]. Some chemicals such as misoprostol, thalidomide, mercury and ethanol induce the generation of reactive oxygen species, leading to the development deficits of cerebellum and brain [12].

Sodium valproate (SV) is a well-known anti-epileptic drug since its first introduction in 1978 [13]. In 1995, it was approved by the Food and Drug Administration and became the mood-stabilizing agent of choice for use in stabilization of mania associated with bipolar disorder [14, 15]. Recently, SV has been also employed in migraine prophylaxis and control of neuropathic pain [4]. As the number of indications for SV usage increases, so does the incidence of both accidental and intentional exposures. It is reported that children exposed to SV in utero can cause fetal-valproate syndrome [16]. Fetal-valproate syndrome shows autism like symptoms such as repetitive behaviors, language and communication deficits, hyper excitability and global delays in behavioral development [1719]. Bairy et al. and Hamza et al. uncovered the reproductive toxicity of SV in male rats, manifesting as the decrease of sperm count and sperm motility [20, 21]. Sivathanu et al. conducted a case report that a female infant presented with global developmental delay and infantile spasms [22]. She was then started on SV but developed encephalopathy after an initial improvement [22]. However, the symptoms of vomiting and seizures were reversal on withdrawal of SV [22]. Mei et al. reported a 42-year-old man received antiepileptic treatment with SV following meningioma surgery [23]. He was diagnosed as liver failure and eventually died [23]. SV was also reported to induce hematologic toxicity including bone marrow failure, macrocytosis, thrombocytopenia and neutropenia [24]. Additionally, SV was also used to induce autism in animal models due to its property of affecting brain neurodevelopment and synaptic integration disturbances [25]. These data indicated that SV-induced toxicity represents an increasing concern for toxicologists. Therefore, it is paramount to investigate the mechanisms responsible for SV-associated toxicity.

Humans and animals evolved in intimate association with microbial communities. Gut microbiota plays a crucial role in organ development, metabolism and immune system [2628]. The microbiota-gut-brain-axis describes the physiological connection to exchange information among the microbiota, the gut and the brain [29]. Gut microbiota represents the greatest density and absolute abundance of microorganisms in the human body. A healthy microbial composition is important to health, as dysbiosis is often observed in irritable bowel syndrome, inflammatory bowel disease, obesity, allergy, depression and of course in ASD [3032]. Researches have demonstrated that gastrointestinal (GI) symptoms such as diarrhoea, constipation and vomiting in children with ASD are fourfold than those of healthy population [33, 34]. Additionally, lower diversities of bacterial communities are found in ASD children, which may be related to the severity of GI symptoms [35]. Meanwhile, the abundance of Clostridia species (spp) was relatively higher in autistic individuals, indicating that it is involved in the pathogenesis of ASD [35]. Interestingly, structural changes in the brain and less sociable behaviors are observed in germ-free mice, which suggests that there may be a functional connection between the microbiota and the brain [36, 37].

In this study, we employed an integrated gas chromatography-mass spectrometry (GC-MS)-based metabolomics approach along with 16S ribosomal RNA (rRNA) sequencing to elucidate potential associations between gut microbiota components and metabolic pathways following exposure to SV. Our findings provide a comprehensive understanding of the toxic effects induced by SV in autism mice model. This study provides a research basis for an in-depth investigation of the mechanisms of ASD caused by SV exposure.

Materials and methods

Chemicals and reagents

Sigma Aldrich (St. Louis, MO, USA) provided sodium valproate (SV), N, O-bis(trimethylsilyl)trifluoroacetamide (containing 1% trimethylchlorosilane) and O-methylhydroxylamine hydrochloride (purity ≥98%). Methanol (chromatographic grade), heptadecanoic acid (purity ≥98%) and pyridine were obtained from Macklin Biochemical (Shanghai, China).

Animals and experimental design

The C57BL/6 mice (10 females and 5 males) weighing between 20 and 25 g and aged 8 weeks were procured from Vital River Laboratory Animal Technology located in Pinghu, China. SV-induced autism mouse model was established as previously depicted [38]. In brief, mice were housed with free water and food under a 12 h light/dark cycle at 20–22°C in plastic cages for acclimatization. One week later, mice were kept in separate cages for mating, with one male and two females in each cage. Female mice were placed in a separate plastic cage when pregnant, and divided into two groups randomly. This day was defined as embryonic day 0 (E0). On the day of E13, one group of pregnant mice (n = 5) was administrated with SV (500 mg/kg dissolved in normal saline) via intraperitoneal injection, while the others (n = 5) were given the equivalent volume of saline. The delivery day was considered as postnatal day 0 (P0). Three weeks later (P21), all mice offspring were weaned and sex-grouped, with 5 mice a cage. Pups from SV-injected mothers were considered as the SV group, while those of mothers treated with saline were used as control mice. Each repeat was performed as a separate, independent observation. All experimental procedures conformed to the Guidelines for the Use of Laboratory Animals, and approved by the Ethical Committee for Animal Experimentation of Jining No.1 People’s Hospital (Approval No. JNRM-2023-DW-034).

Behavioral tests

When the pups grew at the age of 8 weeks, they were used for behavioral tests. The social choice test was performed in mice using a three-chamber apparatus comprising a central compartment and two end chambers. A stimulus mouse was randomly introduced into one end chamber (designated as the social chamber), while the opposite chamber served as the non-social chamber [39, 40]. Two identical transparent plexiglass cylinders with perforations for air exchange were positioned in the end chambers. The experimental mouse was placed in the central compartment and permitted to freely explore all chambers for 10 min. Throughout the test session, the duration of time spent in each chamber was recorded, and baseline locomotor activity was quantified.

To detect repetitive behavior, the test mouse is individually placed in a clean environment with bedding, similar to a home cage. The mouse is first allowed to acclimate for 5 min, after which its activity is recorded for 10 min. The duration of repetitive behaviors (such as grooming or digging) is measured using a stopwatch.

For open field test, an open field box (45 × 45 × 30 cm) was used [41]. The 25 × 25 cm area in the center was defined as the central region. Mice were placed into this area and given 10 min for habituation. The distance moved, average speed and moving time of mice were recorded.

For elevated plus-maze test, an apparatus consisted of two closed arms (50 × 10 × 40 cm) and two open arms (50 × 10 cm) [42]. A 25 × 25 cm central platform connected the arms. The apparatus was raised to 50 cm above the floor. Mice were placed on the center platform for 5 min. The number of open/closed arm entries and the corresponding time spent time were recorded. Ethanol (70%) was used to clean the apparatus following each test.

A novel object recognition test is conducted on the mice. On the first day, the test mouse is individually placed in a white, opaque circular chamber for a 30-minute habituation session. On the second day, the mouse is allowed to freely explore two identical, symmetrically positioned objects in the chamber for 10 min. On the third day, one of the two objects is randomly replaced with a novel object that differs in shape and texture, and the test mouse is allowed to explore both objects freely for 10 min. During the test, Ethovision XT 10.1 (Noldus) is used to analyze and record the sniffing time of the mouse toward the novel object (N) and the familiar object (F). The experiment is conducted over three consecutive days.

Sample collection

After behavioral testing, mice were anesthetized with sodium pentobarbital (50 mg/kg), and blood was collected via orbital extraction. The blood samples were centrifuged at 4°C and 4000 rpm for 10 min and then stored at −80°C for further use. Subsequently, the mice were euthanized via cervical dislocation. The intestine, hippocampus, and cerebral cortex were collected, rapidly frozen in liquid nitrogen, and stored at −80°C for further experiments.

Sample preparation

For serum samples, 100 μL of each sample was mixed with 350 μL of heptadecanoic acid (100 μg/mL) and centrifuged at 14,000 rpm for 15 min at 4°C. The supernatant was then dried with liquid nitrogen at 37°C. Subsequently, O-methylhydroxylamine hydrochloride (15 mg/mL) was added, and the mixture was incubated at 70°C for 90 min. Then, 100 μL of N,O-bis(trimethylsilyl)trifluoroacetamide containing 1% trimethyl chlorosilane was added and incubated at 70°C for one hour. For tissue samples, 50 mg of each sample was homogenized in 1 mL of methanol, followed by the addition of 50 μL of heptadecanoic acid (1 mg/mL). The remaining steps were the same as those for the serum samples. A 10 µL aliquot from both the control and SV groups was pooled to serve as a quality control (QC) sample.

GC-MS analysis

GC-MS analysis was performed using a 7890B GC system and 7000C mass spectrometer (Agilent Technologies, USA) equipped with an HP-5MS fused silica capillary column. Helium was used as the carrier gas at a flow rate of 1 mL/min. A sample volume of 1 μL was injected into the GC-MS system with a split ratio of 50:1. The injection temperature, transfer line temperature, and ion source temperature were set to 280°C, 250°C, and 230°C respectively. Electron collision ionization was set to −70 EV with an acquisition frequency of 20 spectra/s. Mass spectrometry employed electrospray ionization in positive mode with a mass/charge (m/z) full scan range from 50 to 800.

Data processing and multivariate analysis

The raw data obtained from gas chromatography-mass spectrometry (GC-MS) analysis were processed using Agilent MassHunter software (version B.07.00). Metabolites in the quality control (QC) samples with a similarity score greater than 80% were identified using the GC-MS library of the National Institute of Standards and Technology (NIST 14). A reference library containing all QC samples was established for spectral matching of experimental sample metabolites. To minimize deconvolution errors during automated data processing and eliminate misidentifications, manual verification was performed.

The resulting comprehensive data matrix included peak indices (RT-m/z pairs), sample names, and corresponding peak areas. Data normalization was performed using the total peak area normalization method. Further data analysis was conducted using orthogonal partial least squares discriminant analysis (OPLS-DA) in SIMCA-P 14.0 software. A two-tailed Student’s t-test was used to assess differences between the two groups. Compounds with a VIP score greater than 1.0 and a p-value less than 0.05 were considered potential differentially expressed metabolites. Heatmap clustering analysis and pathway enrichment were performed using MetaboAnalyst 5.0 software.

16S rRNA sequencing of gut microbiota

Microbial genomic DNA from colonic contents was extracted using the E.Z.N.A.® Soil DNA Kit and separated by 1% agarose gel electrophoresis, followed by PCR amplification of the V3-V4 region of the bacterial 16S rRNA gene. PCR products were purified and sequenced using the TruSeq™ DNA Sample Preparation Kit. Raw data were obtained following the RS_ReadsOfInsert1 protocol. High-quality sequences were processed using the QIIME software package (version 1.9.1). Operational taxonomic units (OTUs) were clustered at 97% similarity using UPARSE (version 7.1). Taxonomic classification analysis was conducted using the Silva 16S rRNA database with the ribosomal database project classifier and the Bayesian algorithm, applying a confidence threshold of 70%.

Results

Autistic-like social and repetitive behavioral deficits in SV-induced mice

In behavioral tests, SV-induced mice showed significantly reduced social interaction on the three-chamber social approach assay, showing diminished preference for investigating a social stimulus as compared to an object (Figures 1A,B, P < 0.0001). When allowed the opportunity for direct interaction with a novel stimulus mouse, SV-induced mice also spent significantly less time initiating contact in comparison to wild-type controls (Figure 1C, P < 0.001). In a home-cage like environment, SV-induced mice showed significantly increased time engaging in repetitive behaviors such as grooming and digging (Figure 1D, P < 0.001). Additionally, SV-induced mice displayed enhanced recognition memory (novel object-recognition test) (Figure 1E, P < 0.001). Open field test demonstrated that the distance moved and average speed of mice in the SV group were significantly reduced compared to those of mice in the control group (Figure 1F, P < 0.05). But the time spent in the open field between the two groups seemed no significant differences. The elevated plus-maze test indicated that the numbers of open arm entries and closed arm entries of mice in the SV group were dramatically reduced relative to the control group (Figure 1G, P < 0.05), and there were no significant differences between the two groups on the time spent in open or closed arms. Therefore, SV-induced mice show strong autistic-like social and repetitive behavioral deficits.

FIGURE 1

Graphs and charts compare control and SV groups across various behavioral tests. Charts A and B display time, distance moved, and number of entries in phases with cylinders. Charts C, D, and E focus on direct interaction, repetitive behavior, and exploration time, respectively. Chart F shows distance, speed, and movement in an open field test with visual trajectories and heat maps. Chart G depicts performance in a plus maze test, highlighting entries and time spent in open and closed arms. Control and SV are denoted in blue and pink, with statistical significance indicated by asterisks.

Autistic-like social and repetitive behavioral deficits in SV-induced mice. (A–C) The three-chambered social assay, (D) repetitive behaviors, (E) novel object recognition test, (F) open field test, and (G) elevated plus-maze test were performed to assess the anxiety like-behaviors and the learning and memory performance of mice. Error bars are Mean ± SEM. p-values from unpaired Student’s t-tests, *P < 0.05, **P < 0.01, ***P < 0.001.

Metabolomic profiles of SV-exposed mice by GC-MS

To determine the effects of SV on the metabolomic profiles of mice among hippocampus, cortex, intestine and serum, a GC-MS-based untargeted metabolomics approach was applied in this study. The representative GC-MS total ion chromatograms (TICs) of QC from hippocampus, cortex, intestine and serum were shown in Supplementary Figures S1A–D. The results indicated that there were significant differences in TICs among different QC samples. As illustrated in Figures 2A–D and Table 1, OPLS-DA models showed that clear differences were observed between the control and SV groups. In different tissues and serum, the proportion of variance explained by the OPLS-DA model (R2X) was 61.6%, 68.2%, 57.2%, and 71.5%, respectively. Additionally, the intersection points between blue regression line (Q2-point) and vertical axis were all negative values, indicating that the OPLS-DA models and the predication were reliable.

FIGURE 2

Four panels labeled A to D compare control and SV groups across different tissues: hippocampus, cortex, intestine, and serum. Each panel includes a scatter plot on the left showing group separation within an ellipse and a bar graph on the right displaying R2 and Q2 values across 200 permutations with one component. Blue and red indicate control and SV groups, respectively. R2 is marked in green and Q2 in blue.

OPLS-DA score plots and 200 permutation tests. (A) Hippocampus, (B) cortex, (C) intestine and (D) serum.

TABLE 1

Tissues R2X (cum) R2Y (cum) Q2 (cum)
Hippocampus 0.616 0.961 0.791
Cortex 0.682 0.985 0.897
Intestine 0.572 0.995 0.861
Serum 0.715 0.939 0.590

OPLS-DA parameter scores.

R2X: the explanation rate of the X matrices; R2Y: the explanation rate of the Y matrices; Q2: the prediction ability.

VIP value >1 and p-value < 0.05 were considered as the important criteria for potential metabolites. As manifested in Figures 3A–D, the cluster analysis of differentially-expressed metabolites of the control group and SV group were depicted. A total of 11 differential metabolites were identified in hippocampal tissues, including 1-monopalmitin (MG (16:0/0:0/0:0)), formamide, o-phosphoethanolamine, glycine, ethanolamine, gamma-aminobutyric acid, serine, L-threonine, pipecolic acid, glycerol, and myo-inositol, all of which were downregulated following SV treatment. In cortex, 7 downregulated metabolites including ethanolamine, glycine, 2-Aminobenzoic acid, serine, L-glutamic acid, o-phosphoethanolamine and adenosine were identified. Moreover, we observed SV treatment altered a total of 8 metabolites including creatinine, myristic acid, desmosterol, stearic acid, palmitelaidic acid, palmitic acid, pyroglutamic acid and scyllo-inositol in intestine tissues. Among these differentially-expressed metabolites, only creatinine was found to be downregulated after SV treatment. In the serum, 11 upregulated metabolites including pyroglutamic acid, palmitic acid, urea, citric acid, glycine, myo-inositol, L-valine, L-proline, ornithine, L-alanine and serine were identified following SV treatment. Next, the differentially-expressed metabolites were subjected to pathway analysis via Metaboanalyst 6.0 and KEGG database. As shown in Figure 4A and Table 2, we found that SV treatment mainly affected galactose metabolism, glycine, serine and threonine metabolism, and glycerophospholipid metabolism in hippocampus. In cortex, pathways of glutathione metabolism, porphyrin metabolism, and glyoxylate and dicarboxylate metabolism were altered following SV treatment. In intestine, only biosynthesis of unsaturated fatty acids and fatty acid biosynthesis were affected by SV treatment. Notably, SV treatment affected glutathione metabolism, arginine biosynthesis, and alanine, aspartate and glutamate metabolism in serum samples. The detailed metabolic network illustrated in Figure 4B highlights the significant impact of SV treatment. Our analysis revealed that SV-exposed disrupted not only various amino acid metabolic pathways but also glutathione metabolism, which plays a critical role in cellular antioxidant defense, and galactose metabolism, an essential pathway for carbohydrate processing.

FIGURE 3

Four heatmaps (A-D) showing metabolite levels in different tissues. (A) Hippocampus, (B) Cortex, (C) Intestine, (D) Serum. Each map includes a color gradient indicating Z-scores from -2 (blue) to 2 (red) for control and SV groups. Metabolites are listed alongside each heatmap. Hierarchical clustering is shown on the top and side of each map.

Heatmap of differential metabolites in (A) hippocampus, (B) cortex, (C) intestine and (D) serum in the SV groups compared with controls. “Control” represents untreated control mice and “SV” represents mice prenatally exposed to sodium valproate. Rows represent samples and columns represent metabolites. The color legend is located at the top right of the figure, where blue represents the control group and red represents the SV-treated group. Prior to unsupervised hierarchical clustering of the samples (rows), the relative abundance of the compounds was z-score normalized, thus making the down-regulated relative expression abundance appear in blue (−2 to 0) and the up-regulated relative expression abundance appear in red (0–2). Darker colors indicate more upregulated metabolite expression.

FIGURE 4

Diagram illustrating metabolism pathways related to various tissues. Part A displays a bubble plot showing different KEGG pathways for the hippocampus, cerebral cortex, intestines, and serum, with bubble size and color indicating significance. Part B presents a network map outlining interactions between metabolites and pathways, such as alanine, aspartate, and glutamate metabolism, alongside others.

Summary of pathway analysis performed using MetaboAnalyst 6.0. (A) KEGG pathway enrichment analysis of differentially metabolites. (B) Schematic diagram of related metabolic pathways affected by SV in serum and major tissues.

TABLE 2

Samples Pathway names Raw p -log10 p value
Hippocampus Galactose metabolism 7.4E-3 2.128
Glycine, serine and threonine metabolism 0.011 1.958
Glycerophospholipid metabolism 0.013 1.884
Valine, leucine and isoleucine biosynthesis 0.040 1.398
Ascorbate and aldarate metabolism 0.045 1.348
Cortex Glutathione metabolism 4.4E-3 2.359
Porphyrin metabolism 5.4E-3 2.271
Glyoxylate and dicarboxylate metabolism 5.7E-3 2.244
Glycerophospholipid metabolism 7.2E-3 2.143
Nitrogen metabolism 0.023 1.644
Biosynthesis of unsaturated fatty acids 4.9E-3 2.313
Fatty acid biosynthesis 8.2E-3 2.084
Intestine Glutathione metabolism 5.6E-4 3.255
Arginine biosynthesis 3.2E-3 2.499
Serum Alanine, aspartate and glutamate metabolism 0.013 1.901
Glyoxylate and dicarboxylate metabolism 0.016 1.789
Arginine and proline metabolism 0.020 1.691
Valine, leucine and isoleucine biosynthesis 0.0497 1.303
Galactose metabolism 7.4E-3 2.128
Glycine, serine and threonine metabolism 0.011 1.958

Pathway analysis performed using MetaboAnalyst 5.0 software.

The changes of gut microbiota induced by SV treatment

Gut microbiota at OTU level in SV-treated mice were determined. As shown in Figure 5A, a total of 232 OTUs were shared in all the samples, and the results showed a species abundance of the control group > the SV group. Supplementary Figures S2A,B shows that community richness and diversity were altered after the SV-exposed (alpha-diversity indices Chao and Ace). The coverage index representing the sample coverage detected indicated that the sequencing depth was representative of the gut microbiome (Supplementary Figure S2C). Beta-diversity was analyzed using the principal coordinate analysis (PCoA) on OTU level showed a significant separation between the control group and the SV group (Figure 5B). The abscissa represents the first principal component and the ordinate is the second principal component. Furthermore, 66% is the contribution value of the first principal component to the sample difference, while 22.62% is the contribution value of the second principal component to the sample difference. The analyses of bacterial communities showed that at the phylum level, Bacteroidota, Firmicutes, Actinobacteriota were the major components of gut microbiota in the control group, while the abundances of Patescibacteria and Verrucomicrobiota were increased but Proteobacteria and Desulfobacterota were decreased following SV treatment (Figures 5C,D). At the genus level, the gut microbiota in the control group was mainly composed of Muribaculaceae, Lactobacillus, and Alistipes. After induced by SV treatment, the abundances of Dubosiella, Faecalibaculum, and Clostridia_UCG-014 were decreased, while the abundances of Lactobacillus, Alistipes, and Lachnoclostridium were increased (Figures 5E,F).

FIGURE 5

Panel A displays a Venn diagram comparing the Control and SV groups, indicating shared and unique features. Panel B shows a principal component analysis plot with distinct separation between the Control and SV groups. Panel C presents bar graphs of community abundance at the phylum level across different samples. Panel D includes a heatmap illustrating the abundance of various phyla in Control and SV. Panel E offers bar graphs of community abundance at the genus level. Panel F displays a heatmap for genus-level abundance comparisons between both groups.

The changes of gut microbiota induced by SV treatment. (A) OTU distribution Venn diagram. (B) PCoA on OUT level between SV and control groups. (C) Histogram of gut microflora at the phylum level. (D) Heatmap of gut microflora at the phylum level. (E) Histogram of gut microflora at the genus level. (F) Heatmap of gut microflora at the genus level.

Relevance analysis between metabolites and gut microbiota

Spearman’s correlation coefficient was used to assess the relationships between metabolites and gut microbiota at both the phylum and genus levels. At the phylum level, Firmicutes showed negative correlations with pipecolic acid, glycine, and ethanolamine in the hippocampus (Supplementary Figure S3A). Patescibacteria exhibited negative correlations with serine, O-phosphoethanolamine, L-threonine, gamma-aminobutyric acid, and formamide. Proteobacteria exhibited positive correlations with MG (16:0/0:0/0:0) and glycerol. In the cortex (Supplementary Figure S3B), Campilobacterota, Cyanobacteria, Proteobacteria, and unclassified_Bacteria were positively correlated with adenosine. Proteobacteria also showed positive correlations with L-glutamic acid, ethanolamine, and 2-aminobenzoic acid. In the intestine (Supplementary Figure S3C), Cyanobacteria exhibited negative correlations with scyllo-inositol, palmitic acid, palmitelaidic acid, and myristic acid, while Patescibacteria showed positive correlations with these metabolites. Bacteroidota exhibited negative correlations with palmitelaidic acid and myristic acid. In the serum (Supplementary Figure S3D), Patescibacteria showed positive correlations with several metabolites, including citric acid, glycine, and L-alanine. Cyanobacteria and unclassified_Bacteria exhibited negative correlations with serine, pyroglutamic acid, and L-alanine.

At the genus level, several bacterial genera, including Alistipes, Alloprevotella, Bifidobacterium, and Candidatus Saccharimonas, showed strong correlations with various metabolites in the hippocampus (Figure 6A). Specifically, Alloprevotella, Bifidobacterium, and Muribaculum exhibited positive correlations with metabolites such as serine, o-phosphoethanolamine, L-threonine, glycine, gamma-aminobutyric acid, formamide, and ethanolamine. In contrast, Alistipes, Candidatus Saccharimonas, and Lachnospiraceae NK4A136 group were negatively correlated with these same metabolites. In the cortex (Figure 6B), Alloprevotella showed positive correlations with L-glutamic acid, glycine, adenosine, and 2-aminobenzoic acid, while Enterorhabdus exhibited negative correlations with these metabolites. Additionally, Lactobacillus was negatively correlated with adenosine, while the Prevotellaceae NK3B31 group showed a positive correlation. In the intestine (Figure 6C), Alloprevotella, Bifidobacterium, and Muribaculum were negatively correlated with palmitic acid and myristic acid, whereas Candidatus Saccharimonas was positively correlated with these metabolites. In the serum (Figure 6D), Alistipes, Candidatus Saccharimonas, and Lachnoclostridium were found to have positive correlations with several key metabolites, including glycine, L-proline, L-valine, and myo-inositol. Similarly, Bifidobacterium, Dubosiella, Faecalibaculum, and Muribaculum also exhibited positive correlations with these metabolites. Additionally, Bifidobacterium, Dubosiella, Faecalibaculum, and Muribaculum were negatively correlated with urea.

FIGURE 6

Heatmaps show correlations of metabolomics and 16S data in four regions: A) Hippocampus, B) Cortex, C) Intestine, D) Serum. Data points represent correlation strength, using color to indicate positive (red) or negative (blue) significance at various p-values, with circle sizes reflecting correlation magnitude.

Spearman correlation analysis of metabolites and gut microbiota at the genus level. (A) Hippocampus, (B) cortex, (C) intestine, and (D) serum. Red: positive correlations, blue: negative correlations.

Discussion

SV is a widely used anticonvulsant drug that modulates γ-aminobutyric acid (GABA) [4346]. Studies suggest that SV likely improves psychotic symptoms by inhibiting presynaptic GABA transaminase (the enzyme responsible for catalyzing the breakdown of GABA into succinic semialdehyde), stimulating glutamate decarboxylase (the enzyme that catalyzes the synthesis of GABA from glutamate), or indirectly increasing presynaptic GABA levels through negative feedback, thereby inducing its release [47, 48]. SV may also regulate the activity of dopamine, serotonin, and glutamate, thus playing a role in improving symptoms of schizophrenia and behavioral disorders. Although SV was historically regarded as a drug with minimal side effects and was extensively used in clinical practice [4951], however, it has been shown that SV can induce hepatotoxicity, hematotoxicity and neurotoxicity [23, 24, 52]. For example, prenatal exposure to SV increases the risk of a range of fetal disorders, including anencephaly, developmental delay, cognitive impairment, autism, and neural tube defects such as spina bifida [5355]. Studies in rodent models have demonstrated that prenatal SV exposure can result in core symptoms resembling those observed in individuals with autism spectrum disorder (ASD) [56, 57]. Additionally, several studies have reported cognitive impairments in SV-induced autism-like models [5659]. However, studies investigating the specific mechanisms underlying SV-induced autism remain limited. For example, Rahel Feleke et al. [60] found that a group of genes downregulated by valproic acid (VPA) were significantly enriched in pathways related to neurodevelopment and synaptic function, as well as in genetic traits associated with human intelligence, schizophrenia, and bipolar disorder. Park et al. [61] reported that in a VPA-induced ASD mouse model, upregulation of Rnf146 led to dysregulation of the Wnt/β-catenin signaling pathway, resulting in impaired social behaviors in mice. Hyo Sang Go et al. [62] demonstrated that VPA could activate the NF-κB signaling pathway and upregulate Bcl-XL, thereby inhibiting normal apoptosis of neural progenitor cells, which may contribute to the neurodevelopmental defects observed in fetal valproate syndrome. In our study, we used a mouse model of autism induced by SV exposure to further investigate the causes of SV-induced autism pathogenesis through untargeted metabolomics combined with gut microbial 16S RNA sequencing.

As researchers continue to focus on autism, the pathogenesis of autism is becoming clearer. Some studies have shown that autism is associated with abnormalities in the body’s metabolism. A report conducted by Tărlungeanu et al. indicated that disrupted amino acid metabolism is a crucial cause of ASD [63]. A recent study has indicated an association between amino acid synthesis/metabolism and attention deficit/hyperactivity disorder or depression [64]. In a study linking polymorphisms for the risk of autism to protein interaction networks in cortex, amino acid synthesis/metabolism ranks the highest score among the biological pathways [65]. Similarly, Brister et al. also uncovered that the dysfunction of amino acid synthesis/metabolism is observed in participants with ASD, which is related to neurodevelopment and affects ASD-related symptoms [66]. We hypothesized that amino acid synthesis/metabolism might be affected by SV treatment. We found that glycine and serine were reduced in the hippocampus of SV-treated mice. Glycine is a well-known non-essential amino acid that possesses anti-inflammatory capability. Research has shown that as the inflammatory status is an important pathological feature in ASD progression, supplementation with glycine in diet may have therapeutic potential for patients with ASD [67]. Meanwhile, similar to gamma-Aminobutyric acid, glycine acts as an excitatory neurotransmitter to depolarize membrane potentials during the early development of central nervous system [68]. It must shift from excitatory to inhibitory neurotransmitter at birth and maturation, otherwise it may lead to the occurrence of ASD [69]. Serine, derived from glycine, is an effective coagonist of N-methyl-D-aspartate receptor in brain areas and its absence has been reported to be associated with the pathogenesis of neurological disorders [70, 71]. It has been demonstrated that the glycine/serine sites on the N-methyl-D-aspartate receptors can serve as a target for the treatment of ASD [72]. These public data uncovered the great importance of glycine and serine metabolism. These results implied that the decrease of glycine and serine caused by SV-exposed may be related to the development of ASD in mice.

Furthermore, our study suggests that disorders of lipid metabolism occur in autism. Lipids in the brain may influence emotional and perception behaviors, leading to depression and anxiety disorders [73]. The central nervous system (CNS) mainly consists of phospholipids (e.g., phosphatidylcholine), sphingolipids (e.g., ceramides), and sterols (e.g., cholesterol), while the content of neutral glycerides is relatively low under physiological conditions. However, some studies have indicated that certain neutral lipids may also be involved in the regulation of neurological diseases [74]. Neutral lipids can be incorporated into lipid droplets. Lipid droplets establish functional contact sites with other organelles such as mitochondria and the endoplasmic reticulum (ER), playing crucial roles in intracellular signaling and metabolism to promote metabolic cross-talk and regulation. Furthermore, although lipid droplets are typically scarce in a healthy CNS, substantial accumulation of lipid droplets has been observed predominantly in glial cells in the aging brain and in neurodegenerative diseases such as multiple sclerosis, Alzheimer’s disease, Parkinson’s disease, and Huntington’s disease, due to imbalances in lipid uptake, synthesis, and mobilization [74]. This implies that disorders of lipid metabolism might be associated with the development of ASD. In our study, we performed KEGG pathway enrichment analysis of differential metabolites in different tissues of ASD mice, and the results showed that autism caused by SV exposure was associated with the biosynthesis of multiple fatty acids and the glycerophospholipid metabolic pathway. These results suggest that we believe that SV-induced disorders of lipid-phospholipid metabolism may be one of the major causes of ASD development.

The Firmicutes phylum and Bacteroidota phylum are the major components of the gut microbiota. We found that both in the control and SV groups, Firmicutes phylum and Bacteroidota phylum were indeed the major components of gut microbiota. These indicated that SV had few effects on these two microbial communities. Notably, a recent study demonstrated that the abundances of Cyanobacteria were increased with age in typically developing individuals compared to the patients with ASD [75]. Interestingly, in patients with late-life depression, higher abundances of Patescibacteria and Verrucomicrobiota were observed [76]. We therefore speculated that there may be similar alterations in SV-induced autism mice. In this study, we found that the abundances of Patescibacteria and Verrucomicrobiota were increased but Cyanobacteria and unclassified_Bacteria were decreased following SV treatment. There results indicated that SV treatment could result in gut microbiota components changes at the phylum level. At the genus level, the mice treatment with SV showed more altered genera, including the decrease abundances of Dubosiella, Faecalibaculum, Clostridia_UCG-014, and Bifidobacterium and the increase abundances of Lactobacillus, Alistipes, Lachnospiraceae, and Lachnospiraceae_NK4A136. Our results lend credence to some published findings suggesting that the decrease abundances of the Dubosiella and Bifidobacterium were associated with the development of ASD [75, 77]. Several microorganisms can mediate gut–brain signaling through inducing host production of neurotransmitters and metabolites, and generating some neuroactive compounds themselves [78]. For instance, Bifidobacterium and Parabacteroides can produce gamma-Aminobutyric acid, a key neurotransmitter in the brain system [79]. We therefore believed that SV treatment may mediate microbiota-gut-brain-axis to induce neurotoxicity in autism mice through decreasing the abundance of Bifidobacterium in intestine tissues.

Since microorganisms can affect host production of metabolites, we then explored the relationship between gut microbes and metabolic changes in mice with ASD. Currently, numerous studies have demonstrated correlations between changes in the gut microbiome and metabolome. For instance, Wu et al. [80] reported that blood metabolites associated with impaired glycemic control were linked to alterations in the gut microbiota. Fu et al. [81] found that the gut microbiota could influence the development of inflammatory bowel disease (IBD) by modulating bile acid metabolism. Liu et al. [82] also demonstrated a relationship between the gut microbiota and lipid metabolism. Moreover, several studies have shown that alterations in the gut microbiota can impact brain metabolism and subsequently regulate the development of various diseases. For example, Wang et al. [83] found that the gut microbiota regulated insomnia-like behaviors via the gut-brain metabolic axis. Similarly, Xiao et al. [84] demonstrated that fecal microbiota transplantation (FMT) could ameliorate gut dysbiosis, cognitive decline, and depression-like behaviors induced by bilateral common carotid artery occlusion (BCCAO), possibly by increasing the relative abundance of short-chain fatty acid (SCFA)-producing bacteria and enhancing SCFA levels, thereby alleviating chronic cerebral hypoperfusion (CCH) injury.

The results of the analysis show that in hippocampus, Bifidobacterium was positively correlated with serine and glycine, while Alistipes was negatively correlated with them. Meanwhile, in cortex, similar patterns were observed in spite of no significant differences. These results implied that SV challenge may affect glycine, serine and threonine metabolism and further influence the abundances of Bifidobacterium and Alistipes via gut-brain-axis, eventually leading to the development of ASD. Additionally, gut microbiota may have direct or indirect effects on drug metabolism, such as providing a series of additional reactions and modulating drug metabolism in host [85]. Although it has been confirmed that there are some interactions between the microbiota-gut-brain-axis and brain disorders including ASD [86, 87], the potential association between gut microbiome changes and the brain metabolic alterations influenced by SV remains unclear.

Although this study, through multi-omics integration, revealed associations between gut microbiota-metabolite alterations and autism-like behaviors induced by SV exposure, certain limitations remain. While our findings suggest that SV exposure leads to dysbiosis of gut microbiota in mice, further studies utilizing fecal microbiota transplantation or germ-free animal models are necessary to validate the underlying mechanisms. Furthermore, KEGG pathway analysis indicated a potential association between lipid metabolism and the pathogenesis of autism. In future studies, we plan to quantify free fatty acids in intestinal tissues and employ techniques such as Oil Red O staining, intestinal electron microscopy, and targeted lipidomic profiling to investigate impairments in lipid absorption and storage induced by SV exposure. In addition, we intend to extend our research to clinical settings by comparing metabolomic profiles of autism patients with those of SV-induced autism mouse models. By analyzing the shared and distinct metabolic alterations between human patients and animal models, we aim to achieve a more comprehensive understanding of the metabolic disturbances and underlying mechanisms involved in autism spectrum disorder.

Finally, although this study highlights the potential involvement of gut microbiota-metabolite associations in SV-induced autism, it only provides preliminary insights into possible mechanisms underlying SV-induced autism, and does not constitute clinical validation. Therefore, these preclinical findings should not be used to encourage dietary or probiotic interventions without medical supervision.

Conclusion

SV exposure led to toxicity associated with gut microbiota and metabolomic pathways changes in autism mice. SV treatment mainly disrupted lipid metabolism and amino acid synthesis/metabolism in hippocampus and cortex, and more importantly, disturbed glycine, serine and threonine metabolism in hippocampus. Additionally, SV decreased the abundances of Bifidobacterium, and increased Alistipes abundance. These may be an important regulatory mechanism for SV-caused ASD.

Statements

Data Availability Statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material. The metabolomics data presented in the study are deposited in the metabolights repository, accession number MTBLS10013. The 16S rRNA data presented in the study are deposited in the NCBI BioProject repository, accession number PRJNA1103334.

Ethics Statement

The animal study was approved by Ethical Committee for Animal Experimentation of Jining No. 1 People's Hospital. The study was conducted in accordance with the local legislation and institutional requirements. Ethics approval number: JNRM-2023-DW-034.

Author Contributions

Conceptualization, LL and YC; Methodology, YM, XG, QW, WQ, HS, YH, XD, YF; Software, YM, XG, QW, WQ, HS, YH, XD, YF; Validation, LL and YC; Investigation, SZ and XZ; Writing-Original Draft Preparation, SZ and XZ; Writing-Review and Editing, all authors; Visualization, SZ and XZ; Supervision, LL and YC. All authors contributed to the article and approved the submitted version.

Funding

The author(s) declare that financial support was received for the research and/or publication of this article. The experiments in this study were completed with the support of internal institutional funding.

Acknowledgments

We acknowledge the strong support of Jining Maternal and Child Health Family Planning Service Center, Zoucheng People's Hospital, Shandong Yiyang Health Group Baodian Mining Hospital and Jining No.1 People's Hospital.

Conflict of Interest

The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

Generative AI statement

The author(s) declare that no Generative AI was used in the creation of this manuscript.

Supplementary Material

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

SUPPLEMENTARY FIGURE 1

Representative GC–MS TICs of QC from (A) hippocampus, (B) cortex, (C) intestine and (D) serum.

SUPPLEMENTARY FIGURE 2

Alpha-diversity was investigated using the (A) Ace, (B) Chao, and (C) coverage indexes.

SUPPLEMENTARY FIGURE 3

Spearman correlation analysis of metabolites and gut microbiota at the phylum level. (A) Hippocampus, (B) cortex, (C) intestine, and (D) serum. Red: positive correlations, blue: negative correlations.

References

  • 1

    Hyman SL Levy SE Myers SM Council On Children With Disabilities SOD Behavioral P . Identification, evaluation, and management of children with autism spectrum disorder. Pediatrics (2020) 145. 10.1542/peds.2019-3447

  • 2

    Fakhoury M . Autistic spectrum disorders: a review of clinical features, theories and diagnosis. Int J Developmental Neurosci (2015) 43:707. 10.1016/j.ijdevneu.2015.04.003

  • 3

    Nazeer A Ghaziuddin M . Autism spectrum disorders: clinical features and diagnosis. Pediatr Clin North America (2012) 59:1925. ix. 10.1016/j.pcl.2011.10.007

  • 4

    Arafat EA Shabaan DA . The possible neuroprotective role of grape seed extract on the histopathological changes of the cerebellar cortex of rats prenatally exposed to valproic acid: animal model of autism. Acta Histochem (2019) 121:84151. 10.1016/j.acthis.2019.08.002

  • 5

    Bernard S Enayati A Redwood L Roger H Binstock T . Autism: a novel form of mercury poisoning. Med Hypotheses (2001) 56:46271. 10.1054/mehy.2000.1281

  • 6

    Chaste P Leboyer M . Autism risk factors: genes, environment, and gene-environment interactions. Dialogues Clin Neurosci (2012) 14:28192. 10.31887/dcns.2012.14.3/pchaste

  • 7

    Kelly E Meng F Fujita H Morgado F Kazemi Y Rice LC et al Regulation of autism-relevant behaviors by cerebellar-prefrontal cortical circuits. Nat Neurosci (2020) 23:110210. 10.1038/s41593-020-0665-z

  • 8

    Hranilović D Bujas-Petković Z Tomičić M Bordukalo-Nikšić T Blažević S Čičin-Šain L . Hyperserotonemia in autism: activity of 5HT-associated platelet proteins. J Neural Transm (Vienna) (2009) 116:493501. 10.1007/s00702-009-0192-2

  • 9

    Kern JK Jones AM . Evidence of toxicity, oxidative stress, and neuronal insult in autism. J Toxicol Environ Health B (2006) 9:48599. 10.1080/10937400600882079

  • 10

    Dai YC Zhang HF Schon M Bockers TM Han SP Han JS et al Neonatal oxytocin treatment ameliorates autistic-like behaviors and oxytocin deficiency in valproic acid-induced rat model of autism. Front Cell Neurosci (2018) 12:355. 10.3389/fncel.2018.00355

  • 11

    Pham C Symeonides C O’Hely M Sly PD Knibbs LD Thomson S et al Early life environmental factors associated with autism spectrum disorder symptoms in children at age 2 years: a birth cohort study. Autism (2022) 26:186481. 10.1177/13623613211068223

  • 12

    Pragnya B Kameshwari JS Veeresh B . Ameliorating effect of piperine on behavioral abnormalities and oxidative markers in sodium valproate induced autism in BALB/C mice. Behav Brain Res (2014) 270:8694. 10.1016/j.bbr.2014.04.045

  • 13

    Pinkston R Walker LA . Multiorgan system failure caused by valproic acid toxicity. The Am J Emerg Med (1997) 15:5046. 10.1016/s0735-6757(97)90195-9

  • 14

    Davis LL Ryan W Adinoff B Petty F . Comprehensive review of the psychiatric uses of valproate. J Clin Psychopharmacol (2000) 20:1S17S. 10.1097/00004714-200002001-00001

  • 15

    Nemeroff CB . An ever-increasing pharmacopoeia for the management of patients with bipolar disorder. J Clin Psychiatry (2000) 61 Supp 13(Suppl. 13):1925.

  • 16

    Santos de Oliveira R Lajeunie E Arnaud E Renier D . Fetal exposure to sodium valproate associated with baller-gerold syndrome: case report and review of the literature. Childs Nerv Syst (2006) 22:904. 10.1007/s00381-004-1089-x

  • 17

    Markram K Rinaldi T Mendola DL Sandi C Markram H . Abnormal fear conditioning and amygdala processing in an animal model of autism. Neuropsychopharmacology (2008) 33:90112. 10.1038/sj.npp.1301453

  • 18

    Schneider T Przewłocki R . Behavioral alterations in rats prenatally exposed to valproic acid: animal model of autism. Neuropsychopharmacology (2005) 30:809. 10.1038/sj.npp.1300518

  • 19

    Schneider T Labuz D Przewłocki R . Nociceptive changes in rats after prenatal exposure to valproic acid. Pol J Pharmacol (2001) 53:5314.

  • 20

    Bairy L Paul V Rao Y . Reproductive toxicity of sodium valproate in male rats. Indian J Pharmacol (2010) 42:904. 10.4103/0253-7613.64503

  • 21

    Hamza AA Amin A . Apium graveolens modulates sodium valproate-induced reproductive toxicity in rats. J Exp Zoolog A: Ecol Genet Physiol (2007) 307A:199206. 10.1002/jez.357

  • 22

    Sivathanu S Sampath S Veerasamy M Sunderkumar S . Encephalopathy in an infant with infantile spasms: possible role of valproate toxicity. BMJ Case Rep (2014) 2014:bcr2013200895. 10.1136/bcr-2013-200895

  • 23

    Mei X Wu HC Ruan M Cai LR . Acute liver failure with thrombotic microangiopathy due to sodium valproate toxicity: a case report. World J Clin Cases (2021) 9:43107. 10.12998/wjcc.v9.i17.4310

  • 24

    Acharya S Bussel JB . Hematologic toxicity of sodium valproate. J Pediatr Hematology/Oncology (2000) 22:625. 10.1097/00043426-200001000-00012

  • 25

    Rinaldi T Kulangara K Antoniello K Markram H . Elevated NMDA receptor levels and enhanced postsynaptic long-term potentiation induced by prenatal exposure to valproic acid. Proc Natl Acad Sci U S A (2007) 104:135016. 10.1073/pnas.0704391104

  • 26

    Collins J Borojevic R Verdu EF Huizinga JD Ratcliffe EM . Intestinal microbiota influence the early postnatal development of the enteric nervous system. Neurogastroenterology & Motil (2014) 26:98107. 10.1111/nmo.12236

  • 27

    Dabke K Hendrick G Devkota S . The gut microbiome and metabolic syndrome. J Clin Invest (2019) 129:40507. 10.1172/jci129194

  • 28

    Zheng D Liwinski T Elinav E . Interaction between microbiota and immunity in health and disease. Cell Res (2020) 30:492506. 10.1038/s41422-020-0332-7

  • 29

    Santocchi E Guiducci L Fulceri F Billeci L Buzzigoli E Apicella F et al Gut to brain interaction in autism spectrum disorders: a randomized controlled trial on the role of probiotics on clinical, biochemical and neurophysiological parameters. BMC Psychiatry (2016) 16:183. 10.1186/s12888-016-0887-5

  • 30

    Dinan TG Cryan JF . The microbiome-gut-brain axis in health and disease. Gastroenterol Clin North America (2017) 46:7789. 10.1016/j.gtc.2016.09.007

  • 31

    Muszer M Noszczyńska M Kasperkiewicz K Skurnik M . Human microbiome: when a friend becomes an enemy. Arch Immunol Ther Exp (Warsz) (2015) 63:28798. 10.1007/s00005-015-0332-3

  • 32

    Petra AI Panagiotidou S Hatziagelaki E Stewart JM Conti P Theoharides TC . Gut-microbiota-brain axis and its effect on neuropsychiatric disorders with suspected immune dysregulation. Clin Ther (2015) 37:98495. 10.1016/j.clinthera.2015.04.002

  • 33

    Fulceri F Morelli M Santocchi E Cena H Del Bianco T Narzisi A et al Gastrointestinal symptoms and behavioral problems in preschoolers with autism spectrum disorder. Dig Liver Dis (2016) 48:24854. 10.1016/j.dld.2015.11.026

  • 34

    Iovene MR Bombace F Maresca R Sapone A Iardino P Picardi A et al Intestinal dysbiosis and yeast isolation in stool of subjects with autism spectrum disorders. Mycopathologia (2017) 182:34963. 10.1007/s11046-016-0068-6

  • 35

    Kang DW Park JG Ilhan ZE Wallstrom G Labaer J Adams JB et al Reduced incidence of prevotella and other fermenters in intestinal microflora of autistic children. PLoS One (2013) 8:e68322. 10.1371/journal.pone.0068322

  • 36

    Dinan TG Cryan JF . Gut instincts: microbiota as a key regulator of brain development, ageing and neurodegeneration. The J Physiol (2017) 595:489503. 10.1113/jp273106

  • 37

    Vuong HE Yano JM Fung TC Hsiao EY . The microbiome and host behavior. Annu Rev Neurosci (2017) 40:2149. 10.1146/annurev-neuro-072116-031347

  • 38

    Eissa N Azimullah S Jayaprakash P Jayaraj RL Reiner D Ojha SK et al The dual-active histamine H3 receptor antagonist and acetylcholine esterase inhibitor E100 ameliorates stereotyped repetitive behavior and neuroinflammmation in sodium valproate induced autism in mice. Chemico-Biological Interactions (2019) 312:108775. 10.1016/j.cbi.2019.108775

  • 39

    Fairless AH Katz JM Vijayvargiya N Dow HC Kreibich AS Berrettini WH et al Development of home cage social behaviors in BALB/cJ vs. C57BL/6J mice. Behav Brain Res (2013) 237:33847. 10.1016/j.bbr.2012.08.051

  • 40

    Sankoorikal GM Kaercher KA Boon CJ Lee JK Brodkin ES . A mouse model system for genetic analysis of sociability: C57BL/6J versus BALB/cJ inbred mouse strains. Biol Psychiatry (2006) 59:41523. 10.1016/j.biopsych.2005.07.026

  • 41

    Lucchina L Depino AM . Altered peripheral and central inflammatory responses in a mouse model of autism. Autism Res (2014) 7:27389. 10.1002/aur.1338

  • 42

    Fatemi I Delrobaee F Bahmani M Shamsizadeh A Allahtavakoli M . The effect of the anti-diabetic drug metformin on behavioral manifestations associated with ovariectomy in mice. Neurosci Lett (2019) 690:958. 10.1016/j.neulet.2018.10.024

  • 43

    Fukuchi M Nii T Ishimaru N Minamino A Hara D Takasaki I et al Valproic acid induces up- or down-regulation of gene expression responsible for the neuronal excitation and inhibition in rat cortical neurons through its epigenetic actions. Neurosci Res (2009) 65:3543. 10.1016/j.neures.2009.05.002

  • 44

    Janszky J Tényi D Bóné B . Valproate in the treatment of epilepsy and status epilepticus. Ideggyogy Sz (2017) 70:25864. 10.18071/isz.70.0258

  • 45

    Gagnon DJ Fontaine GV Riker RR Fraser GL . Repurposing valproate, enteral clonidine, and phenobarbital for comfort in adult ICU patients: a literature review with practical considerations. Pharmacotherapy (2017) 37:130921. 10.1002/phar.2017

  • 46

    Quinn NJ Hohlfelder B Wanek MR Duggal A Torbic H . Prescribing practices of valproic acid for agitation and delirium in the intensive care unit. Ann Pharmacother (2021) 55:3117. 10.1177/1060028020947173

  • 47

    Löscher W . Basic pharmacology of valproate: a review after 35 years of clinical use for the treatment of epilepsy. CNS Drugs (2002) 16:66994. 10.2165/00023210-200216100-00003

  • 48

    Löscher W . Valproate: a reappraisal of its pharmacodynamic properties and mechanisms of action. Prog Neurobiol (1999) 58:3159. 10.1016/s0301-0082(98)00075-6

  • 49

    Rakitin A Eglit T Kõks S Lember M Haldre S . Comparison of the metabolic syndrome risk in valproate-treated patients with epilepsy and the general population in Estonia. PLoS One (2014) 9:e103856. 10.1371/journal.pone.0103856

  • 50

    Rakitin A Kõks S Reimann E Prans E Haldre S . Changes in the peripheral blood gene expression profile induced by 3 months of valproate treatment in patients with newly diagnosed epilepsy. Front Neurol (2015) 6:188. 10.3389/fneur.2015.00188

  • 51

    Rakitin A Kõks S Haldre S . Metabolic syndrome and anticonvulsants: a comparative study of valproic acid and carbamazepine. Seizure (2016) 38:116. 10.1016/j.seizure.2016.03.008

  • 52

    Triyasakorn K Ubah UDB Roan B Conlin M Aho K Awale PS . The antiepileptic drug and toxic teratogen valproic acid alters microglia in an environmental mouse model of autism. Toxics (2022) 10:379. 10.3390/toxics10070379

  • 53

    Rasalam AD Hailey H Williams JH Moore SJ Turnpenny PD Lloyd DJ et al Characteristics of fetal anticonvulsant syndrome associated autistic disorder. Dev Med Child Neurol (2005) 47:5515. 10.1017/s0012162205001076

  • 54

    Koren G Nava-Ocampo AA Moretti ME Sussman R Nulman I . Major malformations with valproic acid. Can Fam Physician (2006) 52:4417.

  • 55

    Christensen J Grønborg TK Sørensen MJ Schendel D Parner ET Pedersen LH et al Prenatal valproate exposure and risk of autism spectrum disorders and childhood autism. Jama (2013) 309:1696703. 10.1001/jama.2013.2270

  • 56

    Nicolini C Fahnestock M . The valproic acid-induced rodent model of autism. Exp Neurol (2018) 299:21727. 10.1016/j.expneurol.2017.04.017

  • 57

    Manzo J Hernández‐Aguilar ME Toledo‐Cárdenas MR Herrera‐Covarrubias D Coria‐Avila GA . Dysregulation of neural tube vascular development as an aetiological factor in autism spectrum disorder: insights from valproic acid exposure. The J Physiol (2025). 10.1113/jp286899

  • 58

    Zohny SM Habib MZ Mohamad MI Elayat WM Elhossiny RM El-Salam MFA et al Memantine/aripiprazole combination alleviates cognitive dysfunction in valproic acid rat model of autism: hippocampal CREB/BDNF signaling and glutamate homeostasis. Neurotherapeutics (2023) 20:46483. 10.1007/s13311-023-01360-w

  • 59

    Taheri F Joushi S Esmaeilpour K Ebrahimi MN Taherizadeh Z Taheri P et al Transmission of behavioral and cognitive impairments across generations in rats subjected to prenatal valproic acid exposure. Birth Defects Res (2024) 116:e2309. 10.1002/bdr2.2309

  • 60

    Feleke R Jazayeri D Abouzeid M Powell KL Srivastava PK O’Brien TJ et al Integrative genomics reveals pathogenic mediator of valproate-induced neurodevelopmental disability. Brain (2022) 145:383242. 10.1093/brain/awac296

  • 61

    Park G Jang WE Kim S Gonzales EL Ji J Choi S et al Dysregulation of the Wnt/β-catenin signaling pathway via Rnf146 upregulation in a VPA-induced mouse model of autism spectrum disorder. Exp Mol Med (2023) 55:178394. 10.1038/s12276-023-01065-2

  • 62

    Go HS Seo JE Kim KC Han SM Kim P Kang YS et al Valproic acid inhibits neural progenitor cell death by activation of NF-κB signaling pathway and up-regulation of Bcl-XL. J Biomed Sci (2011) 18:48. 10.1186/1423-0127-18-48

  • 63

    Tărlungeanu DC Deliu E Dotter CP Kara M Janiesch PC Scalise M et al Impaired amino acid transport at the blood brain barrier is a cause of autism spectrum disorder. Cell (2016) 167:148194.e18. 10.1016/j.cell.2016.11.013

  • 64

    Yang J Yan B Zhao B Fan Y He X Yang L et al Assessing the causal effects of human serum metabolites on 5 major psychiatric disorders. Schizophrenia Bull (2020) 46:80413. 10.1093/schbul/sbz138

  • 65

    Golovina E Fadason T Lints TJ Walker C Vickers MH O’Sullivan JM . Understanding the impact of SNPs associated with autism spectrum disorder on biological pathways in the human fetal and adult cortex. Sci Rep (2021) 11:15867. 10.1038/s41598-021-95447-z

  • 66

    Brister D Rose S Delhey L Tippett M Jin Y Gu H et al Metabolomic signatures of autism spectrum disorder. J Personalized Med (2022) 12:1727. 10.3390/jpm12101727

  • 67

    van Sadelhoff JHJ Perez Pardo P Wu J Garssen J van Bergenhenegouwen J Hogenkamp A et al The gut-immune-brain axis in autism spectrum disorders; A focus on amino acids. Front Endocrinol (Lausanne) (2019) 10:247. 10.3389/fendo.2019.00247

  • 68

    Kaila K Price TJ Payne JA Puskarjov M Voipio J . Cation-chloride cotransporters in neuronal development, plasticity and disease. Nat Rev Neurosci (2014) 15:63754. 10.1038/nrn3819

  • 69

    Zheng HF Wang WQ Li XM Rauw G Baker GB . Body fluid levels of neuroactive amino acids in autism spectrum disorders: a review of the literature. Amino Acids (2017) 49:5765. 10.1007/s00726-016-2332-y

  • 70

    Ozeki Y Sekine M Fujii K Watanabe T Okayasu H Takano Y et al Phosphoserine phosphatase activity is elevated and correlates negatively with plasma d-serine concentration in patients with schizophrenia. Psychiatry Res (2016) 237:34450. 10.1016/j.psychres.2016.01.010

  • 71

    Deutschenbaur L Beck J Kiyhankhadiv A Muhlhauser M Borgwardt S Walter M et al Role of calcium, glutamate and NMDA in major depression and therapeutic application. Prog Neuro-Psychopharmacology Biol Psychiatry (2016) 64:32533. 10.1016/j.pnpbp.2015.02.015

  • 72

    Basu AC Tsai GE Ma CL Ehmsen JT Mustafa AK Han L et al Targeted disruption of serine racemase affects glutamatergic neurotransmission and behavior. Mol Psychiatry (2009) 14:71927. 10.1038/mp.2008.130

  • 73

    Bozek K Wei Y Yan Z Liu X Xiong J Sugimoto M et al Organization and evolution of brain lipidome revealed by large-scale analysis of human, chimpanzee, macaque, and mouse tissues. Neuron (2015) 85:695702. 10.1016/j.neuron.2015.01.003

  • 74

    Vanherle S Loix M Miron VE Hendriks JJA Bogie JFJ . Lipid metabolism, remodelling and intercellular transfer in the CNS. Nat Rev Neurosci (2025) 26:21431. 10.1038/s41583-025-00908-3

  • 75

    Dan Z Mao X Liu Q Guo M Zhuang Y Liu Z et al Altered gut microbial profile is associated with abnormal metabolism activity of Autism Spectrum disorder. Gut Microbes (2020) 11:124667. 10.1080/19490976.2020.1747329

  • 76

    Tsai CF Chuang CH Wang YP Lin YB Tu PC Liu PY et al Differences in gut microbiota correlate with symptoms and regional brain volumes in patients with late-life depression. Front Aging Neurosci (2022) 14:885393. 10.3389/fnagi.2022.885393

  • 77

    Strati F Cavalieri D Albanese D De Felice C Donati C Hayek J et al New evidences on the altered gut microbiota in autism spectrum disorders. Microbiome (2017) 5:24. 10.1186/s40168-017-0242-1

  • 78

    Morais LH Schreiber HL Mazmanian SK . The gut microbiota-brain axis in behaviour and brain disorders. Nat Rev Microbiol (2021) 19:24155. 10.1038/s41579-020-00460-0

  • 79

    Buckley MM O’Brien R Brosnan E Ross RP Stanton C Buckley JM et al Glucagon-like Peptide-1 secreting L-Cells coupled to sensory nerves translate microbial signals to the host rat nervous system. Front Cell Neurosci (2020) 14:95. 10.3389/fncel.2020.00095

  • 80

    Wu H Lv B Zhi L Shao Y Liu X Mitteregger M et al Microbiome-metabolome dynamics associated with impaired glucose control and responses to lifestyle changes. Nat Med (2025) 31:222231. 10.1038/s41591-025-03642-6

  • 81

    Fu Y Guzior DV Okros M Bridges C Rosset SL González CT et al Balance between bile acid conjugation and hydrolysis activity can alter outcomes of gut inflammation. Nat Commun (2025) 16:3434. 10.1038/s41467-025-58649-x

  • 82

    Liu Y Liu J Ren R Xin Z Luo Y Chen Y et al Short-term and long-term high-fat diet promote metabolic disorder through reprogramming mRNA m(6)A in white adipose tissue by gut microbiota. Microbiome (2025) 13:75. 10.1186/s40168-025-02047-4

  • 83

    Wang Z Wang Z Lu T Yuan G Chen W Jin J et al Gut microbiota regulate insomnia-like behaviors via gut-brain metabolic axis. Mol Psychiatry (2024) 30:2597611. 10.1038/s41380-024-02867-0

  • 84

    Xiao W Su J Gao X Yang H Weng R Ni W et al The microbiota-gut-brain axis participates in chronic cerebral hypoperfusion by disrupting the metabolism of short-chain fatty acids. Microbiome (2022) 10:62. 10.1186/s40168-022-01255-6

  • 85

    Wilson ID Nicholson JK . Gut microbiome interactions with drug metabolism, efficacy, and toxicity. Translational Res (2017) 179:20422. 10.1016/j.trsl.2016.08.002

  • 86

    Maiuolo J Gliozzi M Musolino V Carresi C Scarano F Nucera S et al The contribution of gut microbiota-brain axis in the development of brain disorders. Front Neurosci (2021) 15:616883. 10.3389/fnins.2021.616883

  • 87

    Srikantha P Mohajeri MH . The possible role of the microbiota-gut-brain-axis in autism spectrum disorder. Int J Mol Sci (2019) 20:2115. 10.3390/ijms20092115

Summary

Keywords

gas chromatography-mass spectrometry, 16S ribosomal RNA, gut microbiota, metabolite, autism

Citation

Zhao S, Zhang X, Miao Y, Gao X, Wan Q, Qiu W, Si H, Han Y, Du X, Feng Y, Liu L and Chen Y (2025) Integrated analysis of metabolome and microbiome in a mouse model of sodium valproate-induced autism. Exp. Biol. Med. 250:10452. doi: 10.3389/ebm.2025.10452

Received

29 November 2024

Accepted

06 August 2025

Published

29 August 2025

Volume

250 - 2025

Updates

Copyright

*Correspondence: Lianhua Liu, ; Yuqing Chen,

†These authors have contributed equally to this work

Disclaimer

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

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article