Academic Editor

Article Metrics

  • Fig. 1.

    View in Article
    Full Image
  • Fig. 2.

    View in Article
    Full Image
  • Fig. 3.

    View in Article
    Full Image
  • Fig. 4.

    View in Article
    Full Image
  • Fig. 5.

    View in Article
    Full Image
  • Fig. 6.

    View in Article
    Full Image
  • Information

  • Download

  • Contents

Abstract

Background:

The subjective limitations of neurobehavioral assessment cause a high misdiagnosis rate for disorders of consciousness (DoC). The purpose of this study was to identify the DoC level based on an analysis of multi-dimensional electroencephalogram (EEG) signals to assist with establishing a clinical diagnosis.

Methods:

Sixty-seven patients with DoC [coma, n = 19; vegetative state (VS), n = 23; and minimally conscious state (MCS), n = 25] were included to analyze resting state EEG characteristics. The EEG features were statistically compared among five band powers (delta, theta, alpha, beta, and gamma) and five brain regions (prefrontal, frontal, parietal, temporal, and occipital) by multidimensional analyses, including time-domain analysis, spectral analysis, and functional brain connectivity.

Results:

Amplitude-integrated electroencephalography (aEEG) center amplitude showed significant differences between coma and MCS (p = 0.02688), with no significant differences observed for the other comparison. Spectral analysis revealed that delta and theta power decreased with higher consciousness levels, whereas alpha, beta, and gamma power increased. Relative power differed among groups across specific brain regions (prefrontal, frontal, parietal, temporal, and occipital) and frequency bands. Weighted Phase Lag Index (wPLI) based functional connectivity demonstrated frequency-specific network reorganization with theta band connectivity strongest in VS and alpha/beta/gamma band connectivity enhanced in MCS. Absolute power topographic maps showed expanding high-power regions from coma-to-MCS in high-frequency bands and the left dorsolateral prefrontal cortex (DLPFC) (F3 electrode) exhibited a consistent power gradient of coma < MCS < VS across all bands.

Conclusions:

Multidimensional EEG features have significant value in differentiating the levels of consciousness disorders. aEEG center amplitude discriminated MCS from coma; delta/gamma relative power separated VS from MCS, and alpha/beta relative power separated coma, VS, and MCS. Parieto-occipital connectivity matrix in the theta band distinguishes coma from VS, while absolute power topography of the left DLPFC shows potential for grading levels of impaired consciousness. These electrophysiologic biomarkers complement behavioral assessments, enhancing diagnostic accuracy.

1. Introduction

Disorders of consciousness (DoC) refer to a spectrum of conditions characterized by impaired levels of consciousness due to various underlying causes. DoC are primarily classified into coma, vegetative state (VS), and minimally conscious state (MCS) [1]. Despite advances in clinical practice, accurately determining the specific level of consciousness in DoC patients remains a significant challenge. The precise assessment and differentiation between these states are crucial for guiding diagnosis and treatment but continue to present difficulties in the management of DoC patients.

Several standardized assessment methods have been developed specifically for evaluating DoC patients, among which the Coma Recovery Scale-Revised (CRS-R) is widely recognized for its practicality and effectiveness. The CRS-R is particularly sensitive in detecting the MCS, making the CRS-R a valuable tool in the assessment of DoC patients. However, the reliability of these neurobehavioral measures can be compromised due to the subjective nature of the evaluation process, resulting in variable accuracy of the assessment outcomes [2]. Therefore, while neurobehavioral scales are useful for diagnosing DoC, neurobehavioral scales are not entirely dependable.

Neuroimaging techniques, such as functional Magnetic Resonance Imaging (fMRI) and Positron Emission Tomography (PET), provide objective assessment of DoC by probing underlying brain function, detecting covert awareness, and revealing behaviorally undetectable neural correlates of consciousness [3, 4]. However, routine clinical use of neuroimaging techniques is significantly limited by high costs, restricted accessibility, technical demands (e.g., patient transport and MRI compatibility), and challenges in scanning critically ill patients [5]. In contrast, electroencephalogram (EEG)-based approaches offer substantial advantages, including bedside monitoring of spontaneous brain activity, delivering high temporal resolution, portability, repeatability, and lower cost [6]. Advances in EEG analysis enhance diagnostic accuracy and covert consciousness detection [7, 8], making EEG a critical tool for complementing behavioral and neuroimaging assessments. Within the EEG methodologies, amplitude-integrated EEG (aEEG) presents a streamlined solution. aEEG compresses raw EEG data into min-max amplitude bands on a semi-log scale, enabling rapid visual trend assessment [9]. Although originally developed for neonatal monitoring [9], the artifact tolerance, simplicity, and suitability of aEEG for continuous monitoring facilitate effective translation to adult critical care [10], for which aEEG is gradually being applied.

EEG signals are classified into low-frequency bands, including delta, theta, and alpha, as well as high-frequency bands, such as beta and gamma. The brain activity state is assessed based on the amplitude and fluctuations of EEG waveforms. Under normal awake conditions EEG amplitudes and waveforms exhibit regular patterns. However, the spectral power across these frequency bands is significantly disrupted in patients with DoC [11]. Klimesch [12] reported that variations in alpha and theta reflect cognitive function and memory capabilities. Patients with DoC typically display enhanced delta and suppressed alpha compared to healthy individuals [5]. Naro et al. [13] further noted that patients with MCS exhibit higher power in theta and alpha than patients with VS in whom alpha wave activity is closely associated with consciousness recovery. Specifically, alpha power in the parietal and occipital regions is significantly increased in patients with MCS [13, 14]. In addition, Piarulli et al. observed that beta power is notably higher in patients with MCS compared to patients with VS. However, research on high-frequency bands, such as beta and gamma in DoC, is limited [15].

