Abstract
Differentiation of levels of consciousness in patients with prolonged disorders of consciousness (DOC) remains a major challenge, especially differentiation between vegetative state (VS) and minimally conscious state (MCS). This study was to investigate the alterations of gut microbiota in patients with DOC and to identify potential microbiome biomarkers that can differentiate levels of consciousness. In this study, we collected fecal samples from patients with prolonged DOC, including 19 patients in MCS and 14 patients in VS; 16S-rRNA sequencing was used to investigate the gut microbiome of patients. Gut microbiota diversity, composition, and discriminant bacterial taxa were analyzed to identify potential biomarkers for differentiating levels of consciousness. We found that diversity and composition of gut microbiome were significantly altered in patients with DOC, and decreased alpha diversity was associated with lower levels of consciousness. Specific bacterial taxa including Firmicutes, Escherichia_shigella, Raecalibacterium, Lachnospiraceae, and Ruminococcaceae_UCG_013 were more abundant in patients in MCS, whereas Clostridiales were more abundant in patients in VS. In conclusion, the study results demonstrated that patients with DOC exhibited distinct diversity and composition of gut microbiota. And there was a decreasing trend of alpha diversity of gut microbiota from patients in MCS to patients in VS, which indicates that lower alpha diversity was associated with more severe level of unconsciousness. Specific bacterial taxa may be potential biomarkers to differentiate MCS and VS.
Introduction
The brain injury, including traumatic and nontraumatic causes, is one of the leading causes of death and long-term disability worldwide. 1 The fatality rate of brain-injured patients has decreased significantly due to the advances in clinical treatment and intensive care techniques. However, some patients would fall into varied degrees of disorders of consciousness (DOC) 2 including coma, vegetative state (VS), and minimally conscious state (MCS). Individuals in VS remain unresponsive to external stimulation or only show spontaneous movements, 3 whereas patients in MCS can show fluctuating but discernible behaviors of self or environmental awareness. 4 The differentiation between MCS and VS can be challenging since voluntary and reflexive behaviors can be difficult to distinguish and subtle signs of consciousness may be missed. Therefore, the accurate diagnosis and classification of patients with DOC are of great significance, which can not only determine the therapeutic options but also affect the prognosis of patients with DOC.
At present, the “gold standard” to assess prolonged DOC is the Coma Recovery Scale-Revised (CRS-R), 5 with six subscales including auditory, visual, motor, oromotor, communication, and arousal to evaluate the consciousness state. However, CRS-R score is based on patients’ behavioral response that can be affected by the fluctuations in arousal level, underlying sensory and motor impairments, as well as the subjective judgment of the operator may lead to a high rate of clinical misdiagnosis. 6 Therefore, to develop a more objective and accurate method for diagnosing patients with DOC, a large amount of studies have focused on neuroimaging and electrophysiology techniques in the past two decades, such as functional magnetic resonance imaging (MRI) and electroencephalogram (EEG). For instance, neuroimaging techniques such as MRI, 7,8 18F-fluorodeoxyglucose positron emission tomography (PET), 9 and diffusion tensor imaging 10,11,12 showed some promising value in evaluating patients with DOC. And EEG measures also demonstrated its utility in classifying different states of DOC. 13,14 Nevertheless, neuroimaging techniques are expensive and less practical in patients with DOC, and some task-based experiments required the cooperation of patients. With regard to bedside EEG, the complex paradigms, analysis algorithms, and feature extraction make it difficult to be widely utilized by clinicians. Therefore, to find a cheaper, accurate, and easily performed method to diagnose patients with DOC has always been one of the main focuses of researches in DOC.
Recently, compelling evidence from experimental and clinical studies has indicated the significance of microbiome in human health, and the association between the gut microbiota and neurological diseases has attracted great attention. It is well known that the central nervous system (CNS) and the gastrointestinal system can communicate bi-directionally via various pathways, which is also called brain–gut axis. 15 The brain–gut axis, which integrates neural, hormonal, metabolic, and immunological signals between the central nervous system (CNS) and the gut, provides a pathway for gut microbiota and their metabolites to access the brain. 16 Brain disease can also alter the gastrointestinal and immune function through top-to-bottom pathways. 17 With deeper understanding of the role of gut microbiota, researchers found that it can contribute to the integrity of the intestinal barrier, inflammatory response, and neuroendocrine functions and can produce metabolites with important physiological and pathological significance. 18 Furthermore, the gut microbiota are also involved in the regulation of the physiological processes 19 of the CNS.
A number of studies have reported the important roles of gut microbiota in CNS disorders, including Alzheimer’s disease, 20 Parkinson’s disease, 21 traumatic brain injury (TBI), 22,23 and stroke, 24 –26 making it a promising therapeutic target for these conditions. Our previous study found that gut microbiota changed significantly accompanied by alteration of bile acid metabolism in the acute stage of TBI, with dynamic changes in alpha diversity. 27 However, little information is known about the gut microbiota in patients with prolonged DOC after brain injury, Until recently, limited data existed on the gut microbiota in patients with prolonged disorders of consciousness (DOC) following brain injury. However, a recent study suggested that gut microbiota could differentiate DOC patients, though these findings require further validation. 28 Therefore, we hypothesize that patients with DOC after brain injury exhibited alterations in gut microbiome, which may be associated with different levels of consciousness. In addition, potential biomarkers in gut microbiome can be identified to distinguish MCS from VS, thereby promoting the diagnostic accuracy and recovery of patients with prolonged DOC.
Participants and Methods
Participants and study design
Patients with prolonged disorders of consciousness after brain injury were admitted to the Neurosurgery Department of the First Affiliated Hospital of Zhejiang University School of Medicine and the Rehabilitation Department of Hangzhou Wujing hospital between October 2019 and December 2020. Ethics approval for this study was obtained from the Ethics Committee of the First Affiliated Hospital of Zhejiang University School of Medicine, and informed consent was obtained from patients’ caregivers.
Inclusion criteria were as follows: (1) in line with previous research, patients were at least 3 months post injury and diagnosed as suffering from DOC based on the CRS-R, including MCS and VS 29 –31 ; (2) a history of brain injury, including TBI, cerebral hemorrhage, and ischemic stroke; (3) age 15 years or older; (4) patients with stable vital signs and intracranial conditions, without any indications for emergency surgery such as intracranial hematoma or hydrocephalus; (5) healthy volunteers with age and gender matched were enrolled as the control group; and (6) the consent of patients’ relatives or guardians was taken to enroll in the study. Exclusion criteria were as follows: (1) age 14 years or younger; (2) intravenous or oral antibiotics used within 1 month; (3) a diagnosis of intestinal disease; and (4) coexisted with other serious systemic organ dysfunctions or malignant diseases. In this study, older adolescents (aged from 15 to 18) who were also included for the study showed that gut microbiota in adolescents and adults were relatively stable. 32
Data collection and CRS-R evaluation
Demographic and clinical data were extracted from case records, including gender, age, ethnicity, etiology, past medical history, CRS-R score, diagnosis, and other general physical examination. Prior to study enrollment, the best of three consecutive behavioral assessments using the CRS-R were conducted within a 1-week period. Then the patients were diagnosed as MCS and VS based on CRS-R score. The demographic data of healthy controls (HC) were also collected.
Fecal sample collection and DNA extraction
We collected 2 g stool (about the size of a little finger) from each patient during the first bowel movement of the day. Fecal samples were collected and stored at −80°C, and bacterial DNA was extracted using E.Z.N. Fecal DNA Kit (D4015, Omega, USA), following the manufacturer’s instructions. The total DNA was eluted in 200 μL elution buffer.
16s-rRNA sequencing and quantitative polymerase chain reaction
The gene-specific sequences used in this protocol target the 16S V3 and V4 regions, and amplicons were generated using primers corresponding to the hypervariable regions (forward: CCTACGGRRBGCASCAGKVRVGAAT; reverse: GGACTACNVGGGTWTCTAATCC) for the first polymerase chain reaction (PCR). PCR amplicons were purified in accordance with the manufacturer’s instructions. Each reaction mixture (25 μL) contained 20 ng template DNA, 0.5 μL DNA polymerase, 2.5 μL dNTPs, 2.5 μL buffer solution, 1 μL of each primer, and PCR-grade water to adjust the volume. The cycling conditions were as follows: initial denaturation at 95°C for 3 mins, 24 cycles of denaturation at 94°C for 5 sec, annealing at 57°C for 90 sec, extension at 72°C for 10 sec. The PCR products were quantified to 10 nM by Qubit 3.0 Fluorometer (Invitrogen, USA). The library was sequenced on an Illumina MiSeq (Illumina, USA) platform, and 250 bp/300 bp paired-end reads were generated. DNA sequencing and data analysis were mainly completed by specialized biological sequencing companies.
Data processing and analysis
Two paired-end reads were merged based on overlapping sequences, and merged reads were trimmed and quality filtered using QIIME (version 1.9) and then dereplicated and clustered at 97% identity, generating an operational taxonomic units (OTUs) table. The DADA2 algorithm was used to denoise demultiplexed sequences. In this way, the final feature table and feature sequence were obtained, and subsequent analyses of diversity and species differences were all performed based on these normalized output data. Alpha diversity was applied to analyze complexity of species diversity through five indexes, including Observed-species, Chao1, Shannon, Simpson, and Good-coverage. The beta diversity as depicted in principal coordinate components analysis (PCoA) and clustering analysis (UPGMA) can reveal the diversity of species between different environmental communities. Furthermore, we used the SILVA (Release 132, https://www.arb-silva.de/documentation/release-132/) and the NT-16s database to accurately analyze the species composition and make the species abundance tables at the level of kingdom, phylum, class, order, family, genus, and species. LEfSe (linear discriminant analysis effect size) software was used to compare species differences among groups, and linear discriminant analysis was used to obtain the final differential species (i.e., biomarker). PICRUSt2 was used to predict functional abundances based only on marker gene sequences. The diagrams were implemented using the R package (V3.5.2).
Results
Baseline characteristics
A total of 33 patients with DOC after brain injury were included in this study, including 19 patients in MCS and 14 patients in VS based on CRS-R assessment. The average age of all 33 patients was 47.1 ± 17.6 years, and 21 (63.6%) patients were male. The average duration from brain injury to fecal sample collection was 6.6 months. Thirty HC with age and gender matched were included with an average age of 47.5 ± 16.7 years, and 17 (56.7%) of them were male. There were no significant differences in age, gender, etiology distribution, and time interval since brain injury between DOC and HC groups and between MCS and VS groups. Baseline characteristics of the patients are summarized in Table 1.
Demographic and Clinical Characteristics of Patients with DOC
CRS-R, Coma Recovery Scale-Revised; DOC, disorders of consciousness; MCS, minimally conscious state; TBI, traumatic brain injury; VS, vegetative state.
Decreased alpha diversity of gut microbiota was associated with levels of consciousness in patients with DOC
Alpha diversity analysis was performed to investigate the community richness and coverage within gut microbiome, and it was measured by Chao1, Goods_coverage, observed_otus, Shannon, and Simpson indexes. In this study, we found that Goods_coverage index (p < 0.01) was significantly lower in patients with DOC compared with HC. However, there were no significant differences in Chao1 (p = 0.064), observed_otus (p = 0.064), Shannon (p = 0.84), and Simpson (p = 0.97) indexes between the two groups (Fig. 1A–E). When we further compared these indexes among HC, MCS, and VS groups, we found that Chao1 and observed_otus of the HC group were higher than that of the VS group, while goods_coverage index in the HC group was lower than the VS group. There were no significant differences in the indexes between the MCS and the VS group. However, a descending trend of Chao1 and observed_otus was observed from the HC to the VS group, which indicated decreased alpha diversity was associated with the decreasing level of consciousness (Fig. 2A–E). These results suggested that a lower alpha diversity of gut microbiota may be an indicator of a more severe level of unconsciousness in patients with DOC.

Alterations of gut microbiota diversity in patients with DOC. Goods_coverage index was significantly lower in patients with DOC compared with HC

Alterations of gut microbiota diversity in patients in MCS and patients in VS. Chao1
Altered beta diversity demonstrated different composition of gut microbiota in patients with DOC
Beta diversity of gut microbiota refers to the potential differences in different bacterial communities. In our study, PCA, PCoA, and NMDS were used to analyze the beta diversity. In patients with DOC, a shift in microbial composition was observed in PCA1 and PCoA1, which accounted for 15.47% and 35.08% of the inter-sample variation. NMDS analysis showed that there was a marked difference between the HC and patients with DOC patients (Fig. 1F–H). Then we further compared the beta diversity between MCS, VS, and HC. The results of PCA, PCoA, and NMDS analysis showed that there were significant differences in the composition of gut microbiota among them. These results indicated that patients with DOC with different levels of consciousness exhibited distinguished gut microbial composition (Fig. 2F–H).
Patients with DOC exhibited a different spectrum of abundant bacterial taxa
The top 30 species with the highest abundance in each taxonomic hierarchy are shown in Figure 3. At the phylum level, we observed decrease of abundance in two phyla, including Firmicutes and Actinobacteria and increase in four phyla (Bacteroidetes, Proteobacteria, Verrucomicrobia, and Fusobacteria) in patients with DOC when compared with HC. In further comparison between MCS and VS, the relative abundance of Firmicutes, Proteobacteria, and Fusobacteria was lower in patients in VS than in patients in MCS, while that of Bacteroidetes, Actinobacteria, and Verrucomicrobia was higher in VS than in MCS. At the genus level, Bacteroides, Klebsiella, and Escherichia_shigella were more abundant in patients with DOC, while Bifidobacterium, Faecalibacterium, prevotella_9, and Lachnospira were less abundant in patients with DOC than HC. In further comparison between MCS and VS, the relative abundance of Faecalibacterium and Klebsiella was lower in patients in VS than patients in MCS, while that of Bacteroides and Bifidobacterium was higher in VS than MCS. Due to the limited depth of 16S-RDNA sequencing, we were unable to account for all of the detected sequences at species level.

Taxonomy stacked bar displaying the relative abundance of bacteria taxa in each group at different levels. DOC and HC at phylum level
LefSe analysis identified discriminant taxa in patients with DOC
LefSe analysis was used to identify discriminant bacterial taxa among different groups; we compared the species differences between the patients with DOC and HC, as well as MCS and VS groups. Thirty differentially expressed taxa were obtained in patients with DOC. The top 10 taxa included Bacteroides, Bacteroidaceae, Proteobacteria, Gammaproteobacterial, Enterobacteriales, Enterobacteriaceae, Klebsiella_unclassified, Verrucomicrobia, Verrucomicrobiae, and Akkermansiaceae. Comparing patients in minimally conscious states and vegetative states, we found higher abundances of five taxa—Firmicutes, Escherichia-Shigella, uncultured_Faecalibacterium_sp, Lachnospiraceae_unclassified, and Ruminococcaceae_UCG_013 in MCS patients, while Clostridiales_unclassified abundance was higher in VS patients. (Fig. 4).

Linear discriminant analysis of gut microbiota effect in DOC and HC group
Functional analysis of gut microbiota biomarkers identified potential pathways involved in DOC after brain injury
We further used PICRUSt2 analysis to predict the potential functions of differently abundant genera. And the results revealed that in the HC group, abundant bacterial taxa were associated with adenosine triphosphate synthesis and bacterial cell division, while abundant genera in DOC group were associated with multidrug efflux pump, heat shock protein, and beta lactamase family. In comparison of patients in MCS and patients in VS, we found that abundant genera in the VS group were associated with galactosidase function, while abundant genera in the MCS group were associated with sugar fermentation and protein biosynthesis (Fig. 5). These results indicated that these biological pathways may be involved in DOC after brain injury.

Functional analysis of gut microbiota biomarkers. Functional prediction of abundant genera in DOC and HC group
Discussion
In this study, the results showed that the diversity and composition of gut microbiome were significantly altered in patients with DOC, and the alpha diversity was significantly lower in patients with DOC when compared with the HC. And the alpha diversity of gut microbiota in patients in VS was also lower compared with patients in MCS though not reaching statistical significance. Moreover, we found that specific bacterial taxa including Firmicutes, Escherichia_shigella, Faecalibacterium, Lachnospiraceae, and Ruminococcaceae_UCG_013 were more abundant in patients in MCS, whereas Clostridiales were more abundant in patients in VS. Collectively, our findings demonstrated the alterations of gut microbiome in patients with DOC and identified potential bacterial biomarkers for diagnosis and evaluation of patients with DOC.
It is known that accurate evaluation of levels of consciousness in DOC is the prerequisite for following effective treatment. However, there was still a high rate of misdiagnosis in these patients, and some studies reported that about 30–40% of patients in MCS were misdiagnosed as in VS. 33,34 Electrophysiological and neuroimaging techniques have been proposed to aid the diagnosis of patients with DOC. For example, studies reported that EEG networks had a discriminative capacity (80%) for the differentiation between VS and MCS. 35,36 And resting-state PET has also been used to distinguish MCS from VS. 36,37 In addition, some researchers found that TMS-EEG based on PCI can classify individuals in VS, individuals in MCS, and healthy individuals. 38,39 However, all these existing methods suffered from some shortcomings, such as high expense and inconvenience in clinical use. In recent years, increasing knowledge about microbiota and gut–brain axis shed some insights into the evaluation of patients with DOC. For instance, Lucidi et al. found that Christensenellaceae is a marker of disease progression in bipolar disorder, and decrease in relative abundance of Faecalis was associated with disease severity. 40 Li et al. found that changes in the gut microbiota may occur several years before the onset of dementia, even in the mild cognitive impairment stage. 41 Inspired by these studies, we investigated gut microbiome in patients with DOC and tried to identify some potential biomarkers for evaluating these patients in this study. And encouragingly, we found that patients with had decreased alpha diversity in gut microbiota compared with HC. Further analysis revealed a decreasing trend in alpha diversity (Chao1 index) across groups: highest in healthy controls (HC), followed by MCS patients, and lowest in VS patients. This suggests that reduced alpha diversity may correlate with the severity of consciousness impairment. Therefore, our results indicated that alpha diversity was a potential biomarker in the evaluation of patients with DOC and provided the preliminary evidence of association between gut microbiome alpha diversity and levels of consciousness.
By further analyzing the gut microbiota at different levels, we found that at phylum level, increased abundance of Bacteroides and Proteobacteria was observed in patients with DOC. At the genus level, Bacteroides, Klebsiella, and Escherichia_shigella were more abundant in patients with DOC. While the abundance of Firmicutes, Proteobacteria, and Fusobacteria was lower in patients in VS than in patients in MCS, while that of Bacteroidetes, Actinobacteria, and Verrucomicrobia was higher in VS than MCS. At the genus level, abundance of Faecalibacterium and Klebsiella was lower in patients in VS than in patients in MCS, while that of Bacteroides and Bifidobacterium was higher in VS than MCS, which was similar to previous studies in experimental TBI. 42,43 For instance, Nicholson et al. found that the relative abundance of traditional beneficial bacteria decreased after TBI, especially in the families Trichospiraceae, Mogamiaceae, and Ruminococcaceae within Firmicutes. On the contrary, the number of Bacteroidetes, Proteobacteriaceae, and Enterobacteriaceae increased after TBI. 42 Bacteroidetes were found to be involved in the inflammatory process, also known as opportunistic pathogens. 20 Moreover, increasing brain injury severity was associated with decreased Firmicutes abundance, increased Proteobacteria abundance, and a significant reduction in alpha diversity. 42 Therefore, the results of our study suggested that with the increase of disease severity, the probiotics in gut microbiome decreased while opportunistic pathogens increased. This also suggested that intervention of intestinal flora, such as probiotics, 44 may be a potential method in the future to inhibit the growth of mechanistic pathogenic bacteria, thus reducing secondary injury and improving prognosis of patients with prolonged DOC. In addition, LefSe analysis identified that specific taxa Firmicutes, Escherichia_shigella, Faecalibacterium, Lachnospiraceae, and Ruminococcaceae_UCG_013 were more abundant in patients in MCS, whereas Clostridiales were more abundant in patients in VS. This was partly consistent with the previous study reporting Faecailbacterium, Enterococcus, and Methanobrevibacter as specific bacterial taxa to distinguish between MSC and VS. 28
Another important goal of our study was to explore the biological functions and potential effects of discriminant bacteria taxa in patients with DOC, such as corresponding inflammatory, immune, or metabolic pathways related to these alterations. Therefore, we performed functional analysis of these specific taxa and found that abundant bacterial taxa in the DOC group were associated with multidrug efflux pump, heat shock protein, and beta lactamase family. A comparison of patients in MCS and patients in VS revealed that abundant taxa in VS group were associated with galactosidase function, while abundant taxa in MCS group were associated with sugar fermentation and protein biosynthesis. Previous studies have demonstrated that the gut microbiota are involved in regulating immune and inflammatory responses in both acute and chronic neurological diseases and can further influence patients’ outcome. 45,46 For instance, the abundance of Clostridiales after 24 h TBI was found to be closely related to the intestinal immune response through modulating the migration of TH17 cells, 47,48 While in ischemic brain injury, reduced Proteobacteria levels and elevated Firmicutes levels led to the increase of anti-inflammatory Treg cells as well as the decrease of pro-inflammatory γδ T cells, which will eventually produce certain neuroprotective effects. 49 In addition, some intestinal bacteria influence the brain function through microbial metabolites, such as short-chain fatty acids, 50,51 bile acids, 52,53 and amino acids. 23 One study demonstrated that multiple sclerosis can induce the BA metabolism alterations, and the supplement of bile acids can alleviate the neuroinflammation through inhibiting polarization of astrocytes and microglia toward neurotoxic phenotypes. 53 Our previous study also found that TBI-induced intestinal microbiome dysregulation promoted intestinal inflammation by reducing the level of secondary bile acids. The change of bile acid to lithocholic acid was negatively correlated with Staphylococcus, and tauroursodeoxycholic acid was positively correlated with Lachnospiraceae. 27 Moreover, a retrospective study by our team also found that the serum bile acid level of patients after TBI was significantly lower than that in healthy volunteers, and TBI was an independent factor of the decreased level of serum bile acid. 54 Gut microbiota such as Lachnospiraceae and Firmicutes 55,56 might affect intestinal inflammation and barrier function by changing bile acid metabolism. And peripheral inflammation can stimulate the activation of microglia and astrocytes in the brain, 57 leading to the release of cytokines in the brain, resulting in the death of neurons, the destruction of the integrity of the blood–brain barrier, and even the disorders of consciousness. Taken together, we speculated that serum bile acid level may be an important factor that links gut microbiome and DOC after brain injury. However, further confirmatory studies utilizing fecal transplantation will be required to investigate the potential mechanisms underlying gut microbiome alterations and prognosis of DOC.
While the findings of our study demonstrated the association between DOC and gut microbiome, some limitations should be recognized. First, this is a preliminary study with a relatively small sample size to investigate potential biomarkers in gut microbiome to distinguish levels of consciousness in patients with prolonged DOC. There results based on the differential expression should be interpreted with caution. Second, the 16s-rRNA sequencing we used could only account for changes in the flora at phylum and genus levels but not at species level. So further studies in metagenomic sequencing are necessary to explore the alterations of gut microbiota at species level in patients with DOC. Third, patients with brain injury of varied etiologies were included in this study, and investigation of the impact of etiologies on gut microbiome was not performed due to the relatively small sample size. Nevertheless, previous studies also indicated that it is the brain function status, not the etiology, that played a more important role in modulating gut microbiome. In conclusion, future studies with larger sample size, combined with metagenomic, metabolomics, and fecal transplantation, are needed to further elucidate the value of gut microbiome as biological markers for diagnosis and evaluation of patients with DOC, as well as the possibility of potential therapeutic targets.
Transparency, Rigor, and Reproducibility Statement
This study was conducted to investigate the alterations of gut microbiota in patients with DOC and to identify potential microbiome biomarkers that can differentiate levels of consciousness. This was an exploratory study with a relatively small sample size, so the study was not formally registered. The analysis plan was not formally pre-registered, but the team member with primary responsibility for the analysis certifies that the analysis plan was prespecified. A sample size of 30 subjects per group was planned based on the availability of patients with DOC. As a result, we included 33 patients with DOC and 30 HC in the study. The authors agree to provide the full content of the article on request.
Footnotes
Acknowledgment
The authors thank the participants and their families for their participation in this study.
Authors’ Contributions
W.Y.: Conceptualization, methodology, investigation, data analysis, and writing—original draft. Y.J.L. and A.W.: Sample collection and data analysis. J.C. and L.W.: Language improvement. Y.X.L.: Article revision. X.Y.: Writing—review and editing. Z.W.: Conceptualization, supervision, investigation, and writing—review and editing.
Author Disclosure Statement
The authors declare that they have no competing interests.
Funding Information
The study was supported by research projects (2023YSJYX-ZD-2, 2022-YJRC3935) and Fujian Technology Innovation Joint Project (2023Y9080).
