Next Article in Journal
Cellular Changes in Injured Rat Spinal Cord Following Electrical Brainstem Stimulation
Previous Article in Journal
Characterizing Fatigue-Related White Matter Changes in MS: A Proton Magnetic Resonance Spectroscopy Study
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Brain Sciences
Volume 9
Issue 5
10.3390/brainsci9050123
brainsci-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Feasibility of Non-Gaussian Diffusion Metrics in Chronic Disorders of Consciousness

by
Elena I. Kremneva
1,*,
Liudmila A. Legostaeva
1,
Sofya N. Morozova
1,
Dmitry V. Sergeev
1,
Dmitry O. Sinitsyn
1,
Elizaveta G. Iazeva
1,
Aleksandr S. Suslin
1,
Natalia A. Suponeva
1,
Marina V. Krotenkova
1,
Michael A. Piradov
1 and
Ivan I. Maximov
2,3,*
1
Research Center of Neurology, 80 Volokolamskoe shosse, 125367 Moscow, Russia
2
Department of Psychology, University of Oslo, Forskningsveien 3A, 0373 Oslo, Norway
3
Norwegian Centre for Mental Disorders Research (NORMENT), Norway and Institute of Clinical Medicine, University of Oslo, Oslo Universitetssykehus Bygg 48 Ullevål, 0317 Oslo, Norway
*
Authors to whom correspondence should be addressed.
Brain Sci. 2019, 9(5), 123; https://doi.org/10.3390/brainsci9050123
Submission received: 23 April 2019 / Revised: 23 May 2019 / Accepted: 23 May 2019 / Published: 27 May 2019
(This article belongs to the Section Neurotechnology and Neuroimaging)
Browse Figures
Review Reports Versions Notes

Abstract

:
Diagnostic accuracy of different chronic disorders of consciousness (DOC) can be affected by the false negative errors in up to 40% cases. In the present study, we aimed to investigate the feasibility of a non-Gaussian diffusion approach in chronic DOC and to estimate a sensitivity of diffusion kurtosis imaging (DKI) metrics for the differentiation of vegetative state/unresponsive wakefulness syndrome (VS/UWS) and minimally conscious state (MCS) from a healthy brain state. We acquired diffusion MRI data from 18 patients in chronic DOC (11 VS/UWS, 7 MCS) and 14 healthy controls. A quantitative comparison of the diffusion metrics for grey (GM) and white (WM) matter between the controls and patient group showed a significant (p < 0.05) difference in supratentorial WM and GM for all evaluated diffusion metrics, as well as for brainstem, corpus callosum, and thalamus. An intra-subject VS/UWS and MCS group comparison showed only kurtosis metrics and fractional anisotropy differences using tract-based spatial statistics, owing mainly to macrostructural differences on most severely lesioned hemispheres. As a result, we demonstrated an ability of DKI metrics to localise and detect changes in both WM and GM and showed their capability in order to distinguish patients with a different level of consciousness.

1. Introduction

The term “consciousness” in clinical practice is defined as the presence of two main components: wakefulness and awareness [1]. If wakefulness is present, but signs of awareness are absent, the patient is referred to having vegetative state/unresponsive wakefulness syndrome (VS/UWS). On the contrary, if any evidence of awareness along with wakefulness cycle can be seen, a clinician may diagnose a minimally conscious state (MCS) [2,3]. The mentioned two states are difficult to differentiate unambiguously since clinical evaluation is often subjective and based solely on clinical behavioural scales that can be influenced by sensorimotor impairment, unnoticeable motor activity, pain, apraxia, aphasia, deafness etc. [4]. Even with the use of established diagnostic criteria, the rate of misdiagnosis remains as high as almost 40% [5,6]. Since the correct diagnosis in disorders of consciousness (DOC) is of great importance for a selection of rehabilitation strategy, recovery prognosis and support of patient family, many researchers try to find reliable objective biomarkers to estimate the level of consciousness.
Multimodal neurophysiological and imaging methods such as electroencephalography (EEG), transcranial magnetic stimulation-EEG (TMS-EEG), functional MRI (fMRI), MR-spectroscopy, and positron emission tomography (PET) are widely used for assessments of brain structure, its function and metabolism in terms of residual consciousness and cognitive function [7,8,9]. For example, PET studies showed that DOC patients exhibited a decrease in brain regional metabolism of up to 40% of the normal value, and this decrease differs in VS/UWS versus MCS patients [1].
With the increasing number of similar studies, it had been revealed that there is no a single structural biomarker for consciousness differentiation, but many of them. The stage of DOC is determined by the severity and number of lesions affecting the brain due to diffuse axonal injury, multifocal cortical contusions, laminar necrosis, thalamic damage, and necrosis [10,11]. Moreover, it was suggested that some regions might serve as the main brain “hubs”, i.e., brain areas that are more important for the function of consciousness than the others [12]. For example, changes of the brain metabolism in DOC patients were most prevalent in regions including frontoparietal areas, midline anterior cingulate/mesiofrontal and posterior cingulate/precuneal associative cortices [13,14]. Damage to these brain “hubs” exhibit a strong correlation with resting-state functional MRI data, and significantly lower activity of the default mode network in DOC [7,15]. Thus, within- and between-network functional disconnections regarding large-scale functional networks of the human brain are considered to be a correlate of unconscious states [1,16,17].
Understanding the origin of the microstructural changes of brain tissues underpinning functional disconnections is also crucial for unravelling the mechanisms of DOC [18]. Conventional MRI techniques may detect and estimate the gross pathological processes leading to coma and other DOC states. There were several attempts to evaluate the DOC prognosis using the standard MRI contrasts (T2- and T1-weighted images) based on lesion distribution, pathology, and amount of lesions [19,20]. However, popular MRI modalities failed to detect microstructural changes in the case of DOC patients with minimal brain lesions. Diffusion-weighted MRI (DW-MRI) techniques are used for localisation and estimation of brain changes typically not visible by the conventional T1- or T2-weighted MR imaging. Correlation between functional and microstructural changes in various DOC stages was confirmed earlier in the study that utilised direct comparison of PET and DW-MRI [21]. Thus, the diffusion imaging technique might be a useful technique for differentiation of VS/UWS from MCS patients [22,23,24].
The most widely used diffusion approach in clinical neuroimaging is diffusion tensor imaging (DTI). For example, in DOC patients it has been applied for the prediction of long-term neurological outcomes in acute stage [25], for brain microstructure damage in both acute and chronic stages, including long-term prospective observation [26,27]. DTI assumes that water diffusion can be described by a Gaussian probability distribution in biological tissue [28]. However, water molecules in biological tissue experience complex non-Gaussian processes due to the presence of different barriers such as cell membranes, complex distribution of geometrical constraints etc. Considering this limitation, one can use higher orders of the cumulant expansion of the diffusion propagator in order to partially take into account the complexity of the diffusion processes, for example, by applying diffusional kurtosis imaging (DKI) [29]. In DKI, in addition to the diffusion tensor, deviation from the Gaussian hypothesis is taken into account by the kurtosis of the diffusion displacement probability density function. Thereby, DKI provides a more accurate tissue microstructure characterisation [30] and could be sensitive biomarker in in-vivo studies [31,32,33,34,35,36,37]. Moreover, DKI does not demand a presence of high anisotropic environment such as WM and permits the quantification of the microstructural integrity in general, for example, in the case of GM [29,38,39,40,41].
In the present study, we aimed to investigate the question of whether the diffusion kurtosis metrics provide enriched information and a more sensitive chronic DOC differentiation criterion compared to the DTI metrics regarding tissue changes following anoxia and severe chronic TBI. At the same time, we would like to assess the feasibility of DKI metrics as a source of biomarkers for microstructure changes in grey matter in DOC patients.

2. Materials and Methods

2.1. Patients

