Disruption in structural–functional network repertoire and time-resolved subcortical fronto-temporoparietal connectivity in disorders of consciousness
Abstract
Understanding recovery of consciousness and elucidating its underlying mechanism is believed to be crucial in the field of basic neuroscience and medicine. Ideas such as the global neuronal workspace (GNW) and the mesocircuit theory hypothesize that failure of recovery in conscious states coincide with loss of connectivity between subcortical and frontoparietal areas, a loss of the repertoire of functional networks states and metastable brain activation. We adopted a time-resolved functional connectivity framework to explore these ideas and assessed the repertoire of functional network states as a potential marker of consciousness and its potential ability to tell apart patients in the unresponsive wakefulness syndrome (UWS) and minimally conscious state (MCS). In addition, the prediction of these functional network states by underlying hidden spatial patterns in the anatomical network, that is so-called eigenmodes, was supplemented as potential markers. By analysing time-resolved functional connectivity from functional MRI data, we demonstrated a reduction of metastability and functional network repertoire in UWS compared to MCS patients. This was expressed in terms of diminished dwell times and loss of nonstationarity in the default mode network and subcortical fronto-temporoparietal network in UWS compared to MCS patients. We further demonstrated that these findings co-occurred with a loss of dynamic interplay between structural eigenmodes and emerging time-resolved functional connectivity in UWS. These results are, amongst others, in support of the GNW theory and the mesocircuit hypothesis, underpinning the role of time-resolved thalamo-cortical connections and metastability in the recovery of consciousness.
Editor's evaluation
This is an important paper providing convincing evidence for altered brain dynamics in patients in a minimally conscious state and those with unresponsive wakefulness syndrome relative to healthy control participants. The results indicate reduced metastability and a contracted network repertoire in disorders of consciousness. Overall, the study provides important new information on mechanisms of disorders of consciousness and the functional brain networks involved.
https://doi.org/10.7554/eLife.77462.sa0Introduction
Diagnosis of the level of consciousness after coma due to severe brain injury is a well-known dilemma in the field of neurology and intensive care medicine. Coma after cardiac arrest or after traumatic brain injury may result in sustained altered states of consciousness. These patients with disorders of consciousness (DoC), irrespective of the aetiology, can be grouped into the unresponsive wakefulness syndrome (UWS) (Laureys et al., 2010), characterized by the presence of eye-opening and reflexive behaviours, and the minimally conscious state (MCS), characterized by consistent but fluctuant wilful conscious behaviours, such as command following or visual pursuit (Giacino et al., 2002; Schnakers et al., 2010). Recovery of consciousness is argued to emerge conjointly with restoration of resting-state functional brain networks (Edlow et al., 2021), which refers to patterns of neuronal interactions inferred by indirect (e.g. functional magnetic resonance imaging – fMRI) or direct (e.g. electro- and magneto-encephalography [EEG/MEG]) measurements. Analysis of these resting-state networks could potentially help in the diagnosis of patients with DoC and provide insight into the mechanisms that results in absence of recovery of consciousness in UWS.
Various resting-state networks that play an important role in the recovery of consciousness have been identified, among which the default mode network (DMN), fronto-parietal network (FPN), and the salience network are the most important (Amico et al., 2017; Heine et al., 2012). Recovery of the DMN in combination with recovery of the auditory network could for instance discriminate between MCS and UWS with a very high accuracy (~85%) (Demertzi et al., 2015). The mechanism of resting-state network restoration in DoC is yet unknown, however, thalamic activity and especially thalamo-cortical connectivity may be a driving force behind restorations of cortical network function that sustains conscious states (Fridman et al., 2014; Laureys et al., 2000). Previous work on resting-state networks in DoC have mainly focused on the ‘static’ picture of functional connectivity (Demertzi et al., 2013; Edlow et al., 2021; Giacino et al., 2014; Heine et al., 2012), that is connections are assessed over the entire duration of the (fMRI) recording and fluctuations in connectivity over time are ignored. However, the underlying dynamics of connectivity seem relevant for consciousness (Barttfeld et al., 2014; Luppi et al., 2019) and a static description may therefore be inadequate to provide mechanistic insight into failure of recovery of consciousness in DoC (Demertzi et al., 2019).
The analysis of dynamic or time-resolved functional connectivity, as well as the relationship between the underlying anatomical connections and emergent time-resolved functional connectivity (Avena-Koenigsberger et al., 2018; Suárez et al., 2020), may be clinically relevant in patients with DoC. Previous studies have already explored the role of time-resolved functional connectivity in DoC (Del Pozo et al., 2021; Golkowski et al., 2021; Sanz Perl et al., 2021). A recent study demonstrated that network states with long-distance connections occurred less frequently over time in MCS compared to UWS patients (Demertzi et al., 2019), emphasizing disintegration of interactions across the cortex in unconscious states. However, network states reminiscent of the well-known resting-state networks were not retrieved. Cao et al. used two methods to extract time-varying networks, that is independent component analysis and hidden Markov modelling, and revealed clinically relevant differences in network state durations between patients with DoC patients and healthy subjects (Cao et al., 2019), while lacking comparative analysis between patients in MCS and in UWS. In another fMRI study, the authors focused on the posterior cingulate area and the DMN using a spatiotemporal point process analysis and demonstrated decreased occurrence of DMN-like patterns in UWS. Dynamic connectivity analysis has also recently been applied to EEG data, revealing a loss of network integration and increased network segregation in DoC patients (Rizkallah et al., 2019). Spatiotemporal properties of networks have been explored using whole brain modelling, which shows reduction in stability, heterogeneity, and information flow in loss of consciousness (Escrichs et al., 2021; López-González et al., 2021; Panda et al., 2021). Despite the importance of the previously published work, the role of the well-known resting-state networks and especially thalamo-cortical functional connections (Monti et al., 2015) within the context of time-resolved connectivity and DoC has so far not been fully explored.
Another important aspect in the context of the emergence or restoration of resting-state networks is the underlying structural network, as anatomical connectivity patterns influence the repertoire of possible functional network states (Deco et al., 2013). It is widely assumed that switching between functional network states is achieved by so-called metastability in the brain (Deco and Kringelbach, 2016), that is winnerless competitive dynamics. A promising and robust approach to analyse the relationship between structural and functional network states is the so-called eigenmode approach (Atasoy et al., 2016; Robinson et al., 2016; Tewarie et al., 2020). With this approach, spatial harmonic components or eigenmodes are extracted from the anatomical network. These eigenmodes can be considered as patterns of ‘hidden connectivity’ that come to expression at the level of functional networks. It has been postulated that eigenmodes form elementary building blocks for spatiotemporal brain dynamics (Aqil et al., 2021). There is evidence that the well-known resting-state networks can be explained by the activation of a small set of eigenmodes (Atasoy et al., 2018). It can be hypothesized that switching between functional network states, as can be observed in the metastable brain, is accompanied by fluctuations in the expression of eigenmodes Preti and Van De Ville, 2019; therefore, a potential loss of metastability in DoC could co-occur with loss of modulations in eigenmode expression (Barttfeld et al., 2014).
In this context, the aim of the current study was fourfold. First, we tested whether loss of metastability and resting-state network activity, derived from time-resolved estimates of functional connectivity, could differentiate between MCS and UWS patients, with a potential extraction of a spatiotemporal thalamo-cortical network state. Second, we analysed whether time-resolved connectivity could be explained by modulations in expression of eigenmodes in DoC, and third, whether potential differences in eigenmode expression in DoC patients would co-occur with a loss of metastability. Finally, we conceptually link the findings of altered spatiotemporal dynamics underlying activity observed in the brain of DoC patients with several consciousness theories to increase our understanding of the mechanism behind pathological states of unconsciousness.
Results
We included 34 healthy control subjects (HC, 39 [mean] ± 14 years [standard deviation], 20 males), 30 MCS (41 ± 13 years, 21 males) and 14 UWS patients (48 ± 16 years, 7 males). There was no difference in patients with MCS and UWS in terms of age (p > 0.05), gender (p > 0.05), time since injury (p > 0.05), and aetiology (p > 0.05). There was also no difference in age (p > 0.05) and gender (p > 0.05) between HC and DoC patients. Further details about the patient population are described in the methods and Supplementary file 1.
Metastability and time-resolved functional connectivity in patients with DoC
Time-resolved or dynamic connectivity for all subjects was extracted from the phase information of the fMRI data. We quantified a proxy measure for metastability defined as the standard deviation of the overall phase behaviour over time (i.e. the Kuramoto order parameter). This was followed by the extraction of spatiotemporal patterns using non-negative tensor factorization (NNTF) from phase connectivity data, corresponding to resting-state networks or network states (see Figure 1). Well-known resting-state networks as well as a residual component were used as initial conditions for spatial connectivity patterns for all network states to allow for stable convergence of the algorithm (i.e. DMN, FPN, visual network, sensorimotor network, salience network, and subcortical network; Finn et al., 2015). However, the NNTF algorithm allowed the spatial patterns of these network states to change in order to maximize the explained variance of the data. Temporal statistics from the network states were derived for every network state in terms of excursions from the median proxy for nonstationarity (Zalesky et al., 2014) and state duration (i.e. dwell time).
A reduction of metastability was found in DoC patients compared to HCs (Figure 2A). Lower metastability was observed in UWS patients in comparison to MCS patients (Figure 2A). Reduced metastability is expected to occur with loss of switching between resting-state networks and potentially with dwelling within a more limited subset of resting-state networks in DoC. The output of the NNTF algorithm resulted in spatial topographies of some of the well-known resting-state networks, that is the DMN, a separate posterior DMN around the precuneus, the visual network, the salience network, the FPN, and a network consisting of subcortical fronto-temporoparietal regions (Sub-FTPN) (Figure 2I–N, Supplementary file 1B). Note that these network states were not identical to the initial conditions, for example the subcortical network that was provided as initial condition to NNTF was incorporated with the frontoparietal network (Figure 2N) by the NNTF algorithm. This modulated subcortical fronto-temporoparietal network consists of the following brain regions: bilateral thalamus, caudate, right putamen, bilateral anterior and middle cingulate, inferior and middle frontal areas, supplementary motor cortex, middle and inferior temporal gyrus, right superior temporal, bilateral inferior parietal, and supramarginal gyrus. At the same time, the sensorimotor network that was provided as initial condition disappeared as state. Excursions from the median were lower for most networks (DMN, visual, Salience, Posterior DMN, and Sub-FTPN) in DoC compared to HC (Figure 2C–F, H). Significant loss of nonstationarity was also found in UWS compared to MCS for the DMN, FPN, and Sub-FTPN (Figure 2C, G, H). The NNTF also yielded a residual state, with a lack of spatial structure, accounting for the variance of connectivity data not explained by the resting-state networks. The residual component had longer dwell times for the decreasing levels of consciousness (Figure 2B). In addition, there were lower dwell times in DoC patients for a specific set of resting-state networks (Salience, Posterior DMN, FPN, and Sub-FTPN), and dwell time was shorter in UWS patients compared to MCS patients only in the Sub-FTPN network (see results in Figure 2—figure supplement 1). These findings of very short dwell times in the posterior DMN, FPN, and Sub-FPTN and long dwell time in the residual network can be considered as a contraction of the functional network repertoire in DoC, which is in agreement with a loss in metastability in these patients.
Relationship between structural eigenmodes and time-resolved functional connectivity in DoC
As resting-state network activity can be explained by the activation of structural eigenmodes, we next analyse the role of fluctuations in eigenmode expression over time. In order to put our findings into context, we first analysed the relationship between static functional networks and structural networks, using the Pearson correlation between static functional connectivity and structural connectivity (SC) for the different groups (Figure 3A). These results show that functional connectivity in DoC patients shows more correspondence with the underlying SC as compared to HCs, as the relationship between structural and functional connectivity was stronger for decreasing levels of consciousness (Figure 3A).
We next obtained the eigenmodes from the SC by extracting the eigenvectors of the graph Laplacian. These eigenmodes can be regarded as distinct spatial harmonics within the SC, where the first eigenmodes correspond to patterns with low spatial frequency and subsequent eigenmodes contain patterns with increasingly higher spatial frequencies. Given their spatial configuration, consecutive eigenmodes can be associated with increasing levels of segregation while the first eigenmodes can be linked with network integration. For every time point we predicted the extent to which phase connectivity could be explained by a weighted combination of the eigenmodes (Tewarie et al., 2020). Since phase connectivity can evolve over time, the weighting coefficients for the eigenmodes can modulate as well, resulting in fluctuations in the strength of the expressions of eigenmodes over time. For every eigenmode, we could then quantify the modulation strength (i.e. how much the eigenmode expression varied over time). In addition to the weighting coefficients, we also obtained the goodness-of-fit for the predictions of time-resolved functional connectivity.
The goodness-of-fit for the eigenmode predictions is displayed in Figure 3B, where we show the average correlation between eigenmode predicted FC and empirical FC for the three groups. Results show better predictions for HC and MCS compared to predictions for static FC (median and interquartile range of correlations HC static 0.18 ± 0.04, HC eigenmode 0.39 ± 0.09, Z = −7.1, p < 0.001, MCS static 0.2 ± 0.05, MCS eigenmode 0.35 ± 0.18, Z = −4.8, p < 0.001). In order to test whether these eigenmode predictions of time-varying connectivity could have been obtained by chance, we redid our analysis using surrogate BOLD data (see methods ‘Analysis steps’). Results show that eigenmode predictions for time-resolved connectivity from surrogate data performed significantly worse compared to genuine empirical data (for all comparisons with surrogate data p < 0.001; Figure 3B). We did not test whether the contribution of individual eigenmodes differed between groups since this would come with a serious multiple comparisons problem.
Instead, since SC appeared to correlate stronger with static FC in DoC compared to HC, we expected that eigenmode coefficients in DoC patients would hardly change over time, highlighting the observation of a ‘fixed’ structural–functional network relationship in DoC patients. To analyse this lack of change in the structural–functional network relationship over time in DoC patients, we quantified the modulation strength of the weighting coefficients over time (see methods ‘Analysis steps’). We performed this analysis separately for the dominant (1st to 107th eigenmode, first half) and non-dominant eigenmodes (108th to 214th eigenmode, second half). Results for the dominant eigenmodes show a clear reduction in modulation of the eigenmode weighting in DoC patients compared to HCs (Figure 3C), with a significantly lower modulation of eigenmode expression in UWS compared to MCS patients. This result could not be explained by chance, since the same results could not be obtained from surrogate data (Figure 3C). Note that no between group difference was obtained for non-dominant eigenmodes (Figure 3D).
We have so far shown a reduction in the modulation strength of eigenmode expressions in DoC patients compared to HC subjects, as well as a loss of metastability in DoC patients and dwelling of the brain in fewer network states in DoC patients. This poses the question of whether these two observations are related. In Figure 3E–G, we show that metastability is strongly correlated with modulations in eigenmode expression for the three groups. This underscores the notion that loss of dynamic modulations in functional network patterns due to a loss of metastability could indeed be related to a reduced modulation of eigenmode expression.
Discussion
Differentiation between MCS and UWS is key for adequate diagnosis and prognosis in DoC patients as it is connected to medical-ethical end of life decisions. Use of imaging characteristics allows testing of hypotheses about causes of delayed, or failure of, recovery of consciousness. Here, we used state-of-the-art techniques to quantify time-varying functional connectivity, metastability, and the relationship between the underlying anatomical network and time-resolved functional connections. We demonstrated that these advanced techniques were sensitive to detect clinically relevant differences for the diagnosis of MCS and UWS patients. More specifically, we first demonstrated that UWS patients showed reduced metastability, and spend less time in states outside the natural equilibrium state that would favour cerebral processing in a cooperative and coordinated manner to support consciousness. This is accompanied by shorter state durations that the brain spends in the subcortical fronto-temporoparietal configuration in UWS. A loss of nonstationarity was observed in several resting-state networks (i.e. DMN, frontoparietal, and subcortical fronto-temporoparietal) in UWS compared to MCS patients. Finally, we showed that functional brain networks are more ‘fixed’ to the underlying anatomical connections and are less subject to spatial reconfigurations over time in UWS compared to MCS patients. The extent to which these spatial reconfigurations occurred (i.e. expressed as modulations in eigenmode expression) correlated strongly to metastability.
Our results are in agreement with several hypothesis and theories for the emergence of consciousness, of which most share the importance of thalamo-cortical connectivity for consciousness (Blumenfeld, 2021; Dehaene et al., 2011; Schiff, 2010). The mesocircuit hypothesis states that deafferentation between the frontal cortex and subcortical regions is crucial in explaining the failure of recovery of consciousness (Schiff, 2010). One of the most novel findings in the current work is the generation of the subcortical fronto-temporoparietal network. Although subcortical connections were, among others, used as initial conditions for the decomposition of the time-varying functional connectivity patterns into resting-state networks, incorporation with fronto-temporoparietal connections emerged from the data-driven NNTF algorithm. Another observation confirmed that this NNTF approach extracted DoC-relevant networks, namely that the sensorimotor network disappeared after optimization of spatial network patterns. This latter result is in line with the fact that somatosensory cortices are not directly involved in the emergence of consciousness, based on current theories (Naccache, 2018). In addition, we found that the subcortical fronto-temporoparietal network showed shorter dwell times in DoC patients compared to HC subjects, with even shorter state durations in UWS compared to MCS patients. Finally, this network also demonstrated a loss of nonstationarity in UWS compared to MCS patients. However, it should be noted that the subcortical fronto-temporoparietal network was not the only network with loss of time-resolved network characteristics; other resting-state networks also showed loss of nonstationarity, such as the DMN and FPN. However, a combination of shorter dwell times and loss of nonstationarity was only found for the subcortical fronto-temporoparietal network.
The mesocircuit hypothesis postulates that increased subcortical (specifically thalamic nuclei) inhibition and reduced excitatory output from the thalamus to the cortex (frontal, parietal, temporal, and occipital) leads to weak neural activity and hampers consciousness in DoC patients (Fridman et al., 2014; Giacino et al., 2014). Although our method does not allow to assess excitatory and inhibitory activity, it may be tempting to interpret that the subcortical fronto-temporoparietal network may play a crucial role in orchestrating global network interactions and dwell times. Hence, this subnetwork may be instrumental for the observed loss of nonstationarity in the other subnetworks. Although the importance of functional connections between the thalamus and frontal cortex has been emphasized by the mesocircuit hypothesis, and shown to relate to consciousness in hypothesis-driven functional (e.g., Crone et al., 2018; Fridman et al., 2014; Monti et al., 2015) and structural (e.g., Annen et al., 2016; Weng et al., 2017) neuroimaging studies, this is the first demonstration of the ability of a (semi)data-driven approach to identify this sub-network in the context of time-resolved functional connectivity. Most previous data-driven approaches have been unable to extract such a network (Cao et al., 2019; Demertzi et al., 2019; Rizkallah et al., 2019). Our findings also support the global neuronal workspace (GNW) theory (Dehaene et al., 2011), which emphasizes the importance of long-distance, feedforward and recurrent functional connections, large-scale reverberant networks, and metastable brain states in the emergence and recovery of consciousness (Mashour et al., 2020). Specifically, the GNW postulates that conscious perception should rely on two neural circuits. First, the feedforward circuit encompasses the visual and sensory areas and propagates sensory input to higher order integration areas. Second, in order for sensory input to be accessible to consciousness, recurrent activity in the circuit is required, responsible for ignition in a wide range of high-order functional areas including visual, parieto-temporal, thalamic, and frontal areas (Aru et al., 2020; Baars, 2002; Dehaene and Changeux, 2011; Mashour et al., 2020). The subcortical fronto-temporoparietal network and frontoparietal network obtained in a data-driven fashion conjugates with the latter recurrent circuit proposed to be essential for consciousness by the GNW recurrent circuit. The breakdown of spatio-temporal coordination in the subcortical fronto-temporoparietal network and frontoparietal network was associated to reduced states of consciousness. So far, the importance of metastability has mainly been addressed in the context of recovery of consciousness from anaesthesia (Cavanna et al., 2018). Here, we underscore this finding and demonstrate that a reduction of subcortical–cortical spatiotemporal functional organization into metastable states can even differentiate between UWS and MCS patients aligning with key factors for consciousness as proposed by GNW and mesocircuit hypothesis.
Despite the overlapping cortical circuits important for the mesocircuit model and GNW, differences between the two hypotheses exist. While the occipital cortex is considered part of the same mesocircuit, its role is mostly for feedforward propagation in the GNW. We found a separate visual network that conceptually overlaps more with the feedforward network described by the GNW. In line with GNW predictions, the temporal dynamics in this feedforward network did not distinguish different states of consciousness (Dehaene and Changeux, 2011; Mashour et al., 2020). On the other hand, the involvement of subcortical areas in the subcortical fronto-temporoparietal network is specifically predicted by the mesocircuit model (Fridman et al., 2014; Giacino et al., 2014). Last, the two theories highlight different aspects of brain activity required for consciousness. The GNW proposes specific effective connectivity patterns required for stimuli to reach consciousness. The mesocircuit model, on the other hand, proposes that balanced excitatory and inhibitory activity within the network is a prerequisite for conscious states to be supported. Our resting-state study conceptually thus resonates more with the concepts of the mesocircuit model, yet, we might hypothesize that relative preservation of the recurrent network crucial for the GNW meets the minimal requirement for conscious stimulus processing.
Our observation of altered spatio-temporal dynamics of functional connectivity that are also more restricted to the SC has been observed in pharmacological and pathological loss of consciousness (Demertzi et al., 2019; López-González et al., 2021; Mashour et al., 2021; Uhrig et al., 2018). During unconscious states it has been demonstrated that cerebral regional connectivity leads to an important reconfiguration of the repertoire of functional brain states, that forms simpler connectivity patterns in the brain, which is predominantly shaped by anatomical connectivity and is characterized by the reduction of complex long range connectivity patterns (Barttfeld et al., 2014; Demertzi et al., 2019; Mashour et al., 2021; Uhrig et al., 2018). Previous studies have already demonstrated a strong signature of the underlying structural network in functional networks in DoC. However, these previously described functional networks were mostly obtained using time-averaged estimation methods (Annen et al., 2016; Enciso-Olivera et al., 2021; Kuceyeski et al., 2016; Lant et al., 2016), and have not provided insight into a stronger structural functional coupling in DoC. Using time-resolved functional connectivity and its relationship to eigenmodes, we shed light on this observation. The underlying anatomical connectivity contains so-called inherent ‘hidden patterns’ or eigenmodes with different spatial structures. In a dynamical system such as the brain, these eigenmodes, or a combination of eigenmodes, can be sequentially activated or deactivated (Atasoy et al., 2016; Tewarie et al., 2019), and thereby shape the repertoire of possible functional network states. Since we assume that structural and functional networks share the same eigenmodes (Abdelnour et al., 2018; Naze et al., 2021; Tewarie et al., 2020), the weighting coefficients of the eigenmodes that result in the structural network can be considered as ‘steady state’. Reconfigurations of spatiotemporal networks for the healthy brain are induced by strong deviations from this steady state, demonstrated by strong modulations in eigenmode expression in our work. In DoC however, there is strong tendency towards the ‘steady state’ with minimal deviations of eigenmode expressions over time, and hence leading to a strong signature of structural connections in functional networks in DoC. This inability was even more pronounced for UWS than for MCS.
A few methodological aspects of our retrospective study deserve further discussion. First of all, we did not analyse the contributions of individual eigenmodes for two reasons: (1) although earlier studies have demonstrated that a limited set of eigenmodes could already explain observed functional connectivity pattern (Atasoy et al., 2016), in our case, assessing group differences between MCS and UWS on the basis of individual eigenmodes would impose a multiple comparisons problem; (2) our analytical approach is based on the assumption that all eigenmodes are necessary in the mapping to functional connectivity instead of a statistical selection of eigenmodes. Second, concatenated data from all groups were fed into the NNTF analysis, instead of per group. This assumes that the spatial network structure is similar across groups. Even though this is not necessarily the case, the amount of data to allow for stable NNTF results for individual groups was limited, especially for the UWS group. Future (multicentric) studies with more patients should verify whether the decomposition of the dynamic functional connectivity patterns into the observed subnetworks holds for the separate groups. However, our approach to concatenate data from all groups made group comparison much easier since there was now no need to ‘match’ potentially slightly different networks from the different groups. This also allowed us to focus on networks that were important in DoC. Lastly, we have used a high temporal resolution method to estimate time-resolved connectivity at every time point instead of a sliding window-based method. Previous studies using sliding window approaches have provided novel insights into brain dynamics of loss of consciousness, such as the occurrence and reconfiguration of brain states in DoC patients (Demertzi et al., 2019) and anaesthesia induced loss of consciousness (Barttfeld et al., 2014; Uhrig et al., 2018). However, sliding window approaches have limited sensitivity to nonstationarity in the fMRI BOLD signals (Hindriks et al., 2016). Recent work on time-resolved connectivity shows that brief periods in co-modulation of BOLD signals are an important driving factor for functional connectivity (Zamani Esfahlani et al., 2020; Hindriks et al., 2016), and sliding window approaches may be less sensitive to detect such brief periods in co-modulation. Furthermore, sliding window approaches in DoC have so far demonstrated limited success in retrieving canonical resting-state networks in DoC, while exploration of spatiotemporal aspects of well-known resting-state networks is important to better understand the relation between brain dynamics and loss of consciousness.
Taken together, we have demonstrated that a (semi) data-driven approach has extracted clinically meaningful time-resolved functional brain networks. This unique network-based spatiotemporal characterization accounts for structure–function coupling (i.e. eigenmodes), and shows a relationship with stability of brain dynamics. The measures that differed between UWS and MCS patients most, were the dominant eigenmodes (reflecting structure–function coupling) and time-resolved functional connectivity in the DMN, frontoparietal network, and the subcortical fronto-temporoparietal network. Interestingly, the latter network was generated by the (semi)data-driven approach to better fit the data, and was the sole network to show shorter dwell times in UWS compared to MCS patients. This suggests that the subcortical fronto-temporoparietal network might play a pivotal role in supporting conscious network interactions, as several theoretical and hypothesis-driven studies indicated. Future work will be required to assess to what extent these advanced aspects of connectivity can serve as biomarkers to aid diagnosis and prognosis in DoC.
Methods
Participants
Forty-four adult DoC patients, of whom 30 in MCS (11 females, age range 24–83 years; mean age ± SD, 45 ± 16 years) and 14 with the UWS (6 females, age range 20–74 years; mean age ± SD, 47 ± 16 years) and 34 age and gender matched HC subjects (14 females, age range 19–72 years; mean age ± SD, 40 ± 14 years) without premorbid neurological problems were included.
In this retrospective study, we used the same dataset as in López-González et al., 2021; Panda et al., 2021. The study was approved by the ethic committee of the the University Hospital of Liège (Belgium) (No. 2009-241), where all data were collected. Written informed consents to collect and published the data were obtained from all healthy subjects and the legal representative for DoC patients in accordance with the Declaration of Helsinki.
The diagnosis of the DoC patients was confirmed through two gold standard approaches that is (1) behavioural and (2) fluorodeoxyglucose-positron emission tomography ([18F]FDG-PET), excluding patients for whom these two diagnostic approaches disagreed. (1) Patients were behaviourally diagnosed through the best of at least five Coma Recovery Scale Revised (CRS-R) assessments, evaluating auditory, visual, motor, oromotor function, communication, and arousal (Giacino et al., 2004). (2) Behavioural diagnosis was complemented with the visual assessment of preserved brain metabolism in the frontoparietal network using [18F]FDG-PET as a neurological proxy for consciousness (Stender et al., 2014). Patient-specific clinical information is presented in Supplementary file 1A. We only included patients for whom (1) MRI data were recorded without anaesthesia, (2) diagnosis was based on at least five repetitions of the CRS-R assessment, (3) diagnosed as UWS or MCS, and (4) the [18F]FDG-PET diagnosis was in agreement with the clinical diagnosis. We excluded the patients (1) for whom the patients the structural MRI segmentation was incorrect or (2) if there were excessive head movement artefacts during MR recordings. There were 46 MCS patients in which 16 were discarded due to mismatch of PET and CRS-R diagnosis, 8 for failed segmentation and 4 for head movement artefacts. Amongst the 28 UWS patients, 8 were discarded due to mismatch of the PET and CRS-R diagnosis, 4 for failure of segmentation, and 2 for head movement artefacts.
MRI data acquisition
For the DoC dataset, structural (T1 and DWI) and functional MRI data were acquired using a Siemens 3T Trio scanner. 3D T1-weighted MP-RAGE images (120 transversal slices, repetition time = 2300 ms, voxel size = 1.0 × 1.0 × 1.2 mm3, flip angle = 9°, field of view = 256 × 256 mm2) were acquired prior to the 10 min of BOLD fMRI resting-state (i.e. task free) (EPI, gradient echo, volumes = 300, repetition time = 2000 ms, echo time = 30 ms, flip angle = 78°, voxel size = 3 × 3 × 3 mm3, field of view = 192 × 192 mm2, 32 transversal slices). HC subjects were instructed to keep eyes open and to be in relaxed state during the fMRI data acquisition. Last, diffusion-weighted MRI (DWI) was acquired in 64 directions (b-value = 1000 s/mm2, voxel size = 1.8 × 1.8 × 3.3 mm3, field of view = 230 × 230 mm2, repetition time = 5700 ms, echo time = 87 ms, 45 transverse slices, 128 × 128 voxel matrix) preceded by a single unweighted image (b0).
Resting-state fMRI preprocessing
Preprocessing was performed as in López-González et al., 2021 using MELODIC (Multivariate Exploratory Linear Optimized Decomposition into Independent Components) version 3.14, which is part of FMRIB’s Software Library (FSL, http://fsl.fmrib.ox.ac.uk/fsl). The preprocessing consisted of the following steps: the first five functional images were discarded to reduce scanner inhomogeneity, motion correction was performed using MCFLIRT, non-brain tissue was removed using Bet Extraction Tool (BET), temporal band-pass filtering with sigma 100 s, spatial smoothing was applied using a 5-mm FWHM Gaussian kernel, rigid-body registration was performed, and finally single-session ICA with automatic dimensionality estimation was employed (Griffanti et al., 2014). Then, FIX (FMRIB’s ICA-based X-noiseifier) was applied to remove the noise components and the lesion-driven for each subject. Specifically, FSLeyes in Melodic mode was used to manually identify the single-subject independent components (ICs) into ‘good’ for cerebral signal, ‘bad’ for noise or injury-driven artefacts, and ‘unknown’ for ambiguous components. Each component was evaluated based on the spatial map, the time series, and the temporal power spectrum (Griffanti et al., 2014). Next, for each subject, FIX was applied with default parameters to remove bad and unknown components. Subsequently, the Shen et al., functional atlas (without cerebellum) was applied for brain parcellation to obtain the BOLD time series of the 214 cortical and subcortical brain areas in each individual’s native EPI space (Finn et al., 2015; Shen et al., 2013). The cleaned functional data were co-registered to the T1-weighted structural image by using FLIRT (Jenkinson and Smith, 2001). Then, the T1-weighted image was co-registered to the standard MNI space by using FLIRT (12 DOF) and FNIRT (Andersson et al., 2007). The transformations matrices were inverted and applied to warp the resting-state atlas from MNI space to the single-subject functional data using a nearest-neighbour interpolation method to ensure the preservation of the labels. Finally, the time series for each of the 214 brain areas were extracted using fslmaths and fslmeants.
Probabilistic tractography analysis
A whole-brain SC matrix was computed for each subject as reported in our previous study (López-González et al., 2021). Briefly, the b0 image was co-registered to the T1 structural image using FLIRT and the T1 structural image was co-registered to the MNI space using FLIRT and FNIRT (Andersson et al., 2007; Jenkinson and Smith, 2001). The transformation matrices were inverted and applied to warp the resting-state atlas from MNI space to the native diffusion space using a nearest-neighbour interpolation method. Analysis of diffusion images was applied using the FMRIB’s Diffusion Toolbox (FDT) (https://fsl.fmrib.ox.ac.uk/fsl/fslwiki). BET was computed, eddy current distortions and head motion were corrected using eddy correct tool (Andersson and Sotiropoulos, 2016). Crossing fibres were modelled by using BEDPOSTX, and the probability of multi-fibre orientations was calculated to enhance the sensitivity of non-dominant fibre populations (Behrens et al., 2007; Behrens et al., 2003). Then, probabilistic tractography analysis was calculated in native diffusion space using PROBTRACKX to compute the connectivity probability of each brain region to each of the other 213 brain regions. Subsequently, the value of each brain area was divided by its corresponding number of generated tracts to obtain the structural probability matrix. Finally, the SCpn matrix was then symmetrized by computing their transpose SCnp and averaging both matrices.
Metastability and time-resolved functional connectivity
BOLD time series for every ROI k were filtered within the narrowband of 0.01–0.9 Hz. We estimated both time-resolved functional connectivity and static functional connectivity. Static functional connectivity () was estimated by the Pearson correlation coefficient between pairwise narrowband time courses. Time-resolved connectivity was estimated from the instantaneous phases. The instantaneous phase of the BOLD signal, , was extracted from the analytic signal as obtained from the Hilbert transform. The synchronization between pairs of brain regions was characterized as the difference between their instantaneous phases. At each time point, the phase difference between two regions j and k was used as estimate for instantaneous phase connectivity (Glerean et al., 2012). Evaluation of pairwise connectivity at each time point resulted in a functional connectivity tensor (number of ROIs × number of ROIs × time points).
Furthermore, we investigated how the synchronization between different nodes fluctuates across time using the concept of metastability. In the current sense, metastability quantifies how variable the states of phase configurations are as a function of time. We quantified metastability (i.e. our proxy measure for metastability) in terms of the standard deviation of the Kuramoto order parameter, , where is the number of ROIs and i denotes the imaginary unit (Deco et al., 2017).
Similar to Ponce-Alvarez et al., 2015; Tewarie et al., 2019, we extracted time-evolving networks using NNTF (Bro, 1997). We used the N-way toolbox (version 1.8) for MATLAB for this analysis (Andersson and Bro, 2000). NNTF can be considered as a higher-order principal component analysis. The goal of the approach is to decompose the functional connectivity tensor into components, such that the approximation of the functional connectivity tensor can be written as . Here, and correspond to vectors decoding spatial information for component l, and represents the vector that contains information on temporal fluctuations of component l. Note that the outer product stands for the spatial pattern of functional connectivity of component l. The number of components L can be estimated from the data using established algorithms (Timmerman and Kiers, 2000). However, here, we fixed the number of components based on the number of a priori expected networks in which we were interested. We fed the NNTF algorithm with initial conditions for and based on the spatial components of six expected resting-state networks: salience network, FPN, DMN, subcortical network, sensorimotor network, and visual network. We added a residual network to account for the unexplained variance in the functional connectivity tensor. Note that the spatial initial conditions did not indicate that these spatial components were kept fixed during NNTF calculation; rather, these spatial components were free to be adjusted according to the maximization of the explained variance. Data from all subjects and groups were concatenated to allow convergence to stable results.
Relationship between structural eigenmodes and time-resolved functional connectivity
For every subject, we extracted structural eigenmodes from the graph Laplacian of the SC matrix defined by , where refers to the diagonal degree matrix of SC. We further applied symmetric normalization to obtain a normalized Laplacian . Subsequently, eigenvectors and eigenvalues (together called eigenmodes) were extracted using diagonalization of , resulting in N eigenmodes (N = number of ROIs). Using a recently introduced approach we mapped functional brain networks at each time point from the structural eigenmodes (Tewarie et al., 2022; Tewarie et al., 2020). In other words, we estimated to what extent functional connectivity at each time point could be explained by a linear combination of the eigenmodes as follows:
where is the diagonal node strength matrix of the functional connectivity matrix at time t. The matrix corresponds to the weighting coefficient matrix for the eigenmodes and is obtained after optimization and equal to
Hence, modulations of eigenmode expressions over time are expressed in . Here, is the normalized graph Laplacian of at time t, and is the ith eigenvector of and the ith column of .
Analysis steps
Time-resolved functional connectivity
Individual temporal time courses for expression of component l for every subject were extracted from . A high value of at a certain time point indicates strong expression of this spatial pattern of functional connectivity () at that time point. At each time point, we determined the component with the strongest expression (), and assumed that connectivity at that point in time was dominated by this state or component. The duration that this component retained the strongest expression was considered as state duration or dwell time (see Figure 1). In addition, we also characterized the amount of nonstationarity in (Hindriks et al., 2016; Zalesky et al., 2014), that is excursions from the median. The rationale for using this metric is its sensitivity to detect modulations if the underlying system is indeed dynamic (Hindriks et al., 2016). Mann–Whitney U tests were used to test, for each network separately, differences in state durations and excursion from the median between groups.
Structural vs. functional brain networks
We first estimated the relationship between static functional connectivity and the SC itself (without decomposing SC into eigenmodes). This relationship was estimated in terms of the Pearson correlation coefficient between the SC and , denoted as . We secondly analysed the extent to which time-varying functional networks could be explained by expressions of the structural eigenmodes. We therefore computed the Pearson correlation between the empirical and the eigenmode predicted , denoted as . In order to be able to test whether there was a difference in fluctuations of eigenmode expression over time between groups, we quantified the eigenmode modulation strength defined as , where and correspond to different points in time and to the total duration of the recording. Since the first eigenmodes can be considered as dominant eigenmodes and more important to shape functional brain networks (Atasoy et al., 2016), we evaluated the eigenmode modulation strength for the dominant eigenmodes and non-dominant eigenmodes separately. We further computed the correlation between eigenmode modulation strength and metastability in all groups separately in order to test whether modulations in eigenmode expression related to our proxy measure for metastability. Group differences for all metrics was tested using Mann–Whitney U tests.
In order to test whether eigenmode predictions of functional connectivity could be obtained by chance, we created surrogate data and redid analysis for the surrogate data. Surrogate data for fMRI BOLD time series were obtained using the circular time shifted method (Quian Quiroga et al., 2002). Time-resolved phase connectivity was estimated in the same way as for genuine empirical data. Time-resolved connectivity obtained from surrogate data was subsequently predicted using the eigenmodes.
False discovery rate was used to correct for multiple comparisons for analysis steps 1 and 2: number of metrics (excursions * seven networks, dwell time * seven networks, metastability, , , non-dominant and dominant eigenmodes) * 3 comparisons + surrogate comparison with genuine data 3 * 3, (, non-dominant and dominant eigenmodes) comparisons = 66 tests (Benjamini and Hochberg, 1995).
Data availability
Connectivity matrices have been made available open access through EBRAINS of the Human Brain Project.
-
Human Brain ProjectPhase-interaction matrices of the phase relationships of different pathological states of consciousness.https://doi.org/10.25493/9GES-37K
References
-
The n-way toolbox for matlabChemometrics and Intelligent Laboratory Systems 52:1–4.https://doi.org/10.1016/S0169-7439(00)00071-X
-
Graph neural fields: A framework for spatiotemporal dynamical models on the human connectomePLOS Computational Biology 17:e1008310.https://doi.org/10.1371/journal.pcbi.1008310
-
Cellular mechanisms of conscious processingTrends in Cognitive Sciences 24:814–825.https://doi.org/10.1016/j.tics.2020.07.006
-
Human brain networks function in connectome-specific harmonic wavesNature Communications 7:10340.https://doi.org/10.1038/ncomms10340
-
Communication dynamics in complex brain networksNature Reviews. Neuroscience 19:17–33.https://doi.org/10.1038/nrn.2017.149
-
The conscious access hypothesis: origins and recent evidenceTrends in Cognitive Sciences 6:47–52.https://doi.org/10.1016/s1364-6613(00)01819-2
-
Characterization and propagation of uncertainty in diffusion-weighted MR imagingMagnetic Resonance in Medicine 50:1077–1088.https://doi.org/10.1002/mrm.10609
-
Controlling the false discovery rate: a practical and powerful approach to multiple testingJournal of the Royal Statistical Society 57:289–300.https://doi.org/10.1111/j.2517-6161.1995.tb02031.x
-
PARAFAC: Tutorial and applicationsChemometrics and Intelligent Laboratory Systems 38:149–171.https://doi.org/10.1016/S0169-7439(97)00032-4
-
Restoration of thalamo-cortical connectivity after brain injury: recovery of consciousness, complex behavior, or passage of time?Journal of Neuroscience Research 96:671–687.https://doi.org/10.1002/jnr.24115
-
Resting-state functional connectivity emerges from structurally and dynamically shaped slow linear fluctuationsThe Journal of Neuroscience 33:11239–11252.https://doi.org/10.1523/JNEUROSCI.1091-13.2013
-
The global neuronal workspace model of conscious access: from neuronal architectures to clinical applicationsCharact Conscious From Cogn to Clin 1:55–84.https://doi.org/10.1007/978-3-642-18015-6_4
-
Consciousness supporting networksCurrent Opinion in Neurobiology 23:239–244.https://doi.org/10.1016/j.conb.2012.12.003
-
Intrinsic functional connectivity differentiates minimally conscious from unresponsive patientsBrain: A Journal of Neurology 138:2619–2631.https://doi.org/10.1093/brain/awv169
-
Recovery from disorders of consciousness: mechanisms, prognosis and emerging therapiesNature Reviews. Neurology 17:135–156.https://doi.org/10.1038/s41582-020-00428-x
-
The jfk coma recovery scale-revised: measurement characteristics and diagnostic utilityArchives of Physical Medicine and Rehabilitation 85:2020–2029.https://doi.org/10.1016/j.apmr.2004.02.033
-
Disorders of consciousness after acquired brain injury: the state of the scienceNature Reviews. Neurology 10:99–114.https://doi.org/10.1038/nrneurol.2013.279
-
Dynamic patterns of global brain communication differentiate conscious from unconscious patients after severe brain injuryFrontiers in Systems Neuroscience 15:e919.https://doi.org/10.3389/fnsys.2021.625919
-
Resting state networks and consciousnessFrontiers in Psychology 3:295.https://doi.org/10.3389/fpsyg.2012.00295
-
A global optimisation method for robust affine registration of brain imagesMedical Image Analysis 5:143–156.https://doi.org/10.1016/s1361-8415(01)00036-6
-
Resting-state temporal synchronization networks emerge from connectivity topology and heterogeneityPLOS Computational Biology 11:e1004100.https://doi.org/10.1371/journal.pcbi.1004100
-
Performance of different synchronization measures in real data: A case study on electroencephalographic signalsPhysical Review. E, Statistical, Nonlinear, and Soft Matter Physics 65:041903.https://doi.org/10.1103/PhysRevE.65.041903
-
Recovery of consciousness after brain injury: a mesocircuit hypothesisTrends in Neurosciences 33:1–9.https://doi.org/10.1016/j.tins.2009.11.002
-
Linking structure and function in macroscale brain networksTrends in Cognitive Sciences 24:302–315.https://doi.org/10.1016/j.tics.2020.01.008
-
Three-mode principal components analysis: choosing the numbers of components and sensitivity to local optimaThe British Journal of Mathematical and Statistical Psychology 53 (Pt 1):1–16.https://doi.org/10.1348/000711000159132
-
Abnormal structural connectivity between the basal ganglia, thalamus, and frontal cortex in patients with disorders of consciousnessCortex; a Journal Devoted to the Study of the Nervous System and Behavior 90:71–87.https://doi.org/10.1016/j.cortex.2017.02.011
Article and author information
Author details
Funding
Horizon 2020 Framework Programme (945539)
- Gustavo Deco
- Steven Laureys
Horizon 2020 Framework Programme (785907)
- Steven Laureys
Horizon 2020 Framework Programme (778234)
- Steven Laureys
The Belgian National Funds for Scientific Research (F.R.S-FNRS)
- Rajanikant Panda
- Aurore Thibaut
- Olivia Gosseries
- Steven Laureys
European Space Agency
- Steven Laureys
Swiss National Science Foundation (Sinergia (170873))
- Ane Lopez-Gonzalez
Belgian Federal Science Policy Office (PRODEX Programme)
- Steven Laureys
Fondazione Europea di Ricerca Biomedica
- Steven Laureys
Bial Foundation
- Olivia Gosseries
- Steven Laureys
The fund Generet
- Steven Laureys
King Baudouin Foundation
- Steven Laureys
Leon Fredericq Foundation
- Steven Laureys
- Aurore Thibaut
- Olivia Gosseries
The DoCMA Project (EU-H2020-MSCA–RISE–778234)
- Aurore Thibaut
FLAG-ERA JTC (PCI2018-092891)
- Gustavo Deco
The Spanish Ministry Project PSI2016-75688-P (AEI/FEDER)
- Gustavo Deco
The Catalan Research Group Support 2017 SGR 1
- Gustavo Deco
Ministerio de Ciencia, Innovación y Universidades
- Gustavo Deco
Agencia Estatal de Investigación
- Gustavo Deco
European Regional Development Fund
- Gustavo Deco
Fonds De La Recherche Scientifique - FNRS (PDR project T.0134.21)
- Aurore Thibaut
- Olivia Gosseries
Mind Science Foundation
- Steven Laureys
The Télévie Foundation
- Olivia Gosseries
The Public Utility Foundation ‘Université Européenne du Travail’
- Steven Laureys
The Mind-Care foundation
- Steven Laureys
National Natural Science Foundation of China (Joint Research Project 81471100)
- Steven Laureys
The European Foundation of Biomedical Research FERB Onlus
- Steven Laureys
The ERA-Net FLAG-ERA JTC2021 project ModelDXConsciousness (Human Brain Project Partnering Project)
- Olivia Gosseries
University and University Hospital of Liege
- Aurore Thibaut
- Olivia Gosseries
- Steven Laureys
The funders had no role in study design, data collection, and interpretation, or the decision to submit the work for publication.
Ethics
Written informed consent was obtained from all healthy subjects and the legal representative for DoC patients.
Copyright
© 2022, Panda et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 2,397
- views
-
- 737
- downloads
-
- 25
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Computational and Systems Biology
- Microbiology and Infectious Disease
Timely and effective use of antimicrobial drugs can improve patient outcomes, as well as help safeguard against resistance development. Matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOF MS) is currently routinely used in clinical diagnostics for rapid species identification. Mining additional data from said spectra in the form of antimicrobial resistance (AMR) profiles is, therefore, highly promising. Such AMR profiles could serve as a drop-in solution for drastically improving treatment efficiency, effectiveness, and costs. This study endeavors to develop the first machine learning models capable of predicting AMR profiles for the whole repertoire of species and drugs encountered in clinical microbiology. The resulting models can be interpreted as drug recommender systems for infectious diseases. We find that our dual-branch method delivers considerably higher performance compared to previous approaches. In addition, experiments show that the models can be efficiently fine-tuned to data from other clinical laboratories. MALDI-TOF-based AMR recommender systems can, hence, greatly extend the value of MALDI-TOF MS for clinical diagnostics. All code supporting this study is distributed on PyPI and is packaged at https://github.com/gdewael/maldi-nn.
-
- Computational and Systems Biology
- Genetics and Genomics
Enhancers and promoters are classically considered to be bound by a small set of transcription factors (TFs) in a sequence-specific manner. This assumption has come under increasing skepticism as the datasets of ChIP-seq assays of TFs have expanded. In particular, high-occupancy target (HOT) loci attract hundreds of TFs with often no detectable correlation between ChIP-seq peaks and DNA-binding motif presence. Here, we used a set of 1003 TF ChIP-seq datasets (HepG2, K562, H1) to analyze the patterns of ChIP-seq peak co-occurrence in combination with functional genomics datasets. We identified 43,891 HOT loci forming at the promoter (53%) and enhancer (47%) regions. HOT promoters regulate housekeeping genes, whereas HOT enhancers are involved in tissue-specific process regulation. HOT loci form the foundation of human super-enhancers and evolve under strong negative selection, with some of these loci being located in ultraconserved regions. Sequence-based classification analysis of HOT loci suggested that their formation is driven by the sequence features, and the density of mapped ChIP-seq peaks across TF-bound loci correlates with sequence features and the expression level of flanking genes. Based on the affinities to bind to promoters and enhancers we detected five distinct clusters of TFs that form the core of the HOT loci. We report an abundance of HOT loci in the human genome and a commitment of 51% of all TF ChIP-seq binding events to HOT locus formation thus challenging the classical model of enhancer activity and propose a model of HOT locus formation based on the existence of large transcriptional condensates.