Consciousness-related information in the brain is often complex and multifaceted. The spectral power of EEG and the amplitude changes in specific frequency bands reflect the differences and correlations between various states of consciousness. Numerous studies have examined the relationship between EEG and consciousness states and shown that consciousness is not simply a “present or absent” condition. EEG-based neurofunctional imaging techniques offer an innovative approach to tackling this challenge. Cavinato et al. [16] observed that patients in a VS exhibit impaired cortical integration with significant disruption of the associations between the frontal and parietal lobes. In contrast, Wu et al. [17] reported that patients in a MCS exhibit superior functional connectivity between brain regions with enhanced cortical peripheral connectivity compared to patients with VS [16, 17, 18, 19]. Increasing evidence suggests that resting-state EEG coherence serves as a robust tool for quantifying the relationship between conscious awareness and cerebral information processing and provides more granular insight into differentiating between various levels of consciousness [20, 21]. Patients who regained consciousness demonstrated higher levels of cortical functional connectivity in quantity and strength in a clinical trial assessing the prognosis of patients with DoC over a 3-month period [22]. Schorr et al. [23] proposed that the coherence between frontal and parietal regions serve as a predictive marker for dynamic shifts in consciousness states, such as recovery from VS-to-MCS with delta and theta in the parietal cortex showing high sensitivity and specificity in forecasting recovery in patients with VS.

It is indisputable that significant challenges remain in the diagnosis and management of DoC, primarily due to an insufficient understanding of the neurophysiologic mechanisms underlying consciousness. Therefore, a comprehensive analysis of EEG signals was performed in DoC patients with a focus on the whole-brain power spectral distribution, amplitude fluctuations within distinct frequency bands, and the functional state of brain connectivity in the resting state. By integrating these brain-machine data, a more cohesive model of the consciousness network within the brain was constructed. Furthermore, the most valuable EEG features related to consciousness were identified, thereby advancing the use of resting-state EEG in DoC. This approach will enhance the accuracy of diagnosis, enable better monitoring of therapeutic interventions, and facilitate the prediction of recovery in DoC patients.

2. Materials and Methods
2.1 Subjects

This study recruited a total of 90 patients with DoC from the Department of Rehabilitation Medicine of the Affiliated Hospital of Nantong University (Jiangsu, China). Due to issues, such as poor data quality and changes in clinical condition, 23 patients withdrew during the study. Ultimately, 67 patients were included in the analysis (19 with coma, 23 with VS, and 25 with MCS). The baseline characteristics of the sample are listed in Table 1.

Table 1. Baseline characteristics of the sample.
Variables Coma VS MCS p
n = 19 n = 23 n = 25
Gender, n (%) 0.105
Male 9 (47.4%) 12 (52.2%) 19 (76.0%)
Female 10 (52.6%) 11 (47.8%) 6 (24.0%)
Age, years 0.944
Median (Q𝑙, Q𝑢) 64 (58.0, 71.0) 64 (56.0, 69.0) 64 (48.5, 74.5)
Time since injury, months 0.099
Median (Q𝑙, Q𝑢) 2.0 (1.3, 2.4) 1.4 (1.1, 2.3) 1.3 (1.0, 1.9)
Etiology 0.512
NTBI 10 (52.6%) 11 (47.8%) 9 (36.0%)
TBI 9 (47.4%) 12 (52.2%) 16 (64.0%)
Behavioral scale assessment
Total score of CRS-R 2 (1, 3) 9 (7, 9) 12 (11, 14) <0.001
Total score of GCS 4 (3, 5) 6 (4, 7) 11 (10, 11) <0.001

VS, vegetative state; MCS, minimally conscious state; NTBI, non-traumatic brain injury; TBI, traumatic brain injury; CRS-R, coma recovery scale-revised; GCS, glasgow coma scale.

All enrolled patients met the following inclusion criteria: (1) met the diagnostic criteria for DoC [6, 24]; (2) >18 years of age; (3) DoC caused by various reasons, including both traumatic brain injury (TBI) and non-traumatic brain injury (NTBI). NTBI etiologies encompassed hypoxic-ischemic encephalopathy, cerebrovascular stroke, intracranial infection (e.g., encephalitis), and metabolic/toxic encephalopathy; and (4) informed consent was obtained from the subjects’ legal guardian. The exclusion criteria were as follows: (1) history of neurologic or psychiatric illness; (2) presence of metal implants or implanted electronic brain medical devices, such as pacemakers; (3) history of cranioplasty with unhealed skull flap or severe scalp trauma, inability to place electrodes, or possible infection due to examination; (4) exposure to centrally acting medications within 24 h prior to EEG acquisition or within five drug half-lives, whichever was longer, that are known to significantly alter EEG background activity or functional connectivity, including sedatives/anesthetics (e.g., propofol, benzodiazepines, dexmedetomidine), opioids, antipsychotics, antidepressants, antiepileptic drugs, or pharmacologic arousal agents (e.g., amantadine or zolpidem), as well as continuous infusions or regimens that were newly initiated or dose-adjusted during this period; (5) patients with combined multi-system and -organ dysfunction; (6) unstable vital signs (patients with a systolic BP <90 or >180 mmHg, heart rate <40 or >120 bpm, respiratory rate <8/min or ventilator-dependent, or SpO2 <90% within 24 h pre-EEG were excluded); and (7) excessive motion artifacts (EEG epochs with an amplitude >75 µV, application of independent component analysis (ICA), or visual inspection by two technicians as motion artifacts).

The study was conducted under the Declaration of Helsinki of the World Medical Association and approved by the Ethics Committee of Nantong University Affiliated Hospital on 20 July 2024. The approval number for the study was 2024-K142-01. Before inclusion, the researcher fully informed the legal guardian of each patient about the study protocol and obtained informed consent, especially those in a coma state.

