Your new experience awaits. Try the new design now and help us make it even better

ORIGINAL RESEARCH article

Front. Microbiol., 02 January 2026

Sec. Microorganisms in Vertebrate Digestive Systems

Volume 16 - 2025 | https://doi.org/10.3389/fmicb.2025.1697206

This article is part of the Research TopicGut Microbiome Influences in Mood Disorders: Unveiling the Gut-Brain AxisView all 8 articles

Gut microbiota and traumatic brain injury: insights from an antibiotic-free cohort


Hao Wu&#x;Hao Wu1Yuqi Huo&#x;Yuqi Huo2Wei FangWei Fang3Jinlin GuoJinlin Guo4Xiaoqin WangXiaoqin Wang5Li Su
Li Su5*Gang Cheng
Gang Cheng6*
  • 1Department of Neurosurgery, Shanxi Bethune Hospital, Shanxi Academy of Medical Sciences, Third Hospital of Shanxi Medical University, Tongji Shanxi Hospital, Taiyuan, Shanxi, China
  • 2Department of Surgical Intensive Care Unit, Heping Campus of Shanxi Provincial People's Hospital, Taiyuan, Shanxi, China
  • 3Department of Clinical Laboratory, Heping Campus of Shanxi Provincial People's Hospital, Taiyuan, Shanxi, China
  • 4Department of Pharmacy, Shanxi Provincial People's Hospital, Taiyuan, Shanxi, China
  • 5Department of Neurosurgery, Shenzhen University Affiliated South China Hospital, Shenzhen, Guangdong, China
  • 6Department of Neurosurgery, Shanxi Provincial People's Hospital, Taiyuan, Shanxi, China

Introduction: Traumatic brain injury (TBI) is a major global health concern, leading to persistent neurological deficits and systemic complications. While animal studies have shown the association between TBI and the gut microbiota, human evidence, particularly in the early post-injury period, remains scarce.

Methods: In this study, we profiled the gut microbiota of TBI patients within 8 days of hospitalization using 16S rRNA sequencing, integrating clinical metadata and excluding individuals who had received antibiotics within the preceding month.

Results: Although alpha diversity remained similar between groups, beta diversity was significantly altered in TBI patients relative to healthy controls, accompanied by increased inter-individual variability. Differential abundance analysis revealed a depletion of short-chain fatty acid (SCFA)-producing taxa, e.g., Eubacterium, Agathobacter, and Faecalibacillus intestinalis, alongside an enrichment of dysbiosis-associated genera, including Escherichia and Shigella. Notably, Lactobacillus, a γ-aminobutyric acid-producing genus, was elevated in TBI patients, potentially reflecting a compensatory response to neural injury. Functional prediction suggested reduced SCFA biosynthetic capacity in TBI patients, whereas biofilm formation and several other fitness-related processes were enriched. Integrating gut microbiota with clinical and demographic variables, a machine learning model moderately predicted Glasgow Coma Scale scores, with age, inflammatory markers, and differentially abundant bacterial taxa as major contributors.

Discussion: Collectively, these findings indicate that TBI is associated with dysbiosis of the gut microbial and altered metabolic potential.

Introduction

In 2021, the global incidence of traumatic brain injury (TBI) was estimated at approximately 259 new cases per 100,000 population, translating to over 20 million new TBI cases worldwide (Perin et al., 2022). Individuals with a history of TBI frequently experience significant impairments in lower extremity mobility and functional activities (Schneider et al., 2021). Moreover, survivors of TBI face heightened risks of induced diseases, e.g., post-traumatic epilepsy, neurodegenerative conditions like Alzheimer's and Parkinson's disease, various psychiatric disorders, hormonal imbalances from pituitary dysfunction, and chronic sleep disorders (Schneider et al., 2021).

There are several metabolites in the gut microbiome that have been shown to influence TBI. Short-chain fatty acids (SCFAs), notably butyrate, acetate, and propionate, are key microbial metabolites that influence immune and neurological outcomes (Guo et al., 2022; Lan et al., 2024). In a TBI mouse model, SCFA supplementation restores microglial homeostasis, supports T-cell trafficking into the brain, and promotes neurogenesis, implying a protective role of SCFA-producers in the gut microbiota in TBI (Celorrio et al., 2025). Meanwhile, γ-aminobutyric acid (GABA), synthesized by gut bacteria, e.g., Bifidobacterium adolescentis, acts as a neuroactive postbiotic that modulates neural signaling (Mazzoli and Pessione, 2016; Duranti et al., 2020). Although direct data linking GABA to TBI outcomes remain limited, its broader role in gut-brain communication, including regulation of anxiety, stress, and central neurotransmission, suggests potential relevance in post-injury neurophysiology.

Metabolomic profiling in mouse models has further highlighted the impact of microbial metabolites on TBI (Du et al., 2021; Amaral et al., 2024). Interventions, i.e., probiotic administration and fecal microbiota transplantation (FMT), have shown potential to modulate the microbiota composition and improve neurological outcomes. For instance, in mouse models, FMT restored neurological function by influencing microbial metabolites and signaling pathways (Du et al., 2021; Wang et al., 2021). Similarly, probiotic treatment in another mouse study was associated with shifts in the microbiota community structure and increased markers of neuroprotection (Amaral et al., 2024). Collectively, these findings suggest that the gut microbiome may play a role in TBI recovery.

Consistent with the influence of microbial metabolites on TBI, studies have increasingly highlighted the changes of gut microbiota in the pathophysiology of TBI, with most research conducted using rodent models (Munley et al., 2023; Celorrio et al., 2024). In rat and mouse models, several studies have reported a reduction in alpha diversity of the gut microbiota following TBI, although some studies observed no significant differences (Munley et al., 2023; Celorrio et al., 2024). However, TBI consistently induces significant alterations in beta diversity across different studies, as shown by shifts in community composition, despite variability in the specific bacterial taxa affected. In a rat model of controlled cortical impact or fluid percussion injury, TBI causes acute dysbiosis of the gut microbiota characterized by reduced microbial diversity and notable decreases in butyrate-producing genera, e.g., Agathobacter, Faecalibacterium, and Eubacterium (Taraskina et al., 2022).

Research on the gut microbiota in human TBI is limited. An observational study in moderate-to-severe TBI patients admitted to a neurotrauma intensive care unit found high prevalence of Proteobacteria colonization as well as multidrug-resistant and colistin-resistant strains as early as 48 h post-injury, raising concerns about dysbiosis and infection risk (Mahajan et al., 2023). However, this study did not include taxonomic profiling of the gut microbiota. Evidence from other studies indicates that distinct and persistent alterations in the gut microbiota can persist for years after TBI (Urban et al., 2020; Pyles et al., 2024; Song et al., 2025).

Despite promising insights from animal models and limited human studies, a major barrier to translating gut microbiome research in TBI is the scarcity of human cohort studies, particularly in the early post-injury period. In this study, we characterized the gut microbiota profiles of 12 TBI patients and 13 healthy controls, none of whom had received antibiotics within 1 month prior to sample collection, by integrating high-resolution 16S rRNA sequencing with clinical metadata. Our analyses examined alterations in the diversity of the human gut microbiota associated with TBI, as well as potential functional changes in the gut microbiota of TBI patients within a few days after injury.

Methods

Recruitment of participants

Healthy volunteers without any known diseases and patients with TBI were recruited at Shanxi Provincial People's Hospital between February 2025 and June 2025, with one sample obtained per participant. For the TBI patients, 17 clinical variables were recorded, including antibiotic administration, Glasgow Coma Scale [GCS, an index to quantify severity (Bodien et al., 2021)], surgery, history of hypertension, diet, comorbidities, days since onset, and parameters indicative of inflammation, liver injury, coagulation activity, kidney function, and nutritional status. Except for surgical information, which was obtained during hospitalization, all other clinical data were collected on the same day as fecal sample collection.