We recruited 23 patients who fulfilled the criteria of chronic DOC (VS/UWS, MCS) of traumatic (≥12 months after accident) or anoxic (≥6 months after accident) origin. However, 5 patients were excluded from the study due to severe motion artifacts (4 patients) or excessive brain deformation due to decompressive craniectomy with massive superficial siderosis (1 patient). As a result, we used data obtained from 18 patients (see Figure 1 as an example of patient images).
In Table 1, we summarised demographic and clinical characteristics of the patients. Control group included 18 age-matched healthy volunteers but 4 subjects were excluded due to a lack of image quality (7 males, 7 females; mean age 40.2 ± 12.8 years). In patients, the clinical examination was performed using the standard tool (Coma Recovery Scale—revised (CRS-R)) on the day of scanning. CRS-R comprises 23 hierarchically arranged items associated with brainstem, subcortical and cortical processes [2]. The lowest score on each subscale represents reflexive activity, while the highest item represents cognitively mediated behaviour by addressing to auditory, visual, motor, oromotor, communicative and arousal functions. The CRS-R is established as the most reliable clinical tool for chronic DOC assessment [42].
Whenever it was possible, MRI acquisition was performed without patient sedation. However, in order to avoid severe motor artefacts during the scanning in 8 patients (44%), an anaesthesiologist induced light sedation by dexmedetomidine administration via intravenous infusion at a constant rate of 1 μg/kg·h−1. During the infusion period, the anaesthesiologist monitored the patient’s cuff blood pressure, the electrocardiogram, and the pulse oximetry.
The local ethical committee approved the study. Informed consent was obtained from the legal representatives of patients and from healthy volunteers before any study-related procedures.

2.2. MRI Measurements

MRI data were acquired using a Siemens MAGNETOM Verio 3T scanner (Siemens Medical Systems, Erlangen, Germany) with a standard 32-channel matrix head coil. The patient and the control groups underwent standard axial T2-weighed imaging (TR 4000 ms; TE 118 ms; slice thickness 5.0 mm; in plane resolution 1.5 mm2; duration: 2 min 02 s) in order to evaluate an extend of the brain damage for DOC patients and to exclude any pathology in controls with further diffusion MRI examination using the productive spin-echo echo-planar imaging sequence with three b-values (0, 1000, and 2500 s/mm2) and 64 non-coplanar diffusion directions for each non-zero b-value (TE/TR 115/12600 ms; matrix 100 × 100, resolution 2 × 2 × 2 mm3; number of excitation was equal 1, GRAPPA acceleration factor 2). To improve the signal-to-noise ratio and corrections of the susceptibility-based distortions, we acquired two non-diffusion-weighted images with different phase-encoding directions. In addition, a 3D T1-weighted anatomical image (MPRAGE) was acquired (TE/TR 2.47/1900 ms, TI 900 ms, flip angle 9°, matrix 256 × 256, 176 axial slices to cover the whole brain) with an isotropic voxel resolution of 1.0 × 1.0 × 1.0 mm3.

2.3. Image Analysis

Diffusion image post-processing prior to estimation of the diffusion scalar metrics was done as follows: first, we performed a noise correction for all diffusion images by assuming the non-central χ2 distribution due to employed parallel imaging [43]. Next, the data were corrected for susceptibility and eddy-current induced distortions and subject motions using the topup and eddy utility from the FSL package [44]. A correction of the diffusion directions was performed in accordance with [45]. In order to suppress the Gibbs ringing effect [46] and decrease a possible presence of outliers [31,47] we applied Gaussian smoothing with a kernel of 1.5 × 1.5 × 1.5 mm3. All diffusion metrics were estimated using the ExploreDTI package [48] using weighted linear estimator. As a result, we obtained fractional anisotropy (FA), mean, radial, and axial diffusivities (MD, RD, AD) and mean, radial, and axial kurtosis (MK, RK, AK) metrics [30,49]. For convenience, we have to note that the diffusion scalar metrics estimated in a frame of the kurtosis model might differ from the estimations of the conventional DTI modality [50].
For the region-based analysis, we used the original T1-weighted image for segmentation of the grey and white matter applying the FAST utility from the FSL package [51] with further images binarization and creating GM/WM masks for each subject. Then, each mask was visually inspected in ITK-SNAP viewer [52] with manual correction of the binary masks if necessary. The non-diffusion weighted original image of each subject was aligned and interpolated to its T1-weighted reference image using affine transformation and cubic spline interpolation with a normalised mutual information as a quality criterion. The co-registration procedure was visually inspected and applied for all diffusion-derived scalar metrics. The interpolation step allows one to improve the statistical reliability of the following analysis [31,33].
The brainstem and thalamus regions of interests (ROI) were manually segmented using ITK-SNAP [52]. The brainstem ROI included area from craniovertebral junction to the anterior-posterior commissural line, thereafter the thalamus ROI was outlined. The corpus callosum (CC) ROI included 3 continuous mid-sagittal slices of the CC, outlined on T1-weighted images. An outline of our post-processing approach is summarised in Figure 2.

2.4. Statistical Analysis

We performed statistical analysis based on a group comparison between healthy subjects and three patient subgroups, namely, all patients, MSC patients, and VS/UWS patients, using tract-based spatial statistics (TBSS) [53]. The most representative subject’s structural image was used and aligned to the Montreal Neurological Institute (MNI) 152 standard space using non-linear co-registration utility from FSL. A mean FA image was built up from FA maps of all subjects and thinned in order to produce a mean FA skeleton. We used an FA threshold equal to 0.2 for the FA skeleton to avoid voxels that could be treated as grey matter or as cerebrospinal fluid. Next, all FA images were projected onto the skeleton. At the end, all data were prepared for the voxel-wise statistical analysis. The statistical analysis was performed using permutation-based, voxel-wise non-parametric testing, implemented via the “randomise”-function with 5000 permutations in FSL package. In TBSS analysis we used threshold-free cluster enhancement approach. Resulting p-values were corrected for multiple comparisons and thresholded at p < 0.05. Non-FA images were prepared in a similar way for the statistical tests in according with FSL guidance.
Subsequently, we performed a statistical analysis of the ROI derived diffusion metrics using in-house Matlab scripts (TheMathWorks, Natick, MA, USA) based on the non-parametric Wilcoxon rank sum test.

3. Results

3.1. Patients versus Control Group

The diffusion scalar metrics comparison was performed for supratentorial grey and white matter, and separately for the thalamus, corpus callosum and brainstem in control and patient groups (Figure 3, Table 2). In supratentorial WM, we found a significant difference (p < 0.05) in the DTI metrics (FA, MD, AD, RD) between control and patient groups, as well as in the DKI metrics (MK, AK, RK). The largest effect was detected for FA values in the case of WM (Cohen’s d = 3.28). This result is clearly confirmed by a TBSS analysis (see Figure 4). Due to remarkable brain damage, FA as well as MK, AK, RK in patients was significantly lower, with higher MD, AD, and RD comparing to healthy volunteers. As for GM, the same changes pattern was evident, but without any significant differences for FA (see Table 2). When analysing the thalamus separately, the most reliable difference occurred in MD, AD and RD maps (p < 0.001), with higher values in patient group. Absence of significant FA inter-group difference in the thalamus could be influenced by manual segmentation in the patient group since the exact ROI borders are difficult to outline due to severe brain architecture damages in most patients. Then we repeated the analysis taking a smaller thalamus ROI as 3 continuous circle ROIs (diameter 8 mm) in the middle part of both thalamus sides on the level of intrathalamical adhesion. We obtained the same result as for the whole thalamus analysis, with no significant difference in the FA metric. The brainstem analysis revealed significant differences in all metrics, except for RK, with the FA showing the highest values in control group comparing to patients (Cohen’s d = 4.66). There were significant differences in all diffusion metrics in corpus callosum, which is the most coherent WM tract, with remarkably higher difference for the DKI metrics comparing to the DTI metrics.

3.2. Anoxia versus Trauma Patient Groups