2.2 Clinical Assessment

Patients were systematically evaluated using the CRS-R and Glasgow Coma Scale (GCS) after a preliminary diagnosis of DoC [25, 26, 27]. Two trained therapists independently performed all assessments with discrepancies resolved through consensus. Three separate CRS-R and GCS evaluations were performed to accurately determine the patient’s level of consciousness at least 5 d prior to the EEG assessment.

CRS-R systematically evaluates 6 functional domains (auditory, visual, motor, oromotor/verbal, communication, and arousal) using a 23-point hierarchical scale. Lower subscale scores suggest brainstem-mediated reflexive responses, while higher scores imply cortically-driven cognitive processing. As the gold-standard behavioral assessment for DoC, CRS-R demonstrates superior sensitivity in differentiating unresponsive wakefulness syndrome (UWS) from a MCS by detecting reproducible neurobehavioral markers, such as visual tracking [2].

GCS quantifies acute cerebral dysfunction through triaxial responsiveness, as follows: ocular (1–4); verbal (1–5); and motor (1–6). Total scores stratify injury severity, as follows: 15 (intact); 13–14 (mild); 9–12 (moderate); and 8 (coma/severe dysfunction) [28, 29]. Serial GCS monitoring (hourly when 8) predicts secondary deterioration via a 2-point cumulative decline [29]. The verbal component is annotated as “T” with emphasis on non-verbal consciousness markers (e.g., visual tracking) for intubated patients. Assessments require exclusion of confounders (sedatives and hypoglycemia) [29]. Notably, GCS lacks specificity for chronic DoC subtypes, necessitating CRS-R supplementation to detect covert awareness [30].

2.3 EEG Recordings and Processing

Resting-state EEG data were collected from patients with three types of DoC. A 32-channel EEG (ZhenTec-NT1, ZhenTec Intelligent Technology Co., Ltd., Xi’an, Shaanxi, China) cap from a brain-computer interface (BCI) system was utilized. The electrodes were positioned in accordance with the internationally recognized 10–10 system [30]. The impedance of all electrodes was maintained below 10 kΩ. Each patient underwent a minimum of a 10-min resting-state EEG recording to facilitate the analysis of key EEG characteristics.

The pre-selected EEG data were pre-processed using MATLAB 2021a (The MathWorks Inc, Natick, MA, USA) and the EEGLAB toolbox (version 2022.1; Swartz Center for Computational Neuroscience, University of California San Diego, San Diego, CA, USA). The pre-processing steps were as follows: (1) electrode coordinate file was loaded and mapped to the corresponding electrode positions on the EEG cap (Gelfree-S3, Wuhan Greentek Pty. Ltd., Wuhan, Hubei, China); (2) a notch filter was applied using the ERPLAB filter to eliminate a 50-Hz power line interference; (3) band-pass filtering was performed between 1 and 45 Hz using a sixth-order Butterworth filter [31, 32, 33]; (4) Epochs with evident electromyographic or motion artifacts and non-EEG electrodes (e.g., electrooculogram) were excluded; (5) application of independent component analysis (ICA) to identify and remove the components related to physiologic artifacts, such as eye movements and electromyography; and (6) the average results of 30-channel electroencephalogram data.

2.4 Data Analysis
2.4.1 aEEG Center Amplitude (Median of Upper/Lower Margins)

The aEEG assesses brain electrical activity by real-time monitoring of EEG amplitude characteristics. This method typically involves computing the mean amplitude within specific EEG frequency bands and analyzing waveform trends. Raw EEG was band-pass filtered to the aEEG bandwidth (≈ 2–15 Hz) and time-compressed to generate the amplitude-integrated trace following standard methodology. We extracted the upper envelope (high amplitude) and lower envelope (low amplitude) of the EEG signal. Within each subject, the center amplitude was defined per channel as follows:

center = upper + lower 2

Where, upper and lower denote the time-median values of the aEEG margins over the full recording. We then aggregated across the 30 channels by taking the within-subject median, yielding a single whole-brain aEEG-center value per subject for group comparisons (Coma, VS, MCS) [34, 35, 36]. The primary analysis used Kruskal–Wallis test for independent groups with Dunn post hoc comparisons (all three pairwise contrasts).

2.4.2 Power Spectral Analysis

Power spectral analysis, which is a fundamental methodology for decomposing neural signals into constituent frequency components and quantifying the energy distribution, utilizes the fast Fourier transform (FFT) to convert EEG signals from the temporal domain to the frequency domain. Subsequently, squared amplitude values were calculated to generate power spectra. To enhance visual interpretability, power values are typically normalized using logarithmic decibel (dB) scaling via the following formula: dB = 10log10(power/reference). Spectral profiles systematically reveal energy distribution patterns across canonical frequency bands in electroencephalographic investigations, providing quantitative insights into cerebral oscillatory dynamics. Notably, healthy individuals typically exhibit prominent α-band power (8–13 Hz) over occipital regions during eyes-closed rest. However, individual variations exist. In contrast, consciousness-impaired patients frequently demonstrate marked α-power attenuation [30]. Analytical outcomes are visualized as spectrograms with dual-axis representations, as follows: frequency range (0–45 Hz) on the abscissa; and power intensity (dB) on the ordinate. Comparative analysis of spectral signatures across clinical states enables identification of consciousness-related pathologic patterns, including elevated θ/β power ratios and disrupted cross-frequency coupling profiles [7]. These spectral biomarkers complement multimodal assessments of DoC by objectively characterizing neurophysiologic aberrations in both resting-state and task-evoked conditions.