Inclusion and exclusion criteria

For TBI patients, eligibility required: (1) a confirmed diagnosis of TBI; (2) age between 18 and 80 years; (3) provision of fecal samples during hospitalization; (4) complete clinical records; and (5) written informed consent obtained from the patient or their legal representative. Healthy volunteers were eligible if they: (1) were aged 18–80 years; (2) reported no known diseases; and (3) provided written informed consent. Exclusion criteria included: (1) spontaneous intracerebral hemorrhage; (2) use of antibiotics within 30 days prior to fecal sample collection; (3) presence of gastrointestinal diseases; (4) pregnancy or lactation; and (5) death within 48 h of admission to the intensive care unit. The 30-day antibiotic-free window was selected as the best achievable practice in clinical settings.

Sample collection and sequencing

Fresh stool samples were collected using a Microbiome Sample Collection Kit (Shanghai Personal Biotechnology Co., Ltd., China) and immediately placed in a −80 °C freezer. Samples were shipped on dry ice to Majorbio Co., Ltd. (China) for 16S rRNA sequencing. Briefly, DNA was extracted using the FastPure Stool DNA Isolation Kit (MJYH, Shanghai, China) and amplified with the primer pair 338F (5′-ACTCCTACGGGAGGCAGCAG-3′) and 806R (5′-GGACTACHVGGGTWTCTAAT-3′), targeting the V3-V4 region of the 16S rRNA gene. Sequencing libraries were prepared using the NEXTFLEX Rapid DNA-Seq Kit (Bioo Scientific, USA) and sequenced on the NextSeq 2000 platform with paired-end 300 bp reads.

Pre-treatment of 16S rRNA sequencing data

