Front. Neurosci. Frontiers in Neuroscience Front. Neurosci. 1662-453X Frontiers Media S.A. 10.3389/fnins.2021.766633 Neuroscience Original Research Altered Effective Connectivity Measured by Resting-State Functional Magnetic Resonance Imaging in Posterior Parietal-Frontal-Striatum Circuit in Patients With Disorder of Consciousness Chen Linglong 1 † Rao Bo 1 † Li Sirui 1 Gao Lei 1 Xie Yu 2 Dai Xuan 3 Fu Kai 3 Peng Xu Zhi 2 Xu Haibo 1 * 1Department of Radiology, Zhongnan Hospital of Wuhan University, Wuhan, China 2Department of Neurology, Zhongnan Hospital of Wuhan University, Wuhan, China 3Department of Neurosurgery, Zhongnan Hospital of Wuhan University, Wuhan, China

Edited by: Notger G. Müller, Helmholtz Association of German Research Centers (HZ), Germany

Reviewed by: Lijun Bai, Xi’an Jiaotong University, China; Jun Zhang, Fudan University, China

 *Correspondence: Haibo Xu, xuhaibo@whu.edu.cn

†These authors have contributed equally to this work and share first authorship

This article was submitted to Perception Science, a section of the journal Frontiers in Neuroscience

 20 01 2022 2021 15 766633 29 08 2021 13 12 2021 Copyright © 2022 Chen, Rao, Li, Gao, Xie, Dai, Fu, Peng and Xu. 2022 Chen, Rao, Li, Gao, Xie, Dai, Fu, Peng and Xu

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

 Objective

Disorder of consciousness (DoC) resulting from severe brain injury is characterized by cortical and subcortical dysconnectivity. However, research on seed-based effective connectivity (EC) of DoC might be questioned as to the heterogeneity of prior assumptions.

 Methods

Functional MRI data of 16 DoC patients and 16 demographically matched healthy individuals were analyzed. Revised coma recovery scale (CRS-R) scores of patients were acquired. Seed-based d mapping permutation of subject images (SDM-PSI) of meta-analysis was performed to quantitatively synthesize results from neuroimaging studies that evaluated resting-state functional activity in DoC patients. Spectral dynamic causal modeling (spDCM) was used to assess how EC altered between brain regions in DoC patients compared to healthy individuals.

 Results

We found increased effective connectivity in left striatum and decreased effective connectivity in bilateral precuneus (preCUN)/posterior cingulate cortex (PCC), bilateral midcingulate cortex and left middle frontal gyrus in DoC compared with the healthy controls. The resulting pattern of interaction in DoC indicated disrupted connection and disturbance of posterior parietal-frontal-striatum, and reduced self-inhibition of preCUN/PCC. The strength of self-inhibition of preCUN/PCC was negatively correlated with the total score of CRS-R.

 Conclusion

This impaired EC in DoC may underlie disruption in the posterior parietal-frontal-striatum circuit, particularly damage to the cortico-striatal connection and possible loss of preCUN/PCC function as the main regulatory hub.

 disorders of consciousness seed-based d mapping permutation of subject images spectral dynamic causal modeling corticostriatal connection preCUN/PCC National Natural Science Foundation of China10.13039/501100001809 Introduction

Disorder of consciousness (DoC) includes several neurological conditions ranging from the minimally conscious state to unresponsive wakefulness syndrome or vegetative state, losing both responsiveness and consciousness, as a result of severe brain injuries. Resting-state functional magnetic resonance imaging (fMRI), describing spontaneous brain activity and generating intrinsic brain networks by calculating functional correlations between brain regions based on spontaneous blood oxygenation level dependent (BOLD) signal fluctuations, is one of the most employed techniques in the field of DoC. Several imaging studies have reported abnormality of crucial regions such as precuneus (preCUN)/posterior cingulate cortex (PCC), medial prefrontal cortex (mPFC), thalamus, and striatum related to consciousness dysfunction in DoC patients (Norton et al., 2012; Yao et al., 2015; Rosazza et al., 2016). Aberrant network connectivity has also been robustly and repeatably reported to be affected in DoC and several studies have reported that abnormal functional activity in the default mode network (DMN) is related to consciousness dysfunction and consciousness recovery in DoC patients (Vanhaudenhuyse et al., 2010; Rosazza et al., 2016; Amico et al., 2017). A common finding of studies on DoC is an impairment in the activity of high-order association cortices rather than “low-level” primary cortices.

However, as research continues, current theoretical views on consciousness have indicated widespread functional disconnection sufficient to lose awareness, rather than dysfunction in isolated regions (Dehaene et al., 2006; Seth et al., 2006). The deficits of cortico-subcortical (such as thalamus and striatum) and cortico-cortical connectivity have been recognized as the neuroanatomic basis of losing consciousness. The mesocircuit hypothesis was proposed to explain the vulnerability of the anterior forebrain (frontal/prefrontal cortical-striatopallidal thalamocortical loop systems; Schiff, 2010). Recent research of Adama et al. indicated the key functional role of the mediodorsal (MD) thalamus for the integrative function of the forebrain corticothalamic systems (Rikhye et al., 2012; Peräkylä et al., 2017; Schmitt et al., 2017; Parnaudeau et al., 2018). Monti et al. (2015) also reported impaired thalamo-frontal functional connectivity (FC) based on functional magnetic resonance imaging (fMRI) data acquired when DoC patients were performing a top-down cognitive task. Several other studies have found disrupted thalamocortical FC of DoC patients during resting-state (Tang et al., 2011; Zhou et al., 2011, 2014). Chen et al. (2018) estimated abnormal effective connectivity (EC) of the anterior forebrain regions which may be associated with patient prognosis. Major studies of cortico-cortical disconnection frame debates on the neural correlates of consciousness as “front versus back.” Arguments were mainly centered on contributions of frontal and posterior cortex in maintaining consciousness, with several theories proposed (Tononi et al., 2016; Farrell, 2018; Mashour et al., 2020). A growing body of studies has identified seed regions based on these theories (Tang et al., 2011; Zhou et al., 2011, 2014). Deficits in the frontal cortex, based on lateral fronto-parietal networks, are reported in DoC (Massimini et al., 2012; Crone et al., 2014). The posterior cortex, as DMN is based on midline fronto-parietal networks, is linearly related to the intensity of internal awareness. Silva et al. (2015) reported impaired functional connectivity between the midline preCUN/PCC and anterior cingulate/medial prefrontal cortices, considered as two core midline nodes of the DMN.