2.4.3 Whole-Brain Relative Power Distribution and Absolute Power Topographic Analysis

Relative power analysis quantitatively assesses cerebral electrophysiologic activity by evaluating the proportional contribution of specific frequency bands to total spectral power. Relative power analysis quantifies cerebral electrophysiologic activity by calculating the power spectral density (PSD) for each frequency band via Welch’s periodogram, then normalizing the band power with respect to the total power to obtain relative power ratios. This normalization reduces inter-subject variability due to anatomic and technical factors with results often visualized as cortical heat maps, facilitating reliable cross-cohort spectral comparisons.

Furthermore, absolute power topographic mapping visualizes spatial-frequency abnormalities in DoC via 2D scalp projections derived from 3D data. The computational pipeline includes the following elements: (1) FFT-based PSD extraction; (2) spectral summation within clinically relevant bands (δ-γ); (3) dB conversion (dB = 10log10(P/Pref), where Pref denotes baseline power); (4) spatial interpolation (bilinear or radial basis function algorithms) for continuous scalp projections; and (5) color-coded visualization using standardized jet color maps. These biomarkers provided spatially resolved indices that may reveal DoC-related patterns with the EEGLAB topoplot function. This multimodal approach could enhance sensitivity in consciousness-level stratification and neuromodulation response assessment when integrated with aEEG metrics.

2.4.4 Weighted Phase Lag Index (wPLI)

The wPLI is an extension of the Phase Lag Index (PLI) that integrates a weighting mechanism. This modification provides several potential benefits, such as improved statistical sensitivity for detecting changes in consistent phase lead/lag between signals, enhanced robustness to additional noise sources, and a reduced impact of small phase differences. Consequently, wPLI may provide a more reliable estimate of the coupling strength of neural oscillations in the brain. The calculation method is outlined as follows:

wPLI = | ( X ) | | ( X ) | = | | ( X ) | sign ( ( X ) ) | ( X ) |

where, (X) represents the imaginary component of the cross-spectrum between x(t) and y(t). The value of wPLI ranges from 0–1. Values approaching 1 reflect a stable directional phase lead/lag, whereas values near zero indicate absence of a consistent phase-lagged relationship.

Functional connectivity metrics, including coherence and the wPLI, offer complementary approaches for mapping neural interactions. Coherence captures linear signal synchrony but is confounded by volume conduction artifacts. In contrast, wPLI improves specificity to phase-lagged coupling by weighting phase differences, which effectively suppresses the spurious synchronization arising from volume conduction [37]. In graph-theoretical analyses these metrics are used to define the edge weights of networks where nodes represent brain regions. The resulting connection strengths can be visualized as spatial patterns, aiding the identification of neuropathologic network reorganization. However, interpretation of these findings requires addressing methodologic constraints, including threshold selection biases, signal non-stationarity in wPLI estimation, and integration with structural connectivity data to resolve neurophysiologic ambiguities.

2.4.5 Statistical Analysis

Statistical analyses were performed using SPSS (Version 26.0, IBM Corp., Armonk, NY, USA) and MATLAB (R2021a, The MathWorks Inc., Natick, MA, USA) with an α = 0.05. Categorical variables were summarized as a frequency (percentage). Group comparisons for categorical variables were made using a chi-square test. For continuous variables data with a normal distribution were analyzed by one-way analysis of variance (ANOVA) and expressed as the mean ± standard deviation; otherwise, the Kruskal–Wallis H test was used, with data presented as a median (IQR). If the ANOVA was significant, the Tukey-Kramer post hoc test was performed. Considering that multiple comparisons may increase type I errors, the Benjamini–Hochberg (BH) FDR was used to correct the main effect p-values and report the q-values (significance threshold q <0.05) (Supplementary Tables 1,2).

3. Results
3.1 Comparison of aEEG Center Amplitude and Average EEG Power Spectra

The aEEG system quantifies the functional state of the brain by extracting the upper envelope (high amplitude) and lower envelope (low amplitude) of the EEG signal, then calculating the median values of upper- and lower-margin amplitudes [38]. The Dunn’s post hoc test was performed between the groups (Coma–VS, Coma–MCS, and VS–MCS) when the Kruskal–Wallis test showed a significant result. The Kruskal–Wallis test was significant (p = 0.02228), and Dunn’s post hoc test identified a significant Coma–MCS difference (p = 0.02688), with Coma–VS and VS–MCS remaining non-significant (p = 0.6096 and p = 0.3649, respectively) (Fig. 1A). Collectively, the between-group effect was primarily driven by the Coma vs. MCS contrast, with no significant differences observed for the other pairings. These findings indicated that while aEEG center amplitude can effectively distinguish coma from MCS, it may lack sufficient sensitivity to discern not only the distinctions between coma and VS, but also the more subtle ones between VS and MCS.

Fig. 1.

Power spectrum and aEEG center amplitude of patients with DoC. The yellow, blue, and red columns in (A) correspond to the coma, VS, and MCS groups, respectively. The Dunn’s post hoc test was used for intergroup comparison. **p < 0.01 for intergroup comparison. (B) shows the comparison between the power spectrum analysis results of EEG data among the three groups and the dashed lines represent the boundaries of the standardized frequency bands. Green is for coma, blue is for VS, and red is for MCS. aEEG, amplitude-integrated electroencephalography; DoC, disorders of consciousness.

Systematic analysis of EEG power spectral characteristics in patients with DOC revealed distinct energy distribution patterns across frequency bands (Fig. 1B). The absolute power exhibited a descending gradient (coma > VS > MCS) in the delta band with the coma group showing the highest power and the MCS group showing the lowest power. Conversely, the coma group power decreased to the lowest level, while the MCS group power increased to the highest level in the theta band. The absolute power pattern reversed to MCS > VS > coma in the alpha, beta, and gamma bands.

