- Open Access
Dysbiosis of the intestinal microbiota in neurocritically ill patients and the risk for death
Critical Carevolume 23, Article number: 195 (2019)
Despite the essential functions of the intestinal microbiota in human physiology, little has been reported about the microbiome in neurocritically ill patients. This investigation aimed to evaluate the characteristics of the gut microbiome in neurocritically ill patients and its changes after admission. Furthermore, we investigated whether the characteristics of the gut microbiome at admission were a risk factor for death within 180 days.
This prospective observational cohort study included neurocritically ill patients admitted to the neurological intensive care unit of a large university-affiliated academic hospital in Guangzhou. Faecal samples were collected within 72 h after admission (before antibiotic treatment) and serially each week. Healthy volunteers were recruited from a community in Guangzhou. The gut microbiome was monitored via 16S rRNA gene sequence analysis, and the associations with the clinical outcome were evaluated by a Cox proportional hazards model.
In total, 98 patients and 84 age- and sex-matched healthy subjects were included in the analysis. Compared with healthy subjects, the neurocritically ill patients exhibited significantly different compositions of intestinal microbiota. During hospitalization, the α-diversity and abundance of Ruminococcaceae and Lachnospiraceae decreased significantly over time in patients followed longitudinally. The abundance of Enterobacteriaceae was positively associated with the modified Rankin Scale at discharge. In the multivariate Cox regression analysis, Christensenellaceae and Erysipelotrichaceae were associated with an increased risk of death. The increases in intestinal Enterobacteriales and Enterobacteriaceae during the first week in the neurological intensive care unit were associated with increases of 92% in the risk of 180-day mortality after adjustments.
This analysis of the gut microbiome in 98 neurocritically ill patients indicates that the gut microbiota composition in these patients differs significantly from that in a healthy population and that the magnitude of this dysbiosis increases during hospitalization in a neurological intensive care unit. The gut microbiota characteristics seem to have an impact on patients’ 180-day mortality. Gut microbiota analysis could hopefully predict outcome in the future.
The crosstalk between the host and intestinal flora strongly influences host physiology in health and disease. However, the understanding of the microbiome is still in its infancy. Recently, many reports have revealed that changes in the gut microbiota are related to myasthenia gravis , Parkinson’s disease  and critical illness . Severe insults to the gut have been hypothesized to promote multiple organ dysfunction syndrome as a result of disruption of commensal bacteria . In previous studies, a low microbial diversity has been associated with an increased risk of mortality in critically ill patients [5, 6], and domination by certain pathogens has been identified as an independent risk factor for adverse outcomes [7,8,9]. Neurocritical care in the management of acute neurologic injury has rapidly evolved over the last 30 years . However, little information is available regarding the ways in which the microbiota changes during neurocritical illness and whether the features of the intestinal microbiome are associated with mortality in the setting of neurocritical care.
The implementation of prebiotics , probiotics  or faecal transplant  to restore the “health-promoting” microbiome suggests that gut microbiota modulation can be used to monitor or cure disease [14,15,16]. However, the composition and community structure of microbes that colonize the human gut during critical illness have not been fully understood. Advances in next-generation sequencing provide opportunities to determine the gut microbiome composition in individuals as “microbial fingerprints”. Therefore, additional knowledge of the characteristics of the gut microbiota in neurocritically ill patients is required to better understand the microbial disturbances associated with clinical outcomes to identify patients who would benefit from microbiota-targeted therapies .
In this study, we characterized the phylogenetic composition of the faecal microbiota in 98 patients suffering from neurocritical illness using 16S rRNA gene sequence analysis and performed a subgroup analysis in patients with critically ill stroke. We also examined whether the abundances of specific taxa in the intestinal microbiota were associated with mortality within 180 days of admission.
This study was a prospective observational study conducted in the neurological intensive care unit (neuroICU) of Nanfang Hospital, a university-affiliated academic hospital in Guangzhou, China, from December 2016 to February 2018. Adult patients admitted to the neuroICU with a Glasgow Coma Scale (GCS) of < 11 and an expected length of intensive care unit (ICU) stay of > 48 h were eligible for inclusion. Exclusion criteria included (1) transfer from other ICUs or hospitals; (2) treatment with antibiotics, prebiotics or probiotics before the first sample was obtained; (3) inability to obtain the first faecal sample within 72 h after admission; (4) use of antibiotics or admission to hospitals within 6 months prior to participation in the current study; and concomitant (5) gastrointestinal disease, (6) malignant cancer or (7) pregnancy. GCS, Acute Physiology and Chronic Health Evaluation (APACHE-II) and Sequential Organ Failure Assessment (SOFA) scores were recorded by specialized neurocritical care physicians at admission. The modified Rankin Scale (mRS) was recorded by specialized neurocritical care physicians at discharge. The healthy subjects were recruited from a community (the Bureau of Reclamation in Guangzhou) between November 2016 and January 2017. The volunteers who met the following inclusion criteria were recruited: (1) self-reported as healthy; (2) had no cross-border tourism recently; (3) had not taken antibiotics, prebiotics or probiotics within 1 year; and (4) had not been admitted to hospitals within 1 year prior to participation in the current study. Age- and sex-matching processes were further performed to select comparable controls. Written informed consent was obtained from all healthy subjects and patients or their legal representatives. Ethical approval for both the patients and healthy subjects was received from the Medical Ethics Committee of Nanfang Hospital (No. NFEC-2018-034), and all studies were conducted in accordance with the Declaration of Helsinki.
Faecal sample collection and DNA extraction
The faecal samples were obtained from the patients within 72 h after admission and serially every week and were collected once from individuals in the control group. Fresh stool samples were frozen at − 80 °C within 3 h after voiding, and 0.2 g of each sample was aliquoted for DNA extraction. Bacterial DNA was extracted with a stool DNA extraction kit using a magnetic bead-based method (Shenzhen Bioeasy Biotechnology Co., Ltd., China) according to the manufacturer’s instructions .
Polymerase chain reaction amplification of bacterial 16S rRNA genes
The V4 region of the bacterial 16S rRNA gene was amplified by quantitative real-time polymerase chain reaction (qPCR) with the bar-coded primers V4F (5′-GTGTGYCAGCMGCCGCGGTAA-3′) and V4R (5′-CCGGACTACNVGGGTWTCTAAT-3′) using a LightCycler 480 II real-time fluorescence quantitative PCR system (Roche Diagnostics Ltd., Switzerland). The qPCR protocol was as follows: (1) initial denaturation (94 °C, 2 min); (2) PCR amplification (32 cycles, 94 °C, 30 s; 52 °C, 30 s; followed by 72 °C, 30 s); (3) melting (95 °C, 5 s and 60 °C, 1 min; followed by 95 °C, continuous); and (4) cooling (37 °C, 30 s). Agarose gel electrophoresis (1%) was used to detect the PCR products. Samples that produced a visible product 290–310 bp in length were used for further experiments. GeneTools analysis software (version 4.03.05.0, SynGene) was used. The PCR products were mixed in equimolar ratios and purified using an EZNA Gel Extraction Kit (Omega, USA). Sequencing libraries were established using an NEBNext Ultra™ DNA Library Prep Kit for Illumina (New England Biolabs, USA) according to the manufacturer’s recommendations, and the index codes were added. The library quality was assessed using the Qubit@ 2.0 Fluorometer (Thermo Scientific) and Agilent Bioanalyzer 2100 systems. Finally, the library was sequenced on an Illumina HiSeq 2500 platform, and 250-bp paired-end reads were generated.
Sequencing and microbial analysis
Sequences longer than 200 bp were trimmed to 200 bp, and those shorter than 200 bp were removed. Sequencing 200 bp at one end of each fragment had the advantages of a good overlap effect, minimal use of computing resources and improved sequence quality [18, 19]. Depending on the overlap between the two paired-end sequences, we then used SeqPrep to merge the paired-end sequences and assessed the quality of the result using the open-source software Quantitative Insights into Microbial Ecology (QIIME, version 1.9.1) . A QIIME workflow script, split_libraries_fastq.py, was used to check the quality of the sequences, and sequences with a Phred score of ≥ Q20 were considered qualified sequences. Then, we used a home-brewed script to split FASTA files according to the paired-end barcode information, which met the following criteria: a 100% match between the barcode and the primer and more than 200 bp remaining after the removal of the barcode and primer. Finally, we used a QIIME workflow script, pick_closed_reference_otus.py, to remove chimaeras, perform reference-based operational taxonomic unit (OTU) clustering and generate a BIOM file. We used Greengenes 13.8 as the reference database and SortMeRNA for clustering and classification with an identity cut-off value of 0.97. All samples were normalized to the same level to avoid possible errors stemming from the use of different sequencing depths. Each sample was normalized to 7000 sequences for subsequent analysis.
The α-diversity (the complexity within a community) was estimated by four indexes and calculated by QIIME : (a) Chao1, based on the richness of the sample; (b) observed species, a direct measure of the species number; (c) Shannon, intended to represent the species abundance and evenness; and (d) phylogenetic diversity (PD)–whole tree, an alternate method that accounts for the phylogenetic differences between the species that comprise the community. The β-diversity (which estimates differences between microbial communities) was analysed using the Bray–Curtis distance and unweighted UniFrac distance approaches [21, 22]. To determine the significantly different taxa between two groups, linear discriminant analysis (LDA) coupled with effect size measurement (LEfSe) was performed using an online utility (http://huttenhower.sph.harvard.edu/galaxy) . Significantly different bacteria with LDA scores of ≥ 2.0 or ≥ 3.5 were diagrammed on taxonomic bar plots.
The continuous nonparametric data are presented as medians (interquartile ranges, IQRs) and were analysed using Mann–Whitney U or Wilcoxon tests. The continuous parametric data are presented as the means (standard deviations, SDs) and were analysed with Student’s t test. The categorical data are presented as numbers (percentages, %) and were analysed using chi-squared tests. For microbial analysis, QIIME was additionally performed using the Adonis test as previously described . A univariate Cox proportional hazards model was first performed, and the candidate variables with a p value of < 0.05 were further included in the multivariate Cox regression model for adjustment. The candidate variables included demographic data (age and sex), comorbidities (diabetes, hypertension, heart failure, intracranial infection, intracranial hypertension, and pneumonia), serum markers (white blood cell count, creatinine, blood urea nitrogen, cystatin C, brain natriuretic peptide, and C-reactive protein), indicators of critical illness (GCS, APACHE-II and SOFA scores and the length of ICU stay), medical treatments (mechanical ventilation and enteral nutrition), the biodiversity of the microbiota (Shannon, PD–whole tree, Chao 1, observed species and Simpson indexes) and the z-scores of specific taxa abundances or their changes. The abundance changes (marked as “Δtaxon”, e.g., “ΔEnterobacteriaceae”) were calculated by subtracting the abundance of specific taxa in the first sample from the abundance of specific taxa in the second sample. The association between abundances and 180-day mortality was analysed based on the increase in the abundances per SD. Correlations between variables were determined with the Spearman’s rank correlation test. SPSS version 20 (Statistical Package for Social Sciences, Chicago, IL, USA) was used for statistical analysis. Two-tailed p values of < 0.05 were considered statistically significant. The figures were generated using R version 3.4.3 (https://www.r-project.org/).
The first faecal samples were collected from 98 neurocritically ill patients (median age 58.5 years; 62.2% male) within 72 h after neuroICU admission. The flow diagram of the patient selection process is shown in Fig. 1. Of these 98 patients, 38 were admitted for ischaemic stroke, 20 for intracerebral haemorrhage, 13 for intracranial infection and 22 for seizure, hypoxic-ischaemic encephalopathy or another conditions (Additional file 3). Table 1 provides a summary of the characteristics of the study patients. In total, 206 faecal samples were collected for microbial analysis. These samples included 98 collected within 72 h after admission, 50 in the second week, 17 in the third week, 12 in the fourth week, 6 in the fifth week, 6 in the sixth week, 5 in the seventh week, 4 in the eighth week, 3 in the ninth week and 5 in the tenth to fourteenth weeks (Additional file 4). Eighty-four healthy subjects served as the healthy controls (HCs) and had faecal samples collected one time. The neuroICU group and HCs did not significantly differ in age (Mann–Whitney U test, p > 0.05) or sex (chi-squared test, p > 0.05) and did not differ in history of smoking, hypertension or diabetes (chi-squared test, p > 0.05) (Additional file 2: Table S1). The principal coordinate analysis (PCoA) plot did not show a significant difference in β-diversity (Bray–Curtis distance) between patients with different primary diagnoses, indicating that regardless of the primary diagnosis at admission, the gut microbiome of patients in the neuroICU differed substantially from that of healthy subjects (Additional file 2: Table S2; Additional file 1: Figure S1).
Comparison of the gut microbiota characteristics in neurocritically ill patients and healthy subjects
We examined the first 98 faecal samples to explore the characteristics of the gut microbiota in the neuroICU patients (Additional file 1: Figure S2). The β-diversity and α-diversity between neuroICU patients and HCs were significantly different. As indicated by taxonomic summary and LEfSe analysis, the relative abundances of the phyla Proteobacteria, Deferribacteres and Verrucomicrobia were higher in the neuroICU group than in the HC group. At the family level, Enterobacteriaceae, Porphyromonadaceae, Enterococcaceae, Verrucomicrobiaceae, Rikenellaceae and Lactobacillaceae were enriched in the neuroICU group. To show the heterogeneity among individual gut microbiota compositions, we generated area plots from the first sample of each individual. At the phylum level, a high interindividual diversity in the faecal microbiota composition was observed in the neuroICU group (Fig. 2a). The faecal samples from patients P40 revealed a high abundance of Verrucomicrobia (> 30%), similar to the results reported by Lankelma et al. . The bacterial composition analysis at the family level revealed an even greater interpatient diversity than that seen at the phylum level (Fig. 2b). Enterococcaceae comprised more than 70% of the microbiota in P31, and Enterobacteriaceae comprised more than 70% in two other patients (P53 and P64). However, none of the microbiomes in the HCs were dominated by any single pathogen to an extent of greater than 70%.
Changes in the gut microbiota composition during neuroICU hospitalization
To examine the gut microbiota dynamics throughout the neuroICU stay, we next studied the community β-diversity in 50 patients from whom paired stool samples were obtained within 72 h after admission and at discharge. The other 48 patients were excluded, as they had only one available stool sample during their hospital stay. The bacterial community composition in all of these patients exhibited a substantial compositional divergence at discharge relative to that at admission (Fig. 3a). The samples collected at admission exhibited significantly higher community composition diversity (as indicated by the Shannon index) than the samples collected at discharge (Wilcoxon test, p < 0.001) (Fig. 3b). The microbiota changed radically after neuroICU admission. The α-diversities of the samples, which included only samples collected more than once from a single patient, were correlated inversely with the number of days spent in the neuroICU (r = − 0.211, p = 0.008; r = − 0.196, p = 0.015, respectively) (Fig. 3c, d). The distances between each sample in the neuroICU and HC groups were calculated by the Bray–Curtis and unweighted UniFrac metrics, which were plotted against the length of hospital stay. The Bray–Curtis distances and unweighted UniFrac were correlated positively with the number of days spent in the neuroICU (r = 0.207, p = 0.010; r = 0.220, p = 0.006, respectively) (Fig. 3e, f). The α-diversities changed greatly in the serial samples collected from the same patients several weeks apart (Additional file 1: Figure S3). In some subjects, the α-diversity recovered by the end of hospitalization (P01, P05 and P26), which was similar to a previous report .
The correlation analysis between six high-abundance taxa and sample time is shown in Fig. 4. At the phylum level, Bacteroidetes and Proteobacteria were not correlated with the length of hospital stay (r = 0.001, p = 0.990; r = 0.115, p = 0.127, respectively). In addition, Firmicutes was negatively associated with the sample collection time (r = − 0.217, p = 0.007). Ruminococcaceae and Lachnospiraceae were significantly correlated with length of stay (r = − 0.280, p < 0.001; r = − 0.357, p < 0.001, respectively). At the family level, the correlation between Enterobacteriaceae and the length of stay exhibited a trend towards significance (r = 0.153, p = 0.059).
To explore the correlation between health status and intestinal microbiota, we further performed Spearman correlation analyses in patients who survived at discharge (n = 79), which showed that the increases in Proteobacteria (r = 0.254, p = 0.012), Gammaproteobacteria (r = 0.203, p = 0.045) and Enterobacteriaceae (r = 0.240, p = 0.017) were paralleled with the increase in mRS scores assessed at discharge.
Associations between the abundance of intestinal microbiota and 180-day mortality
During the follow-up period, 25 (25.5%) patients died within 180 days of admission. A Cox proportional hazards regression model was constructed to examine the associations between abundances of specific taxa and 180-day mortality after ICU admission. The results of the univariate model showed that Christensenellaceae and Erysipelotrichaceae were associated with 180-day mortality and were independently associated with 180-day mortality in the multivariate Cox regression analysis adjusted for age, APACHE-II, white blood cell count and serum creatinine (adjusted hazard ratio (aHR) = 1.545; 95% confidence interval (CI) 1.163–2.053, p = 0.003; and aHR = 1.493; 95% CI, 1.094–2.038, p = 0.012, respectively) (Table 2).
Then, the association between the microbiota abundance changes and 180-day mortality after ICU admission was analysed by a Cox regression model. Fifty patients with paired samples both in the first 72 h and in the second week were included. The other 48 patients were excluded as they had only one available stool sample during their hospital stay. The results of the univariate model showed that the abundance changes in Fusobacteria, Fusobacteriales, Enterobacteriales, Veillonellaceae, Fusobacteriaceae and Enterobacteriaceae were associated with the 180-day mortality (Table 3). However, only increases in Enterobacteriales and Enterobacteriaceae were found to be independently associated with 180-day mortality in the multivariate Cox regression analysis adjusted for APACHE-II, respiratory failure, intracranial hypertension and serum lactate (Table 3). The increases in Enterobacteriales and Enterobacteriaceae during the first week in the neuroICU were both associated with increases of 92% in the risk of 180-day mortality after adjustments.
To eliminate confounding from the primary diagnosis and the differences in pathophysiology, we performed a subgroup analysis in patients with critically ill stroke, including 58 patients diagnosed with ischaemic stroke and intracerebral haemorrhage (median age 63 years; 60.3% male; median NIHSS score 12) (Table 4). Fifty-eight age- and sex-matched healthy subjects were selected as HCs (Additional file 2: Table S3). The characteristics of the gut microbiota in patients with stroke were significantly different from those in healthy subjects (Additional file 1: Figure S4), as indicated by the β-diversity, α-diversity, taxonomic summary and LEfSe analysis.
In this study, an in-depth faecal microbiota analysis was performed in neurocritically ill patients. We documented both changes in the diversities and large interindividual diversities in these patients relative to those characteristics in HCs. An overgrowth of opportunistic pathogens defined dysbiosis in patients with neurocritical illness. We also found that the bacterial community composition in the intestine at discharge significantly diverged from that observed at admission. In addition, the Cox regression analyses demonstrated that the abundances of Christensenellaceae and Erysipelotrichaceae at admission, as well as the abundance increase in Enterobacteriaceae during the ICU stay, were significantly associated with patients’ mortality within 180 days. To the best of our knowledge, this prospective observational study represents the largest cohort study of neurocritically ill patients to perform 16S rRNA gene sequence-based analyses of the gut microbiota and explore the association between the abundances of specific bacteria and survival status during the following 180 days in the setting of neurocritical illness.
Intestinal microbiota dysbiosis can be defined as an altered dialogue between the host’s cells and enteric bacteria due to disruption of microbial diversity usually manifested as the dominance of a given taxon. In recent decades, basic and translational studies have elucidated the functions of and architectural shifts in the intestinal microbiota in specific conditions and processes such as metabolic disease , central nervous system function [28, 29], systemic inflammation [30, 31] and immune function . More recently, convincing evidence has emerged to report the appearance of pathogens, disappearance of commensals and loss of microbial diversity in the pathogenesis of critical illness [3, 9, 25, 26, 33,34,35]. We previously reported significant dysbiosis of the gut microbiota in patients with stroke and transient ischaemic attack  or chronic kidney disease . Very recently, our research showed that participants at a high risk of stroke were characterized by the enrichment of opportunistic pathogens, low abundance of butyrate-producing bacteria and reduced concentrations of faecal butyrate . As a whole, these investigations could contribute to an increased awareness among clinicians that disturbances in the gut microbiota might be associated with poor prognosis. Dysbiosis of the microbiota can be reasonably presumed to increase the risk of adverse events, including infection, undernutrition and unconsciousness. Consistent with the findings of prior studies, our results confirmed that the microbiome in patients in neurocritical care differed substantially from that in healthy subjects. For example, the samples from the neuroICU group were characterized by gut microbiota dysbiosis and a loss of health-promoting commensal microbes. These differences were further characterized using LEfSe, which showed that neuroICU patients tended to have lower relative abundances of Firmicutes and Bacteroidetes and an increased level of Proteobacteria than healthy subjects. This difference is unlikely to be a study effect, as the samples from patients were rigorously processed and run in parallel with those from healthy subjects. Moreover, we collected our samples in the same city to partially avoid geographical effects . In addition, the disruption of the microbial community appeared to be greater at discharge than during the ICU hospitalization, with the decreased abundance of the butyrate-producing bacteria Lachnospiraceae and Ruminococcaceae. This ICU-acquired dysbiosis likely resulted from the pressures imposed by various patient-related factors and modern intensive care therapy. The use of antibiotics , vasoactive agents, artificial nutrition, proton pump inhibitors, analgesics or sedative agents that impair intestinal motility  may affect the gut microbiota. Antibiotic treatment precluded distinguishing the effect of critical illness from the effect of antibiotic pressure. However, the current study failed to address the impact of antibiotics on microbiota since most patients enrolled in this cohort received antimicrobial treatments, which is a common clinical practice in critical care. Because of the ubiquitous use of antibiotics worldwide , antibiotic pressure should be considered both a medical treatment and an inevitable pathophysiological insult that could exert negative effects on beneficial organisms.
Our findings that the abundances of Christensenellaceae and Erysipelotrichaceae are potential risk indicators of mortality within 180 days are consistent with previous reports. Christensenellaceae and Erysipelotrichaceae have been found to be associated with metabolic disorders and gastrointestinal diseases [41,42,43]. Christensenellaceae was reported to be enriched in patients with multiple sclerosis . Some studies have shown that HIV-infected individuals  and individuals with non-alcoholic fatty liver disease  and Rett syndrome  exhibited a high abundance of Erysipelotrichaceae. Moreover, our results showed that the abundance increase in Enterobacteriaceae is associated with death within 180 days. As previously reported, Enterobacteriaceae is one of the most detrimental pathogens to humans, especially for patients in ICUs [48, 49]. Our results provide evidence potentially suggesting that high abundances of some specific pathogens or their abundance changes may contribute to poor clinical outcomes. However, due to the limited sample size, lack of confounder control and regional variations, the results from the Cox regression analysis are insufficient for worldwide generalization. Nonetheless, the relationship between the gut microbiota and poor prognosis might serve as a platform for identifying potential targets for mortality prediction in the future. Freedberg et al.  provided novel insights into the possibility of predicting death or subsequent infection in critically ill patients by the characteristics of the gut microbiota upon ICU admission. Both colonization with vancomycin-resistant enterococci and the dominance of Enterococcus species were associated with death or all-cause infection adjusted by illness severity. Some studies have suggested that critical illness can lead to translocation of gut microbes into the bloodstream [50, 51], potentially explaining the association between specific pathogens and mortality. Identifying the critically ill patients most at risk of mortality is crucial for medical treatment strategies and doctor-patient communication. An intriguing possibility is analysing stool samples by 16S rRNA gene sequence as a non-invasive method for pretreatment assessment at admission to describe the role of the faecal “fingerprint” in predicting poor outcomes.
As the gut is hypothesized to play key roles in the progression of critical illness and multiple organ failure [4, 52], the restoration of commensal “healthy microbes” or eradication of pathogens might exert beneficial effects in critically ill patients. Microbiome intervention therapy includes antibiotic treatment, probiotic administration and faecal microbiota transplantation. Two case reports have discussed the successful application of faecal microbiota transplantation in patients with diarrhoea and therapy-resistant sepsis [53, 54]. However, these strategies should be implemented with personalized knowledge of the bacterial colonization patterns in the patients, as well as the selection of appropriately high-risk patients prior to implementing any intervention. Prior investigations of the ICU microbiome have generally been restricted by the limitations of culture-based approaches, which can detect only certain species of bacteria and cannot provide a complete picture immediately during acute illness. The 16S rRNA gene sequence technique has allowed clinicians to perform sequence-based investigations of microbes previously thought to be inaccessible . This achievement may inspire the development of simple and effective microbiota-targeted therapies.
Our study has several strengths. This prospective observational study used 16S rRNA sequence analysis to identify members of the intestinal microbiota, exploring early samples after admission and longitudinal samples obtained weekly to provide dynamic observations. This research was the first to investigate the microbiota in neurocritically ill patients, with a correlation analysis between the abundances of taxa and clinical severity scores, as well as survival analyses based on the abundances of or changes in specific taxa and the rigorously adjudicated primary outcome of death. However, this research also has some inherent limitations that warrant further discussion. First, the window of our first sample enrolment was 72 h after admission, and the underlying changes in microbial diversity may occur during the period of the first sample collection. Although rectal swabs could obtain samples immediately and have been shown to provide similar efficacies as stool , faecal sampling was used in this study to avoid harmful effects on the patients. Therefore, the sample collection process was restricted to patients’ defecation cycle, and selection bias was present because faecal samples could be collected only from patients who defecated spontaneously. Second, the limited sample size and the substantial heterogeneity of the selected patients make the interpretation of the relationship between the microbiota composition and prognosis challenging. Third, the 16S rRNA gene sequence method used in this study could only evaluate the composition percentages but not the absolute numbers of microorganisms in the samples. An advanced technique is needed to define the gut microbiota more precisely. Finally, we could not continue to analyse faecal samples after discharge, and the recovery of the microbiome might be prolonged. We would like to evaluate post-discharge microbiota dynamics in future studies.
In summary, we identified significant alterations in the gut microbiota composition between neurocritically ill patients and healthy subjects as well as marked dynamic changes in the patients’ gut microbiomes during hospitalization. The abundance and its changes were important risk factors for 180-day mortality, indicating that a gut microbiome examination could be a useful tool for identifying high-risk patients. High-quality prospective investigations with substantial samples collected on the day of admission and, ideally, with concurrent experimental models, are warranted to confirm the findings of this research and evaluate the net effect of microbiota dysbiosis on the progression and outcome of neurocritical illness.
Availability of data and materials
The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.
Adjusted hazard ratio
Acute Physiology and Chronic Health Evaluation-II
Glasgow Coma Scale
Intensive care unit
Linear discriminant analysis
Linear discriminant analysis effect size
Modified Rankin Scale
Neurological intensive care unit
Operational taxonomic unit
Principal coordinate analysis
Quantitative Insights into Microbial Ecology
Quantitative real-time polymerase chain reaction
Sequential Organ Failure Assessment
Moris G, Arboleya S, Mancabelli L, Milani C, Ventura M, de Los Reyes-Gavilan CG, Gueimonde M. Fecal microbiota profile in a group of myasthenia gravis patients. Sci Rep. 2018;8(1):14384.
Barichella M, Severgnini M, Cilia R, Cassani E, Bolliri C, Caronni S, Ferri V, Cancello R, Ceccarani C, Faierman S, et al. Unraveling gut microbiota in Parkinson’s disease and atypical parkinsonism. Mov Disord. 2019;34(3):396–405.
McDonald D, Ackermann G, Khailova L, Baird C, Heyland D, Kozar R, Lemieux M, Derenski K, King J, Vis-Kampen C, et al. Extreme dysbiosis of the microbiome in critical illness. mSphere. 2016;1(4):e00199–16.
Clark JA, Coopersmith CM. Intestinal crosstalk: a new paradigm for understanding the gut as the “motor” of critical illness. Shock. 2007;28(4):384–93.
Taur Y, Jenq RR, Perales MA, Littmann ER, Morjaria S, Ling L, No D, Gobourne A, Viale A, Dahi PB, et al. The effects of intestinal tract bacterial diversity on mortality following allogeneic hematopoietic stem cell transplantation. Blood. 2014;124(7):1174–82.
Shimizu K, Ogura H, Hamasaki T, Goto M, Tasaki O, Asahara T, Nomoto K, Morotomi M, Matsushima A, Kuwagata Y, et al. Altered gut flora are associated with septic complications and death in critically ill patients with systemic inflammatory response syndrome. Digest Dis Sci. 2011;56(4):1171–7.
Jung E, Byun S, Lee H, Moon SY, Lee H. Vancomycin-resistant Enterococcus colonization in the intensive care unit: clinical outcomes and attributable costs of hospitalization. Am J Infect Control. 2014;42(10):1062–6.
Taur Y, Xavier JB, Lipuma L, Ubeda C, Goldberg J, Gobourne A, Lee YJ, Dubin KA, Socci ND, Viale A, et al. Intestinal domination and the risk of bacteremia in patients undergoing allogeneic hematopoietic stem cell transplantation. Clin Infect Dis. 2012;55(7):905–14.
Freedberg DE, Zhou MJ, Cohen ME, Annavajhala MK, Khan S, Moscoso DI, Brooks C, Whittier S, Chong DH, Uhlemann AC, et al. Pathogen colonization of the gastrointestinal microbiome at intensive care unit admission and risk for subsequent death or infection. Intensive Care Med. 2018;44(8):1203–11.
Korbakis G, Bleck T. The evolution of neurocritical care. Crit Care Clin. 2014;30(4):657–71.
Majid HA, Cole J, Emery PW, Whelan K. Additional oligofructose/inulin does not increase faecal bifidobacteria in critically ill patients receiving enteral nutrition: a randomised controlled trial. Clin Nutr. 2014;33(6):966–72.
Petrof EO, Dhaliwal R, Manzanares W, Johnstone J, Cook D, Heyland DK. Probiotics in the critically ill: a systematic review of the randomized trial evidence. Crit Care Med. 2012;40(12):3290–302.
van Nood E, Vrieze A, Nieuwdorp M, Fuentes S, Zoetendal EG, de Vos WM, Visser CE, Kuijper EJ, Bartelsman JF, Tijssen JG, et al. Duodenal infusion of donor feces for recurrent Clostridium difficile. N Engl J Med. 2013;368(5):407–15.
Wischmeyer PE, McDonald D, Knight R. Role of the microbiome, probiotics, and ‘dysbiosis therapy’ in critical illness. Curr Opin Crit Care. 2016;22(4):347–53.
Debast SB, Bauer MP, Kuijper EJ. European Society of Clinical Microbiology and Infectious Diseases. European Society of Clinical Microbiology and Infectious Diseases: update of the treatment guidance document for Clostridium difficile infection. Clin Microbiol Infect. 2014;20(Suppl 2):1–26.
Shimizu K, Yamada T, Ogura H, Mohri T, Kiguchi T, Fujimi S, Asahara T, Yamada T, Ojima M, Ikeda M, et al. Synbiotics modulate gut microbiota and reduce enteritis and ventilator-associated pneumonia in patients with sepsis: a randomized controlled trial. Crit Care. 2018;22(1):239.
Peng X, Yu KQ, Deng GH, Jiang YX, Wang Y, Zhang GX, Zhou HW. Comparison of direct boiling method with commercial kits for extracting fecal microbiome DNA by Illumina sequencing of 16S rRNA tags. J Microbiol Methods. 2013;95(3):455–62.
He Y, Wu W, Wu S, Zheng HM, Li P, Sheng HF, Chen MX, Chen ZH, Ji GY, Zheng ZD, et al. Linking gut microbiota, metabolic syndrome and economic status based on a population-level analysis. Microbiome. 2018;6(1):172.
He Y, Wu W, Zheng HM, Li P, McDonald D, Sheng HF, Chen MX, Chen ZH, Ji GY, Zheng ZD, et al. Regional variation limits applications of healthy gut microbiome reference ranges and disease models. Nat Med. 2018;24(10):1532–5.
Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, Fierer N, Pena AG, Goodrich JK, Gordon JI, et al. QIIME allows analysis of high-throughput community sequencing data. Nat Methods. 2010;7(5):335–6.
Lozupone C, Hamady M, Knight R. UniFrac--an online tool for comparing microbial community diversity in a phylogenetic context. BMC Bioinformatics. 2006;7:371.
Lozupone C, Knight R. UniFrac: a new phylogenetic method for comparing microbial communities. Appl Environ Microbiol. 2005;71(12):8228–35.
Segata N, Izard J, Waldron L, Gevers D, Miropolsky L, Garrett WS, Huttenhower C. Metagenomic biomarker discovery and explanation. Genome Biol. 2011;12(6):R60.
Yin J, Liao SX, He Y, Wang S, Xia GH, Liu FT, Zhu JJ, You C, Chen Q, Zhou L, et al. Dysbiosis of gut microbiota with reduced trimethylamine-N-oxide level in patients with large-artery atherosclerotic stroke or transient ischemic attack. J Am Heart Assoc. 2015;4(11):e002699.
Lankelma JM, van Vught LA, Belzer C, Schultz MJ, van der Poll T, de Vos WM, Wiersinga WJ. Critically ill patients demonstrate large interpersonal variation in intestinal microbiota dysregulation: a pilot study. Intensive Care Med. 2017;43(1):59–68.
Rogers MB, Firek B, Shi M, Yeh A, Brower-Sinning R, Aveson V, Kohl BL, Fabio A, Carcillo JA, Morowitz MJ. Disruption of the microbiota across multiple body sites in critically ill children. Microbiome. 2016;4(1):66.
Tremaroli V, Backhed F. Functional interactions between the gut microbiota and host metabolism. Nature. 2012;489(7415):242–9.
Mayer EA, Tillisch K, Gupta A. Gut/brain axis and the microbiota. J Clin Invest. 2015;125(3):926–38.
Victor DW, Quigley EM. Hepatic encephalopathy involves interactions among the microbiota, gut, brain. Clin Gastroenterol Hepatol. 2014;12(6):1009–11.
Trompette A, Gollwitzer ES, Yadava K, Sichelstiel AK, Sprenger N, Ngom-Bru C, Blanchard C, Junt T, Nicod LP, Harris NL, et al. Gut microbiota metabolism of dietary fiber influences allergic airway disease and hematopoiesis. Nat Med. 2014;20(2):159–66.
Han YW, Wang X. Mobile microbiome: oral bacteria in extra-oral infections and inflammation. J Dent Res. 2013;92(6):485–91.
Surana NK, Kasper DL. Deciphering the tete-a-tete between the microbiota and the immune system. J Clin Invest. 2014;124(10):4197–203.
Zaborin A, Smith D, Garfield K, Quensen J, Shakhsheer B, Kade M, Tirrell M, Tiedje J, Gilbert JA, Zaborina O, et al. Membership and behavior of ultra-low-diversity pathogen communities present in the gut of humans during prolonged critical illness. MBio. 2014;5(5):e01361–14.
Panzer AR, Lynch SV, Langelier C, Christie JD, McCauley K, Nelson M, Cheung CK, Benowitz NL, Cohen MJ, Calfee CS. Lung microbiota is related to smoking status and to development of acute respiratory distress syndrome in critically ill trauma patients. Am J Respir Crit Care Med. 2018;197(5):621–31.
Wan YD, Zhu RX, Wu ZQ, Lyu SY, Zhao LX, Du ZJ, Pan XT. Gut microbiota disruption in septic shock patients: a pilot study. Med Sci Monit. 2018;24:8639–46.
Xu KY, Xia GH, Lu JQ, Chen MX, Zhen X, Wang S, You C, Nie J, Zhou HW, Yin J. Impaired renal function and dysbiosis of gut microbiota contribute to increased trimethylamine-N-oxide in chronic kidney disease patients. Sci Rep. 2017;7(1):1445.
Zeng X, Gao X, Peng Y, Wu Q, Zhu J, Tan C, Xia G, You C, Xu R, Pan S, et al. Higher risk of stroke is correlated with increased opportunistic pathogen load and reduced levels of butyrate-producing bacteria in the gut. Front Cell Infect Microbiol. 2019;9:4.
Iapichino G, Callegari ML, Marzorati S, Cigada M, Corbella D, Ferrari S, Morelli L. Impact of antibiotics on the gut microbiota of critically ill patients. J Med Microbiol. 2008;57(Pt 8:1007–14.
Rhee SH, Pothoulakis C, Mayer EA. Principles and clinical implications of the brain-gut-enteric microbiota axis. Nat Rev Gastroenterol Hepatol. 2009;6(5):306–14.
Vincent JL, Rello J, Marshall J, Silva E, Anzueto A, Martin CD, Moreno R, Lipman J, Gomersall C, Sakr Y, et al. International study of the prevalence and outcomes of infection in intensive care units. JAMA. 2009;302(21):2323–9.
Kvasnovsky CL, Leong LEX, Choo JM, Abell GCJ, Papagrigoriadis S, Bruce KD, Rogers GB. Clinical and symptom scores are significantly correlated with fecal microbiota features in patients with symptomatic uncomplicated diverticular disease: a pilot study. Eur J Gastroenterol Hepatol. 2018;30(1):107–12.
Kaakoush NO. Insights into the role of Erysipelotrichaceae in the human host. Front Cell Infect Microbiol. 2015;5:84.
Stenman LK, Burcelin R, Lahtinen S. Establishing a causal link between gut microbes, body weight gain and glucose metabolism in humans - towards treatment with probiotics. Benef Microbes. 2016;7(1):11–22.
Tremlett H, Fadrosh DW, Faruqi AA, Zhu F, Hart J, Roalstad S, Graves J, Lynch S, Waubant E, US Network of Pediatric MS Centers. Gut microbiota in early pediatric multiple sclerosis: a case-control study. Eur J Neurol. 2016;23(8):1308–21.
Gootenberg DB, Paer JM, Luevano JM, Kwon DS. HIV-associated changes in the enteric microbial community: potential role in loss of homeostasis and development of systemic inflammation. Curr Opin Infect Dis. 2017;30(1):31–43.
Shen F, Zheng RD, Sun XQ, Ding WJ, Wang XY, Fan JG. Gut microbiota dysbiosis in patients with non-alcoholic fatty liver disease. Hepatobiliary Pancreat Dis Int. 2017;16(4):375–81.
Strati F, Cavalieri D, Albanese D, De Felice C, Donati C, Hayek J, Jousson O, Leoncini S, Pindo M, Renzi D, et al. Altered gut microbiota in Rett syndrome. Microbiome. 2016;4(1):41.
Barbier F, Pommier C, Essaied W, Garrouste-Orgeas M, Schwebel C, Ruckly S, Dumenil AS, Lemiale V, Mourvillier B, Clec'h C, et al. Colonization and infection with extended-spectrum beta-lactamase-producing Enterobacteriaceae in ICU patients: what impact on outcomes and carbapenem exposure? J Antimicrob Chemother. 2016;71(4):1088–97.
Ferrer M, Difrancesco LF, Liapikou A, Rinaudo M, Carbonara M, Li Bassi G, Gabarrus A, Torres A. Polymicrobial intensive care unit-acquired pneumonia: prevalence, microbiology and outcome. Crit Care. 2015;19:450.
Kale IT, Kuzu MA, Berkem H, Berkem R, Acar N. The presence of hemorrhagic shock increases the rate of bacterial translocation in blunt abdominal trauma. J Trauma. 1998;44(1):171–4.
Yao YM, Sheng ZY, Tian HM, Wang YP, Yu Y, Fu XB, Lu LR, Wang DW, Ma YY. Gut-derived endotoxemia and multiple system organ failure following gunshot wounds combined with hemorrhagic shock: an experimental study in the dog. J Trauma. 1995;38(5):742–6.
Mittal R, Coopersmith CM. Redefining the gut as the motor of critical illness. Trends Mol Med. 2014;20(4):214–23.
Li Q, Wang C, Tang C, He Q, Zhao X, Li N, Li J. Therapeutic modulation and reestablishment of the intestinal microbiota with fecal microbiota transplantation resolves sepsis and diarrhea in a patient. Am J Gastroenterol. 2014;109(11):1832–4.
Li Q, Wang C, Tang C, He Q, Zhao X, Li N, Li J. Successful treatment of severe sepsis and diarrhea after vagotomy utilizing fecal microbiota transplantation: a case report. Crit Care. 2015;19:37.
Arumugam M, Raes J, Pelletier E, Le Paslier D, Yamada T, Mende DR, Fernandes GR, Tap J, Bruls T, Batto JM, et al. Enterotypes of the human gut microbiome. Nature. 2011;473(7346):174–80.
Bansal S, Nguyen JP, Leligdowicz A, Zhang Y, Kain KC, Ricciuto DR, Coburn B. Rectal and naris swabs: practical and informative samples for analyzing the microbiota of critically ill patients. mSphere. 2018;3(3):e00219–8.
The authors thank all of the study participants and the clinical and research staff from Nanfang Hospital and Zhujiang Hospital for their contributions to this study. We thank Dr. Hao Wang (Department of General Surgery, Nanfang Hospital, Southern Medical University, Guangzhou, China) for his assistance with the statistical analysis.
This study was supported by the National Natural Science Foundation of China (NSFC81671171, NSFC81870936, NSFC81800746), the Guangdong Natural Science Foundation (2017A030313821), China Postdoctoral Science Foundation (2018M633077) and the Clinical Research Startup Program of Southern Medical University by High-level University Construction Funding of the Guangdong Provincial Department of Education (LC2016PY025).
Ethics approval and consent to participate
This study was approved by the Medical Ethics Committee of Nanfang Hospital (No. NFEC-2018-034). Informed consent was obtained from all healthy subjects and patients or their legal representatives. This research was conducted in accordance with the Declaration of Helsinki.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Figure S1. PCoA plot illustrating the grouping patterns of the samples collected from patients with different primary diagnoses at admission. Figure S2. The gut microbiota composition of neuroICU patients was significantly different from that of 84 healthy subjects. Figure S3. Dynamic changes in α-diversities of samples from seven patients with neuroICU stay lengths of greater than six weeks. Figure S4. The gut microbiota composition of 58 patients with critically ill stroke was significantly different from that of 58 healthy subjects. (DOCX 1152 kb)
Table S1. Baseline characteristics of the patients from the neuroICU and the healthy subjects. Table S2. Comparison of β-diversity between patients with different primary diagnoses and healthy controls. Table S3. Baseline characteristics of patients with critically ill stroke and the healthy subjects. (DOCX 25 kb)
Clinical information of all patients in the neuroICU group in this study. (XLSX 25 kb)
α-Diversities and taxa of all samples collected from neuroICU and HCs in this study. (XLSX 922 kb)