The DOC intra-group comparison based on anoxic versus trauma brain lesion aetiology showed a significant difference only for AK in GM, with anoxic AK values being higher than the traumatic ones (0.55 versus 0.46, Cohen’s d = 1.13) (see Figure 5a). The TBSS analysis did not reveal any significant differences in the case of Anoxia versus Trauma group comparison for all metrics.

3.3. VS/UWS versus MCS Patients Groups

There was no significant difference between the groups when we compared the metrics using ROI-based analysis (see Figure 5b), but we found a clear tract-specific difference in the TBSS analysis (Figure 6). TBSS analysis showed that DOC intra-group comparison possesses a significant difference (p < 0.05) in FA and kurtosis metrics in all types of fibres projection, commissural, associative (Figure 6a), predominantly in the left hemisphere. An increased significance level (p < 0.001) reduced a number of useful diffusion metrics and only the kurtosis metrics (MK and AK) localised the difference in the left hemisphere (Figure 6b).

4. Discussion

To the best of our knowledge, the present study is a first attempt to apply the DKI approach in order to investigate non-Gaussian diffusion patterns in the brain and associated tissue degradation in chronic disorders of consciousness. We performed a group analysis based on “patients versus controls”, “traumatic versus anoxic patients”, and “VS/UWS versus MCS patients” comparisons. Firstly, we assessed microstructure changes using diffusion kurtosis metrics in DOC patients in contrast to the healthy subjects. In turn, in different subgroups of DOC patients, we investigated a question whether advanced diffusion imaging method works in the frame of severe brain lesions. Therefore, these group comparisons have a methodological implication allowed us to differentiate microstructure alternations for inter- and intra-DOC groups. Moreover, clinically related difficulties in distinguishing different levels of consciousness inspired us to search more reliable biomarkers among diffusion metrics. Taking into account a different aetiology of the patient subgroups, we aimed to assess in vivo variations in neuropathology of traumatic brain injuries and anoxia that might improve a prognostication and to help in treatment planning.
We demonstrated the ability of DKI metrics to detect and localise the changes in both WM and GM and showed their feasibility to distinguish the patients with different forms of chronic DOC. Our results replicated previous studies utilising DTI metrics only [22,23,24,54] and allowed us to assess and differentiate white matter damages in DOC. The results of TBSS analyses suggested that microstructure of almost all WM tracts is dramatically affected in DOC. As a result, the most important structural connectivity disruption patterns were intricate and problematic for predicting the level of consciousness. Along with WM changes, extensive cortex, thalamus and brainstem lesions of different aetiology and their combinations contributed to the consciousness decrease, in accordance with [55]. The use of advanced diffusion approach such as DKI allowed us to evaluate not only WM tract degradation but to perform a GM assessment as well. For example, significant AK differences were obtained for supratentorial GM in DOC subjects of different aetiology. Such diffusion approach seems to be a more robust and reliable method taking into account severe brain damage and deformation of brain structures in chronic DOC patients. It allows us to avoid additional mistakes in fibre tract delineation and its reconstruction, in particular, recalling group analysis methods with brain normalisation algorithms.
The analysis based on kurtosis as well as diffusion tensor metrics showed here a good impact in distinguishing DOC patients from healthy subjects since not so many diffusion kurtosis research has been performed in clinical settings so far [30,56,57,58,59], including only mild TBI kurtosis studies in acute period where only kurtosis metrics revealed alterations in some brain structures that correlated with cognitive impairment [60,61,62,63]. In our study we found dramatic diffusion metric changes in all regions of interest such as increased MD, AD, RD and decreased kurtosis values in patients comparing to the healthy controls (for example, see Figure 4 and Figure 5). The found tissue degradation with confirmed extensive brain damages in chronic DOC patients affected both WM and GM (see Figure 2 and Figure 3). Notably, all used diffusion metrics represent integrative values and do not reflect the unique microstructure changes and its causes. One can expect the mutually adjusted diffusion changes (for example, MD increase and FA decrease) correlated with more severe brain damage and lower level of consciousness in chronic DOC patients [22,25,26]. We should emphasise that our results are in agreement with the previous findings based on DTI metrics only [22,23,24,55], and provided additional information based on the kurtosis metrics.
FA is one of the important and simple summaries allowing one to probe microstructural integrity of coherent fibres. In the case of DOC patients, the FA shows the highest difference in WM and is associated with severe disruption. Nevertheless, FA does not work for GM (in both cortex and deep grey matter), since there is no high anisotropic tissue environment detectable by DTI at that scale. In opposite, all kurtosis metrics of supratentorial GM were higher in healthy controls as a measure of heterogeneity and complexity of the GM microstructure [36]. At the same time, increased values of diffusivity (MD, AD, and RD) in patient’s GM reflected significant GM degradation. Kurtosis metrics did not reveal any difference in thalamus in contrast to the detected changes of MD, AD and RD values. This discrepancy can be explained by a phenomenological background of the DTI and DKI metrics coming from the cumulant expansion of the diffusion propagator. Thus, an absence of underlying biophysical models based on DTI or DKI allows one to search empiric correlations associated with tissue changes only.
It is well known that thalamus is severely affected in DOC patients [64]. We found that changes of increased MD, AD and RD values could be a criterion for thalamus degradation and are consistent with the data obtained in [23]. However, severely damaged tissue rarefication due to cells and fibres destruction can be partly compensated by reactive astrogliosis and scar tissue formation with collagen and reticulin fibres laid own, the whole appearing as a glio-mesodermal reaction [65]. These changes might lead to complex microstructure with no significant difference for FA and kurtosis metrics between healthy controls and patients similar to our findings for kurtosis metrics. From this point of view, the reduced MK values can be associated with a loss of cellular structures [60]. In turn, various nuclei of the thalami were altered differently due to cells loss or cells atrophy in association with reactive changes in microglia, macrophages, and astrocytes [66]. Thus, we assume that a ROI-specified analysis of different thalamus regions might be more accurate but tricky because of severe thalamus deformation and size changes. The same explanation can be applied for the inter-group difference for MD, as a possible inverse measure of the membrane density since MD metric could be more influenced by tissue gliosis and reorganisation of membrane barriers.
Another interesting finding was that AD was increased in all ROIs in DOC patients. Decreased AD metrics in clinical neuroradiology are often considered to be associated with axonal damages and disruption of white matter bundles [67,68], which certainly has place in chronic DOC patients. However, other studies, including simulations and synthetic models, strongly discourage researchers from interpreting changes of the “axial” and “radial” diffusivities on the basis of the underlying tissue structure, unless accompanied by a thorough investigation of their mathematical and geometrical properties in each studied dataset [69]. The direction of the “axial” diffusivity would not be always preserved in pathological tissue and aligned with the underlying real tissue architecture [70]. Thus, explanation could include an additional interpretation such as axonal loss in measured area with interstitial fluid increase, i.e., by cerebrospinal fluid area [41]. As a result, all diffusion eigenvalues are increased with leading increase of AD, RD and MD metrics.
The highest inter-group difference was observed in corpus callosum where the kurtosis metrics showed more reliable changes (Cohen’s d = 7.52 for RK). The radial diffusion metrics (RD, RK) changes in CC experienced prevalence over axial ones. It could be associated with a profound demyelination. Thus, RK may serve as a new method for examining the integrity of the cell membrane and the surrounding myelin sheaths [36,41]. On the other hand, the analysis did not reveal a significant difference in RK for the brainstem between healthy and DOC patients as we found for RD. Morphological changes included not only demyelination and the presence of myelin debris, but axonal damage leading to decrease in the axonal diameter and reactive gliosis. Individual components of this process and the combination of morphological components have different effects on the behaviour of water molecule diffusion. It can lead to intricate diffusion processes in both residual healthy and severe damaged tissue [41]. At the same time, the presence of residual short white matter tracts following other directions than superior-inferior or inferior-superior, and so on, could influence the result and partially explain the absence of RK differences of the brainstem metrics in contrast to CC one.
ROI-based analysis between patients with different aetiology or stage of consciousness did not reveal any significant differences. However, we should mention that possible difficulties in the detection of such difference might originate from the fact that a different combination of brain lesion patterns (in TBI—mainly diffuse axonal injury and brainstem lesions, in anoxia—cortical lesions and selective necrosis of thalamus) may lead to the same clinical presentation of DOC [55]. For example, VS/UWS can be the result of either severe neocortical diffuse damage in asphyxia as well as of severe bilateral thalamic damage due to cardiopulmonary arrest with much of the neocortex spared [11]. Recent DKI studies in traumatic brain injury (TBI) performed either on patients or on animal models, showed that kurtosis metrics correspond to active processes involving reactive astrocytes, in particular, in GM, which cannot be revealed with other MR imaging techniques [61,71]. Diffuse damages of subcortical white matter are common to all chronic DOC patients. In such a way, despite structurally intact cortex the brain function is disabled due to the impairment of connections between different cortical areas via the thalamic nuclei, as DOC is associated with extensive damages to afferent and efferent cerebral connections [64]. In turn, the TBSS analysis did not find any differences based on aetiology. However, it detected WM changes between VS/UWS and MCS groups (see Figure 6). Moreover, the results of TBSS comparison in the case of kurtosis metrics demonstrated a higher sensitivity to structural changes in these subgroups, in particular, with increasing significance level (Figure 6b). The difference is seen mostly in the left hemisphere, and this side-predominance could originate from the fact that the most severe brain lesions were localised in the left hemisphere in 5 VS/UWS patient with the right side being predominantly affected in 2 MCS patients. Thus, a group concerned limitation based on non-uniform brain lesion pattern can affect the TBSS results. However, a found differentiation between MCS and VS/UWS patients offers a great opportunity for non-Gaussian models to become a sensitive biomarker of level of consciousness in chronic DOC in doubtful cases.
The major limitation of this study is the small sample size. A recruiting of a larger patient cohort would certainly improve a statistical power and diminish spurious or unreliable findings in some anatomical regions. However, we have to note that DOC patients demand a special treatment and operation. It presents some difficulties in recruiting and measurements. We hope to increase the sample size in further ongoing studies. Another problem is an image manipulation and data processing related to the brains with massive lesions and severe deformations [72]. Currently, there is no reliable method of DKI/DTI post-processing similar to reference [73], which allows one to validate and verify a robustness of numerical algorithms in the case of heavy brain deformation/deflection. Thus, it might introduce a larger variability in the results. Moreover, severe brain atrophy with subsequent cerebrospinal fluid enlargement increases bulk motion due to periodic tissue motion caused by cardiac pulsatility, especially in the brainstem. As a result, it could influence the final diffusion metrics [74]. There is also not enough of histological confirmation of the diffusion imaging findings in DOC patients, since the present studies are based on animal models [75] or report a very few number of human observations [76]. As a consequence, subsequent longitudinal studies are strongly required, with DKI/DTI multiple time points measurement, from an acute stage of brain lesion to the chronic DOC with morphology confirmation.