3.2 Comparison of Relative Power Data of Each Brain Region

Relative power, a core quantitative indicator in EEG analysis, untangles the distribution characteristics of different band oscillatory energy within the total power spectrum [39] (Fig. 2; Supplementary Tables 1,2). Group differences in each frequency band were analyzed to clarify the frequency- and region-specific differences across DoC states, as follows. Significant differences were observed in the delta band between the VS and MCS groups across the entire brain. In addition, the coma and MCS groups exhibited significant differences within the frontal lobe. Significant differences were identified in the theta band between the coma and MCS groups in the parietal and temporal regions with differences between the VS and MCS groups (specifically in the parietal region). Significant differences were identified in the alpha band between the coma and MCS groups and between the VS and MCS groups in the prefrontal cortex. Furthermore, all pairwise comparisons among the three groups (coma, VS, and MCS) had significant differences across all other brain regions. Significant differences were observed in the beta band between the coma and MCS groups and between the VS and MCS groups within the occipital lobe. Moreover, all pairwise comparisons among the coma, VS, and MCS groups demonstrated significant differences across all other brain regions. Significant differences were noted in the gamma band between the VS and MCS groups in the parietal lobe. Significant differences were observed in the occipital lobe between the coma and MCS groups and between the VS and MCS groups. Moreover, all pairwise comparisons among the three groups demonstrated significant differences across the prefrontal, frontal, and temporal cortical regions.

Fig. 2.

Average relative EEG power of patients with DoC. The average data of three groups of patients with respect to different brain regions and frequency bands were compared. In the histogram, yellow, blue, and red represented coma, VS, and MCS groups, respectively. The error line represents the standard deviation within the group. The Tukey–Kramer test was performed on the data between each two groups when the results of the ANOVA were significant. *p < 0.05 for intergroup comparison. ANOVA, one-way analysis of variance.

3.3 Connectivity Analysis of Brain Network Based on wPLI

Functional connectivity of the brain refers to the synergistic relationship between different brain regions, reflecting collaboration in processing information [40]. The wPLI was used to reveal functional interaction patterns between different brain regions and systematically analyze functional brain network connectivity patterns in different frequency bands in three groups of patients to quantitatively characterize functional brain network connectivity [41, 42] (Fig. 3).

Fig. 3.

wPLI-based functional connectivity matrices of brain regions across five frequency bands. Reflecting to the functional connectivity between different brain compartments, the nodes (Nodes) usually represent electrode positions arranged in a 10–10 system and the lines connecting the nodes are edges (Edges), which indicate the strength of connectivity as calculated by the wPLI. The deeper the red color, the stronger the intensity; the deeper the blue color, the weaker the intensity. The connection lines for intensities greater than the threshold were drawn using the 90th percentile as the threshold and marking it at the top of the image. wPLI, Weighted Phase Lag Index.

All three groups showed dense connectivity in the delta band within the occipital lobe. Dense connectivity was observed in the theta band within the temporal, frontal, and parietal lobes across the three groups. The connectivity strength in the coma and VS groups was significantly higher than the MCS group with the VS group exhibiting the highest connectivity strength among the three. Dense connectivity was present in the alpha band within the bilateral parietal and temporal lobes across the three groups. In addition, the coma and MCS groups showed dense connectivity between the bilateral parietal, and temporal lobes and the right prefrontal cortex. The three groups demonstrated a radial connectivity pattern centered on the central brain region in the beta band. A diffuse connectivity pattern across multiple brain regions was observed in the gamma band.

3.4 Connectivity Matrix Analysis of Different Brain Regions Based on wPLI Connections

The connectivity matrix of wPLI brain network connectivity indices for different brain regions reflects the connectivity characteristics and strengths of the patients under different bands (Fig. 4) [43]. The connection characteristics of each band were systematically compared to describe the network recombination of specific frequencies [44]. In addition, statistical matrices of brain networks were plotted for patients with DoC and different consciousness levels across each frequency band (Fig. 5).

Fig. 4.

Connectivity matrix of wPLI brain connections in different brain regions of patients with DoC. The connectivity matrix of wPLI brain connections reflects the connection characteristics and intensity in different frequency bands and the closer the value is to 1 (the closer the color is to red), the stronger the connection coherence between the two channels. Rows and columns represent brain regions in turn (prefrontal, frontal, temporal, parietal, and occipital).

Fig. 5.

Statistical analysis matrix of brain networks connectivity of patients with DoC. Edge-wise comparisons of functional connectivity indices were performed within each frequency band across three independent groups (coma, VS, and MCS). Connectivity strengths of the three groups were pooled (missing values were excluded) for each channel pair (upper triangle of the 30 × 30 matrix; 435 edges total) with each edge treated as an independent unit for intergroup testing (diagonal self-connections excluded). A one-way ANOVA was first applied to each edge to test the main group effect. Tukey–Kramer pairwise comparisons (coma–MCS, coma–VS, and VS–MCS) generated p-value matrix plots with a uniform color scale, overlaid brain region partition boundaries, and significant units marked by asterisks (*: 0.01 <p < 0.05; **: p < 0.01).