According to the research mentioned above, studies on global brain spontaneous activity using independent component analysis (Norton et al., 2012), amplitude of low-frequency fluctuation or FC strength (Wu et al., 2019) revealed altered brain intrinsic function within regions rather than interaction between regions. Abnormal effective connectivity of isolated brain regions discussed at circuit-level should be the goal of future studies. However, the above-mentioned fMRI research revealed abnormal functional connectivity in cortico-subcortical and cortico-cortical regions using a seed-based method based on a priori assumptions which may introduce bias, indicating a need for an unbiased strategy of selecting brain regions in further research.

To explore the crucial cortical and subcortical regions and verify contributions of these regions in DoC, we performed a meta-analysis on resting-state fMRI data extracted from an online database. Given that data are gathered from isolated brain regions using different methods, quantitative meta-analysis allows heterogeneous results of individual studies to be pooled and analyzed using a rigorous statistical framework that can identify regions of vulnerability associated with different levels of conscious awareness. Seed-based d mapping permutation of subject images (SDM-PSI) initially developed by Albajes-Eizagirre et al. (2019) is helpful to summarize the data of interest and provide brain regions of interest for further analysis. This voxel-based meta-analytic method allows a more exhaustive and unbiased inclusion of multiple voxel-based neuroimaging studies of normal brain functions and brain abnormalities in neuropsychiatric disorders, as well as more accurate estimations. Selecting the spatially extensive set of seed regions of interest (ROIs) provided in the meta-analysis, we aimed to organize the results into a coherent model of large-scale reciprocal connections between circuit-level regions.

The present study was designed to investigate the localization of the brain regions displaying altered activity in patients with DoC with a quantitative meta-analysis, summarizing the resting-state fMRI literature available to date. We compared EC among these regions and their directionality between DoC patients and healthy controls using spectral dynamic causal modeling (spDCM). Finally, we investigated the relationship between DoC and altered ECs between the regions.

 Materials and Methods Participants

In the current study, eye-closed and awake resting-state fMRI data were obtained in a total of 19 DoC patients and 17 age- and gender-matched healthy controls (HC). Patients’ diagnosis based on Coma Recovery-Scale-Revised (CRS-R) assessment (Giacino et al., 2004) was made by two specialists. Data from subjects with head motion more than 3 mm translation, 3° of rotation or large cerebral deformation were excluded from further analyses. Careful visual image quality inspection was used to exclude unsuccessful realignment and segmentation, and data from a total of 3 DoC patients and one control subject were excluded on this basis. Data from sixteen DoC patients (13M/3F, mean age 46 years, SD 15 years, 6MCS/10VS) and 16 demographically matched HCs (12M/4F, mean age 45 years, SD 16 years) were therefore included in further analysis. Demographic and clinical data of DoC patients and HCs are shown in Tables 1, 2. Informed consent was obtained directly from all healthy participants and the legal representative of all patients. The study was approved by the ethics committee of Wuhan University Zhongnan Hospital and was conducted in accordance with the tenets of the Declaration of Helsinki.

Demographic characteristics of the patients with disorders of consciousness (DOC) and the healthy controls (HC).

Patient index Gender Age (years) Time of fMRI (days after insult) Etiology CRS-R scores MCS/VS
P01 M 75 15 Hypoxic brain injury 4 VS
P02 M 63 23 Hypoxic brain injury 7 VS
P03 F 26 48 Intracranial surgery 4 VS
P04 M 53 81 Hypoxic brain injury 7 VS
P05 M 75 592 Traumatic brain injury 9 MCS
P06 M 55 991 Intracerebral hemorrhage 9 MCS
P07 M 65 57 Intracerebral hemorrhage 4 VS
P08 M 63 172 Intracerebral hemorrhage 6 VS
P09 M 38 125 Traumatic brain injury 12 MCS
P10 F 48 96 Traumatic brain injury 6 VS
P11 F 26 165 Intracranial surgery 11 MCS
P12 M 32 310 Intracranial surgery 4 VS
P13 M 47 45 Hypoxic brain injury 7 VS
P14 M 73 64 Hypoxic brain injury 6 VS
P15 M 19 149 Hypoxic brain injury 9 MCS
P16 M 19 179 Hypoxic brain injury 11 MCS

Demographic and clinical data of DoC patients.

Characteristics DoC HC P-value
Gender (male/female) 13/3 12/4 P > 0.05
Age (years) 46 ± 15 45 ± 16 P > 0.05
Handedness (R/L) 15/1 16/0 P > 0.05
CRS-R scores 7.23 ± 2.67 – –

M, male; F, female.

 Behavioral Assessment

The CRS-R score was acquired using largely preserved brainstem reflexes and preserved sleep–wake-cycle checked by neurological examination.

 MRI Data Acquisition

MRI data of all subjects were acquired on a 3.0 T Siemens Prisma scanner (Siemens AG, Healthcare Sector) using a 64-channel head coil. Each subject was required to keep in the supine position by a belt and foam pads during rest and awake condition with eyes closed. fMRI data was scanned using a gradient-echo echo-planer imaging (EPI) sequence of 240 volumes in an ascending interleaved order using the following protocols: repetition time (TR)/echo time (TE) = 2,000/40 ms, flip angle = 90°, field of view (FOV) = 240 mm × 240 mm, slice thickness = 4.0 mm, inplane resolution = 64 × 64, 32 axial slices with a slice gap of 1 mm. Then, high-resolution brain structural images were collected with a T1-weighted 3D magnetization-prepared rapid gradient-echo (MPRAGE) sequence (TR/TE = 1,900 ms/2.26 ms, matrix = 240 × 256, FOV = 215 mm × 230 mm, slice thickness/gap = 1.0/0 mm, 176 sagittal slices).

 Selection of Regions of Interest (ROI) Search Strategy and Selection Criteria

