Gut microbiota and metabolite variations in a migraine mouse model

1 Introduction

Migraine, a common neurological disorder, is characterized by recurrent unilateral throbbing headaches accompanied by symptoms such as nausea, vomiting, photophobia, and phonophobia (Schiffman et al., 2013). Despite its high prevalence, the diagnosis, pathophysiology, and treatment of migraine remain contentious. Migraine has been widely acknowledged as a distinct disorder associated with neurological disturbances. However, the precise mechanisms underlying its pathogenesis are not fully understood, and there is a lack of unified explanations for patient symptoms in both clinical studies and basic experiments.

The hypothesis that the activation of the trigeminovascular system can trigger the release of various neurotransmitters, such as calcitonin gene–related peptide and substance P leading to headaches, has gained wide acceptance (Vass et al., 2001; Vass et al., 2004). Another study has suggested that migraine results from neurohormonal and metabolic integrity loss (Dzugan and Dzugan, 2015). In recent years, the human gut microbiota has emerged as a focal point of research due to its role as the largest reservoir of human microbiota. It forms a symbiotic ecosystem with the host, influencing human health through metabolites, and interactions with the host’s immune system. Previously, the central nervous system was thought to influence the gut microbiota. However, with the rapid development of microbial sequencing technologies, it has become evident that the gut microbiota functions like an endocrine organ, producing various biologically active substances that enter the bloodstream and modulate distant organ functions (Tang et al., 2017). The gut microbiota exerts considerable influence on the central nervous system through multiple pathways. This realization has led to the concept of the “gut–brain axis” (Cryan et al., 2019). Abundant research has demonstrated the pivotal role of the gut microbiota in regulating various central nervous system disorders, such as anxiety disorders, depression, and Parkinson’s disease (Ma et al., 2019; Li et al., 2023). However, the mechanisms through which gut microbes influence the nervous system and their roles in disease pathogenesis remain unclear. Potential pathways through which this interaction occurs include the vagus nerve pathway, which connects the brain and the gut; the neuroendocrine pathway, especially the hypothalamic–pituitary–adrenal axis; metabolic pathways influenced by gut microbiota metabolites; and immune system pathways involving various immune factors (Décarie-Spain et al., 2023). Within the neuroendocrine pathway, the gut microbiota can affect the hypothalamic–pituitary–adrenal axis by modulating the secretion of 5-hydroxytryptamine(5-HT) from enteroendocrine cells into the peripheral bloodstream, thereby affecting the central nervous system (Kasarello et al., 2023). The gut microbiota produces specific metabolites, including amines (such as histamine, tyramine, putrescine, and cadaverine). For example, amines can be absorbed in situ and may affect the host’s behavior through unknown mechanisms, potentially contributing to central nervous system disorders (Gabastou et al., 1996). Alterations in the gut microbiota have been associated with various conditions, including metabolic disorders, obesity, hematological malignancies, neurological or behavioral disorders, and migraines. Metabolites produced by the gut microbiome have been shown to influence the gut–brain axis. The use of probiotics as dietary supplements may reduce the frequency and severity of migraine attacks, For instance, a study found that an IgG elimination diet combined with probiotics reduced migraine attack frequency in individuals with migraine and irritable bowel syndrome (Xie et al., 2019). Fecal microbiota transplantation represents a promising approach to treating migraine (Kappéter et al., 2023).

Therefore, this study established a mouse model of migraine to investigate the differences in gut microbiota and gut metabolites between migraine-afflicted mice and normal mice and to explore the mechanisms underlying the pathogenesis of migraine and potential therapeutic avenues.

2 Materials and methods

2.1 Experimental design and modeling

2.1.1 Animals

C57BL/6J WT mice (C57BL/6-Htr1Dem1Smoc, aged 6–8 weeks) obtained from Nanmo Bio were used as experimental subjects and kept under a normal diet and regular 12 h light/dark cycles. All animal procedures were performed by the experimental protocol approved by the animal ethics committee of the medical school at Fudan University in China. This research involved a cohort of 40 mice, with 18 designated for the primary experimental protocols and 22 utilized for procedural training sessions.

2.1.2 Experiment protocol

To augment the reliability of our experimental outcomes, the mice were divided into 3 groups, the gender ratio in each group of mice was 1:1, and we implemented two control groups: one subjected to no interventions (CON: n = 6) and another administered intraperitoneal injections of physiological saline solution(NS: n = 6), the migraine model group (CM: n = 6). The CM group was intraperitoneally injected with nitroglycerin (10 mg/kg) every other day for nine consecutive days (Demartini et al., 2019). The NS group received intraperitoneal injections of an equivalent volume of 0.9% saline. The CON group (CON) did not receive any intervention.