5. Conclusions

In conclusion, we found that the scalar kurtosis metrics represent additional biomarkers, in particular, for detection of micro- and macrostructural brain changes in the case of chronic DOC. AK was found to be sensitive for DOC aetiology differentiation and MK and AK—for VS/UWS versus MCS comparison. The application of non-Gaussian diffusion metrics yield enriched information about status and perspectives of chronic DOC patients. Since the DKI dataset holds DTI metrics as a part of its total diffusion information, it might be very effective to use both types of metrics. Thus, it might help clinicians to improve the diagnostic approaches in DOC and to estimate the prognosis with higher accuracy.

Author Contributions

Conceptualization, I.I.M. and E.I.K.; Data Curation, S.N.M., E.I.K., L.A.L. and E.G.I.; Methodology, I.I.M. and E.I.K.; CRS-R rating, L.A.L., E.G.I., D.V.S.; Formal Analysis, E.I.K., A.S.S., I.I.M. and D.O.S.; Software, I.I.M., A.S.S.; Investigation, S.N.M., E.I.K., L.A.L., E.G.I., D.V.S.; Validation, E.I.K., D.O.S., D.V.S., N.A.S.; Resources, N.A.S., and M.V.K.; Writing—Original Draft Preparation, E.I.K. and I.I.M.; Writing—Review & Editing, D.V.S., D.O.S., N.A.S., M.A.P.; Supervision, N.A.S., M.V.K. and M.A.P.; Project Administration, M.A.P.; Funding Acquisition, M.A.P.

Funding