Raw sequencing data were processed through a quality control pipeline involving sequence trimming with fasterq-dump, paired-end merging using FLASH, and removal of human-derived reads via Bowtie2. The resulting reads were aligned to the Greengenes2 database (https://greengenes2.ucsd.edu/) for taxonomic classification (McDonald et al., 2024). Species-level annotations were used where available. To minimize noise, taxa with relative abundances of at least 0.1% (or 0.01%) present in a minimum of 5% (or 15%) of samples were retained, as previously described (Wang et al., 2023), resulting in a total of 280 species-level taxa. All samples were sequenced to a depth exceeding 34,000 reads, and therefore, none were excluded due to insufficient quality.

Diversity analysis

The 16S rRNA feature table was rarefied to a minimum sequencing depth of 34,000 reads per sample. Alpha diversity metrics, including Shannon index, observed taxa richness, and evenness, were calculated, with group comparisons conducted using two-sided Mann–Whitney U-tests. Beta diversity differences were visualized by Non-metric Multidimensional Scaling (NMDS) plots based on Bray-Curtis distance matrices and statistically tested using permutational multivariate analysis of variance (PERMANOVA) via the adonis2 function from the vegan R package (Dixon, 2003). The adonis2 function illustrates the effects of variables on the Bray-Curtis distance matrix or, in other words, the composition of the gut microbiota.

Differential abundance analysis

Taxonomic differential abundance was evaluated using LefSe analysis implemented through the run_lefse function within the microbiomeMarker R package, which identified taxa that were statistically different between groups and estimated their effect sizes.

Functional enrichment of KEGG orthologs

Based on the 16S rRNA profiles of the gut microbiota, PICRUSt2 (Douglas et al., 2020) was used to infer the abundance of KEGG orthologs (K numbers) by converting the 16S rRNA feature table into a KEGG ortholog feature table. Differential abundance analysis of K numbers between healthy and TBI patients was performed using DESeq2, a tool to analyze count-based data by modeling counts with negative binomial distributions, performing normalization, estimating dispersion, and identifying differentially expressed features between conditions (Love et al., 2014). Subsequently, K numbers were mapped to KEGG pathways according to the KEGG database annotations (Kanehisa et al., 2017). Enrichment of KEGG pathways was assessed through over-representation analysis utilizing the phyper function in R. The phyper function checked whether the number of genes found in a KEGG pathway was higher than what would be expected just by chance.

Modeling of GCS scores using machine learning

Three missing entries in the diet variable were labeled as “unknown.” One missing value in the white blood cell variable was imputed by the mice function in R. The imputation procedure is widely used in machine learning modeling which fills in missing data by creating possible values based on the patterns in the rest of the dataset. All available data, including clinical and demographic features as well as 16S rRNA-based gut microbiota profiles, were standardized using the scale function in R. The scale function standardizes data by centering the values around the mean and scaling them by the standard deviation, so that each feature is on the same scale for fair comparison. The GCS score served as the response variable, with all other factors included as predictors. Models were built with the train function from the caret package (Kuhn, 2008) and evaluated via two-fold cross-validation repeated 20 times. Variable importance within each model, measuring how much each variable contributes to a model's predictions, was assessed using the varImp function in R. Model performance was summarized by the root mean square error, R2, mean absolute error, and their corresponding standard deviations. Nine algorithms were explored: tree-based eXtreme Gradient Boosting, k-nearest neighbors, random forest, boosted generalized linear model, conditional inference random forest, support vector machine with linear kernel, linear eXtreme Gradient Boosting, generalized linear model, and simple linear regression.

Correlation analysis

Microbial abundances in the gut microbiota were normalized using the OTU table generated by the run_lefse function in R. Spearman's correlation, with p-values adjusted via the Benjamini–Hochberg method, was used to examine associations among demographic and clinical characteristics of TBI patients, as well as between these characteristics and abundances of the gut microbes.

Results

Profiles of the dataset

To investigate the role of the gut microbiota in TBI, fecal samples were collected from 13 healthy volunteers and 12 patients with TBI. No significant differences in age or sex were observed between the healthy and TBI groups (Table 1). The characteristics of these clinical variables are summarized in Table 2, and details are provided in Supplementary Datasets 1, 2. The average time of sample collection was 4.92 days after the initial visit (Table 2).

Table 1
www.frontiersin.org

Table 1. Demographic characteristics of participants in the cohort.

Table 2
www.frontiersin.org

Table 2. Clinical characteristics of TBI patients in the cohort.

Because antibiotics exert a strong influence on the composition of the human microbiota (Patangia et al., 2022), none of the participants in this study had received antibiotic treatment within 1 month prior to sample collection (Table 2). Therefore, potential effects attributable to antibiotic use were minimized.

Gut microbiota profiles from the 25 fecal samples were characterized using 16S rRNA gene sequencing. The minimum sequencing depth across all samples exceeded 34,000 reads (Supplementary Figure S1a). Alpha rarefaction analysis showed no significant difference in the number of observed bacterial taxa when rarefied to 34,000 vs. 32,000 reads (Supplementary Figure S1b). Consequently, all samples were rarefied to 34,000 reads for subsequent diversity analyses.

Association between diversity of the gut microbiota and TBI

In rodent models, multiple investigations have shown that TBI can lead to decreased alpha diversity in the gut microbiota, whereas other studies have found no notable changes (Munley et al., 2023; Celorrio et al., 2024). In the present human cohort, no significant differences in alpha diversity, i.e., Shannon index, observed taxa richness, and evenness, were observed between healthy controls and TBI patients (Figures 1ac). In contrast, consistent with findings from animal studies (Munley et al., 2023; Celorrio et al., 2024), the Adonis test revealed a significant association between beta diversity of the gut microbiota and TBI status (Figure 1d). Specifically, samples from the healthy group exhibited tighter clustering on the NMDS plot compared to those from the TBI group (Figure 1d), and the within-group similarity quantified by the Bray-Curtis distance was significantly lower among healthy controls (Figure 1e). These results imply that TBI could disrupt the homogeneity of the gut microbiota, resulting in increased inter-individual variability in the composition of the gut microbiota.

Figure 1
Violin plots and a scatter plot comparing Healthy and TBI groups. Plots a-c show Shannon index, observed taxa, and evenness with non-significant p-values. Plot d is an NMDS plot with a significant p-value (0.044) indicating group separation. Plot e shows within-group distance with a borderline significant p-value (0.045). Healthy and TBI groups are color-coded.

Figure 1. Gut microbiota diversity and compositional differences associated with TBI. Alpha diversity metrics, including Shannon index (a), observed taxa (b), and evenness (c), are shown for healthy and TBI groups, with statistical comparisons performed using two-sided Mann–Whitney U-tests. (d) Differences in overall microbiota composition are illustrated using a NMDS plot based on Bray-Curtis distances and evaluated with the Adonis test. (e) Within-group Bray-Curtis distances are compared between healthy and TBI individuals to assess within-group variability, using the two-sided Mann–Whitney U-test.

This increased dispersion in beta diversity among the TBI patients may be influenced by variability in treatment regimens. To investigate this, an Adonis test incorporating multiple clinical variables using a marginal model was conducted on the 12 TBI samples. The effect of each variable on the gut microbiota was tested marginally, i.e., controlling for all the other variables included in the model. As indicated by the significance “Pr(>F),” none of the tested factors, i.e., sex, age, surgery, hypertension history, diet, or comorbidity, had a significant effect on beta diversity (Supplementary Table S1). Thus, the tested variables may have had only a limited influence on the composition of the gut microbiota in this study, likely due to the small number of cases.

Differential abundance of gut microbes in TBI patients

Differential abundance analysis was performed using the Lefse method (Segata et al., 2011), revealing distinct differences in the abundance of bacterial taxa between TBI patients and healthy controls (Figure 2). As introduced above, SCFA is a key metabolite that maintains neurological and psychiatric health (Guo et al., 2022; Lan et al., 2024) and SCFA supplementation has a potential protective role mediated via the gut-brain axis in a TBI mouse model (Taraskina et al., 2022; Celorrio et al., 2025). Consistently, taxa enriched in healthy individuals included members of the Lachnospiraceae family, e.g., Eubacterium, Agathobacter, and Faecalibacillus intestinalis, which are recognized for their production of SCFAs like butyrate and propionate (Guo et al., 2022; Lan et al., 2024). Thus, the reduced abundance of SCFA-producing taxa in TBI patients may impair gut-brain axis homeostasis, thereby potentially aggravating neuroinflammation and neurological dysfunction associated with TBI.

Figure 2
Bubble chart displaying LDA scores for different microbial taxa, indicating enrichment in either healthy (turquoise) or TBI (purple) groups. Bubble size represents the -Log10(p-value), with sizes ranging from 1.5 to 3.0.

Figure 2. Differential gut microbiota composition associated with TBI. LefSe analysis was used to identify bacterial taxa differentially enriched between healthy controls and TBI patients. The LDA score reflects the effect size of each taxon, indicating its contribution to group differences. Dot size corresponds to statistical significance as computed by LefSe, and dot color indicates the group in which each taxon is enriched.

In contrast, Escherichia and Shigella species were among the most enriched taxa in the TBI group (Figure 2). These genera are generally considered opportunistic and are frequently associated with gut dysbiosis and intestinal inflammation (Baldelli et al., 2021; Li et al., 2022). Such microbial imbalances can exacerbate systemic inflammation and contribute to adverse neurological outcomes (Munley et al., 2023). Collectively, these findings indicate that TBI induces dysbiosis within the human gut microbiota.

Additionally, among taxa enriched in the TBI group, species in the genus of Lactobacillus are well-known producers of GABA, a key inhibitory neurotransmitter involved in neuroimmune regulation and gut-brain axis signaling (Figure 2) (Ma et al., 2019; Bin-Khattaf et al., 2022; Conn et al., 2024). Because neuroprotective effects of Lactobacillus have been demonstrated in mouse models of TBI through the modulation of inflammation, the increased abundance of Lactobacillus in TBI patients may reflect a compensatory microbiome response to neural injury and stress, potentially aimed at modulating gut-brain axis communication. However, contrasting findings show that TBI can also lead to a transient reduction in Lactobacillus, particularly L. gasseri, during the early post-injury phase in murine models (Treangen et al., 2018).

Together, these findings suggest that TBI resulted in gut dysbiosis characterized by the depletion of SCFA production-associated microbes and the overgrowth of dysbiosis-associated taxa, e.g., Escherichia. However, elevated Lactobacillus in TBI patients may serve as a compensatory mechanism. Whether shifts in the gut microbiota, particularly in taxa involved in SCFA and GABA pathways, impact the gut-brain axis during head injury and recovery remains to be elucidated.

Functional enrichment of KEGG orthologs and pathways in the gut microbiota associated with TBI

To further explore the metabolic potential of the gut microbiota, functional profiles were inferred from 16S rRNA sequencing data by predicting the abundance of KEGG Orthologs (K numbers) using the PICRUSt2 software (Douglas et al., 2020). Differential enrichment of specific K numbers was assessed using the DESeq2 software (Love et al., 2014). Subsequently, each enriched K number was mapped to corresponding KEGG pathways, and pathway-level enrichment was evaluated through over-representation analysis (Wieder et al., 2021).

There were 43 and 211 K numbers significantly enriched in the gut microbiota of healthy controls and TBI patients, respectively (Figure 3a and Supplementary Dataset 3). Consistent with the differential abundance analysis described above, several K numbers (annotated as matched gene names in Figure 3a) encode key enzymes catalyzing rate-limiting or critical steps in metabolic pathways leading to SCFA production, including fucA (Cheong et al., 2022), accA (Cheng et al., 2023), celS (Liao et al., 2024), iolN/iolS (Yoshida et al., 2008), and kdpgA (Reher et al., 2010), all enriched in the healthy group. For example, fucA, involved in fucose metabolism, has been shown to promote SCFA and 1,2-propanediol production (Cheong et al., 2022). The accA gene encodes the α-subunit of acetyl-CoA carboxylase, which catalyzes the conversion of acetyl-CoA to malonyl-CoA, the first committed step in fatty acid synthesis (Soriano et al., 2006). In microbial community studies, increased abundance of accA has been positively correlated with enhanced SCFA production (Cheng et al., 2023). Furthermore, several key genes directly responsible for SCFA production, such as pccA and mcmA, showed reduced but nonsignificant abundance changes in the TBI group, while others, e.g., ackA and pta (Frolova et al., 2022), showed no significant difference between groups. Overall, the enrichment of multiple upstream and regulatory SCFA-associated genes in healthy participants suggests a greater potential for SCFA biosynthesis capacity in this group compared to TBI patients.

Figure 3
Scatter plot and bar graph showing bacterial enrichment in healthy and TBI individuals. In the scatter plot, genes such as kdpG and ifhB are highlighted, showing fold change and FDR. The bar graph depicts KEGG pathways like peptidoglycan biosynthesis with significant enrichment in TBI.

Figure 3. Functional enrichment of KEGG orthologs (K numbers) and pathways in the gut microbiota associated with TBI. (a) K numbers, representing KEGG orthologs (functionally similar genes across species), were quantified using PICRUSt2. Differential abundance analysis of K numbers between healthy and TBI groups was performed using DESeq2. The volcano plot illustrates both the magnitude of change (log2 fold change) and statistical significance (FDR) for each K number. Selected significantly altered K numbers are labeled with their corresponding gene names. More details are shown in Supplementary Dataset 3. (b) Over-representation analysis was conducted to identify KEGG pathways significantly enriched in the gut microbiota of either healthy controls or TBI patients based on differentially abundant K numbers. More details are shown in Supplementary Dataset 4.

In contrast, many K numbers with matched genes associated with biofilm formation, a key virulence factor in the gut microbiota (Palandurkar and Kumar, 2023; Singh et al., 2025), were enriched in the TBI group (Figure 3a). These included genes for exopolysaccharide biosynthesis, such as pelB, pelC, and pelG (Amyx-Sherer and Reichhardt, 2024), as well as ompU, acfD, vpsE, and gspN (Teschler et al., 2022). Genes involved in flagellar motility, e.g., fliF and fliC, which facilitate initial surface attachment (Guttenplan and Kearns, 2013), and LPS core biosynthesis genes waaR and waaW, which influence biofilm stability (Aquilini et al., 2010), were also more abundant in the TBI group. Furthermore, global regulators ihfA and ihfB (Fan et al., 2023) and transcriptional regulators rob, marA, soxS, and ramA (Holden and Webber, 2020; Middlemiss et al., 2023), which modulate biofilm-related genes and stress responses, were enriched. Such enrichment implies a microbial community in TBI patients with an enhanced capacity for adhesion, matrix production, and persistence within host environments, potentially contributing to chronic inflammation or impaired recovery.

At the pathway level, none of the KEGG pathways remained significant after multiple-testing correction (Supplementary Dataset 4). However, nominal enrichment (p ≤ 0.05) revealed that three pathways, i.e., secondary metabolism, peptidoglycan biosynthesis, and other glycan degradation, were enriched in the TBI group (Figure 3b), consistent with a specific set of differentially abundant genes (FDR ≤ 0.05, Supplementary Dataset 3) identified at the KEGG Ortholog level. Genes involved in secondary metabolism, including thnK, thnN, and pltB, are responsible for the biosynthesis of naturally produced antimicrobial compounds such as thienamycin (Rodríguez et al., 2008) and pyoluteorin (Nowak-Thompson et al., 1999). Their enrichment implies an increased production of bioactive molecules that may enable microbes to outcompete others or adapt to altered host immune responses. The enrichment of the peptidoglycan biosynthesis pathway, marked by genes like murE, bacA, mraY, and pbpA, indicates elevated synthesis of peptidoglycan, a major component of bacterial cell walls especially in both Gram-positive and Gram-negative bacteria (Garde et al., 2021). Additionally, the upregulation of glycan degradation genes, including hex, manB, and afcA, indicates an adaptive microbial response to exploit host-derived glycans as alternative nutrient sources, potentially triggered by changes in nutrient availability in TBI patients. A shared feature in secondary metabolism and glycan degradation pathways is their upregulation in response to nutrient limitation or stress environment. This pattern may correspond with reduced food intake and altered nutrient availability in the gastrointestinal tract of TBI patients.

Predicting GCS scores from gut microbiota profiles and clinical-demographic factors

A GCS score is essential to identify the severity of TBI. To further understand the role of gut microbes in TBI, all the data involved in this study for the 12 TBI patients were applied to establish predictive modeling of GCS scores using a 10-times repeated two-fold cross-validation strategy. Nine different algorithms were applied and, based on the root mean square error values, the tree-based eXtreme Gradient Boosting (XGB tree) algorithm had the best performance (Figure 4a). More specifically, this XGB tree model achieved a root mean square error of 3.734, indicating that, on average, its predictions deviated from the observed values by about 3.7 units. The mean absolute error was 2.864, suggesting that the average magnitude of prediction errors, without considering direction, was approximately 2.86 units. The R2 value of 0.346 indicated that the model explained roughly 34.6% of the variance in the observed GCS scores, implying moderate predictive performance with substantial unexplained variability. Therefore, it appears that the primary factors accounting for the remaining 65% of the variance in GCS scores were not captured in the dataset, potentially including the degree of injury at the time of hospitalization and biomarkers that directly reflect the status of the nervous, cardiovascular, respiratory, immune, and other related systems.

Figure 4
Data visualization comprising various analyses on clinical variables and their impact on the Glasgow Coma Scale (GCS). Panel (a) presents a table ranking model performance metrics like RMSE and Rsquared for different statistical models. Panel (b) lists overall importance scores for clinical variables such as age and D-dimer. Panels (c) to (k) depict scatter plots detailing the relationship between GCS and different variables, including age, diet, and serum markers like D-dimer, white blood cells, and C-reactive protein, with trend lines and statistical values shown.

Figure 4. Predicting GCS scores from tut microbiota profiles and clinical-demographic factors. (a) Nine machine learning algorithms were employed to predict GCS scores, a clinical measure of consciousness level following brain injury. Model performance was assessed through two-fold cross-validation repeated 20 times. Performance metrics included root mean square error (RMSE), R-squared (R2), mean absolute error (MAE), and their respective standard deviations (RMSESD, RsquaredSD, MAESD). Algorithms used include: tree-based eXtreme Gradient Boosting (XGB_tree), k-nearest neighbors (KNN), random forest (RF), boosted generalized linear model (GLMB), conditional inference random forest (CFOREST), support vector machine with linear kernel (SVM), linear eXtreme Gradient Boosting (XGB_linear), generalized linear model (GLM), and simple linear model (LM). (b) The variable importance ranking from the XGB_tree model illustrates the relative influence of each feature on predictive performance. (c–k) Correlations among clinical or demographic variables were quantified by Spearman's rank correlation (see details in Supplementary Dataset 6). p-Values were adjusted for multiple testing using the Benjamini–Hochberg method to control the false discovery rate (FDR). Associations between GCS scores and each clinical or demographic variable included in the XGB_linear model are shown.

It is worth noting that most decision tree-based algorithms, e.g., XGB tree and random forest, demonstrated better performance than linear models, like linear eXtreme Gradient Boosting, generalized linear model, and simple linear regression (Figure 4a). Hence, the relationships between the predictors and GCS scores were likely nonlinear and tree-based methods were inherently better suited to capture than linear approaches.

The contribution of each factor in the XGB tree model was evaluated using feature importance, which quantifies how much each predictor improves model performance, measured by gain, cover, and frequency, across all trees in the ensemble (Figure 4b). Consistent with the differential abundance analysis (Figure 2), all the bacterial taxa contributing to the XGB tree model are the gut microbes that had significant differential abundances between healthy controls and TBI patients. Particularly, Clostridium innocuum and the genus of Escherichia had increased abundance in TBI patients, but all other taxa were enriched in healthy participants (Figures 2, 4b). Consistent with these observations, a rat model demonstrated that head injury led to a dramatic increase in Escherichia coli populations (Waligora-Dupriet et al., 2018), and C. innocuum has been associated with gut inflammation and disease (Cherny et al., 2021).

Spearman's correlation analysis was used to evaluate the relationships between clinical-demographic factors and GCS scores (Supplementary Dataset 5). Consistent with previous studies showing that older patients tend to experience more severe TBI (Hukkelhoven et al., 2003), age had the strongest influence on GCS scores among the tested factors and was negatively, though not significantly, correlated with GCS (Figure 4c). Other clinical variables included in the XGB tree model, i.e., D-dimer, total parenteral nutrition, days since injury onset, white blood cell count, neutrophil percentage, procalcitonin, and C-reactive protein, were also associated with greater TBI severity (Figures 4d–j). Conversely, albumin levels were higher in patients with better GCS scores (Figure 4k), although many of these correlations did not reach statistical significance. These clinical markers reflect inflammatory and immune responses, as well as metabolic and nutritional status. Elevated D-dimer levels have been linked to increased coagulation and potential microvascular injury (Weitz et al., 2017), while higher procalcitonin and C-reactive protein, together with lower albumin levels, indicate systemic inflammation and possible infection (Karaca et al., 2023). Alterations in white blood cell counts and neutrophil percentages further imply immune activation (Abramson and Melton, 2000; Buonacera et al., 2022). Moreover, the need for total parenteral nutrition and longer time since injury may reflect greater physiological stress and impaired recovery. Collectively, these findings highlight how TBI induces widespread systemic effects, including alterations in the gut microbiota, that extend beyond the central nervous system.

Additionally, Spearman's correlations between clinical-demographic characteristics of the TBI patients and the abundances of the gut microbial taxa were evaluated (Supplementary Dataset 6). The only significant outcome after multiple correction of the p-values was the positive correlation between the species of Wujia chipingensis and albumin (Supplementary Figure S2). Aligned with this observation, W. chipingensis was enriched in healthy participants (Figure 2) and albumin was associated with lower levels of severity of TBI (Figure 4k). Additionally, W. chipingensis has been associated with lower levels of C-reactive protein levels, indicating its potential involvement in modulating inflammatory responses (Zhu et al., 2025). Consistent with the neuroprotective effects of Lactobacillus in mouse models (Ma et al., 2019; Bin-Khattaf et al., 2022; Conn et al., 2024), Lactobacillus was significantly negatively correlated with C-reactive protein, a well-established inflammatory biomarker, prior to p-value adjustment (Supplementary Dataset 6).

Discussion

In this study, our findings showed that TBI in humans was associated with alterations in gut microbial composition and potential biological function, characterized by depletion of SCFA-producing taxa and enrichment of taxa commonly associated with dysbiosis. SCFAs are well-characterized mediators of gut barrier integrity, immune regulation, and gut-brain signaling (Guo et al., 2022; Lan et al., 2024; Celorrio et al., 2025). The reduced SCFA biosynthetic potential in TBI patients could therefore contribute to increased gut permeability, systemic inflammation, and exacerbation of neuroinflammatory cascades after TBI. Functionally, the enrichment of KEGG orthologs linked to biofilm formation, secondary metabolism, and glycan degradation in the TBI group further implies that post-injury microbial communities may be more adhesive, stress-tolerant, and capable of exploiting host-derived glycans under nutrient-limited conditions, features that can promote persistence and inflammation in injured hosts.

GABA, an inhibitory neurotransmitter involved in neuroimmune regulation and gut-brain communication, plays a protective role in TBI (Mazzoli and Pessione, 2016; Duranti et al., 2020). The observed enrichment of Lactobacillus in TBI patients suggests a potential increase in gut GABA, which could represent a compensatory response to injury-induced dysbiosis or a transient bloom influenced by diet, hospitalization, or other treatments (Ma et al., 2019; Bin-Khattaf et al., 2022; Conn et al., 2024). Supporting this idea, previous research has reported an elevated GABA-to-glutamate ratio in the gut microbiota as a metabolic signature in mild autism spectrum disorder, linked to overrepresentation of Escherichia (Wang et al., 2025) and our study found Escherichia to be enriched in the TBI group. Furthermore, Lactobacillus is enriched in schizophrenia, accompanied by changes in the glutamate-glutamine-GABA cycle (Nuncio-Mora et al., 2023). Together, these findings highlight the need for further investigation into GABA production by the gut microbes and their potential role in TBI recovery.

Linking microbial signatures and clinical variables to the severity of TBI, our predictive modeling (XGB tree) explained a modest fraction of GCS variance (R2 ≈ 0.35), indicating that gut microbiota and the clinical-demographic factors measured here capture part but not the majority of determinants of consciousness after injury. Several clinical markers associated with inflammation, coagulation, and nutritional status, e.g., D-dimer, procalcitonin, CRP, albumin, correlated with GCS, consistent with a model in which systemic inflammation and metabolic derangement after TBI influence gut microbial composition. Notably, albumin was positively associated with W. chipingensis, a taxon enriched in healthy participants, and with better GCS scores, implying that preserved nutritional status may mitigate dysbiosis of the gut microbiota.

This study has several important limitations. The sample size is small and includes only a single time point, which limits the statistical power of our analyses. Consequently, the detection of correlations between bacterial taxa and clinical variables and detecting pathway-level changes were constrained, particularly after correcting for multiple comparisons. Detection of alpha diversity associated with TBI could potentially be limited by the small sample size. Additionally, despite efforts to exclude recent antibiotic exposure, the cohort remains heterogeneous in terms of injury timing, clinical management (including nutrition and surgery), disease severity, and potentially probiotic and dietary factors, which may limit the generalizability of these findings to broader populations.

Taxonomic resolution from 16S rRNA sequencing is limited at the strain level and precludes definitive assignments of pathogenic vs. beneficial roles for taxa within polyphyletic genera, e.g., Clostridium. Achieving strain-level resolution requires metagenomic sequencing, and functional inferences based on predicted KEGG orthologs likewise need validation through metabolomic measurements or other quantification methods (e.g., direct quantification of SCFAs, GABA, and biofilm-related products).

It was challenging to collect fecal samples from patients with severe TBI who had not received antibiotics, due to several factors, e.g., constipation, critical condition, fasting, reduced stool output associated with total parenteral nutrition, and prophylactic antibiotic use in patients with open or penetrating TBI. Consequently, seven of the 12 TBI patients included in this study presented with only mild impairment of consciousness (GCS score ≥ 13, Table 2), which may introduce bias in sample collection with respect to disease severity.

In conclusion, our results indicate that human TBI is associated with a shift in gut microbial composition away from SCFA-producing, homeostatic taxa toward communities enriched in opportunistic, adhesive, and potentially pro-inflammatory organisms. Elucidating causal relationships and therapeutic opportunities will require larger, higher-resolution, and longitudinal investigations.

Data availability statement

Raw 16S rRNA sequencing data are available in the Genome Sequence Archive database belonging to China National Center for Bioinformation (accession code: PRJCA045145).

Ethics statement

The study was approved by the Ethics Committee of Shanxi Provincial People's Hospital (Approval No. 2025-198). The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.

Author contributions

HW: Writing – original draft, Software, Visualization, Formal analysis, Methodology, Writing – review & editing. YH: Resources, Writing – review & editing, Data curation. WF: Data curation, Resources, Writing – review & editing. JG: Data curation, Writing – review & editing, Resources. XW: Formal analysis, Visualization, Writing – review & editing, Software. LS: Writing – review & editing, Conceptualization, Supervision. GC: Supervision, Conceptualization, Project administration, Funding acquisition, Writing – review & editing, Writing – original draft.

Funding

The author(s) declared that financial support was received for this work and/or its publication. This work was funded by the Special Research Project on Science and Technology Strategy of Shanxi Province (GC, ID: 202404030401186).

Acknowledgments

We would like to thank all the participants and their families for their support during the case collection process. We thank Dr. Jia Lu of the Department of Neurosurgery, Shanxi Provincial People's Hospital, for his creative suggestions on the work.

Conflict of interest

The author(s) declared that this work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declared that generative AI was not used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

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.

Supplementary material

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

References

Abramson, N., and Melton, B. (2000). Leukocytosis: basics of clinical assessment. Am. Fam. Physician 62, 2053–2060.

PubMed Abstract | Google Scholar

Amaral, W. Z., Kokroko, N., Treangen, T. J., Villapol, S., and Gomez-Pinilla, F. (2024). Probiotic therapy modulates the brain-gut-liver microbiota axis in a mouse model of traumatic brain injury. Biochim. Biophys. Acta BBA-Mol. Basis Dis. 1870:167483. doi: 10.1016/j.bbadis.2024.167483

PubMed Abstract | Crossref Full Text | Google Scholar

Amyx-Sherer, K., and Reichhardt, C. (2024). Challenges and opportunities in elucidating the structures of biofilm exopolysaccharides: a case study of the Pseudomonas aeruginosa exopolysaccharide called Pel. Magn. Reson. Chem. 62, 361–369. doi: 10.1002/mrc.5405

PubMed Abstract | Crossref Full Text | Google Scholar

Aquilini, E., Azevedo, J., Jimenez, N., Bouamama, L., Tomás, J. M., Regué, M., et al. (2010). Functional identification of the Proteus mirabilis core lipopolysaccharide biosynthesis genes. J. Bacteriol. 192, 4413–4424. doi: 10.1128/JB.00494-10

PubMed Abstract | Crossref Full Text | Google Scholar

Baldelli, V., Scaldaferri, F., Putignani, L., and Del Chierico, F. (2021). The role of enterobacteriaceae in gut microbiota dysbiosis in inflammatory bowel diseases. Microorganisms 9:697. doi: 10.3390/microorganisms9040697

PubMed Abstract | Crossref Full Text | Google Scholar

Bin-Khattaf, R. M., Alonazi, M. A., Al-Dbass, A. M., Almnaizel, A. T., Aloudah, H. S., Soliman, D. A., et al. (2022). Probiotic ameliorating effects of altered GABA/glutamate signaling in a rodent model of autism. Metabolites 12:720. doi: 10.3390/metabo12080720

PubMed Abstract | Crossref Full Text | Google Scholar

Bodien, Y. G., Barra, A., Temkin, N. R., Barber, J., Foreman, B., Vassar, M., et al. (2021). Diagnosing level of consciousness: the limits of the glasgow coma scale total score. J. Neurotrauma 38, 3295–3305. doi: 10.1089/neu.2021.0199

PubMed Abstract | Crossref Full Text | Google Scholar

Buonacera, A., Stancanelli, B., Colaci, M., and Malatino, L. (2022). Neutrophil to lymphocyte ratio: an emerging marker of the relationships between the immune system and diseases. Int. J. Mol. Sci. 23:3636. doi: 10.3390/ijms23073636

PubMed Abstract | Crossref Full Text | Google Scholar

Celorrio, M., Shumilov, K., and Friess, S. H. (2024). Gut microbial regulation of innate and adaptive immunity after traumatic brain injury. Neural Regen. Res. 19, 272–276. doi: 10.4103/1673-5374.379014

PubMed Abstract | Crossref Full Text | Google Scholar

Celorrio, M., Shumilov, K., Ni, A., Ayerra, L., Self, W. K., de Francisca, N. V., et al. (2025). Short-chain fatty acids are a key mediator of gut microbial regulation of T cell trafficking and differentiation after traumatic brain injury. Exp. Neurol. 392:115349. doi: 10.1016/j.expneurol.2025.115349

PubMed Abstract | Crossref Full Text | Google Scholar

Cheng, B., Zhang, D., Lin, Q., Xi Ma, Ma, S., Zan, J. F., et al. (2023). Short-chain fatty acid production and phosphorous recovery from waste activated sludge via anaerobic fermentation: a comparison of in-situ and ex-situ thiosulfate-assisted Fe(2+)/persulfate pretreatment. Sci. Total Environ. 871:162172. doi: 10.1016/j.scitotenv.2023.162172

PubMed Abstract | Crossref Full Text | Google Scholar

Cheong, Y. E., Kim, J., Jin, Y-. S., and Kim, K. H. (2022). Elucidation of the fucose metabolism of probiotic Lactobacillus rhamnosus GG by metabolomic and flux balance analyses. J. Biotechnol. 360, 110–116. doi: 10.1016/j.jbiotec.2022.11.002

PubMed Abstract | Crossref Full Text | Google Scholar

Cherny, K. E., Muscat, E. B., Reyna, M. E., and Kociolek, L. K. (2021). Clostridium innocuum: microbiological and clinical characteristics of a potential emerging pathogen. Anaerobe 71:102418. doi: 10.1016/j.anaerobe.2021.102418

PubMed Abstract | Crossref Full Text | Google Scholar

Conn, K. A., Borsom, E. M., and Cope, E. K. (2024). Implications of microbe-derived ɤ-aminobutyric acid (GABA) in gut and brain barrier integrity and GABAergic signaling in Alzheimer's disease. Gut Microbes 16:2371950. doi: 10.1080/19490976.2024.2371950

PubMed Abstract | Crossref Full Text | Google Scholar

Dixon, P. (2003). VEGAN, a package of R functions for community ecology. J. Veg. Sci. 14, 927–930. doi: 10.1111/j.1654-1103.2003.tb02228.x

Crossref Full Text | Google Scholar

Douglas, G. M., Maffei, V. J., Zaneveld, J. R., Yurgel, S. N., Brown, J. R., Taylor, C. M., et al. (2020). PICRUSt2 for prediction of metagenome functions. Nat. Biotechnol. 38, 685–688. doi: 10.1038/s41587-020-0548-6

PubMed Abstract | Crossref Full Text | Google Scholar

Du, D., Tang, W., Zhou, C., Sun, X., Wei, Z., Zhong, J., et al. (2021). Fecal microbiota transplantation is a promising method to restore gut microbiota dysbiosis and relieve neurological deficits after traumatic brain injury. Oxid. Med. Cell. Longev. 2021:5816837. doi: 10.1155/2021/5816837

PubMed Abstract | Crossref Full Text | Google Scholar

Duranti, S., Ruiz, L., Lugli, G. A., Tames, H., Milani, C., Mancabelli, L., et al. (2020). Bifidobacterium adolescentis as a key member of the human gut microbiota in the production of GABA. Sci. Rep. 10:14112. doi: 10.1038/s41598-020-70986-z

PubMed Abstract | Crossref Full Text | Google Scholar

Fan, Z., Fu, T., Li, Z., Du, B., Cui, X., Zhang, R., et al. (2023). The role of integration host factor in biofilm and virulence of high-alcohol-producing Klebsiella pneumoniae. Microbiol. Spectr. 11:e0117023. doi: 10.1128/spectrum.01170-23

PubMed Abstract | Crossref Full Text | Google Scholar

Frolova, M. S., Suvorova, I. A., Iablokov, S. N., Petrov, S. N., and Rodionov, D. A. (2022). Genomic reconstruction of short-chain fatty acid production by the human gut microbiota. Front. Mol. Biosci. 9:949563. doi: 10.3389/fmolb.2022.949563

PubMed Abstract | Crossref Full Text | Google Scholar

Garde, S., Chodisetti, P. K., and Reddy, M. (2021). Peptidoglycan: structure, synthesis, and regulation. EcoSal Plus 9. doi: 10.1128/ecosalplus.ESP-0010-2020

PubMed Abstract | Crossref Full Text | Google Scholar

Guo, C., Huo, Y-. J., Li, Y., Han, Y., and Zhou, D. (2022). Gut-brain axis: focus on gut metabolites short-chain fatty acids. World J. Clin. Cases 10:1754. doi: 10.12998/wjcc.v10.i6.1754

PubMed Abstract | Crossref Full Text | Google Scholar

Guttenplan, S. B., and Kearns, D. B. (2013). Regulation of flagellar motility during biofilm formation. FEMS Microbiol. Rev. 37, 849–871. doi: 10.1111/1574-6976.12018

PubMed Abstract | Crossref Full Text | Google Scholar

Holden, E. R., and Webber, M. A. (2020). MarA, RamA, and SoxS as mediators of the stress response: survival at a cost. Front. Microbiol. 11:828. doi: 10.3389/fmicb.2020.00828

PubMed Abstract | Crossref Full Text | Google Scholar

Hukkelhoven, C. W. P. M., Steyerberg, E. W., Rampen, A. J. J., Farace, E., Habbema, J. D. F., Marshall, L. F., et al. (2003). Patient age and outcome following severe traumatic brain injury: an analysis of 5600 patients. J. Neurosurg. 99, 666–673. doi: 10.3171/jns.2003.99.4.0666

PubMed Abstract | Crossref Full Text | Google Scholar

Kanehisa, M., Furumichi, M., Tanabe, M., Sato, Y., and Morishima, K. (2017). KEGG: new perspectives on genomes, pathways, diseases and drugs. Nucleic Acids Res. 45, D353–D361. doi: 10.1093/nar/gkw1092

PubMed Abstract | Crossref Full Text | Google Scholar

Karaca, B., Esin, F., Tiryaki, M. M., Akkan, G., and Kiris, T. (2023). Combining C-reactive protein, procalcitonin, and serum albumin to predict long-term mortality in patients with infective endocarditis. J. Int. Med. Res. 51:3000605231208910. doi: 10.1177/03000605231208910

PubMed Abstract | Crossref Full Text | Google Scholar

Kuhn, M. (2008). Building predictive models in R using the caret package. J. Stat. Softw. 28, 1–26. doi: 10.18637/jss.v028.i05

Crossref Full Text | Google Scholar

Lan, Z., Tang, X., Lu, M., Hu, Z., and Tang, Z. (2024). The role of short-chain fatty acids in central nervous system diseases: a bibliometric and visualized analysis with future directions. Heliyon 10:e26377. doi: 10.1016/j.heliyon.2024.e26377

PubMed Abstract | Crossref Full Text | Google Scholar

Li, S., Guo, J., Liu, R., Zhang, F., Wen, S., Liu, Y., et al. (2022). Predominance of Escherichia-Shigella in gut microbiome and its potential correlation with elevated level of plasma tumor necrosis factor alpha in patients with tuberculous meningitis. Microbiol. Spectr. 10:e01926-22. doi: 10.1128/spectrum.01926-22

PubMed Abstract | Crossref Full Text | Google Scholar

Liao, Y., Wu, S., Zhou, G., Mei, S., Yang, Z., Li, S., et al. (2024). Cellulolytic Bacillus cereus produces a variety of short-chain fatty acids and has potential as a probiotic. Microbiol. Spectr. 12:e0326723. doi: 10.1128/spectrum.03267-23

PubMed Abstract | Crossref Full Text | Google Scholar

Love, M. I., Huber, W., and Anders, S. (2014). Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 15:550. doi: 10.1186/s13059-014-0550-8

PubMed Abstract | Crossref Full Text | Google Scholar

Ma, Y., Liu, T., Fu, J., Fu, S., Hu, C., Sun, B., et al. (2019). Lactobacillus acidophilus exerts neuroprotective effects in mice with traumatic brain injury. J. Nutr. 149, 1543–1552. doi: 10.1093/jn/nxz105

PubMed Abstract | Crossref Full Text | Google Scholar

Mahajan, C., Khurana, S., Kapoor, I., Sokhal, S., Kumar, S., Prabhakar, H., et al. (2023). Characteristics of gut microbiome after traumatic brain injury. J. Neurosurg. Anesthesiol. 35, 86–90. doi: 10.1097/ANA.0000000000000789

PubMed Abstract | Crossref Full Text | Google Scholar

Mazzoli, R., and Pessione, E. (2016). The neuro-endocrinological role of microbial glutamate and GABA signaling. Front. Microbiol. 7:1934. doi: 10.3389/fmicb.2016.01934

PubMed Abstract | Crossref Full Text | Google Scholar

McDonald, D., Jiang, Y., Balaban, M., Cantrell, K., Zhu, Q., Gonzalez, A., et al. (2024). Greengenes2 unifies microbial data in a single reference tree. Nat. Biotechnol. 42, 715–718. doi: 10.1038/s41587-023-01845-1

PubMed Abstract | Crossref Full Text | Google Scholar

Middlemiss, A. D., Haycocks, J. R. J., Stringer, A. M., Piddock, L. J. V., Wade, J. T., Grainger, D. C., et al. (2023). Mapping direct and indirect MarA/SoxS/Rob/RamA regulons in Salmonella typhimurium reveals repression of csgD and biofilm formation. Microbiol. Read. Engl. 169:001330. doi: 10.1099/mic.0.001330

PubMed Abstract | Crossref Full Text | Google Scholar

Munley, J. A., Kirkpatrick, S. L., Gillies, G. S., Bible, L. E., Efron, P. A., Nagpal, R., et al. (2023). The intestinal microbiome after traumatic injury. Microorganisms 11:1990. doi: 10.3390/microorganisms11081990

PubMed Abstract | Crossref Full Text | Google Scholar

Nowak-Thompson, B., Chaney, N., Wing, J. S., Gould, S. J., and Loper, J. E. (1999). Characterization of the pyoluteorin biosynthetic gene cluster of Pseudomonas fluorescens Pf-5. J. Bacteriol. 181, 2166–2174. doi: 10.1128/JB.181.7.2166-2174.1999

PubMed Abstract | Crossref Full Text | Google Scholar

Nuncio-Mora, L., Lanzagorta, N., Nicolini, H., Sarmiento, E., Ortiz, G., Sosa, F., et al. (2023). The role of the microbiome in first episode of psychosis. Biomedicines 11:1770. doi: 10.3390/biomedicines11061770

PubMed Abstract | Crossref Full Text | Google Scholar

Palandurkar, G. S., and Kumar, S. (2023). Biofilm's impact on inflammatory bowel diseases. Cureus 15:e45510. doi: 10.7759/cureus.45510

PubMed Abstract | Crossref Full Text | Google Scholar

Patangia, D. V., Anthony Ryan, C., Dempsey, E., Paul Ross, R., and Stanton, C. (2022). Impact of antibiotics on the human microbiome and consequences for host health. Microbiologyopen 11:e1260. doi: 10.1002/mbo3.1260

PubMed Abstract | Crossref Full Text | Google Scholar

Perin, J., Mulick, A., Yeung, D., Villavicencio, F., Lopez, G., Strong, K. L., et al. (2022). Global, regional, and national causes of under-5 mortality in 2000–19: an updated systematic analysis with implications for the sustainable development goals. Lancet Child Adolesc. Health 6, 106–115. doi: 10.1016/S2352-4642(21)00311-4

PubMed Abstract | Crossref Full Text | Google Scholar

Pyles, R. B., Miller, A. L., Urban, R. J., Sheffield-Moore, M., Wright, T. J., and Maxwell, C. A.. (2024). The altered TBI fecal microbiome is stable and functionally distinct. Front. Mol. Neurosci. 17:1341808. doi: 10.3389/fnmol.2024.1341808

PubMed Abstract | Crossref Full Text | Google Scholar

Reher, M., Fuhrer, T., Bott, M., and Schönheit, P. (2010). The nonphosphorylative Entner-Doudoroff pathway in the thermoacidophilic euryarchaeon Picrophilus torridus involves a novel 2-keto-3-deoxygluconate- specific aldolase. J. Bacteriol. 192, 964–974. doi: 10.1128/JB.01281-09

PubMed Abstract | Crossref Full Text | Google Scholar

Rodríguez, M., Núñez, L. E., Braña, A. F., Méndez, C., Salas, J. A., Blanco, G., et al. (2008). Identification of transcriptional activators for thienamycin and cephamycin C biosynthetic genes within the thienamycin gene cluster from Streptomyces cattleya. Mol. Microbiol. 69, 633–645. doi: 10.1111/j.1365-2958.2008.06312.x

PubMed Abstract | Crossref Full Text | Google Scholar

Schneider, A. L., Wang, D., Gottesman, R. F., and Selvin, E. (2021). Prevalence of disability associated with head injury with loss of consciousness in adults in the United States: a population-based study. Neurology 97, e124–e135. doi: 10.1212/WNL.0000000000012148

PubMed Abstract | Crossref Full Text | Google Scholar

Segata, N., Izard, J., Waldron, L., Gevers, D., Miropolsky, L., Garrett, W. S., et al. (2011). Metagenomic biomarker discovery and explanation. Genome Biol. 12:R60. doi: 10.1186/gb-2011-12-6-r60

PubMed Abstract | Crossref Full Text | Google Scholar

Singh, A. A., Khan, F., and Song, M. (2025). Biofilm-associated amyloid proteins linked with the progression of neurodegenerative diseases. Int. J. Mol. Sci. 26:2695. doi: 10.3390/ijms26062695

PubMed Abstract | Crossref Full Text | Google Scholar

Song, B., Qiu, Y., Wang, Z., Tao, Y., Wang, M., Duan, A., et al. (2025). The causal relationship between gut microbiomes, inflammatory mediators, and traumatic brain injury in Europeans: evidence from genetic correlation and functional mapping annotation analyses. Biomedicines 13:753. doi: 10.3390/biomedicines13030753

PubMed Abstract | Crossref Full Text | Google Scholar

Soriano, A., Radice, A. D., Herbitter, A. H., Langsdorf, E. F., Stafford, J. M., Chan, S., et al. (2006). Escherichia coli acetyl-coenzyme a carboxylase: characterization and development of a high-throughput assay. Anal. Biochem. 349, 268–276. doi: 10.1016/j.ab.2005.10.044

PubMed Abstract | Crossref Full Text | Google Scholar

Taraskina, A., Ignatyeva, O., Lisovaya, D., Ivanov, M., Ivanova, L., Golovicheva, V., et al. (2022). Effects of traumatic brain injury on the gut microbiota composition and serum amino acid profile in rats. Cells 11:1409. doi: 10.3390/cells11091409

PubMed Abstract | Crossref Full Text | Google Scholar

Teschler, J. K., Nadell, C. D., Drescher, K., and Yildiz, F. H. (2022). Mechanisms underlying vibrio cholerae Biofilm formation and dispersion. Annu. Rev. Microbiol. 76, 503–532. doi: 10.1146/annurev-micro-111021-053553

PubMed Abstract | Crossref Full Text | Google Scholar

Treangen, T. J., Wagner, J., Burns, M. P., and Villapol, S. (2018). Traumatic brain onjury in mice induces acute bacterial dysbiosis within the fecal microbiome. Front. Immunol. 9:2757. doi: 10.3389/fimmu.2018.02757

Crossref Full Text | Google Scholar

Urban, R. J., Pyles, R. B., Stewart, C. J., Ajami, N., Randolph, K. M., Durham, W. J., et al. (2020). Altered fecal microbiome years after traumatic brain injury. J. Neurotrauma 37, 1037–1051. doi: 10.1089/neu.2019.6688

PubMed Abstract | Crossref Full Text | Google Scholar

Waligora-Dupriet, A-. J., Lafleur, S., Charrueau, C., Choisy, C., Cynober, L., Butel, M-. J., et al. (2018). Head injury profoundly affects gut microbiota homeostasis: results of a pilot study. Nutrition. 45, 104–107. doi: 10.1016/j.nut.2017.06.026

PubMed Abstract | Crossref Full Text | Google Scholar

Wang, D., Jiang, Y., Jiang, J., Pan, Y., Yang, Y., Fang, X., et al. (2025). Gut microbial GABA imbalance emerges as a metabolic signature in mild autism spectrum disorder linked to overrepresented Escherichia. Cell Rep. Med. 6:101919. doi: 10.1016/j.xcrm.2024.101919

PubMed Abstract | Crossref Full Text | Google Scholar

Wang, S., Zhu, K., Hou, X., and Hou, L. (2021). The association of traumatic brain injury, gut microbiota and the corresponding metabolites in mice. Brain Res. 1762:147450. doi: 10.1016/j.brainres.2021.147450

PubMed Abstract | Crossref Full Text | Google Scholar

Wang, Y., Thompson, K. N., Yan, Y., Short, M. I., Zhang, Y., Franzosa, E. A., et al. (2023). RNA-based amplicon sequencing is ineffective in measuring metabolic activity in environmental microbial communities. Microbiome 11:131. doi: 10.1186/s40168-022-01449-y

PubMed Abstract | Crossref Full Text | Google Scholar

Weitz, J. I., Fredenburgh, J. C., and Eikelboom, J. W. (2017). A test in context: D-Dimer. J. Am. Coll. Cardiol. 70, 2411–2420. doi: 10.1016/j.jacc.2017.09.024

PubMed Abstract | Crossref Full Text | Google Scholar

Wieder, C., Frainay, C., Poupin, N., Rodríguez-Mier, P., Vinson, F., Cooke, J., et al. (2021). Pathway analysis in metabolomics: recommendations for the use of over-representation analysis. PLoS Comput. Biol. 17:e1009105. doi: 10.1371/journal.pcbi.1009105

PubMed Abstract | Crossref Full Text | Google Scholar

Yoshida, K., Yamaguchi, M., Morinaga, T., Kinehara, M., Ikeuchi, M., Ashida, H., et al. (2008). myo-Inositol catabolism in Bacillus subtilis. J. Biol. Chem. 283, 10415–10424. doi: 10.1074/jbc.M708043200

PubMed Abstract | Crossref Full Text | Google Scholar

Zhu, B., Liang, L., Chen, S., Li, H., Huang, Y., Wang, W., et al. (2025). Multi-kingdom microbial changes and their associations with the clinical characteristics in schizophrenia patients. Transl. Psychiatry 15:228. doi: 10.1038/s41398-025-03449-6

PubMed Abstract | Crossref Full Text | Google Scholar

Keywords: traumatic brain injury, gut microbiota, short-chain fatty acid, functional prediction, biofilm, machine learning modeling

Citation: Wu H, Huo Y, Fang W, Guo J, Wang X, Su L and Cheng G (2026) Gut microbiota and traumatic brain injury: insights from an antibiotic-free cohort. Front. Microbiol. 16:1697206. doi: 10.3389/fmicb.2025.1697206

Received: 02 September 2025; Revised: 06 November 2025;
Accepted: 10 November 2025; Published: 02 January 2026.

Edited by:

Bhagavathi Sundaram Sivamaruthi, Chiang Mai University, Thailand

Reviewed by:

Li Shan Chiu, Ear Science Institute Australia, Australia
Yajun Bi, Wuhan Pulmonary Hospital, China
Elias Barbosa Da Silva Junior, Federal University of Rio de Janeiro, Brazil

Copyright © 2026 Wu, Huo, Fang, Guo, Wang, Su and Cheng. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Gang Cheng, Z2FuZ2NoZW4wODA4MDZAZ21haWwuY29t; Li Su, c3VsaTI5N0AxNjMuY29t

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.