2.1.3 Pain threshold detection

Using von Frey mechanical stimulation filaments (Stoelting), changes in mechanical stimulation sensitivity thresholds around the eyes and hind paws of the mice were measured on days 1, 3, 5, 7, and 9 of modeling and days 1, 3, 5, 7, and 9 after modeling. A series of von Frey filaments with incrementally increasing stimulus forces (0.07, 0.16, 0.40, 0.60, 1.00, 1.40, and 2.00 g) were used. To acclimate themselves to the testing environment, the mice were placed in it for 1 h/day for two to three consecutive days before testing. The temperature and humidity of the experimental environment were kept stable. The sensitivity thresholds were recorded.

2.2 16S rRNA gene sequencing of the gut microbiome

Total genome DNA from samples was extracted using the cetyltrimethylammonium bromide method. DNA concentration and purity was monitored on 1% agarose gels. According to the concentration, DNA wasdiluted to 1ng/L using sterile water.16S rRNA genes of distinct regions (16S V4/16S V3/16S V3-V4/16SV4-V5) were amplified used the specific primer (16SV4: 515F- 806R), with the barcode. All PCR reactions were carried out with 15 L of Phusion® High -Fidelity PCR Master Mix (New England Biolabs); mixture PCR products was purified with Qiagen Gel Extraction Kit (Qiagen, Germany). Sequencing libraries were generated usingTruSeq® DNA PCR-Free Sample Preparation Kit(Illumina, USA) following manufacturer’s recommendations and index codes were added.The library quality was assessed on the Qubit@ 2.0 Fluorometer (Thermo Scientific) and Agilent Bioanalyzer 2100 system. At last, the library was sequenced on an Illumina NovaSeq platform and 250 bp paired-end reads were generated.

In order to study phylogenetic relationship of different OTUs, and dominant species in different samples (groups), multiple sequence alignment were conducted using the MAFFT (v7.490, https://mafft.cbrc.jp/alignment/software/) (Caporaso et al., 2010). OTUs abundance information were normalized using a standard of sequence number corresponding to the sample with the least sequences. Subsequent analysis of alpha diversity and beta diversity were all performed basing on this output normalized data. Alpha diversity is applied in analyzing complexity of species diversity for a sample through 6 indices, including Observed-species, Chao1, Shannon, Simpson, ACE, Good-coverage. All this indices in our samples were calculated with QIIME and displayed with R software (Version 4.1.2). Beta diversity analysis was used to evaluate differences of samples in species complexity, Beta diversity on both weighted and unweighted unifrac were calculated by QIIME software. Cluster analysis was preceded by principal component analysis (PCA), which was applied to reduce the dimension of the original variables using the stats package and ggplot2 package in R software. Principal Coordinate Analysis (PCoA) was performed to get principal coordinates and visualize from complex, multidimensional data. A distance matrix of weighted or unweighted unifrac among samples obtained before was transformed to a new set of orthogonal axes, by which the maximum variation factor is demonstrated by first principal coordinate, and the second maximum one by the second principal coordinate, and so on. PCoA analysis was displayed by stats package and ggplot2 package in R software. Unweighted Pair-group Method with Arithmetic Means (UPGMA) Clustering was performed as a type of hierarchical clustering method to interpret the distance matrix using average linkage and was conducted by QIIME software.

2.3 Gut metabolomics analysis

The original data file acquisited by LC-MS was converted into mzML format by ProteoWizard software. Peak extraction, peak alignment and retention time correction were respectively performed by XCMS program. The “SVR” method was used to correct the peak area. The peaks with detetion rate lower than 50% in each group of samples were discarded. After that, metabolic identification information was obtained by searching the laboratory’s self-built database, integrated public database, AI database and metDNA.

Unsupervised PCA (principal component analysis) was performed by statistics function prcomp within R (www.r-project.org). The data was unit variance scaled before unsupervised PCA.The HCA (hierarchical cluster analysis) results of samples and metabolites were presented as heatmaps with dendrograms, while pearson correlation coefficients (PCC) between samples were caculated by the cor function in R and presented as only heatmaps. Both HCA and PCC were carried out by R package ComplexHeatmap. For HCA, normalized signal intensities of metabolites (unit variance scaling) are visualized as a color spectrum.

For two-group analysis, differential metabolites were determined by VIP (VIP > 1) and P-value (P-value < 0.05, Student’s t test). For multi-group analysis, differential metabolites were determined by VIP (VIP > 1) and P-value (P-value < 0.05, ANOVA). VIP values were extracted from OPLS-DA result, which also contain score plots and permutation plots, was generated using R package MetaboAnalystR. The data was log transform (log22) and mean centering before OPLS-DA. In order to avoid overfitting, a permutation test (200 permutations) was performed.

Identified metabolites were annotated using KEGG Compound database (http://www.kegg.jp/kegg/compound/), annotated metabolites were then mapped to KEGG Pathway database (http://www.kegg.jp/kegg/pathway.html). Significantly enriched pathways are identified with a hypergeometric test’s P-value for a given list of metabolites.

2.4 Statistical analysis

The Wilcoxon signed-rank test or Student’s t-test was used for comparisons between multiple groups. Spearman’s correlation coefficients were used to determine the correlations between gut microbiota and metabolites. P-values of <0.05 were considered statistically significant. Figures and charts were created using R software (Version 4.1.2).

3 Results

3.1 Pain threshold

The pain thresholds of mice in the CM group remained consistently lower than those in the NS group (p < 0.05; Figure 1). This confirmed the successful establishment of the migraine model, and no adverse or toxic effects were observed in the model.

Figure 1

Figure 1 This line chart demonstrates the mechanical stimulus sensitivity thresholds of the NS and CM groups assessed on days 1, 3, 5, 7, and 9 after modeling. The pain threshold of the CM group was notably lower than that of the NS group. *P<0.05, the differences are statistically significant.

3.2 Mouse gut microbiota

3.2.1 Microbial diversity analysis

The Shannon index (Figure 2A) and the richness estimator Chao1 (Figure 2B) revealed no significant differences in gut microbiota diversity between the three groups.However, the unweighted_unifrac WilcoxTest analysis (Figure 2C) revealed significant differences in beta diversity between the NS and CM groups (p = 0.01). Moreover, ANOSIM (Figure 2D) revealed significant differences between the three groups (r = 0.51, p = 0.001).

Figure 2

Figure 2 Alpha diversity of gut bacteria: Box plots representing (A) community diversity and (B) richness were used to estimate the alpha diversity in the CM (blue), CON (red), and NS (green) groups. The boxes show the medians, interquartile ranges, and minimum and maximum values. Beta diversity of gut microbiota: (C) Box plot based on OTUs’ unweighted UniFrac beta diversity. (D) ANOSIM (r= 0.51, p = 0.001) revealed significant differences between the NS and CM groups.

3.2.2 Species difference and contribution analysis

Fecal 16S rRNA sequencing and t-tests revealed significant intergroup differences in species based on OTUs. Notably, the abundance of Akkermansia muciniphila was markedly higher in the NS group than in the CM group (Figure 3A). Akkermansia muciniphila exhibited the greatest dissimilarity between the NS and CM groups, with a substantially greater abundance in the former than in the latter (Figure 3B).

Figure 3

Figure 3 (A) Statistical analysis of gut microbiota in the CM and NS groups based on Student’s t-test. The abundance of Akkermansia muciniphila was significantly higher in the NS group than in the CM group. (B) SIMPER analysis showing the top 10 species contributing to the differences between the NS and CM groups. Akkermansia muciniphila had the greatest contribution.

3.2.3 Comparison of population abundances

The results of the UPGMA clustering analysis (Figure 4A) indicated differential relative abundances between the groups at the phylum level. The linear discriminant analysis (LDA) value distribution histogram obtained from the LEfSe analysis (Figure 4B) also revealed microbes that were differentially abundant in the NS and CM groups. The NS group was characterized by a high abundance of Akkermansia muciniphila, while the CM group was characterized by a high abundance of unidentified Clostridiaceae. A phylogenetic tree (Figure 4C) showed that Akkermansiaceae constituted the most abundant taxon in the NS group, while Lachnospiraceae constituted the most prevalent group in the CM group.

Figure 4

Figure 4 (A) The UPGMA clustering analysis indicates the differences between the three groups at the phylum level. he brightness of each point corresponds to the magnitude of its LDA effect size (B). The LDA scores are combined with the effect size measurements in the CM and NS groups. Abundant taxonomic groups in the NS (green) and CM (red) groups are displayed along with the LDA scores. Only taxonomic groups with LDA scores exceeding the threshold of 4.0 are displayed (B). (C) Cladogram representing the gut microbiota in the CM and NS groups constructed using 16S rRNA sequencing. Abundant taxonomic groups are indicated for both the NS (green) and CM (red) groups.

3.3 Mouse gut metabolomics

3.3.1 Metabolite PCA

A PCA plot of metabolites was conducted for the NS and CM groups (Figure 5A) and for the VM, NS, and CON groups (Figure 5C). The results showed distinct intergroup differences in the gut metabolomes of the mice. The orthogonal partial least squares discriminant analysis (OPLS-DA) scores indicated clear discrimination between the NS, CM, and CON groups (Figures 5A, C). The goodness of fit values and predictive ability values (NS vs. CM: R2X = 0.511, R2Y = 0.980, Q2 = 0.833, p < 0.05; CON vs. CM and CON: R2X = 0.320, R2Y = 0.961, Q2 = 0.790, p < 0.05) showed that the OPLS-DA model exhibited satisfactory fitting and great predictive ability (Figures 5B, D).

Figure 5

Figure 5 OPLS-DA analysis was performed to show the metabolite profile (NS: n = 6, CON: n = 6, CM: n = 6). The OPLS-DA score and model validation plots showed clear discrimination between the NS and CM groups, with R2X = 0.511, R2Y = 0.980, Q2 = 0.833, and p < 0.05 (A, B). They also showed clear discrimination between the NS, CM, and CON groups, with R2X = 0.320, R2Y = 0.961, Q2 = 0.790, and p < 0.05 (C, D).

3.3.2 Differential metabolites

Based on the OPLS-DA results, we used variable importance in projection (VIP) as a threshold to further screen for metabolite differences between the groups, considering metabolites with a VIP value of ≥1.0, a p-value of <0.05, and a fold change value of ≥2 or ≤0.5 to be differentially regulated. A total of 633 differential metabolites were detected in the NS and CM groups. Of those, 145 showed significantly higher levels in the CM group than in the NS group, while 488 exhibited significantly lower levels (Figure 6).

Figure 6

Figure 6 Volcano plot showing differential metabolites between the NS and CM groups (A). Linear discriminant analysis (LDA) effect size (LEfSe) results showed the top 20 differential metabolites between the NS and CM groups (B). Hierarchical cluster analysis of differential metabolites between the NS and CMgroups (C). Each column represents a sample, and each row stands for a metabolite. Each column represents a sample, and each row represents a metabolite. The differential metabolites were determined using the following criteria: VIP ≥ 1.0, p < 0.05, and fold change ≥ 2 or ≤ 0.5.

3.3.3 Metabolic pathway analysis

To identify the involved metabolic pathways, the differential metabolites were examined against the KEGG database. The results revealed the involvement of various metabolic pathways (Figure 7). Some of the main pathways were bile secretion, steroid hormone biosynthesis, biosynthesis of neomycin and kanamycin, glutathione metabolism, and biotin metabolism.

Figure 7

Figure 7 Top 20 pathways with the most pronounced differences in gut metabolites between the NS and CM groups. The size of each origin corresponds to the number of differential metabolites, while the colors indicate different levels of statistical significance. The differential metabolites were primarily associated with pathways such as bile secretion, steroid hormone biosynthesis, biosynthesis of neomycin and kanamycin, glutathione metabolism, and biotin metabolism.

3.4 Relationship between differential microbes and differential metabolites

We calculated Spearman’s correlation coefficients to determine the functional correlations between variations in gut microbiota and metabolite profiles in the groups, with r ≥ 0.8 and p < 0.05 considered to indicate significant correlations. The results (Figure 8) indicated the following correlations between the abundances of specific microbes and metabolites: Akkermansia muciniphila–glycycoumarin (positive), Lachnospiraceae bacterium GAM79S–dehydrolithocholic acid and carnitine C14:3 (positive), Akkermansia muciniphila–Ala-Lys-Asn-Glu (negative), and S-Lachnospiraceae bacterium 10-1–withaferin A and cysteine (negative).

Figure 8

Figure 8 Scatterplots illustrating significant correlations between the altered relative abundances of gut bacterial genera and the LC-MS spectral intensities of differential metabolites. (A): the correlations between the abundances of S-Akkermansia muciniphila and glycycoumarin is positive, (B) and (C): the correlations between the abundances of Lachnospiraceae bacterium GAM79S and dehydrolithocholic acid and carnitine C14:3 is positive. (D) and (E): the correlations between the abundances of Lachnospiraceae bacterium 10-1 and withaferin A and cysteine is negative, (F) the correlations between the abundances of S-Akkermansia muciniphila and Ala-Lys-Asn-Glu is negative.

4 Conclusion

This study provides evidence of differences in gut microbiota composition and metabolites between a migraine mouse model and normal mice, which showed that Akkermansiaceae constituted the most abundant taxon in the NS group, while Lachnospiraceae constituted the most prevalent group in the CM group. The associations between the abundances of Akkermansia muciniphila and Lachnospiraceae bacteria and metabolites suggested their potential roles in the pathogenesis of migraine. The altered abundance of Lachnospiraceae observed in migraine-afflicted mice and its correlations with changes in metabolites suggest that it may affect the host’s health. Overall, this study provides insights into the potential roles of the gut microbiota and metabolites in migraines, indicating potential research avenues and therapeutic strategies.

5 Discussion

In this study, we observed a decrease in the abundance of Akkermansia muciniphila in the guts of migraine-afflicted mice. A combined analysis of differential metabolites and differential bacteria revealed a positive correlation between the abundance of Akkermansia muciniphila and glycycoumarin. Glycycoumarin, a representative coumarin compound, exhibits significant pharmacological activities, including antitumor, hepatoprotective, antispasmodic, antibacterial, and antiviral effects (Tang et al., 2023). We also found a negative correlation between the abundance of Akkermansia muciniphila and the molecular substances named Ala-Lys-Asn-Glu. The biological effects of this peptide segment remain unexplored. However, they may be related to amino acid metabolism.

Akkermansia muciniphila, a Gram-negative anaerobic mucinophilic bacterium, is a unique representative of the Verrucomicrobia phylum. It thrives in the mucus layer of the intestinal epithelium by utilizing mucus as a source of carbon and nitrogen. It has been associated with improving host metabolic functions and immune responses, which makes it a promising candidate for probiotic therapy (Derrien et al., 2004; Zhang et al., 2019; Zhou and Zhang, 2019). A decrease in the abundance of Akkermansia muciniphila in the guts of migraine-afflicted mice and its correlation with specific metabolites, such as glycycoumarin, indicate potential protective action against migraine (Cani et al., 2022). Recent studies have indicated that Akkermansia in the gut can regulate serum 5-HT concentrations and influence central 5-HT levels through the gut–brain axis, suggesting a novel therapeutic strategy for treating 5-HT-related disorders (Yaghoubfar et al., 2020; Yaghoubfar et al., 2021; Sun et al., 2023). Certainly, some studies indicate exercise might reduce migraine frequency, however, the optimal exercise and intensity lack consensus. Moreover, evidence suggests vigorous exercise might provoke migraines (Amin et al., 2018). Nevertheless, contradictory findings show reduced Akkermansia levels in the gut microbiota of elite athletes (Hintikka et al., 2022).Perhaps high-intensity exercise could lead to a reduction in Akkermansia abundance, which might have a detrimental effect on migraine control.

In our study, the abundance of Lachnospiraceae in the guts of migraine-afflicted mice was greater than in those of normal mice. A combined analysis of differential metabolites and differential bacteria revealed a positive correlation between carnitine C14:3 and the abundance of Lachnospiraceae bacterium GAM79S. Carnitine plays a crucial role in energy metabolism. However, many of its functions remain unclear (Borum, 1983). We also found a negative correlation between cyasterone and S-Lachnospiraceae bacterium 10-1. Cyasterone can promote mesenchymal stem cell migration and osteogenesis, accelerating bone fracture healing (Zhu et al., 2021). Furthermore, we found a positive correlation between dehydrolithocholic acid and Lachnospiraceae bacterium GAM79S. We also found a negative correlation between withaferin A and S-Lachnospiraceae bacterium 10-1. Withaferin A is a C28 steroidal lactone derived from the plant Withania somnifera, commonly known as ashwagandha. It has attracted considerable attention due to its anticancer properties, which have been observed in various cancer cell lines from different sources, as well as its anti-inflammatory, metabolic, and proapoptotic activities (Atteeq, 2022). Our findings suggest that changes in the abundance of Lachnospiraceae can lead to alterations in gut metabolites and that different bacteria have different effects on metabolites, which, in turn, have different effects on the host’s health. Based on our findings, we can speculate that S-Lachnospiraceae bacterium 10-1 has adverse effects on the host’s health, while Lachnospiraceae bacterium GAM79S plays a beneficial role.

Different taxonomic groups within the Lachnospiraceae family have been associated with various diseases, and their impacts on host physiology often vary (Vacca et al., 2020). According to one analysis, the diversity of Lachnospiraceae may affect the ability of isolates to influence the host’s health (Sorbara et al., 2020). Schizophrenia is a debilitating psychiatric disorder with undefined underlying molecular mechanisms. The gut microbiota can regulate brain function and behavior through the microbiota–gut–brain axis. A previous study found a lower gut microbiota alpha diversity index and significantly greater gut microbial dysbiosis in both medicated and medication-naive SCZ patients than in healthy controls. The same study reported that several unique bacterial taxa, such as Veillonellaceae and Lachnospiraceae, were associated with the severity of SCZ (Zheng et al., 2019). Fusimonas intestini, a commensal species of the Lachnospiraceae family, highly colonizes obese and hyperglycemic humans and mice, producing long-chain fatty acids, such as elaidic acid, thus promoting diet-induced obesity (Takeuchi et al., 2023). However, another study reported that the abundances of unclassified Mogibacteriaceae, Lachnospira, and Slackia exhibited significant negative correlations with dehydrolithocholic acid (Liu et al., 2022).

In this study, we found significant differences in metabolites between the migraine mouse model (CM group) and the normal mice (NS group). The pathways most closely associated with these differential metabolites included bile secretion, steroid hormone biosynthesis, biosynthesis of neomycin and kanamycin, glutathione metabolism, and biotin metabolism. Among them, bile secretion and steroid hormone biosynthesis were the most severely affected, suggesting that abnormalities in bile and steroid metabolism may be associated with the pathogenesis of migraine. Regrettably, we did not establish a baseline comparison in animals before modeling. Our analysis of the gut microbiota in mice occurred post-modeling. Furthermore, our study lacked testing at various intervals, precluding precise determination of the timeline for changes in the gut microbiota of the animal model. Overall, this study provides insights into the potential roles of the gut microbiota and metabolites in migraines, indicating potential research avenues and therapeutic strategies, For instance, a study found that an IgG elimination diet combined with probiotics reduced migraine attack frequency in individuals with migraine and irritable bowel syndrome (Xie et al., 2019). However, given the current ambiguity surrounding the relationship between migraines and changes in gut microbiota, additional investigations are warranted. Consequently, evaluating its efficacy in mitigating migraines remains a complex task at present.Firstly, we need to ascertain the probiotics’ effectiveness timeline in the model, conducting pain threshold assessments at different intervals will provide insights into the onset of probiotic efficacy based on the pain threshold variations observed in mice. Secondly, For human subjects, prospective cohort studies are imperative to scrutinize the gut microbiota at different disease stages within the same patient, allowing inferences about the timing of microbiota alterations. Concerning the efficacy of probiotics, evaluating migraine symptom severity in patients post-probiotic treatment will facilitate discernment of the drug’s onset of action.

The study’s limitations include its small sample size and a lack of clinical validation. Further studies with larger samples and clinical investigations are necessary to confirm our findings. Furthermore, this study lacks an analysis of changes in pain thresholds in mice models post-fecal transplantation. This aspect will be addressed in our subsequent research endeavors. Moreover, this study focused on the acute stage of migraine. Future studies should conduct longitudinal observations and investigate the potential dynamic changes in the gut microbiota occurring during the development and treatment of migraine. Blood samples should be concurrently analyzed to validate the functional effects of differential metabolites on remote organs.

Data availability statement

The accession number of our sequencing data is PRJNA1061318 in the SRA database and has been released successfully.

Ethics statement

The animal study was approved by The ethical committee of the EENT Hospital. The study was conducted in accordance with the local legislation and institutional requirements.

Author contributions

DW: Formal analysis, Methodology, Validation, Writing – review & editing. XL: Data curation, Formal analysis, Writing – original draft, Writing – review & editing. SS: Data curation, Writing – review & editing, Investigation. TR: Data curation, Investigation, Writing – review & editing. WW: Funding acquisition, Project administration, Supervision, Validation, Writing – review & editing.

Funding

The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This study was supported by the Natural Science Foundation of Shanghai (No. 20ZR1409600).

Conflict of interest

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

Publisher’s note

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

References

Amin, F. M., Aristeidou, S., Baraldi, C., Czapinska-Ciepiela, E. K., Ariadni, D. D., Di Lenola, D., et al. (2018). The association between migraine and physical exercise. J. Headache Pain. 19 (1), 83. doi: 10.1186/s10194-018-0902-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Atteeq, M. (2022). Evaluating anticancer properties of Withaferin A-a potent phytochemical. Front. Pharmacol. 13, 975320. doi: 10.3389/fphar.2022.975320

PubMed Abstract | CrossRef Full Text | Google Scholar

Borum, P. R. (1983). Carnitine. Annu. Rev. Nutr. 3, 233–259. doi: 10.1146/annurev.nu.03.070183.001313

PubMed Abstract | CrossRef Full Text | Google Scholar

Cani, P. D., Depommier, C., Derrien, M., Everard, A., de Vos, W. M. (2022). Akkermansia muciniphila: paradigm for next-generation beneficial microorganisms. Nat. Rev. Gastroenterol. Hepatol. 19 (10), 625–637. doi: 10.1038/s41575-022-00631-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Caporaso, J. G., Kuczynski, J., Stombaugh, J., Bittinger, K., Bushman, F. D., Costello, E. K., et al. (2010). QIIME allows analysis of high-throughput community sequencing data. Nat. Methods 7 (5), 335–336. doi: 10.1038/nmeth.f.303

PubMed Abstract | CrossRef Full Text | Google Scholar

Cryan, J. F., O'Riordan, K. J., Cowan, C. S. M., Sandhu, K. V., Bastiaanssen, T. F. S., Boehme, M., et al. (2019). The microbiota-gut-brain axis. Physiol. Rev. 99 (4), 1877–2013. doi: 10.1152/physrev.00018.2018

PubMed Abstract | CrossRef Full Text | Google Scholar

Décarie-Spain, L., Hayes, A. M. R., Lauer, L. T., Kanoski, S. E. (2023). The gut-brain axis and cognitive control: A role for the vagus nerve. Semin. Cell Dev. Biol. 156, 201–209. doi: 10.1016/j.semcdb.2023.02.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Demartini, C., Greco, R., Zanaboni, A. M., Sances, G., De Icco, R., Borsook, D., et al. (2019). Nitroglycerin as a comparative experimental model of migraine pain: From animal to human and back. Prog. Neurobiol. 177, 15–32. doi: 10.1016/j.pneurobio.2019.02.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Derrien, M., Vaughan, E. E., Plugge, C. M., de Vos, W. M. (2004). Akkermansia muciniphila gen. nov., sp. nov., a human intestinal mucin-degrading bacterium. Int. J. Syst. Evol. Microbiol. 54 (Pt 5), 1469–1476. doi: 10.1099/ijs.0.02873-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Dzugan, S. A., Dzugan, K. S. (2015). Is migraine a consequence of a loss of neurohormonal and metabolic integrity? A new hypothesis. Neuro Endocrinol. Lett. 36 (5), 421–429.

PubMed Abstract | Google Scholar

Gabastou, J. M., Nugon-Baudon, L., Robert, Y., Manuel, C., Vaissade, P., Bourgeon, E., et al. (1996). Digestive amines of bacterial origin and behavior disorders. Apropos of a case. Pathol. Biol. (Paris). 44 (4), 275–281.

PubMed Abstract | Google Scholar

Hintikka, J. E., Munukka, E., Valtonen, M., Luoto, R., Ihalainen, J. K., Kallonen, T., et al. (2022). Gut microbiota and serum metabolome in elite cross-country skiers: A controlled study. Metabolites 12 (4), 335. doi: 10.3390/metabo12040335

PubMed Abstract | CrossRef Full Text | Google Scholar

Kappéter, Á., Sipos, D., Varga, A., Vigvári, S., Halda-Kiss, B., Péterfi, Z. (2023). Migraine as a disease associated with dysbiosis and possible therapy with Fecal microbiota transplantation. Microorganisms 11 (8), 2083. doi: 10.3390/microorganisms11082083

PubMed Abstract | CrossRef Full Text | Google Scholar

Kasarello, K., Cudnoch-Jedrzejewska, A., Czarzasta, K. (2023). Communication of gut microbiota and brain via immune and neuroendocrine signaling. Front. Microbiol. 14, 1118529. doi: 10.3389/fmicb.2023.1118529

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, Y., Chen, Y., Fan, Y., Chen, Y., Chen, Y. (2023). Dynamic network modeling of gut microbiota during Alzheimer's disease progression in mice. Gut Microbes 15 (1), 2172672. doi: 10.1080/19490976.2023.2172672

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, Y., Azad, M. A. K., Zhang, W., Xiong, L., Blachier, F., Yu, Z., et al. (2022). Intrauterine growth retardation affects liver bile acid metabolism in growing pigs: effects associated with the changes of colonic bile acid derivatives. J. Anim. Sci. Biotechnol. 13 (1), 117. doi: 10.1186/s40104-022-00772-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Ma, Q., Xing, C., Long, W., Wang, H. Y., Liu, Q., Wang, R. F. (2019). Impact of microbiota on central nervous system and neurological diseases: the gut-brain axis. J. Neuroinflamm. 16 (1), 53. doi: 10.1186/s12974-019-1434-3

CrossRef Full Text | Google Scholar

Schiffman, E., Ohrbach, R., Truelove, E., Look, J., Anderson, G., Goulet, J-P., et al. (2013). The International Classification of Headache Disorders. 3rd edition. (Cephalalgia). An International Journal of Headache. 33 (9), 629–808. doi: 10.1177/0333102413485658

CrossRef Full Text | Google Scholar

Sorbara, M. T., Littmann, E. R., Fontana, E., Moody, T. U., Kohout, C. E., Gjonbalaj, M., et al. (2020). Functional and genomic variation between human-derived isolates of Lachnospiraceae reveals inter- and intra-species diversity. Cell Host Microbe 28 (1), 134–46.e4. doi: 10.1016/j.chom.2020.05.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Sun, Y., Zhu, H., Cheng, R., Tang, Z., Zhang, M. (2023). Outer membrane protein Amuc_1100 of Akkermansia muciniphila alleviates antibiotic-induced anxiety and depression-like behavior in mice. Physiol. Behav. 258, 114023. doi: 10.1016/j.physbeh.2022.114023

PubMed Abstract | CrossRef Full Text | Google Scholar

Takeuchi, T., Kameyama, K., Miyauchi, E., Nakanishi, Y., Kanaya, T., Fujii, T., et al. (2023). Fatty acid overproduction by gut commensal microbiota exacerbates obesity. Cell Metab. 35 (2), 361–75.e9. doi: 10.1016/j.cmet.2022.12.013

PubMed Abstract | CrossRef Full Text | Google Scholar

Tang, W. H., Kitai, T., Hazen, S. L. (2017). Gut microbiota in cardiovascular health and disease. Circ. Res. 120 (7), 1183–1196. doi: 10.1161/CIRCRESAHA.117.309715

PubMed Abstract | CrossRef Full Text | Google Scholar

Tang, Y., Ou, S., Ye, L., Wang, S. (2023). Pharmacological activities and pharmacokinetics of glycycoumarin. Rev. Bras. Farmacogn. 33 (3), 471–483. doi: 10.1007/s43450-022-00342-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Vacca, M., Celano, G., Calabrese, F. M., Portincasa, P., Gobbetti, M., De Angelis, M. (2020). The controversial role of human gut Lachnospiraceae. Microorganisms 8 (4), 573. doi: 10.3390/microorganisms8040573

PubMed Abstract | CrossRef Full Text | Google Scholar

Vass, Z., Dai, C. F., Steyger, P. S., Jancsó, G., Trune, D. R., Nuttall, A. L. (2004). Co-localization of the vanilloid capsaicin receptor and substance P in sensory nerve fibers innervating cochlear and vertebro-basilar arteries. Neuroscience 124 (4), 919–927. doi: 10.1016/j.neuroscience.2003.12.030

PubMed Abstract | CrossRef Full Text | Google Scholar

Vass, Z., Steyger, P. S., Hordichok, A. J., Trune, D. R., Jancsó, G., Nuttall, A. L. (2001). Capsaicin stimulation of the cochlea and electric stimulation of the trigeminal ganglion mediate vascular permeability in cochlear and vertebro-basilar arteries: a potential cause of inner ear dysfunction in headache. Neuroscience 103 (1), 189–201. doi: 10.1016/S0306-4522(00)00521-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Xie, Y., Zhou, G., Xu, Y., He, B., Wang, Y., Ma, R., et al. (2019). Effects of diet based on IgG elimination combined with probiotics on migraine plus irritable bowel syndrome. Pain Res. Manage. 2019, 7890461. doi: 10.1155/2019/7890461

CrossRef Full Text | Google Scholar

Yaghoubfar, R., Behrouzi, A., Ashrafian, F., Shahryari, A., Moradi, H. R., Choopani, S., et al. (2020). Modulation of serotonin signaling/metabolism by Akkermansia muciniphila and its extracellular vesicles through the gut-brain axis in mice. Sci. Rep. 10 (1), 22119. doi: 10.1038/s41598-020-79171-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Yaghoubfar, R., Behrouzi, A., Zare Banadkoki, E., Ashrafian, F., Lari, A., Vaziri, F., et al. (2021). Effect of Akkermansia muciniphila, Faecalibacterium prausnitzii, and Their Extracellular Vesicles on the Serotonin System in Intestinal Epithelial Cells. Probiotics Antimicrob. Proteins. 13 (6), 1546–1556. doi: 10.1007/s12602-021-09786-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, T., Li, Q., Cheng, L., Buch, H., Zhang, F. (2019). Akkermansia muciniphila is a promising probiotic. Microb. Biotechnol. 12 (6), 1109–1125. doi: 10.1111/1751-7915.13410

PubMed Abstract | CrossRef Full Text | Google Scholar

Zheng, P., Zeng, B., Liu, M., Chen, J., Pan, J., Han, Y., et al. (2019). The gut microbiome from patients with schizophrenia modulates the glutamate-glutamine-GABA cycle and schizophrenia-relevant behaviors in mice. Sci. Adv. 5 (2), eaau8317. doi: 10.1126/sciadv.aau8317

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhou, J. C., Zhang, X. W. (2019). Akkermansia muciniphila: a promising target for the therapy of metabolic syndrome and related diseases. Chin. J. Nat. Med. 17 (11), 835–841. doi: 10.1016/S1875-5364(19)30101-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhu, J., Liu, Y., Chen, C., Chen, H., Huang, J., Luo, Y., et al. (2021). Cyasterone accelerates fracture healing by promoting MSCs migration and osteogenesis. J. Orthop. Translat. 28, 28–38. doi: 10.1016/j.jot.2020.11.004

PubMed Abstract | CrossRef Full Text | Google Scholar