This study was funded by Russian Science Foundation grant No. 16-15-00274.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Laureys, S. The neural correlate of (un)awareness: Lessons from the vegetative state. Trends Cogn. Sci. 2005, 9, 556–559. [Google Scholar] [CrossRef]
  2. Giacino, J.T.; Ashwal, S.; Childs, N.; Cranford, R.; Jennett, B.; Katz, D.I.; Kelly, J.P.; Rosenberg, J.H.; Whyte, J.; Zafonte, R.D.; et al. The minimally conscious state: Definition and diagnostic criteria. Neurology 2002, 58, 349–353. [Google Scholar] [CrossRef]
  3. Bernat, J.L. Chronic disorders of consciousness. Lancet 2006, 367, 1181–1192. [Google Scholar] [CrossRef]
  4. Schiff, N.D.; Fins, J.J. Brain death and disorders of consciousness. Curr. Biol. 2016, 26, R572–R576. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Van Erp, W.; Larvrijsen, J.; Vos, P.; Bor, H.; Laureys, S.; Koopmans, R. The vegetative state: Prevalence, misdiagnosis and treatment limitations. JAMDA 2015, 85, 9–14. [Google Scholar] [CrossRef]
  6. Wannez, S.; Heine, L.; Thonnard, M.; Gosseries, O.; Laureys, S. The repetition of behavioral assessments in diagnosis of disorders of consciousness. Ann. Neurol. 2017, 81, 883–889. [Google Scholar] [CrossRef] [Green Version]
  7. Di Perri, C.; Thibaut, A.; Heine, L.; Soddu, A.; Demertzi, A.; Laureys, S. Measuring consciousness in coma and related states. World J. Radiol. 2014, 6, 589–597. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Bodart, O.; Gosseries, O.; Wannez, S.; Thibaut, A.; Annen, J.; Boly, M.; Rosanova, M.; Casali, A.G.; Casarotto, S.; Tononi, G.; et al. Measures of metabolism and complexity in the brain of patients with disorders of consciousness. NeuroImage Clin. 2017, 14, 354–362. [Google Scholar] [CrossRef] [PubMed]
  9. Stender, J.; Gosseries, O.; Bruno, M.-A.; Charland-Verville, V.; Vanhaudenhuyse, A.; Demertzi, A.; Chatelle, C.; Thonnard, M.; Thibaut, A.; Heine, L.; et al. Diagnostic precision of PET imaging and functional MRI in disorders of consciousness: A clinical validation study. Lancet 2014, 384, 514–522. [Google Scholar] [CrossRef]
  10. Jellinger, K.A. Neuropathology of prolonged unresponsive wakefulness syndrome after blunt head injury: Review of 100 post-mortem cases. Brain Inj. 2013, 27, 917–923. [Google Scholar] [CrossRef]
  11. Kinney, H.C.; Korein, J.; Panigraphy, A.; Dikkes, P.; Goode, R. Neuropathological findings in the brain of Karen Ann Quinlan. The role of thalamus in the persistent vegetative state. N. Engl. J. Med. 1994, 330, 1469–1475. [Google Scholar] [CrossRef] [PubMed]
  12. Boly, M.; Tshibanda, L.; Vanhaudenhuyse, A.; Noirhomme, Q.; Schnakers, C.; Ledoux, D.; Boveroux, P.; Garweg, C.; Lambermont, B.; Phillips, C. Functional connectivity in the default network during resting state is preserved in a vegetative but not in a brain dead patient. Hum. Brain Mapp. 2009, 30, 2393–2400. [Google Scholar] [CrossRef] [Green Version]
  13. Nakayama, N.; Okumura, A.; Shinoda, J.; Nakashima, T.; Iwama, T. Relationship between regional cerebral metabolism and consciousness disturbance in traumatic diffuse brain injury without large focal lesions: An FDG-PET study with statistical parametric mapping analysis. J. Neurol. Neurosurg. Psychiatry 2006, 77, 856–862. [Google Scholar] [CrossRef] [PubMed]
  14. Silva, S.; Alacoque, X.; Fourcade, O.; Samii, K.; Marque, P.; Woods, R.; Mazziotta, J.; Chollet, F.; Loubinoux, I. Wakefulness and loss of awareness: Brain and brainstem interaction in the vegetative state. Neurology 2010, 74, 313–320. [Google Scholar] [CrossRef]
  15. Di Perri, C.; Bahri, M.A.; Amico, E.; Thibaut, A.; Heine, L.; Antonopoulos, G.; Charland-Verville, V.; Wannez, S.; Gomez, F.; Hustinx, R.; et al. Neural correlates of consciousness in patients who have emerged from a minimally conscious state: A cross-sectional multimodal imaging study. Lancet Neurol. 2016, 15, 830–842. [Google Scholar] [CrossRef]
  16. Laureys, S.; Tononi, G. The Neurology of Consciousness; Academic Press-Elsevier: London, UK, 2008. [Google Scholar]
  17. Mashour, G.A.; Hudetz, A.G. Neural correlates of unconsciousness in large-scale brain networks trends. Neuroscience 2018, 41, 150–160. [Google Scholar] [CrossRef]
  18. Laouchedi, M.; Galanaud, D.; Delmaire, C.; Fernandez-Vidal, S.; Messé, A.; Mesmoudi, S.; Oulebsir Boumgharb, F.; Pélégrini-Issaca, M.; Puybasseta, L.; Benali, H.; et al. Deafferentation in thalamic and pontine areas in severe traumatic brain injury. J. Neuroradiol. 2015, 42, 202–211. [Google Scholar] [CrossRef]
  19. Kampfl, A.; Schmutzhard, E.; Franz, G.; Pfausler, B.; Haring, H.P.; Ulmer, H.; Felber, S.; Golaszewski, S.; Aichner, F. Prediction of recovery from post-traumatic vegetative state with cerebral magnetic-resonance imaging. Lancet 1998, 351, 1763–1767. [Google Scholar] [CrossRef]
  20. Hoelper, B.M.; Soldner, F.; Choné, L.; Wallenfang, T. Effect of intracerebral lesions detected in early MRI on outcome after acute brain injury. Acta Neurochir. 2000, 76, 265–267. [Google Scholar]
  21. Annen, J.; Heine, L.; Ziegler, E.; Frasso, G.; Bahri, M.; Di Perri, C.; Stender, J.; Martial, C.; Wannez, S.; D’ostilio, K.; et al. Function–structure connectivity in patients with severe brain injury as measured by MRI-DWI and FDG-PET. Hum. Brain Mapp. 2016, 37, 3707–3720. [Google Scholar] [CrossRef]
  22. Fernandez-Espejo, D.; Bekinschtein, T.; Monti, M.M.; Pickard, J.D.; Junque, C.; Coleman, M.R.; Owen, A.M. Diffusion weighted imaging distinguishes the vegetative state from the minimally conscious state. NeuroImage 2011, 54, 103–112. [Google Scholar] [CrossRef]
  23. Newcombe, V.F.; Williams, G.B.; Scoffings, D.; Cross, J.; Carpenter, T.A.; Pickard, J.D.; Menon, D.K. Aetiological differences in neuroanatomy of the vegetative state: Insights from diffusion tensor imaging and functional implications. J. Neurol. Neurosurg. Psychiatry 2010, 81, 552–561. [Google Scholar] [CrossRef]
  24. Zheng, Z.S.; Reggente, N.; Lutkenhoff, E.; Owen, A.M.; Monti, M.M. Disentangling disorders of consciousness: Insights from diffusion tensor imaging and machine learning. Hum. Brain Mapp. 2017, 38, 431–443. [Google Scholar] [CrossRef]
  25. Velly, L.; Perlbarg, V.; Boulier, T.; Adam, N.; Delphine, S.; Luyt, C.-E.; Battisti, V.; Torkomian, G.; Arbelot, C.; Chabanne, R.; et al. Use of brain diffusion tensor imaging for the prediction of long-term neurological outcomes in patients after cardiac arrest: A multicentre, international, prospective, observational, cohort study. Lancet Neurol. 2018, 17, 317–326. [Google Scholar] [CrossRef]
  26. Van der Eerden, A.W.; Khalilzadeh, O.; Perlbarg, V.; Dinkel, J.; Sanchez, P.; Vos, P.E.; Luyt, C.E.; Stevens, R.D.; Menjot de Champfleur, N.; Delmaire, C.; et al. White Matter Changes in Comatose Survivors of Anoxic Ischemic Encephalopathy and Traumatic Brain Injury: Comparative Diffusion-Tensor Imaging Study. Radiology 2014, 270, 506–516. [Google Scholar] [CrossRef]
  27. Dinkel, J.; Drier, A.; Khalilzadeh, O.; Perlbarg, V.; Czernecki, V.; Gupta, R.; Gomas, F.; Sanchez, P.; Dormont, D.; Galanaud, D.; et al. Long-Term White Matter Changes after Severe Traumatic Brain Injury: A 5-Year Prospective Cohort. Am. J. Neuroradiol. 2013, 35, 23–29. [Google Scholar] [CrossRef] [Green Version]
  28. Basser, P.J.; Jones, D.K. Diffusion-tensor MRI: Theory, experimental design and data analysis—A technical review. NMR Biomed. 2002, 15, 456–467. [Google Scholar] [CrossRef]
  29. Jensen, J.H.; Helpern, J.A.; Ramani, A.; Lu, H.; Kaczynski, K. Diffusional kurtosis imaging: The quantification of non-gaussian water diffusion by means of magnetic resonance imaging. Magn. Reson. Med. 2005, 53, 1432–1440. [Google Scholar] [CrossRef]
  30. Jensen, J.H.; Helpern, J.A. MRI quantification of non-Gaussian water diffusion by kurtosis analysis. NMR Biomed. 2010, 23, 698–710. [Google Scholar] [CrossRef]
  31. Vellmer, S.; Tonoyan, A.S.; Suter, D.; Pronin, I.N.; Maximov, I.I. Validation of DWI pre-processing procedures for reliable differentiation between human brain gliomas. Z. Med. Phys. 2018, 28, 14–24. [Google Scholar] [CrossRef]
  32. De Santis, S.; Bastiani, M.; Droby, A.; Kolber, P.; Zipp, F.; Pracht, E.; Stoecker, T.; Groppa, S.; Roebroeck, A. Characterizing Microstructural Tissue Properties in Multiple Sclerosis with Diffusion MRI at 7 T and 3 T: The Impact of the Experimental Design. Neuroscience 2018, 403, 17–26. [Google Scholar] [CrossRef]
  33. Maximov, I.I.; Tonoyan, A.S.; Pronin, I.N. Differentiation of glioma malignancy grade using diffusion MRI. Phys. Med. 2017, 40, 24–32. [Google Scholar] [CrossRef] [Green Version]
  34. Grinberg, F.; Maximov, I.I.; Farrher, E.; Neuner, I.; Amort, L.; Thönneßen, H.; Oberwelland, E.; Konrad, K.; Shah, N.J. Diffusion kurtosis metrics as biomarkers of microstructural development: A comparative study of a group of children and a group of adults. NeuroImage 2017, 44, 12–22. [Google Scholar] [CrossRef]
  35. Porcari, P.; Hegi, M.E.; Lei, H.; Hamou, M.F.; Vassallo, I.; Capuani, S.; Gruetter, R.; Mlynarik, V. Early detection of human glioma sphere xenografts in mouse brain using diffusion MRI at 14.1 T: Early detection of glioma xenografts using diffusion MRI at 14.1 T. NMR Biomed. 2016, 29, 1577–1589. [Google Scholar] [CrossRef]
  36. Steven, A.J.; Zhuo, J.; Melhem, E.R. Diffusion kurtosis imaging: An emerging technique for evaluating the microstructural environment of the brain. Am. J. Roentgenol. 2014, 202, 26–33. [Google Scholar] [CrossRef]
  37. Cheung, J.S.; Wang, E.; Lo, E.H.; Sun, P.Z. Stratification of heterogeneous diffusion MRI ischemic lesion with kurtosis imaging. Stroke 2012, 43, 2252–2254. [Google Scholar] [CrossRef]
  38. Surova, Y.; Nilsson, M.; Lampinen, B.; Lätt, J.; Hall, S.; Widner, H.; van Westen, D.; Hansson, O. Alteration of putaminal fractional anisotropy in Parkinson’s disease: A longitudinal diffusion kurtosis imaging study. Neuroradiology 2018, 60, 247–254. [Google Scholar] [CrossRef]
  39. Kamagata, K.; Zalesky, A.; Hatano, T.; Ueda, R.; Di Biase, M.A.; Okuzumi, A.; Shimoji, K.; Hori, M.; Caeyenberghs, K.; Pantelis, C.; et al. Gray Matter Abnormalities in Idiopathic Parkinson’s Disease: Evaluation by Diffusional Kurtosis Imaging and Neurite Orientation Dispersion and Density Imaging. Hum. Brain Mapp. 2017, 38, 3704–3722. [Google Scholar] [CrossRef]
  40. Van Cauter, S.; Veraart, J.; Sijbers, J.; Peeters, R.R.; Himmelreich, U.; De Keyzer, F.; van Gool, S.W.; van Calenbergh, F.; de Vleeschouwer, S.; van Hecke, W.; et al. Gliomas: Diffusion kurtosis MR imaging in grading. Radiology 2012, 263, 492–501. [Google Scholar] [CrossRef]
  41. Falangola, M.F.; Guilfoyle, D.N.; Tabesh, A.; Hui, E.S.; Nie, X.; Jensen, J.H.; Gerum, S.V.; Hu, C.; La Francois, J.; Collins, H.R.; et al. Histological Correlation of Diffusional Kurtosis and White Matter Modeling Metrics in the Cuprizone-Induced Corpus Callosum Demyelination. NMR Biomed. 2014, 27, 948–957. [Google Scholar] [CrossRef]
  42. Seel, R.T.; Sherer, M.; Whyte, J.; Katz, D.I.; Giacino, J.T.; Rosenbaum, A.M.; Hammond, F.M.; Kalmar, K.; Pape, T.L.; Zafonte, R.; et al. Assessment scales for disorders of consciousness: Evidence-based recommendations for clinical practice and research. Arch. Phys. Med. Rehabil. 2010, 91, 1795–1813. [Google Scholar] [CrossRef]
  43. Aja-Fernández, S.; Tristán-Vega, A.; Hoge, S. Statistical noise analysis in GRAPPA using a parametrized non-central chi approximation model. Magn. Reson. Med. 2011, 65, 1195–1206. [Google Scholar] [CrossRef]
  44. Andersson, J.; Sotiropoulos, S. An integrated approach to correction for offresonance effects and subject movement in diffusion MR imaging. Neuroimage 2016, 125, 1063–1078. [Google Scholar] [CrossRef] [PubMed]
  45. Leemans, A.; Jones, D. The b-matrix must be rotated when correcting for subject motion in DTI data. Magn. Reson. Med. 2009, 62, 1336–1349. [Google Scholar] [CrossRef] [PubMed]
  46. Veraart, J.; Fieremans, E.; Jelescu, I.; Knoll, F.; Novikov, D.S. Gibbs ringing in diffusion MRI. Magn. Reson. Med. 2016, 76, 301–314. [Google Scholar] [CrossRef]
  47. Maximov, I.I.; Thönneßen, H.; Konrad, K.; Amort, L.; Neuner, I.; Shah, N.J. Statistical instability of TBSS analysis based on DTI fitting algorithm. J. Neuroimaging 2015, 25, 883–891. [Google Scholar] [CrossRef]
  48. Leemans, A.; Jeurissen, B.; Sijbers, J.; Jones, D. ExploreDTI: A graphical toolbox for processing, analyzing, and visualizing diffusion MR data. Proc. Intl. Soc. Magn. Reson. Med. 2009, 17, 3537. [Google Scholar]
  49. Poot, D.H.; den Dekker, A.J.; Achten, E.; Verhoye, M.; Sijbers, J. Optimal experimental design for diffusion kurtosis imaging. IEEE Trans. Med. Imaging 2010, 29, 819–829. [Google Scholar] [CrossRef]
  50. Lanzafame, S.; Giannelli, M.; Garaci, F.; Floris, R.; Duggento, A.; Guerrisi, M.; Toschi, N. Differences in Gaussian diffusion tensor imaging and non-Gaussian diffusion kurtosis imaging model-based estimates of diffusion tensor invariants in the human brain. Med. Phys. 2016, 43, 2464. [Google Scholar] [CrossRef]
  51. Zhang, Y.; Brady, M.; Smith, S. Segmentation of brain MR images through a hidden Markov random field model and the expectation-maximization algorithm. IEEE Trans. Med. Imaging 2001, 20, 45–57. [Google Scholar] [CrossRef]
  52. Yushkevich, P.A.; Piven, J.; Hazlett, H.C.; Smith, R.G.; Ho, S.; Gee, J.C.; Gerig, G. User-guided 3D active contour segmentation of anatomical structures: Significantly improved efficiency and reliability. Neuroimage 2006, 31, 1116–1128. [Google Scholar] [CrossRef] [Green Version]
  53. Smith, S.M.; Jenkinson, M.; Johansen-Berg, H.; Rueckert, D.; Nichols, T.E.; Mackay, C.E.; Watkins, K.E.; Ciccarelli, O.; Cader, M.Z.; Matthews, P.M.; et al. Tract-based spatial statistics: Voxelwise analysis of multi-subject diffusion data. Neuroimage 2006, 31, 1487–1505. [Google Scholar] [CrossRef]
  54. Bodart, O.; Amico, E.; Gómez, F.; Casali, A.G.; Wannez, S.; Heine, L.; Thibaut, A.; Annen, J.; Boly, M.; Casarotto, S.; et al. Global structural integrity and effective connectivity in patients with disorders of consciousness. Brain Stimul. 2017, 11, 358–365. [Google Scholar] [CrossRef]
  55. Kinney, H.C.; Samuels, M.A. Neuropathology of the persistent vegetative state. A review. J. Neuropath. Exp. Neurol. 1994, 53, 548–558. [Google Scholar] [CrossRef]
  56. Helpern, J.A.; Adisetiyo, V.; Falangola, M.F.; Hu, C.; Di Martino, A.; Williams, K.; Castellanos, F.X.; Jensen, J.H. Preliminary evidence of altered gray and white matter microstructural development in the frontal lobe of adolescents with attention-deficit hyperactivity disorder: A diffusional kurtosis imaging study. J. Magn. Reson. Imaging 2011, 33, 17–23. [Google Scholar] [CrossRef]
  57. Lu, H.; Jensen, J.H.; Ramani, A.; Helpern, J.A. Three-dimensional characterization of non-gaussian water diffusion in humans using diffusion kurtosis imaging. NMR Biomed. 2006, 19, 236–247. [Google Scholar] [CrossRef]
  58. Raab, P.; Hattingen, E.; Franz, K.; Zanella, F.E.; Lanfermann, H. Cerebral gliomas: Diffusional kurtosis imaging analysis of microstructural differences. Radiology 2010, 254, 876–881. [Google Scholar] [CrossRef]
  59. Ramani, A.; Jensen, J.H.; Szulc, K.U.; Ali, O.; Hu, C.; Lu, H.; Brodle, J.D.; Helpern, J.A. Assessment of Abnormalities in the Cerebral Microstructure of Schizophrenia Patients: A Diffusional Kurtosis Imaging Study. Proc. Int. Soc. Magn. Reson. Med. 2007, 15, 648. [Google Scholar]
  60. Grossman, E.J.; Ge, Y.; Jensen, J.H.; Babb, J.S.; Miles, L.; Reaume, J.; Silver, J.M.; Grossman, R.I.; Inglese, M. Thalamus and Cognitive Impairment in Mild Traumatic Brain Injury: A Diffusional Kurtosis Imaging Study. J. Neurotrauma 2012, 29, 2318–2327. [Google Scholar] [CrossRef]
  61. Yu, F.; Shukla, D.K.; Armstrong, R.C.; Marion, C.M.; Radomski, K.L.; Selwyn, R.G.; Dardzinski, B.J. Repetitive Model of Mild Traumatic Brain Injury Produces Cortical Abnormalities Detectable by Magnetic Resonance Diffusion Imaging, Histopathology, and Behavior. J. Neurotrauma 2017, 34, 1364–1381. [Google Scholar] [CrossRef] [Green Version]
  62. Sours, C.; Raghavan, P.; Medina, A.E.; Roys, S.; Jiang, L.; Zhuo, J.; Gullapalli, R.P. Structural and Functional Integrity of the Intraparietal Sulcus in Moderate and Severe Traumatic Brain Injury. J. Neurotrauma 2017, 34, 1473–1481. [Google Scholar] [CrossRef] [Green Version]
  63. Stokum, J.A.; Sours, C.; Zhuo, J.; Kane, R.; Shanmuganathan, K.; Gullapalli, R.P. A longitudinal evaluation of diffusion kurtosis imaging in patients with mild traumatic brain injury. Brain Inj. 2015, 29, 47–57. [Google Scholar] [CrossRef]
  64. Adams, J.H.; Graham, D.I.; Jennett, B. The neuropathology of the vegetative state after an acute brain insult. Brain 2000, 123, 1327–1338. [Google Scholar] [CrossRef] [Green Version]
  65. Graham, D.I. Pathology of hypoxic brain damage in man. J. Clin. Path. 1977, 30, 170–180. [Google Scholar] [CrossRef]
  66. Maxwell, W.L.; Dhillon, K.; Harper, L.; Espin, J.; MacIntosh, T.K.; Smith, D.H.; Graham, D.I. There is differential loss of pyramidal cells from the human hippocampus with survival after blunt head injury. J. Neuropathol. Exp. Neurol. 2003, 62, 272–279. [Google Scholar] [CrossRef]
  67. Feldman, H.M.; Yeatman, J.D.; Lee, E.S.; Barde, L.H.; Gaman-Bean, S. Diffusion Tensor Imaging: A Review for Pediatric Researchers and Clinicians. J. Dev. Behav. Pediatr. 2010, 31, 346–356. [Google Scholar] [CrossRef] [Green Version]
  68. Alexander, A.L.; Hurley, S.A.; Samsonov, A.A.; Adluru, N.; Hosseinbor, A.P.; Mossahebi, P.; Tromp do, P.M.; Zakszewski, E.; Field, A.S. Characterization of Cerebral White Matter Properties Using Quantitative Magnetic Resonance Imaging Stains. Brain Connect. 2011, 1, 423–446. [Google Scholar] [CrossRef]
  69. Wheeler–Kingshott, C.A.M.; Cercignani, M. About “axial” and “radial” diffusivities. Magn. Reson. Med. 2009, 61, 1255–1260. [Google Scholar] [CrossRef]
  70. Field, A.S.; Alexander, A.L.; Wu, Y.C.; Hasan, K.M.; Witwer, B.; Badie, B. Diffusion tensor eigenvector directional color imaging patterns in the evaluation of cerebral white matter tracts altered by tumor. J. Magn. Reson. Imaging 2004, 20, 555–562. [Google Scholar] [CrossRef]
  71. Zhuo, J.; Xu, S.; Proctor, J.L.; Mullins, R.J.; Simon, J.Z.; Fiskum, G.; Gullapalli, R.P. Diffusion kurtosis as an in vivo imaging marker for reactive astrogliosis in traumatic brain injury. Neuroimage 2012, 59, 467–477. [Google Scholar] [CrossRef] [Green Version]
  72. Cavaliere, C.; Aiello, M.; Di Perri, C.; Fernandez-Espejo, D.; Owen, A.M.; Soddu, A. Diffusion tensor imaging and white matter abnormalities in patients with disorders of consciousness. Front. Hum. Neurosci. 2015, 8, 1028. [Google Scholar] [CrossRef] [PubMed]
  73. Maximov, I.I.; Alnaes, D.; Westlye, L.T. Towards an Optimised Processing Pipeline for Diffusion MRI Data: Effects of Artefact Corrections on Diffusion Metrics and Their Age Associations in UK Biobank. Available online: https://www.biorxiv.org/content/10.1101/511964v1.abstract (accessed on 26 May 2019).
  74. Bammer, R.; Holdsworth, S.J.; Veldhuis, W.B.; Skare, S.T. New methods in Diffusion Weighted and Diffusion Tensor Imaging. Magn. Reson. Imaging Clin. N. Am. 2009, 17, 175–204. [Google Scholar] [CrossRef] [PubMed]
  75. Mac Donald, C.L.; Dikranian, K.; Bayly, P.; Holtzman, D.; Brody, D. Diffusion Tensor Imaging Reliably Detects Experimental Traumatic Axonal Injury and Indicates Approximate Time of Injury. J. Neurosci. 2007, 27, 11869–11876. [Google Scholar] [CrossRef] [PubMed]
  76. Edlow, B.L.; Haynes, R.L.; Takahashi, E.; Klein, J.P.; Cummings, P.; Benner, T.; Greer, D.M.; Greenberg, S.M.; Wu, O.; Kinney, H.C.; et al. Disconnection of the Ascending Arousal System in Traumatic Coma. J. Neuropathol. Exp. Neurol. 2013, 72, 505–523. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Brainsci 09 00123 g001 550