The literature search was conducted on whole-brain fMRI studies in DoC using the PubMed1 and Web of Science2 databases without publication date filter and combinations of the following keywords: “disorder of consciousness/DoC,” “unresponsive wakefulness syndrome,” “minimally conscious state,” “vegetative state” plus “resting-state fMRI” or “neuroimaging.” Due to the limited number of publications, we included any cause of DoC, including traumatic brain injury, anoxic brain injury, cardiac/cerebral vascular events, intoxication, and hypoxic ischemic encephalopathy, etc. After carefully screening the titles and abstracts of related studies, 57 potential studies were selected from 205 original studies for further scrutiny. Based on the rules for neuroimaging meta-analysis (Muller et al., 2018; Tahmasian et al., 2019), higher-quality criteria for study inclusion were as following: images covering the whole brain, inclusion of more than 10 subjects, whole-brain analyses, match for age and gender, use of standardized categorical or dimensional measures, definitive inclusion/exclusion criteria and description of software use and statistical methods. Studies were excluded if they were reviews, unrelated issues, without reporting the Montreal Neurological Institute (MNI)/Talairach coordinates (Lancaster et al., 2000) and T/Z value of clusters, ROI-based analysis, small-volume correction, studies of independent component analysis focused on only one resting-state network, or repeating data in two studies. Need of special note is that only studies in which DoC patients were combined rather than reporting separate information for each condition (e.g., VS and MCS) were included in the analysis. A total of 8 studies were finally selected for the following research (see Figure 1 and Table 3).

Flowchart of literature screening.

Subject demographic in included studies.

Study N Analysis method Gender (M/F) Age in years (SD) Space Contrast
Huang et al., 2019 20 fALFF 13/7 53.30 (NA) MNI Doc<HC
Wu et al., 2019 15 FCS & fALFF 12/3 41.40 (13.22) MNI Doc<HC
He et al., 2014 12 fALFF 8/4 44.7 (17.5) MNI Doc<HC
Soddu et al., 2016 11 total neuronal activity 5/6 50 (14) MNI Doc<HC
Demertzi et al., 2014 53 ICA 33/20 50 (18) MNI Doc<HC
Huang et al., 2014 11 ALFF 10/1 44 (20) MNI Doc<HC
Yao et al., 2015 11 ALFF 6/5 43.1 (15.6) MNI Doc<HC
Salvato et al., 2020 14 fALFF 5/9 63.6 (8.5) MNI Doc<HC

N, number; M, male; F, female; SD, Standard deviation; ReHo, Regional Homogeneity; fALFF, Fractional Amplitude of Low Frequency Fluctuations; FCS, functional connectivity strength; ALFF, Amplitude of Low Frequency Fluctuations; ICA, Independent Component Analysis; MNI, Montreal Neurological Institute; DoC, disorder of consciousness; HC, healthy control.

 SDM Meta-Analysis

Quantitative evaluation of meta-analysis was performed using the SDM software3 with MNI coordinates (Talairach coordinates were first converted into MNI coordinates) on the brain activation patterns for DoC and HC. The meta-analysis was performed in the following steps which have been described previously (Radua and Mataix-Cols, 2009; Radua et al., 2012). Firstly, an effect-size signed map and an effect-size variance map of given peak coordinates and their between-group t-values were extracted from individual studies using an anisotropic unnormalized Gaussian kernel with 20 mm full width at half maximum (FWHM) to control false-positive results. Then, the SDM generated a mean map by voxel-wise calculation of the mean of the dataset maps with a random-effect linear model, weighted by the sample size of each study, estimating intrastudy variability and inter-study heterogeneity. The result of heterogeneity analysis using a random-effects model with Q statistics represents between-study variability in the results. The above analyses were performed with a statistical significance of the following combination of thresholds: voxel probability threshold p < 0.005, peak height Z ≥ 1 and cluster extent threshold >10 voxels, achieving an optimal balance between sensitivity and specificity (Radua et al., 2012).

 Functional Characterization of Identified SDM Meta-Analysis Result

To facilitate the functional interpretability of the identified clusters resulting from SDM meta-analysis, we used the Neurosynth database4 for data-driven characterization. Neurosynth is a database for large-scale, automated synthesis of functional magnetic resonance imaging (fMRI) data from >14,000 functional MRI studies and can be searched for the functional decoding of chosen clusters (Yarkoni et al., 2011). We selected the psychological terms from the association list with meta-analysis maps of the identified cluster in the Neurosynth database.

 Acquisition and Preprocessing of Functional Magnetic Resonance Imaging Data

After converting the DICOM format of the functional data to NIFTI format, all preprocessing was completed using the Data Processing and Analysis of Brain Imaging (DPABI v4.0)5 toolbox (Yan et al., 2016) in a MATLAB (2013b, MathWorks Inc., Natick, MA, United States) platform. The standard preprocessing procedures based on the following steps: discarding the first ten volumes of each data set to avoid magnetization instability, slice-timing correction, head motion correction, spatial normalization to the Montreal Neurological Institute (MNI) space based on high-resolution T1WI registration, resampling to a voxel size of 3 × 3 × 3 mm3 and smoothed with a Gaussian kernel of 6 × 6 × 6 mm3 depending on the full width at half maximum (FWHM). Besides, the criterion of head motion was limited within 3 mm translations and 3° rotations.

 Dynamic Casual Modeling

The spDCM analyses were performed using the DCM12 module in SPM12 (Wellcome Department of Cognitive Neurology, London, United Kingdom6). Based on results of meta-analysis, the averaged time series of each voxel in the left striatum, middle frontal gyrus (MFG), midcingulate cortex (MCC) and preCUN/PCC were extracted from each subject. The first eigenvectors were then extracted after the influences of head motion and low-frequency drift were removed using a generalized linear model (GLM) implemented in SPM12. A fully connected mode was specified meaning that bi-directional connections were identified between every pair of ROIs in each subject. No exogenous input was provided to this model of resting state fMRI data. For the given four ROIs, 24 free parameters performing the effective connections among the ROIs. Spectral DCM obtains a prediction of data feature based on the Fourier transform of the cross-correlation of the time series, previously described in detail (Yan et al., 2016).

To determine the best fitting model for each group through all possible dynamic causal models, a post hoc model selection routine fitting the full model of all free parameters was applied (Friston and Penny, 2011). As there are 24 free parameters in the current estimation, a “greedy search” was efficiently implemented to result in 28 = 256 reduced models (Rosa et al., 2012). A post hoc model optimization routine was then used to determine the best fitting model for each group based on posterior probability. The fully connected model has the highest evidence among distribution of log probability over the 256 reduced models. After the optimized model was selected, one sample T-test and two sample T-test were performed to compare the parameter inference of endogenous connections (DCM.Ep.A) within and between the groups.

Figure 2 shows the experiment designing of the current study.

Flow diagram shows the experiment designing of the current study.

 Results Result of Meta-Analyses

One hundred and forty-seven subjects with DoC (92 males) were included from eight studies, details of which are shown in Table 3. Figure 3 and Table 4 show the SDM of altered activity in patients with DoC at resting-state compared to healthy subjects. Result on the entire sample of DoC patients was constituted of 4 separate regions highlighting the involvement of deactivated cortical areas (bilateral preCUN/PCC, left MCC, bilateral MFG) and activation in left striatum.

Results of SDM-PSI meta-analysis. Yellow represents the brain region with increased activity. Blue represents the brain region with decreased activity. preCUN, precuneus; PCC, posterior cingulate cortex; MCC, middle cingulate cortex; MFG, middle frontal gyrus.

Results of meta-analysis were selected for following spDCM analyses.

Regions Hemisphere MNI coordinates SDM-Z value P value voxel
Striatum L −20, 0, −6 6.142 <0.05 16
PreCUN/PCC B −2, −66, 26 −8.057 <0.0001 594
MCC B 2, 8, 44 −6.529 <0.001 131
MFG L −26, 18, 48 −6.658 <0.005 43

MNI, Montreal Neurological Institute; L, left; R, right; preCUN, precuneus; PCC, posterior cingulate cortex; MCC, midcingulate cortex; MFG, middle frontal gyrus.

 Functional Characterization of the Identified Clusters

Figure 4 depicts the functional profiles of preCUN/PCC, MFG, MCC and left striatum defined by NeuroSynth. Several highly correlated terms with similar functional meanings were merged, for example, memory retrieval was the combination of memory retrieval and retrieval. The exact r values of extracted psychological terms mean value of Pearson correlation between functional terms and the selected brain region.

The left column showed brain regions showing positive (blue) or negative (yellow) correlation with DoC in the meta-analysis. Clusters were displayed at voxel probability threshold P < 0.005, peak height Z ≥ 1 and cluster extent threshold >10 voxels. The right column showed functional characterization of identified clusters. The statistics in the radar charts are r values from Pearson’s correlations. In the current Neurosynth framework, the r values reflect the correlation between two maps. preCUN, precuneus; PCC, posterior cingulate cortex; MCC, middle cingulate cortex; MFG, middle frontal gyrus.

 Bayesian Model Selection and Dynamic Causal Modeling

The significant differences in effective connectivity strength between the DoC patients and controls are shown in Figures 5, 6. Strength of positive connectivity between both MFG and MCC and the left striatum was significantly lower in DoC patients than controls (P < 0.05). The strength of the connection between the left striatum and preCUN/PCC in HCs was not reported in DoC (while reported as increased mathematically) (P < 0.05). We observed that the strength of negative connectivity between MFG and MCC, disappeared negative connectivity strength between MFG and preCUN/PCC and self-connection of preCUN/PCC was observed as positive in DoC patients (P < 0.05).

The result of DCM of HC (A) and DoC (B). L, left; R, right; Stria, striatum; preCUN, precuneus; PCC, posterior cingulate cortex; MCC, middle cingulate cortex; MFG, middle frontal gyrus.

Significant changes in effective connections (Hz) in the DOC patients. Red arrows, connections with significant increases; blue arrows, connections with significant decreases (P < 0.05) (A). Mean effective connections of significantly changed pairs of regions (B). The strength of self-connection of preCUN/PCC is negatively correlated with the score of CRS-R (C). L, left; R, right; Stria, striatum; preCUN, precuneus; PCC, posterior cingulate cortex; MCC, middle cingulate cortex; MFG, middle frontal gyrus.

To further explore the relationship between clinical scale and properties of effective connectivity in patients, Pearson’s correlation was calculated between the CRS-R scores and those values, showing statistically significant differences between the two groups. The strength of self-connection of preCUN/PCC was negatively correlated with the CRS-R score in DoC patients (P < 0.05) (Figure 6C).

 Discussion

Our meta-analysis results show diminished activity in preCUN/PCC, MCC and MFG (higher order areas of cognition) in DoC, in agreement with previous research (Laureys, 2005). In DCM analysis of resting-state fMRI, values of connectivity strength made sense of rate constants of neuronal responses, interpreted as reciprocal connectivity. We identified the specific ECs and directionality among these regions at resting-state in DoC patients compared with controls. In particular, our results of the functional integrity of cortical-striatum and posterior parietal-frontal connectivity indicated the abnormal functional connectivity between network of intrinsic awareness and extrinsic awareness. Importantly, the self-inhibition of preCUN/PCC appeared to be related to the degree of consciousness impairment.

Differences in the Cortical-Cortical Connectivity Between the Two Groups

PreCUN/PCC, as a component of the DMN, is considered crucial to maintain consciousness. At resting state, preCUN/PCC plays the role of a central and structural hub of information input from multiple brain regions (Hagmann et al., 2008; Deshpande et al., 2011). Previous studies of effective connectivity using DCM of resting-state fMRI in healthy individuals have highlighted the role of the posterior cingulate cortex (PCC) as a main connector hub between distinct networks in healthy individuals (Hagmann et al., 2008; van den Heuvel and Sporns, 2011; Li et al., 2012). We demonstrated that the disruption of self-inhibition and neuronal oscillations in preCUN/PCC is negatively correlated with the total score of CRS-R. This finding is generally consistent with previous research (Crone et al., 2015). The long afferent and efferent axonal fibers of preCUN/PCC connecting to distant cortical and subcortical targets, as well as the high resting state metabolic activity in this region make it vulnerable to damage (Hannawi et al., 2015). Afrasiabi et al. (2021) showed the contribution of parietal cortex in detecting changes in consciousness and in the integration of cortico-striatal-thalamic activity. Our results suggest that preCUN/PCC functions as the regulatory hub among cortical regions maintaining equilibrium between excitatory and inhibitory connectivity may be lost in DoC patients.

The regions MFG and MCC play a role in advanced cognitive functions such as aspects of memory, mental state and self-reference (Figure 1; Jin et al., 2018). Our results indicated that DoC patients have disrupted connectivity between the posterior parietal (preCUN/PCC) and frontal (MFG) cortex, related to two distinct and negatively correlated networks of consciousness. It has been proposed that loss of consciousness may be associated with a deterioration of the functional dynamics between mid-line regions of the frontal and parietal cortices. Evidence has accumulated from studies showing widespread disconnection between frontal and parietal regions in patients with impaired consciousness (Boly et al., 2009; Threlkeld et al., 2018). Laureys and Schiff proposed a model of recovery of consciousness which emphasized the contribution of network connectivity between frontal and parietal regions (Laureys and Schiff, 2012). Bonini et al. (2016) found changes in synchrony in frontal and parietal regions associated with loss of consciousness in frontal lobe epilepsy. Similar breakdowns in intracortical connectivity of frontal regions (ventral attention network) and parietal regions (DMN) have also been identified in sevoflurane-induced unconsciousness (Palanca et al., 2017). Our current result also indicated the altered effective function connectivity from MFG to preCUN/PCC may related with maintenance of consciousness.

 Impaired Effective Connections of Corticostriatal Circuitry in Disorder of Consciousness Patients

The central role of the striatum in consciousness and the necessary link between the disrupted corticostriatal connections and loss of consciousness are highlighted by the present results. Previous studies consider the striatum to be an important subcortical structure with several motor and cognitive functions (Liu et al., 2020) rather than a brain area contributing to consciousness (Liu et al., 2020). However, a growing body of literature has suggested altered striatum activity in hallucinogen-induced altered states of consciousness (Slagter et al., 2017), general anesthesia (Mhuircheartaigh et al., 2010) and the sleep-wake cycle (Braun et al., 1997). The striatum also contributes to consciousness through reciprocal connections with the thalamus and cortical regions according to the mesocircuit theory (Schiff, 2010). Lacey et al. (2007) indicated a sharply reduced output of medium spiny neurons in the striatum of diffuse brain injury subjects which lead to a reduction of direct excitatory input from the central thalamus and down-regulation of corticostriatal input (Schiff, 2010). These existing models have explained the mechanism of impaired metabolic activity of the striatum during changing states of consciousness.

In support of the complex and extensive corticostriatal circuitry theory (Haber, 2016), we found reduced MFG-striatum, MCC-striatum and preCUN/PCC-striatum connectivity in DoC patients compared with controls. As the main source of input to the basal ganglia, the striatum is connected with the cerebral cortex through the ganglia-cortical loop (Albin et al., 1989). In healthy volunteers, DCM has demonstrated a negative effective connection between cortical nodes with preCUN/PCC as the main hub (Soch et al., 2017), with the striatum receiving excitatory input from all other cortical nodes. Both input and output of cortical regions show a trend toward reduced or absent activity in DoC. Several current models indicate that the direct and indirect efferent pathways from the striatum, influencing basal ganglia output nuclei, play important roles in the regulation of thalamocortical and brainstem motor circuits (Smith et al., 1998; Kravitz et al., 2010; Freeze et al., 2013), thought to be fundamental to basic motor functions of survival. The frontal-striatal circuit is also well described in healthy populations, implicated in behavioral alterations in diseases of motor dysfunction (Quan et al., 2013; Chen et al., 2018). Alexander et al found reduced or absent effective connectivity between cortex and striatum in anesthetized rats compared with free behavior rats (Nakhnikian et al., 2014). Striatal activity is also observed to be strongly reduced when epilepsy patients lose consciousness (David et al., 2008; Moeller et al., 2008). Taking these studies together, the disappeared MFG-striatum may explain movement disorders in DoC patients. However, Nicholas et al reported disrupted structural integrity of preCUN and bilateral striatum in DoC patients, perhaps suggesting a loss of striatal output and dysregulation between the preCUN and the striatum (David et al., 2008; Moeller et al., 2008). Another DCM study by Crone et al. (2017) indicated a breakdown of basal ganglia-cortical (including PCC) connectivity in a loss of consciousness. Considering the sensory inputs that the striatum receives from multiple cortical sources, our result likely reflects the necessity of integrated corticostriatal circuitry to maintain consciousness.

 Limitations

The present study has some limitations. Firstly, we included any cause of DoC when recruiting patients and in the meta-analysis. We excluded some patients before data acquisition and discarded several others based on relatively strict procedures of quality control, reducing the sample size. A larger cohort is needed to study the different types of DoC caused by different etiologies and this will be a focus in our future work in this area. Secondly, our study has limitations inherent to meta-analyses. The meta-analysis study design involves data extracted from published studies rather than raw data and statistics, increasing the variance in the pooled data and possibly reducing the accuracy of results. Thirdly, although a number of DOC patients were recruited for data acquisition, we were only able to utilize 16 patients for the analysis. There was an agonizing trade-off between subject quantity and image quality. The extensive brain lesion in some patients would have produced poor normalization results and overlarge head-motions would have introduced artifacts. Because of these reasons, the sample size of this study is small and it is difficult to divide DoC patients into MCS and VS again, which is also what we need to make efforts in the future.

 Conclusion

Using meta-analysis, we identified decreased activity of preCUN/PCC, MFG and MCC and increased activity of the left striatum in DoC patients. The resulting pattern of interaction in DoC identified disrupted connection, disturbance of the posterior parietal-frontal-striatum circuit, and reduced self-inhibition of preCUN/PCC. These impairments may be principally due to the disruption of mechanisms underlying damage of corticostriatal connection and possible loss of function of preCUN/PCC, a significant regulatory hub.

 Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

 Ethics Statement

The studies involving human participants were reviewed and approved by the Ethics Committee of Zhongnan Hospital of Wuhan University. The participants and family member of patients provided their written informed consent to participate in this study.

 Author Contributions

All authors listed have made a substantial, direct, and intellectual contribution to the work, and approved it for publication.

 Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

 Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

 Funding

This work was supported by grants from the National Natural Science Foundation of China (NSFC) (82001799 and 81771819).

 References Afrasiabi M. Redinbaugh M. J. Phillips J. M. Kambi N. A. Mohanta S. Raz A. (2021). Consciousness depends on integration between parietal cortex, striatum, and thalamus. Cell Syst. 12 363–373.e11. 10.1016/j.cels.2021.02.003 33730543 Albajes-Eizagirre A. Solanes A. Vieta E. Radua J. (2019). Voxel-based meta-analysis via permutation of subject images (PSI): theory and implementation for SDM. Neuroimage 186 174–184. 10.1016/j.neuroimage.2018.10.077 30389629 Albin R. L. Young A. B. Penney J. B. (1989). The functional anatomy of basal ganglia disorders. Trends Neurosci. 12 366–375. 10.1016/0166-2236(89)90074-x Amico E. Marinazzo D. Di Perri C. Heine L. Annen J. Martial C. (2017). Mapping the functional connectome traits of levels of consciousness. Neuroimage 148 201–211. 10.1016/j.neuroimage.2017.01.020 28093358 Boly M. Tshibanda L. Vanhaudenhuyse A. Noirhomme Q. Schnakers C. Ledoux D. (2009). Functional connectivity in the default network during resting state is preserved in a vegetative but not in a brain dead patient. Hum. Brain Mapp. 30 2393–2400. 10.1002/hbm.20672 19350563 Bonini F. Lambert I. Wendling F. Mcgonigal A. Bartolomei F. (2016). Altered synchrony and loss of consciousness during frontal lobe seizures. Clin. Neurophysiol. 127 1170–1175. 10.1016/j.clinph.2015.04.050 25912335 Braun A. R. Balkin T. J. Wesenten N. J. Carson R. E. Varga M. Baldwin P. (1997). Regional cerebral blood flow throughout the sleep-wake cycle. An H2(15)O PET study. Brain 120 1173–1197. 10.1093/brain/120.7.1173 9236630 Chen P. Xie Q. Wu X. Huang H. Lv W. Chen L. (2018). Abnormal effective connectivity of the anterior forebrain regions in disorders of consciousness. Neurosci. Bull. 34 647–658. 10.1007/s12264-018-0250-6 29959668 Crone J. S. Lutkenhoff E. S. Bio B. J. Laureys S. Monti M. M. (2017). Testing proposed neuronal models of effective connectivity within the cortico-basal ganglia-thalamo-cortical loop during loss of consciousness. Cereb. Cortex 27, 2727–2738. 10.1093/cercor/bhw112 27114177 Crone J. S. Schurz M. Holler Y. Bergmann J. Monti M. Schmid E. (2015). Impaired consciousness is linked to changes in effective connectivity of the posterior cingulate cortex within the default mode network. Neuroimage 110 101–109. 10.1016/j.neuroimage.2015.01.037 25620493 Crone J. S. Soddu A. Holler Y. Vanhaudenhuyse A. Schurz M. Bergmann J. (2014). Altered network properties of the fronto-parietal network and the thalamus in impaired consciousness. Neuroimage Clin. 4 240–248. 10.1016/j.nicl.2013.12.005 24455474 David O. Guillemain I. Saillet S. Reyt S. Deransart C. Segebarth C. (2008). Identifying neural drivers with functional MRI: an electrophysiological validation. PLoS Biol. 6:e315. 10.1371/journal.pbio.0060315 19108604 Dehaene S. Changeux J. P. Naccache L. Sackur J. Sergent C. (2006). Conscious, preconscious, and subliminal processing: a testable taxonomy. Trends Cogn. Sci. 10 204–211. 10.1016/j.tics.2006.03.007 16603406 Demertzi A. Gomez F. Crone J. S. Vanhaudenhuyse A. Tshibanda L. Noirhomme Q. (2014). Multiple fMRI system-level baseline connectivity is disrupted in patients with consciousness alterations. Cortex 52 35–46. 10.1016/j.cortex.2013.11.005 24480455 Deshpande G. Santhanam P. Hu X. (2011). Instantaneous and causal connectivity in resting state brain networks derived from functional MRI data. Neuroimage 54 1043–1052. 10.1016/j.neuroimage.2010.09.024 20850549 Farrell J. (2018). Higher-order theories of consciousness and what-it-is-like-ness. Philos. Stud. 175 2743–2761. 10.1007/s11098-017-0980-8 30393397 Freeze B. S. Kravitz A. V. Hammack N. Berke J. D. Kreitzer A. C. (2013). Control of basal ganglia output by direct and indirect pathway projection neurons. J. Neurosci. 33 18531–18539. 10.1523/JNEUROSCI.1278-13.2013 24259575 Friston K. Penny W. (2011). Post hoc Bayesian model selection. Neuroimage 56 2089–2099. 10.1016/j.neuroimage.2011.03.062 21459150 Giacino J. T. Kalmar K. Whyte J. (2004). The JFK coma recovery scale-revised: measurement characteristics and diagnostic utility. Arch. Phys. Med. Rehabil. 85 2020–2029. 10.1016/j.apmr.2004.02.033 15605342 Haber S. N. (2016). Corticostriatal circuitry. Dialog. Clin. Neurosci. 18 7–21. Hagmann P. Cammoun L. Gigandet X. Meuli R. Honey C. J. Wedeen V. J. (2008). Mapping the structural core of human cerebral cortex. PLoS Biol. 6:e159. 10.1371/journal.pbio.0060159 18597554 Hannawi Y. Lindquist M. A. Caffo B. S. Sair H. I. Stevens R. D. (2015). Resting brain activity in disorders of consciousness: a systematic review and meta-analysis. Neurology 84 1272–1280. 10.1212/WNL.0000000000001404 25713001 He J. H. Yang Y. Zhang Y. Qiu S. Y. Zhou Z. Y. Dang Y. Y. (2014). Hyperactive external awareness against hypoactive internal awareness in disorders of consciousness using resting-state functional MRI: highlighting the involvement of visuo-motor modulation. NMR Biomed. 27 880–886. 10.1002/nbm.3130 24820617 Huang L. Zheng Y. Zeng Z. Li M. Zhang L. Gao Y. (2019). Fractional amplitude of low-frequency fluctuations and functional connectivity in comatose patients subjected to resting-state functional magnetic resonance imaging. Ann. Indian Acad. Neurol. 22 203–209. 10.4103/aian.AIAN_420_17 33692645 Huang Z. Dai R. Wu X. Yang Z. Liu D. Hu J. (2014). The self and its resting state in consciousness: an investigation of the vegetative state. Hum. Brain Mapp. 35 1997–2008. 10.1002/hbm.22308 23818102 Jin F. Zheng P. Liu H. Guo H. Sun Z. (2018). Functional and anatomical connectivity-based parcellation of human cingulate cortex. Brain Behav. 8:e01070. 10.1002/brb3.1070 30039643 Kravitz A. V. Freeze B. S. Parker P. R. Kay K. Thwin M. T. Deisseroth K. (2010). Regulation of parkinsonian motor behaviours by optogenetic control of basal ganglia circuitry. Nature 466 622–626. 10.1038/nature09159 20613723 Lacey C. J. Bolam J. P. Magill P. J. (2007). Novel and distinct operational principles of intralaminar thalamic neurons and their striatal projections. J. Neurosci. 27 4374–4384. 10.1523/JNEUROSCI.5519-06.2007 17442822 Lancaster J. L. Woldorff M. G. Parsons L. M. Liotti M. Freitas C. S. Rainey L. (2000). Automated Talairach atlas labels for functional brain mapping. Hum. Brain Mapp. 10 120–131. 10.1002/1097-0193(200007)10:3<120::aid-hbm30<3.0.co;2-8 Laureys S. (2005). The neural correlate of (un)awareness: lessons from the vegetative state. Trends Cogn. Sci. 9 556–559. 10.1016/j.tics.2005.10.010 16271507 Laureys S. Schiff N. D. (2012). Coma and consciousness: paradigms (re)framed by neuroimaging. Neuroimage 61 478–491. 10.1016/j.neuroimage.2011.12.041 22227888 Li B. Wang X. Yao S. Hu D. Friston K. (2012). Task-dependent modulation of effective connectivity within the default mode network. Front. Psychol. 3:206. 10.3389/fpsyg.2012.00206 22737141 Liu X. Eickhoff S. B. Hoffstaedter F. Genon S. Caspers S. Reetz K. (2020). Joint multi-modal parcellation of the human striatum: functions and clinical relevance. Neurosci. Bull. 36 1123–1136. 10.1007/s12264-020-00543-1 32700142 Mashour G. A. Roelfsema P. Changeux J. P. Dehaene S. (2020). Conscious processing and the global neuronal workspace hypothesis. Neuron 105 776–798. 10.1016/j.neuron.2020.01.026 32135090 Massimini M. Ferrarelli F. Sarasso S. Tononi G. (2012). Cortical mechanisms of loss of consciousness: insight from TMS/EEG studies. Arch. Ital. Biol. 150 44–55. 10.4449/aib.v150i2.1361 23165870 Mhuircheartaigh R. N. Rosenorn-Lanng D. Wise R. Jbabdi S. Rogers R. Tracey I. (2010). Cortical and subcortical connectivity changes during decreasing levels of consciousness in humans: a functional magnetic resonance imaging study using propofol. J. Neurosci. 30 9095–9102. 10.1523/JNEUROSCI.5516-09.2010 20610743 Moeller F. Siebner H. R. Wolff S. Muhle H. Granert O. Jansen O. (2008). Simultaneous EEG-fMRI in drug-naive children with newly diagnosed absence epilepsy. Epilepsia 49 1510–1519. 10.1111/j.1528-1167.2008.01626.x 18435752 Monti M. M. Rosenberg M. Finoia P. Kamau E. Pickard J. D. Owen A. M. (2015). Thalamo-frontal connectivity mediates top-down cognitive functions in disorders of consciousness. Neurology 84 167–173. 10.1212/WNL.0000000000001123 25480912 Muller V. I. Cieslik E. C. Laird A. R. Fox P. T. Radua J. Mataix-Cols D. (2018). Ten simple rules for neuroimaging meta-analysis. Neurosci. Biobehav. Rev. 84 151–161. 10.1016/j.neubiorev.2017.11.012 29180258 Nakhnikian A. Rebec G. V. Grasse L. M. Dwiel L. L. Shimono M. Beggs J. M. (2014). Behavior modulates effective connectivity between cortex and striatum. PLoS One 9:e89443. 10.1371/journal.pone.0089443 24618981 Norton L. Hutchison R. M. Young G. B. Lee D. H. Sharpe M. D. Mirsattari S. M. (2012). Disruptions of functional connectivity in the default mode network of comatose patients. Neurology 78 175–181. 10.1212/WNL.0b013e31823fcd61 22218274 Palanca B. Avidan M. S. Mashour G. A. (2017). Human neural correlates of sevoflurane-induced unconsciousness. Br. J. Anaesth. 119 573–582. 10.1093/bja/aex244 29121298 Parnaudeau S. Bolkan S. S. Kellendonk C. (2018). The mediodorsal thalamus: an essential partner of the prefrontal cortex for cognition. Biol. psychiatry 83 648–656. 10.1016/j.biopsych.2017.11.008 29275841 Peräkylä J. Sun L. Lehtimäki K. Peltola J. Öhman J. Möttönen T. (2017). Causal evidence from humans for the role of mediodorsal nucleus of the thalamus in working memory. J. Cogn. Neurosci. 29 2090–2102. 10.1162/jocn_a_01176 Quan M. Lee S. H. Kubicki M. Kikinis Z. Rathi Y. Seidman L. J. (2013). White matter tract abnormalities between rostral middle frontal gyrus, inferior frontal gyrus and striatum in first-episode schizophrenia. Schizophr. Res. 145 1–10. 10.1016/j.schres.2012.11.028 23415471 Radua J. Mataix-Cols D. (2009). Voxel-wise meta-analysis of grey matter changes in obsessive-compulsive disorder. Br. J. Psychiatry 195 393–402. 10.1192/bjp.bp.108.055046 19880927 Radua J. Mataix-Cols D. Phillips M. L. El-Hage W. Kronhaus D. M. Cardoner N. (2012). A new meta-analytic method for neuroimaging studies that combines reported peak coordinates and statistical parametric maps. Eur. Psychiatry 27 605–611. 10.1016/j.eurpsy.2011.04.001 21658917 Rikhye R. V. Gilra A. Halassa M. M. (2012). Thalamic regulation of switching between cortical representations enables cognitive flexibility. Nat. Neurosci. 21 1753–1763. 10.1038/s41593-018-0269-z 30455456 Rosa M. J. Friston K. Penny W. (2012). Post-hoc selection of dynamic causal models. J. Neurosci. Methods 208 66–78. 10.1016/j.jneumeth.2012.04.013 22561579 Rosazza C. Andronache A. Sattin D. Bruzzone M. G. Marotta G. Nigri A. (2016). Multimodal study of default-mode network integrity in disorders of consciousness. Ann. Neurol. 79 841–853. 10.1002/ana.24634 26970235 Salvato G. Berlingeri M. De Maio G. Curto F. Chieregato A. Magnani F. G. (2020). Autonomic responses to emotional linguistic stimuli and amplitude of low-frequency fluctuations predict outcome after severe brain injury. Neuroimage Clin. 28:102356. 10.1016/j.nicl.2020.102356 32750635 Schiff N. D. (2010). Recovery of consciousness after brain injury: a mesocircuit hypothesis. Trends Neurosci. 33 1–9. 10.1016/j.tins.2009.11.002 19954851 Schmitt L. I. Wimmer R. D. Nakajima M. Happ M. Mofakham S. Halassa M. M. (2017). Thalamic amplification of cortical connectivity sustains attentional control. Nature 545 219–223. 10.1038/nature22073 28467827 Seth A. K. Izhikevich E. Reeke G. N. Edelman G. M. (2006). Theories and measures of consciousness: an extended framework. Proc. Natl. Acad. Sci. U.S.A. 103 10799–10804. 10.1073/pnas.0604347103 16818879 Silva S. de Pasquale F. Vuillaume C. Riu B. Loubinoux I. Geeraerts T. (2015). Disruption of posteromedial large-scale neural communication predicts recovery from coma. Neurology. 85 2036–2044. 10.1212/WNL.0000000000002196 26561296 Slagter H. A. Mazaheri A. Reteig L. C. Smolders R. Figee M. Mantione M. (2017). Contributions of the ventral striatum to conscious perception: an intracranial EEG study of the attentional blink. J. Neurosci. 37 1081–1089. 10.1523/JNEUROSCI.2282-16.2016 27986925 Smith Y. Bevan M. D. Shink E. Bolam J. P. (1998). Microcircuitry of the direct and indirect pathways of the basal ganglia. Neuroscience 86 353–387. 10.1016/s0306-4522(98)00004-9 Soch J. Deserno L. Assmann A. Barman A. Walter H. Richardson-Klavehn A. (2017). Inhibition of information flow to the default mode network during self-reference versus reference to others. Cereb. Cortex 27 3930–3942. 10.1093/cercor/bhw206 27405334 Soddu A. Gomez F. Heine L. Di Perri C. Bahri M. A. Voss H. U. (2016). Correlation between resting state fMRI total neuronal activity and PET metabolism in healthy controls and patients with disorders of consciousness. Brain Behav. 6:e00424. 10.1002/brb3.424 27110443 Tahmasian M. Sepehry A. A. Samea F. Khodadadifar T. Soltaninejad Z. Javaheripour N. (2019). Practical recommendations to conduct a neuroimaging meta-analysis for neuropsychiatric disorders. Hum. Brain Mapp. 40 5142–5154. 10.1002/hbm.24746 31379049 Tang L. Ge Y. Sodickson D. K. Miles L. Zhou Y. Reaume J. (2011). Thalamic resting-state functional networks: disruption in patients with mild traumatic brain injury. Radiology 260 831–840. 10.1148/radiol.11110014 21775670 Threlkeld Z. D. Bodien Y. G. Rosenthal E. S. Giacino J. T. Nieto-Castanon A. Wu O. (2018). Functional networks reemerge during recovery of consciousness after acute severe traumatic brain injury. Cortex 106 299–308. 10.1016/j.cortex.2018.05.004 29871771 Tononi G. Boly M. Massimini M. Koch C. (2016). Integrated information theory: from consciousness to its physical substrate. Nat. Rev. Neurosci. 17 450–461. 10.1038/nrn.2016.44 27225071 van den Heuvel M. P. Sporns O. (2011). Rich-club organization of the human connectome. J. Neurosci. 31 15775–15786. 10.1523/JNEUROSCI.3539-11.2011 22049421 Vanhaudenhuyse A. Noirhomme Q. Tshibanda L. J. Bruno M. A. Boveroux P. Schnakers C. (2010). Default network connectivity reflects the level of consciousness in non-communicative brain-damaged patients. Brain 133 161–171. 10.1093/brain/awp313 20034928 Wu X. Xie Q. Liu X. Huang H. Ma Q. Wang J. (2019). Spatially overlapping regions show abnormal thalamo-frontal circuit and abnormal precuneus in disorders of consciousness. Brain Topogr. 32 445–460. 10.1007/s10548-018-0693-0 30707390 Yan C. G. Wang X. D. Zuo X. N. Zang Y. F. (2016). DPABI: data processing & analysis for (resting-state) brain imaging. Neuroinformatics 14 339–351. 10.1007/s12021-016-9299-4 27075850 Yao S. Song J. Gao L. Yan Y. Huang C. Ding H. (2015). Thalamocortical sensorimotor circuit damage associated with disorders of consciousness for diffuse axonal injury patients. J. Neurol. Sci. 356 168–174. 10.1016/j.jns.2015.06.044 26165776 Yarkoni T. Poldrack R. A. Nichols T. E. Van Essen D. C. Wager T. D. (2011). Large-scale automated synthesis of human functional neuroimaging data. Nat. Methods 8 665–670. 10.1038/nmeth.1635 21706013 Zhou J. Liu X. Song W. Yang Y. Zhao Z. Ling F. (2011). Specific and nonspecific thalamocortical functional connectivity in normal and vegetative states. Conscious. Cogn. 20 257–268. 10.1016/j.concog.2010.08.003 21078562 Zhou Y. Lui Y. W. Zuo X. N. Milham M. P. Reaume J. Grossman R. I. (2014). Characterization of thalamo-cortical association using amplitude and connectivity of functional MRI in mild traumatic brain injury. J. Magn. Reson. Imaging 39 1558–1568. 10.1002/jmri.24310 24014176

https://pubmed.ncbi.nlm.nih.gov/

https://www.webofknowledge.com/

https://www.sdmproject.com/

http://www.neurosynth.org

http://www.rfmri.org/dpabi

http://www.fil.ion.ucl.ac.uk/spm/