The overall connectivity strength among the three groups showed a gradually decreasing trend in the delta band (coma > VS > MCS) and there was a dense connectivity in the occipital lobe region. Dense connectivity was observed in the bilateral parietal lobes, bilateral occipital lobes, and between the parieto-occipital region and other brain regions across the three groups in the theta band with the VS group showing the strongest connectivity. Significant differences were detected between the coma and VS groups in the bilateral parietal lobes, bilateral occipital lobes, and between the parieto-occipital region and other brain regions as illustrated in Fig. 6. In addition, significant differences were identified between the coma and MCS groups in the connections between temporal parietal lobes, the parietal lobe and the prefrontal, frontal, and temporal lobes. The three groups exhibited a connection pattern in the alpha band, which was centered on the frontal lobe and the right prefrontal lobe with radiation to other brain regions. The three groups shared similar connectivity patterns in the beta band with prominent connectivity in the bilateral parietal lobes. The MCS group exhibited the highest overall connectivity strength among the three groups and showed significant differences from the coma group in the parietal lobe. Interestingly, the VS group demonstrated the highest overall connectivity strength in the gamma band.

Fig. 6.

Absolute power topographic maps of five rhythms in patients with DoC. Nodes typically represent electrode positions arranged in a 10–10 system. The deeper the red color, the higher the absolute power and the deeper the blue color, the lower the absolute power.

3.5 Distribution Characteristics of Absolute Power Topographic Maps in Different Frequency Bands

The absolute power topographic map utilizes spatial visualization techniques to untangle the regional distribution characteristics of EEG energy intensity within specific frequency bands [45]. The color gradient mapping intuitively reflects the absolute power of neural oscillations in different brain regions. Based on this method, the current study compared the whole-brain and regional absolute power topographic maps (Fig. 6) across various bands for the three groups of patients with DoC. Distinct patterns of absolute power were observed among the three patient groups across the frequency bands.

The absolute power topographies of different DoC states across five frequency bands revealed the following. There was no significant difference among the three groups in all five frequency bands. The high power of the three groups was mainly concentrated in the delta frequency band within the frontal lobe region. The high power of the coma group in the theta band was primarily localized to the prefrontal, bilateral frontal, and bilateral temporal regions. The VS and MCS groups demonstrated a generalized increase in absolute power compared to the coma group, which was accompanied by a more widespread distribution of these high-power areas. The elevated power of the coma group remained concentrated in the alpha band within the prefrontal area. The VS group showed enhanced activity extending to the prefrontal, bilateral frontal, bilateral temporal, and partial occipital regions. The MCS group, however, displayed a global increase in absolute power across the entire brain with the most extensive distribution of high-power regions. The distribution of high-power areas showed a gradually expanding pattern in the beta and gamma frequency bands from the coma, VS, and MCS groups. A pattern of increasing absolute power at the F3 electrode (left dorsolateral prefrontal cortex) was observed across multiple frequency bands, following the descriptive sequence: coma < MCS < VS.

4. Discussion

Recent studies have demonstrated that multidimensional EEG indicators can effectively distinguish patients’ levels of consciousness and provide an important basis for clinical assessment and prognosis. Our analysis revealed characteristic electrophysiological patterns across different states of impaired consciousness. A significant difference in aEEG center values was observed between coma and MCS patients. In terms of relative spectral power, significant differences were found in the δ and γ bands between VS and MCS patients, while relative power in the α and β bands showed significant differences among coma, VS, and MCS groups. Regarding functional connectivity, parieto-occipital connectivity matrices in the theta band exhibited more pronounced differences between VS and coma patients. Notably, in absolute power topography, activity over the left DLPFC (electrode F3) remained consistently elevated across all groups, following a gradient of coma < MCS < VS.

The results of the current study aligned closely with the existing literature and further refined the spectrum of EEG features associated with different levels of consciousness. Significant intergroup differences in quantitative scores were detected. The low amplitude range in the coma group indicated a significant decrease in brain activity levels, which reflected a severe disruption of brain network integration [46]. The “Theory of Consciousness Index” proposed by Casali et al. [47] suggests that changes in conscious states may be accompanied by reorganization of some neural networks. This reorganization may manifest as the partial preservation of cortical-thalamic connections. Notably, while pediatric EEG studies typically use limited-channel setups due to cranial size constraints, this adult-focused investigation utilized high-density 30-channel EEG [48]. This approach compensated for the spatial resolution limitations inherent in low-channel configurations [49], thereby enabling precise characterization of the dynamic reorganization of the prefrontal-parietal network [50]. Even if overall functional connectivity weakens, connections in some key areas may persist. The medium amplitude observed in the VS group may reflect this partial preservation of cortical-thalamic connections, which align with the Casali perspective [47]. The highest amplitude observed in the MCS group was associated with the enhanced functional connectivity of the default mode network (DMN) [30]. Patients in the MCS state may exhibit limited consciousness or responsiveness, demonstrating some degree of reaction to certain external stimuli, such as sound or touch [51]. The aEEG data from the three groups in the current study showed significant differences and exhibited a change corresponding to the level of consciousness. This finding indicated that aEEG may be a useful quantitative indicator for assessing levels of consciousness. The thalamo-cortical network functional integrity model (ABCD model) is an EEG-based power spectral density classification system that assesses the functional integrity of the thalamo-cortical network through EEG dynamic characteristics [52]. The ABCD model is used to describe neural activity patterns in different states of consciousness, as proposed by Curley et al. [52] in a 2022 study. Type A is characterized by dominance in the delta band, reflecting a complete lack of cortical input, and is commonly observed in patients with severe brain injuries that disrupt thalamocortical connections [53]. The dominance of delta band power in coma patients aligned with the “A-type power spectrum” described in the literature, which is characterized by a loss of integrity and reduced excitability, reflecting severe inhibition of brain function [52, 54]. In contrast, the preservation and enhancement of high-frequency activity (alpha/beta/gamma) in the MCS group supported the theory in which restoration of consciousness depended on high-frequency oscillations [7, 55]. In addition, the elevated power in the gamma frequency band within the MCS group may be linked to higher-order cognitive processing [56]. The aEEG and power spectrum characteristics aligned with the classic pattern of low-frequency dominance and high-frequency attenuation observed in patients with DoC.

The close correlation between the power distribution of different frequency bands in various brain regions and levels of consciousness was highly consistent with recent studies on the neural mechanisms of consciousness [7]. The coma group exhibited the highest relative power in the delta band across all brain regions among the three patient groups. This finding supported the classic conclusion that low-frequency oscillations serve as markers of conscious inhibition, suggesting that the pathomechanism may be related to the generalization of slow-wave activity due to disintegration of the thalamo-cortical loop. The difference in theta power in the parietal lobe suggests that this area has a crucial role in maintaining consciousness [57]. Engemann DA et al. [58] reported that alpha band power was one of the most prominently performing indicators for distinguishing MCS from UWS among more than 100 EEG recordings. Stefan et al. [14] reported the neural characteristics of enhanced alpha band activity and elevated beta coherence in patients with MCS through multidimensional EEG analysis, including power spectra, connectivity, and information entropy, which provided a reliable quantitative index for clinically differentiating between UWS and MCS. Based on our analysis, patients in the MCS group demonstrated significantly higher spectral power in the α and β bands across nearly all brain regions compared to the other groups. Elevated power in these frequency bands reflected the preserved cortical synchronization, functional network integration, and neurometabolic coupling capacity of MCS patients, making these features key biomarkers for differentiating states of consciousness [32]. Gamma oscillations are associated with higher-order cognitive processing and information binding across different brain regions [59]. The widespread enhancement of activities in these frequency bands in MCS patients may underpin the residual awareness and cognitive functions, such as the ability to intermittently follow commands. This finding aligns with previous research suggesting that high-frequency gamma activity is a neural correlate of conscious perception [30].

The brain functional connectivity patterns of three patient groups across different frequency bands was systematically analyzed based on wPLI. The five-band wPLI brain connectivity matrix (Fig. 5) for different brain regions revealed the dynamic reorganization characteristics of patients with varying levels of consciousness in large-scale brain network integration and dissociation. The connection strength in the delta frequency band increased as the level of consciousness decreased, which may be related to inhibition of the cortical-thalamus circuit and enhancement of slow-wave activity [60]. All three groups of patients showed strong posterior regional connectivity in the delta frequency band. Delta activity is a hallmark of deep sleep and severe brain injury. The predominant distribution in the occipital cortex may reflect deafferentation of the visual pathways and a generalized suppression of cortical function [61]. This pattern was consistently observed across the three groups, indicating that severe cortical dysfunction is a common pathologic basis for patients with DoC. A particularly intriguing finding emerged in the theta band. The VS group exhibited the strongest frontal-temporal-parietal connectivity in the theta frequency, surpassing that of the MCS group. This result challenges the simplistic linear model that “stronger connectivity equals higher level of consciousness” and supports the “over-activation of thalamocortical circuits” hypothesis. The highest theta connectivity in the VS group may reflect enhanced thalamic drive leading to pathologic synchronization between cortical regions, whereas the relatively lower theta connectivity in MCS patients suggests that the thalamocortical system is beginning to break free from this rigid synchrony, transitioning toward more flexible and complex dynamic activity [62]. This finding is consistent with several fMRI studies, which have reported that thalamic functional connectivity abnormalities in VS patients are even more severe than in coma patients [63]. The alpha-band connectivity pattern, which is centered on the frontal-right prefrontal regions, may underlie the network foundation that sustains a basic level of consciousness [43]. The dense connections in the beta band among temporal, parietal, and frontal areas observed in the VS and MCS groups could correspond to residual cognitive processing and self-awareness capacities [64]. These findings echo fMRI studies that report damage to the DMN and executive control networks in DoC [65]. Gamma-band activity is considered the neural basis for consciousness integration and the “binding problem” [66]. Only diffuse, disorganized gamma connectivity was observed among all patients with DoC, strongly suggesting that the breakdown of high-frequency connections may be a key neurophysiologic substrate for the loss of conscious content.

The analysis of EEG absolute power topography across the entire brain and different brain regions in the three patient groups under five EEG rhythms is essential. This study compared EEG characteristics across frequency bands in patients with coma, VS, and MCS using absolute power topographic maps. The ubiquitous elevation of delta power in the prefrontal regions is consistent with previous literature, potentially stemming from the desynchronization of thalamocortical projections and the hyper-synchronous firing of cortical neurons [30, 67]. VS and MCS patients demonstrated stronger power and more widespread brain involvement in the theta band compared to the coma group, which aligns with the role of theta oscillations in attention, memory, and network integration [68, 69]. Some studies have indicated that the preservation of theta activity is associated with the recovery of consciousness. Theta activity enhancement in MCS patients may reflect the partial functional preservation of thalamocortical circuits and the default mode network [70]. Chennu et al. [30] reported that alpha power is positively correlated with the retention of default network function, supporting use of alpha power as a predictor of recovery of consciousness. De Koninck et al. [71] noted that transcranial alternating current stimulation improves the state of consciousness by enhancing alpha synchronization, suggesting that alpha activity in patients with MCS may reflect residual neuroplasticity. Prefrontal beta activity correlated with dorsal attentional network function. The VS and MCS groups exhibited oscillatory activity in the beta band that was concentrated in the frontal regions with significantly greater power than the coma group. Neural oscillations in the beta frequency band had a key role in frontal-specific activation associated with motor intention. The increase in parietal beta power in patients with MCS may reflect a capacity for sensorimotor integration [72]. Gamma power was low in all three groups but prefrontal localization was active in the MCS and VS groups, which is consistent with the general observation of high-frequency oscillatory suppression in DoC [73, 74]. Notably, the high-power regions across all frequency bands involved the prefrontal lobes. Stender et al. [75] concluded that for every 10% increase in prefrontal metabolic rate, the probability of patients with MCS converting to conscious wakefulness increased 2.3-fold in the PET study. Specifically, the left DLPFC (corresponding to the F3 electrode) functions as a critical hub between the executive control network (ECN) and the DMN. The absolute power at the F3 electrode exhibits a consistent trend of coma  <  MCS  <  VS across all frequency bands in the three groups, which is in agreement with the findings herein that the higher power observed in the VS group may arise from network dysregulation leading to excessive synchronization [20].

EEG has made significant progress in the diagnosis and prognostic assessment of DoC, particularly in distinguishing between UWS/ VS and MCS. Toplutaş et al. [76] successfully distinguished between patients with VS and MCS by analyzing the duration, transition probabilities, and spatial topologic features of resting EEG microstates [76, 77, 78, 79]. The EEG spatial-spectral gradient model developed by Colombo et al. [58] in combination with support vector machines, achieved a classification accuracy of 92% in cross-center validation, outperforming the traditional CRS-R, which has an accuracy of approximately 75%. In addition, Pan et al. [80] successfully achieved online consciousness detection for patients with DoC using a P300 brain-computer interface network, attaining a response accuracy of 68%. Portable EEG systems have made significant advances in ICUs and rehabilitation wards, including real-time microstate tracking through mobile-based algorithms, multimodal stimulation paradigms integrating auditory, pain, and visual stimuli to enhance covert consciousness detection by 15%, and the integration of EEG features with clinical data to develop automated prognostic prediction models using elastic net regression [78]. The current study used a systematic analysis of multi-dimensional EEG features to uncover the differences in neural mechanisms among patients with DoC at varying levels of consciousness. The aEEG exhibited an increase corresponding to levels of consciousness and the power of high-frequency bands significantly increased with recovery of consciousness, further elucidating the association between high-frequency activity and consciousness recovery [30]. The elevated gamma-band power in MCS aligns with the existing literature. Gamma oscillations are associated with feature binding, attention, and conscious perception, and preservation of gamma oscillations in MCS may indicate residual cognitive capacity [33]. The increase in gamma power observed herein might reflect partially preserved thalamocortical and frontoparietal circuits, which are critical for consciousness. While previous studies have employed aEEG threshold criteria similar to those used in this research [10]—thereby providing a methodological precedent—it should be noted that standardized adult-specific threshold criteria remain to be established. Against this backdrop, the current study represents a methodologic breakthrough by integrating multi-dimensional EEG features, analyzing mechanisms across frequency bands, and validating findings with networks [81]. This approach transcends the limitations of single biomarkers and constructs a multi-scale model for consciousness recovery. In the future combining real-time neurofeedback with closed-loop regulation technology is expected to advance the assessment of DoC toward precision medicine.

However, although this study has made some advancements in methodology and existing evidence indicates that the level of consciousness serves as the primary determinant of neurophysiologic changes, it is still necessary to acknowledge the existence of certain limitations. First, its sample size was relatively limited, restricting the generalizability of findings. Second, it only focused on EEG features without integrating other neuroimaging modalities or comprehensive clinical variables. Third, the single-center design may introduce selection bias, and the follow-up period was not detailed. Fourth, diverse etiologies of DoC were not fully controlle. Future studies should increase the sample size, integrate multimodal neuroimaging and clinical data, extend the follow-up period, and conduct more explicit subgroup analyses in order to further validate and optimize the proposed model.

5. Conclusions

In summary, multidimensional EEG contributed to distinguishing different degrees of DoC. The aEEG center amplitude effectively distinguishes MCS from coma patients; relative power in the δ and γ bands shows discriminatory significance between VS and MCS patients; relative power in the α and β bands enables differentiation among coma, VS, and MCS patients. Furthermore, the brain network connectivity matrix in the θ band over the parieto-occipital region exhibits more distinct features when distinguishing coma from VS patients. The absolute power topography of the left DLPFC also demonstrates potential for discriminating levels of impaired consciousness, further supporting the utility of multidimensional EEG features in clinical differentiation.

Multidimensional EEG biomarkers extend prior single-domain EEG approaches by offering a concise, multiparametric framework that complements behavioral assessments and refines diagnostic classification in DoC. Future integration into bedside monitoring systems and automated analysis could enhance diagnostic accuracy, prognostication, and treatment monitoring.

Availability of Data and Materials

The data presented in this study are available on request from the corresponding author.

Author Contributions

SL and JG designed the study. YHS, SYL, JWC, MYZ collected the data. WJY and JZ analyzed the data and prepared the figures. YHS, SYL, JZ, WJY, JWC, MYZ drafted the manuscript. SL and JG brought major revisions in significant proportions of the manuscript. All authors read and approved the final manuscript. All authors have participated sufficiently in the work and agreed to be accountable for all aspects of the work.

Ethics Approval and Consent to Participate

The study was conducted under the Declaration of Helsinki of the World Medical Association and approved by the Ethics Committee of the Affiliated Hospital of Nantong University. The approval number for the study was 2024-K142-01. Before inclusion, the researcher fully informed the legal guardian of each patient about the study protocol and obtained informed consent.

Acknowledgment

We would like to express our gratitude to all those who acted as peer reviewers, providing us with their valuable opinions and suggestions.

Funding

This work was supported by the Jiangsu Province Research Hospital, China [grant numbers YJXYY202204], Natural Science Foundation of Jiangsu Province, China [grant number BK20241839] and Jiangsu Commission of Health, China [grant number M2022052].

Conflict of Interest

The authors declare no conflict of interest.

Supplementary Material

Supplementary material associated with this article can be found, in the online version, at https://doi.org/10.31083/JIN44233.

References

Publisher’s Note: IMR Press stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.