Figure 1. MRI (axial T2-weighted image) data present different chronic disorders of consciousness (DOC) patients from the study. Upper row—trauma patients (T), bottom row—anoxic patients (A).
Figure 1. MRI (axial T2-weighted image) data present different chronic disorders of consciousness (DOC) patients from the study. Upper row—trauma patients (T), bottom row—anoxic patients (A).
Brainsci 09 00123 g001
Brainsci 09 00123 g002 550
Figure 2. Schematic representation of the diffusion post-processing procedure.
Figure 2. Schematic representation of the diffusion post-processing procedure.
Brainsci 09 00123 g002
Brainsci 09 00123 g003 550
Figure 3. Scatter plots of estimated diffusion metrics for grey/white matter between control/patient groups. Points on the scatter plot are encoded by the colour in respect to the axial slice position and control/patient groups. All diffusion metrics were aligned to T1-weighted individual images. Diffusion units for mean, axial, and radial diffusivities (MD, AD, and RD) are in (µm2/ms). The dashed line presents the linear dependence.
Figure 3. Scatter plots of estimated diffusion metrics for grey/white matter between control/patient groups. Points on the scatter plot are encoded by the colour in respect to the axial slice position and control/patient groups. All diffusion metrics were aligned to T1-weighted individual images. Diffusion units for mean, axial, and radial diffusivities (MD, AD, and RD) are in (µm2/ms). The dashed line presents the linear dependence.
Brainsci 09 00123 g003
Brainsci 09 00123 g004 550
Figure 4. Tract-based spatial statistics (TBSS) results demonstrating locations of metric changes in patients (P) comparing to healthy volunteers (C) (p < 0.05). The analysed white matter skeleton is shown in green with the colour overlays of significant differences. Yellow/red colour scheme represents the areas where the diffusion metrics of the patients were significantly lower than in the healthy control group; blue/light blue colour scheme is the areas where the diffusion metrics of the patients were significantly higher, respectively.
Figure 4. Tract-based spatial statistics (TBSS) results demonstrating locations of metric changes in patients (P) comparing to healthy volunteers (C) (p < 0.05). The analysed white matter skeleton is shown in green with the colour overlays of significant differences. Yellow/red colour scheme represents the areas where the diffusion metrics of the patients were significantly lower than in the healthy control group; blue/light blue colour scheme is the areas where the diffusion metrics of the patients were significantly higher, respectively.
Brainsci 09 00123 g004
Brainsci 09 00123 g005 550
Figure 5. Scatter plots of estimated diffusion metrics for grey/white matter between Trauma/Anox groups (a) and unresponsive wakefulness syndrome/ minimally conscious state (UWS/MCS) (b). Points on scatter plot are encoded by the colour in respect to the axial slice position and control/patient groups. All diffusion metrics were aligned to T1-weighted individual images. Diffusion units for MD, AD, and RD are in (µm2/ms). The dashed line presents the linear dependence.
Figure 5. Scatter plots of estimated diffusion metrics for grey/white matter between Trauma/Anox groups (a) and unresponsive wakefulness syndrome/ minimally conscious state (UWS/MCS) (b). Points on scatter plot are encoded by the colour in respect to the axial slice position and control/patient groups. All diffusion metrics were aligned to T1-weighted individual images. Diffusion units for MD, AD, and RD are in (µm2/ms). The dashed line presents the linear dependence.
Brainsci 09 00123 g005
Brainsci 09 00123 g006 550
Figure 6. TBSS results demonstrating locations of diffusion metric changes in vegetative state/unresponsive wakefulness syndrome (VS/UWS) versus MCS patients. The analysed white matter skeleton is shown in green with the colour overlay of significant differences. Yellow/red colour scheme represents the areas where the diffusion metrics of MCS patients were significantly higher than VS/UWS. (a) the significance level, p < 0.05, (b) significance level, p < 0.001.
Figure 6. TBSS results demonstrating locations of diffusion metric changes in vegetative state/unresponsive wakefulness syndrome (VS/UWS) versus MCS patients. The analysed white matter skeleton is shown in green with the colour overlay of significant differences. Yellow/red colour scheme represents the areas where the diffusion metrics of MCS patients were significantly higher than VS/UWS. (a) the significance level, p < 0.05, (b) significance level, p < 0.001.
Brainsci 09 00123 g006
Table
Table 1. Demographic and clinical data of patients.
Table 1. Demographic and clinical data of patients.
IDSexAgeMPICRS-RAetiologyLesion SiteDiagnosis
1M201521TBMCS
2M48137AL > RVS/UWS
3F47265AL > RVS/UWS
4F31518AL > RMCS
5M2356ABVS/UWS
6M26512TBMCS
7F24227TBVS/UWS
8M231212TL > RMCS
9M47912ABMCS
10M22125TBVS/UWS
11F25127TBVS/UWS
12M6776ABVS/UWS
13M47306ABVS/UWS
15F43713AR > LMCS
15M51165ABVS/UWS
16F24734ABVS/UWS
17M331012AR > LMCS
18M4236AL > RVS/UWS
Abbreviations: MPI, months post-ictus; A, anoxic; T, traumatic; M, male; F, female; MCS, minimally conscious state, VS/UWS—vegetative state/unresponsive wakefulness syndrome. Lesion site—predominant lesion hemisphere site; B, both hemispheres, L, predominant left hemisphere, R, predominant right hemisphere.
Table
Table 2. Control (C) versus patients (P) groups comparison of mean/standard deviation diffusion scalar metrics.
Table 2. Control (C) versus patients (P) groups comparison of mean/standard deviation diffusion scalar metrics.
FAMD (µm2/ms)AD (µm2/ms)RD (µm2/ms)MKAKRK
GM C < PC < PC < PC > PC > PC > P
Control 1.5/0.101.6/0.111.4/0.100.66/0.050.63/0.030.66/0.04
Patient 2.1/0.452.3/0.462.0/0.450.58/0.250.52/0.090.59/0.23
Cohen’s d 1.972.051.91−0.4−1.6−0.39
WMC > PC < PC < PC < PC > PC > PC > P
Control0.28/0.011.0/0.081.3/0.080.9/0.090.87/0.050.76/0.040.70/0.07
Patient0.17/0.031.6/0.401.9/0.381.5/0.410.68/0.220.64/0.140.98/0.24
Cohen’s d−3.281.991.862.07−1.13−1.19−1.51
BrainstemC > PC < PC < PC < PC > PC > P
Control0.32/0.101.4/0.191.8/0.231.1/0.150.85/0.130.76/0.12
Patient0.22/0.041.9/0.252.3/0.311.7/0.260.75/0.090.68/0.06
Cohen’s d−4.663.033.063.06−1.5−2.05
Thalamus C < PC < PC < P
Control 1.2/0.201.4/0.231.0/0.20
Patient 1.7/0.332.0/0.331.5/0.33
Cohen’s d 2.312.72.17
CCC > PC < PC < PC < PC > PC > PC > P
Control0.44/0.121.5/0.172.2/0.311.1/0.180.94/0.130.62/0.041.28/0.25
Patient0.20/0.082.2/0.422.6/0.362.0/0.460.58/0.120.52/0.050.60/0.14
Cohen’s d−4.462.211.422.5−4.17−2.47−7.52
Diffusion metrics demonstrated the significant difference (p < 0.025) using Wilcoxon rank test are marked as C > (<) P. Wilcoxon one-sided test for the Patient-Control group is used; results are presented for a significant difference (p < 0.025). In Table 2 the mean values over the regions are presented. White matter (WM), grey matter (GM), fractional anisotropy (FA), mean diffusivity (MD), axial diffusivity (AD), radial diffusivity (RD), mean kurtosis (MK), axial kurtosis (AK), and radial kurtosis (RK).

Share and Cite

MDPI and ACS Style

Kremneva, E.I.; Legostaeva, L.A.; Morozova, S.N.; Sergeev, D.V.; Sinitsyn, D.O.; Iazeva, E.G.; Suslin, A.S.; Suponeva, N.A.; Krotenkova, M.V.; Piradov, M.A.; et al. Feasibility of Non-Gaussian Diffusion Metrics in Chronic Disorders of Consciousness. Brain Sci. 2019, 9, 123. https://doi.org/10.3390/brainsci9050123

AMA Style

Kremneva EI, Legostaeva LA, Morozova SN, Sergeev DV, Sinitsyn DO, Iazeva EG, Suslin AS, Suponeva NA, Krotenkova MV, Piradov MA, et al. Feasibility of Non-Gaussian Diffusion Metrics in Chronic Disorders of Consciousness. Brain Sciences. 2019; 9(5):123. https://doi.org/10.3390/brainsci9050123

Chicago/Turabian Style

Kremneva, Elena I., Liudmila A. Legostaeva, Sofya N. Morozova, Dmitry V. Sergeev, Dmitry O. Sinitsyn, Elizaveta G. Iazeva, Aleksandr S. Suslin, Natalia A. Suponeva, Marina V. Krotenkova, Michael A. Piradov, and et al. 2019. "Feasibility of Non-Gaussian Diffusion Metrics in Chronic Disorders of Consciousness" Brain Sciences 9, no. 5: 123. https://doi.org/10.3390/brainsci9050123

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop