REKINDLE: Robust extraction of kurtosis INDices with linear estimation

ROI-based analysis of diffusion indices in healthy subjects and subjects with deficit or non-deficit syndrome schizophrenia

We analyzed DTI data involving 22 healthy subjects (HC), 15 patients with deficit syndrome schizophrenia (DSZ), and 25 patients with non-deficit syndrome schizophrenia (NDSZ). We used a 1.5-T MRI scanner to collect diffusion-weighted images and T1 images, which were employed to correct distortions and deformations within the diffusion-weighted images. For 156 regions of interest (ROI), we calculated the average fractional anisotropy (FA), mean diffusion (MD), and radial diffusion (RD). Each ROI underwent a group-wise comparison using permutation F-test, followed by post hoc pairwise comparisons with Bonferroni correction. In general, we observed lower FA in both schizophrenia groups compared to HC (i.e., HC>(DSZ=NDSZ)), while MD and RD showed the opposite pattern. Notably, specific ROIs with reduced FA in schizophrenia patients included bilateral nucleus accumbens, left fusiform area, brain stem, anterior corpus callosum, left rostral and caudal anterior cingulate, right posterior cingulate, left thalamus, left hippocampus, left inferior temporal cortex, right superior temporal cortex, left pars triangularis and right lingual gyrus. Significantly, the right cuneus exhibited lower FA in the DSZ group compared to other groups ((HC=NDSZ)>DSZ), without affecting MD and RD. These results indicate that compromised neural integrity in the cuneus may contribute to the pathophysiological distinctions between DSZ and NDSZ.

Outliers in diffusion-weighted MRI: Exploring detection models and mitigation strategies

Diffusion-weighted MRI (dMRI) is a medical imaging method that can be used to investigate the brain microstructure and structural connections between different brain regions. The method, however, requires relatively complex data processing frameworks and analysis pipelines. Many of these approaches are vulnerable to signal dropout artefacts that can originate from subjects moving their head during the scan. To combat these artefacts and eliminate such outliers, researchers have proposed two approaches: to replace outliers or to downweight outliers during modelling and analysis. With the rising interest in dMRI for clinical research, these types of corrections are increasingly important. Therefore, we set out to investigate the differences between outlier replacement and weighting approaches to help the dMRI community to select the best tool for their data processing pipelines. We evaluated dMRI motion correction registration and single tensor model fit pipelines using Gaussian Process and Spherical Harmonic based replacement approaches and outlier downweighting using highly realistic whole-brain simulations. As a proof of concept, we applied these approaches to dMRI infant data sets that contained varying numbers of dropout artefacts. Based on our results, we concluded that the Gaussian Process based outlier replacement provided similar tensor fit results to Gaussian Process based outlier detection and downweighting. Therefore, if only the least-squares estimate of the single tensor model is of interest, our recommendation is to use outlier replacement. However, outlier downweighting can potentially provide a more accurate estimate of the model precision which could be relevant for applications such as probabilistic tractoraphy.

Altered corpus callosum structure in adolescents with cerebral palsy: connection to gait and balance

Cerebral palsy (CP) is the most common motor disorder in childhood. Recent studies in children with CP have associated weakened sensorimotor performance with impairments in the major brain white-matter (WM) structure, corpus callosum (CC). However, the relationship between CC structure and lower extremity performance, specifically gait and balance, remains unknown. This study investigated the transcallosal WM structure and lower limb motor stability performance in adolescents aged 10-18 years with spastic hemiplegic (n = 18) or diplegic (n = 13) CP and in their age-matched controls (n = 34). The modern diffusion-weighted MRI analysis included the diffusivity properties of seven CC subparts and the transcallosal lower limb sensorimotor tract of the dominant hemisphere. Children with CP had comprehensive impairments in the cross-sectional area, fractional anisotropy, and mean diffusivity of the CC and sensorimotor tract. Additionally, the extent of WM alterations varied between hemiplegic and diplegic subgroups, which was seen especially in the fractional anisotropy values along the sensorimotor tract. The diffusion properties of transcallosal WM were further associated with static stability in all groups, and with dynamic stability in healthy controls. Our novel results clarify the mechanistic role of the corpus callosum in adolescents with and without CP offering valuable insight into the complex interplay between the brain's WM organization and motor performance. A better understanding of the brain basis of weakened stability performance could, in addition, improve the specificity of clinical diagnosis and targeted rehabilitation in CP.

The effects of cytomegalovirus on brain structure following sport-related concussion

The neurotrophic herpes virus cytomegalovirus is a known cause of neuropathology in utero and in immunocompromised populations. Cytomegalovirus is reactivated by stress and inflammation, possibly explaining the emerging evidence linking it to subtle brain changes in the context of more minor disturbances of immune function. Even mild forms of traumatic brain injury, including sport-related concussion, are major physiological stressors that produce neuroinflammation. In theory, concussion could predispose to the reactivation of cytomegalovirus and amplify the effects of physical injury on brain structure. However, to our knowledge this hypothesis remains untested. This study evaluated the effect of cytomegalovirus serostatus on white and grey matter structure in a prospective study of athletes with concussion and matched contact-sport controls. Athletes who sustained concussion (n = 88) completed MRI at 1, 8, 15 and 45 days post-injury; matched uninjured athletes (n = 73) completed similar visits. Cytomegalovirus serostatus was determined by measuring serum IgG antibodies (n = 30 concussed athletes and n = 21 controls were seropositive). Inverse probability of treatment weighting was used to adjust for confounding factors between athletes with and without cytomegalovirus. White matter microstructure was assessed using diffusion kurtosis imaging metrics in regions previously shown to be sensitive to concussion. T1-weighted images were used to quantify mean cortical thickness and total surface area. Concussion-related symptoms, psychological distress, and serum concentration of C-reactive protein at 1 day post-injury were included as exploratory outcomes. Planned contrasts compared the effects of cytomegalovirus seropositivity in athletes with concussion and controls, separately. There was a significant effect of cytomegalovirus on axial and radial kurtosis in athletes with concussion but not controls. Cytomegalovirus positive athletes with concussion showed greater axial (P = 0.007, d = 0.44) and radial (P = 0.010, d = 0.41) kurtosis than cytomegalovirus negative athletes with concussion. Similarly, there was a significant association of cytomegalovirus with cortical thickness in athletes with concussion but not controls. Cytomegalovirus positive athletes with concussion had reduced mean cortical thickness of the right hemisphere (P = 0.009, d = 0.42) compared with cytomegalovirus negative athletes with concussion and showed a similar trend for the left hemisphere (P = 0.036, d = 0.33). There was no significant effect of cytomegalovirus on kurtosis fractional anisotropy, surface area, symptoms and C-reactive protein. The results raise the possibility that cytomegalovirus infection contributes to structural brain abnormalities in the aftermath of concussion perhaps via an amplification of concussion-associated neuroinflammation. More work is needed to identify the biological pathways underlying this process and to clarify the clinical relevance of this putative viral effect.

Limb-specific thalamocortical tracts are impaired differently in hemiplegic and diplegic subtypes of cerebral palsy

Thalamocortical pathways are considered crucial in the sensorimotor functioning of children with cerebral palsy (CP). However, previous research has been limited by non-specific tractography seeding and the lack of comparison between different CP subtypes. We compared limb-specific thalamocortical tracts between children with hemiplegic (HP, N = 15) or diplegic (DP, N = 10) CP and typically developed peers (N = 19). The cortical seed-points for the upper and lower extremities were selected (i) manually based on anatomical landmarks or (ii) using functional magnetic resonance imaging (fMRI) activations following proprioceptive-limb stimulation. Correlations were investigated between tract structure (mean diffusivity, MD; fractional anisotropy, FA; apparent fiber density, AFD) and sensorimotor performance (hand skill and postural stability). Compared to controls, our results revealed increased MD in both upper and lower limb thalamocortical tracts in the non-dominant hemisphere in HP and bilaterally in DP subgroup. MD was strongly lateralized in participants with hemiplegia, while AFD seemed lateralized only in controls. fMRI-based tractography results were comparable. The correlation analysis indicated an association between the white matter structure and sensorimotor performance. These findings suggest distinct impairment of functionally relevant thalamocortical pathways in HP and DP subtypes. Thus, the organization of thalamocortical white matter tracts may offer valuable guidance for targeted, life-long rehabilitation in children with CP.

Mean kurtosis-Curve (MK-Curve) correction improves the test-retest reproducibility of diffusion kurtosis imaging at 3 T

Diffusion kurtosis imaging (DKI) is applied to gain insights into the microstructural organization of brain tissues. However, the reproducibility of DKI outside brain white matter, particularly in combination with advanced estimation to remedy its noise sensitivity, remains poorly characterized. Therefore, in this study, we investigated the variability and reliability of DKI metrics while correcting implausible values with a fit method called mean kurtosis (MK)-Curve. A total of 10 volunteers (four women; age: 41.4 ± 9.6 years) were included and underwent two MRI examinations of the brain. The images were acquired on a clinical 3-T scanner and included a T1-weighted image and a diffusion sequence with multiple diffusion weightings suitable for DKI. Region of interest analysis of common kurtosis and tensor metrics derived with the MK-Curve DKI fit was performed, including intraclass correlation (ICC) and Bland-Altman (BA) plot statistics. A p value of less than 0.05 was considered statistically significant. The analyses showed good to excellent agreement of both kurtosis tensor- and diffusion tensor-derived MK-Curve-corrected metrics (ICC values: 0.77-0.98 and 0.87-0.98, respectively), with the exception of two DKI-derived metrics (axial kurtosis in the cortex: ICC = 0.68, and radial kurtosis in deep gray matter: ICC = 0.544). Non-MK-Curve-corrected kurtosis tensor-derived metrics ranged from 0.01 to 0.52 and diffusion tensor-derived metrics from 0.06 to 0.66, indicating poor to moderate reliability. No structural bias was observed in the BA plots for any of the diffusion metrics. In conclusion, MK-Curve-corrected DKI metrics of the human brain can be reliably acquired in white and gray matter at 3 T and DKI metrics have good to excellent agreement in a test-retest setting.

Structural filtering of functional data offered discriminative features for autism spectrum disorder

This study attempted to answer the question, "Can filtering the functional data through the frequency bands of the structural graph provide data with valuable features which are not valuable in unfiltered data"?. The valuable features discriminate between autism spectrum disorder (ASD) and typically control (TC) groups. The resting-state fMRI data was passed through the structural graph's low, middle, and high-frequency band (LFB, MFB, and HFB) filters to answer the posed question. The structural graph was computed using the diffusion tensor imaging data. Then, the global metrics of functional graphs and metrics of functional triadic interactions were computed for filtered and unfiltered rfMRI data. Compared to TCs, ASDs had significantly higher clustering coefficients in the MFB, higher efficiencies and strengths in the MFB and HFB, and lower small-world propensity in the HFB. These results show over-connectivity, more global integration, and decreased local specialization in ASDs compared to TCs. Triadic analysis showed that the numbers of unbalanced triads were significantly lower for ASDs in the MFB. This finding may indicate the reason for restricted and repetitive behavior in ASDs. Also, in the MFB and HFB, the numbers of balanced triads and the energies of triadic interactions were significantly higher and lower for ASDs, respectively. These findings may reflect the disruption of the optimum balance between functional integration and specialization. There was no significant difference between ASDs and TCs when using the unfiltered data. All of these results demonstrated that significant differences between ASDs and TCs existed in the MFB and HFB of the structural graph when analyzing the global metrics of the functional graph and triadic interaction metrics. Also, these results demonstrated that frequency bands of the structural graph could offer significant findings which were not found in the unfiltered data. In conclusion, the results demonstrated the promising perspective of using structural graph frequency bands for attaining discriminative features and new knowledge, especially in the case of ASD.

Topology of diffusion changes in corpus callosum in Alzheimer's disease: An exploratory case-control study

Aim

This study aims to assess the integrity of white matter in various segments of the corpus callosum in Alzheimer's disease (AD) by using metrics derived from diffusion tensor imaging (DTI), diffusion kurtosis imaging (DKI) and white matter tract integrity model (WMTI) and compare these findings to healthy controls (HC).

Methods

The study was approved by the institutional ethics board. 12 AD patients and 12 HC formed the study population. All AD patients were recruited from a tertiary neurology memory clinic. A standardized battery of neuropsychological assessments was administered to the study participants by a trained rater. MRI scans were performed with a Philips Ingenia 3.0T scanner equipped with a 32-channel head coil. The protocol included a T1-weighted sequence, FLAIR and a dMRI acquisition. The dMRI scan included a total of 71 volumes, 8 at b = 0 s/mm2, 15 at b = 1,000 s/mm2 and 48 at b = 2,000 s/mm2. Diffusion data fit was performed using DKI REKINDLE and WMTI models.

Results and discussion

We detected changes suggesting demyelination and axonal degeneration throughout the corpus callosum of patients with AD, most prominent in the mid-anterior and mid-posterior segments of CC. Axial kurtosis was the most significantly altered metric, being reduced in AD patients in almost all segments of corpus callosum. Reduced axial kurtosis in the CC segments correlated with poor cognition scores in AD patients in the visuospatial, language and attention domains.

Sensory processing sensitivity and axonal microarchitecture: identifying brain structural characteristics for behavior

Previous research using functional MRI identified brain regions associated with sensory processing sensitivity (SPS), a proposed normal phenotype trait. To further validate SPS, to characterize it anatomically, and to test the usefulness in psychology of methodologies that assess axonal properties, the present study correlated SPS proxy questionnaire scores (adjusted for neuroticism) with diffusion tensor imaging (DTI) measures. Participants (n = 408) from the Human Connectome Project were studied. Voxelwise analysis showed that mean- and radial diffusivity correlated positively with SPS scores in the right and left subcallosal and anterior-ventral cingulum bundle, and the right forceps minor of the corpus callosum, all frontal cortex areas generally underlying emotion, motivation, and cognition. Further analyses showed correlations throughout medial frontal cortical regions in the right and left ventromedial prefrontal cortex, including the superior longitudinal fasciculus, inferior fronto-occipital fasciculus, uncinate, and arcuate fasciculus. Fractional anisotropy was negatively correlated with SPS scores in white matter (WM) of the right premotor/motor/somatosensory/supramarginal gyrus regions. Region of interest (ROI) analysis showed small effect sizes (- 0.165 to 0.148) in WM of the precuneus and inferior frontal gyrus. Other ROI effects were found in the dorsal-, ventral visual pathways and primary auditory cortex. The results reveal that in a large group of participants, axonal microarchitectural differences can be identified with SPS traits that are subtle and in the range of typical behavior. The results suggest that the heightened sensory processing in people who show that SPS may be influenced by the microstructure of WM in specific cortical regions. Although previous fMRI studies had identified most of these areas, the DTI results put a new focus on brain areas related to attention and cognitive flexibility, empathy, emotion, and first levels of sensory processing, as in primary auditory cortex. Psychological trait characterization may benefit from DTI methodology by identifying influential brain systems for traits.

Topological analysis of brain dynamics in autism based on graph and persistent homology

Autism spectrum disorder (ASD) is a heterogeneous disorder with a rapidly growing prevalence. In recent years, the dynamic functional connectivity (DFC) technique has been used to reveal the transient connectivity behavior of ASDs' brains by clustering connectivity matrices in different states. However, the states of DFC have not been yet studied from a topological point of view. In this paper, this study was performed using global metrics of the graph and persistent homology (PH) and resting-state functional magnetic resonance imaging (fMRI) data. The PH has been recently developed in topological data analysis and deals with persistent structures of data. The structural connectivity (SC) and static FC (SFC) were also studied to know which one of the SC, SFC, and DFC could provide more discriminative topological features when comparing ASDs with typical controls (TCs). Significant discriminative features were only found in states of DFC. Moreover, the best classification performance was offered by persistent homology-based metrics and in two out of four states. In these two states, some networks of ASDs compared to TCs were more segregated and isolated (showing the disruption of network integration in ASDs). The results of this study demonstrated that topological analysis of DFC states could offer discriminative features which were not discriminative in SFC and SC. Also, PH metrics can provide a promising perspective for studying ASD and finding candidate biomarkers.

Advanced spinal cord MRI in multiple sclerosis: Current techniques and future directions

Spinal cord magnetic resonance imaging (MRI) has a central role in multiple sclerosis (MS) clinical practice for diagnosis and disease monitoring. Advanced MRI sequences capable of visualizing and quantifying tissue macro- and microstructure and reflecting different pathological disease processes have been used in MS research; however, the spinal cord remains under-explored, partly due to technical obstacles inherent to imaging this structure. We propose that the study of the spinal cord merits equal ambition in overcoming technical challenges, and that there is much information to be exploited to make valuable contributions to our understanding of MS. We present a narrative review on the latest progress in advanced spinal cord MRI in MS, covering in the first part structural, functional, metabolic and vascular imaging methods. We focus on recent studies of MS and those making significant technical steps, noting the challenges that remain to be addressed and what stands to be gained from such advances. Throughout we also refer to other works that presend more in-depth review on specific themes. In the second part, we present several topics that, in our view, hold particular potential. The need for better imaging of gray matter is discussed. We stress the importance of developing imaging beyond the cervical spinal cord, and explore the use of ultra-high field MRI. Finally, some recommendations are given for future research, from study design to newer developments in analysis, and the need for harmonization of sequences and methods within the field. This review is aimed at researchers and clinicians with an interest in gaining an overview of the current state of advanced MRI research in this field and what is primed to be the future of spinal cord imaging in MS research.

Improved sensitivity and precision in multicentre diffusion MRI network analysis using thresholding and harmonization

PURPOSE

To investigate if network thresholding and raw data harmonization improve consistency of diffusion MRI (dMRI)-based brain networks while also increasing precision and sensitivity to detect disease effects in multicentre datasets.

METHODS

Brain networks were reconstructed from dMRI of five samples with cerebral small vessel disease (SVD; 629 patients, 166 controls), as a clinically relevant exemplar condition for studies on network integrity. We evaluated consistency of network architecture in age-matched controls, by calculating cross-site differences in connection probability and fractional anisotropy (FA). Subsequently we evaluated precision and sensitivity to disease effects by identifying connections with low FA in sporadic SVD patients relative to controls, using more severely affected patients with a pure form of genetically defined SVD as reference.

RESULTS

In controls, thresholding and harmonization improved consistency of network architecture, minimizing cross-site differences in connection probability and FA. In patients relative to controls, thresholding improved precision to detect disrupted connections by removing false positive connections (precision, before: 0.09-0.19; after: 0.38-0.70). Before harmonization, sensitivity was low within individual sites, with few connections surviving multiple testing correction (k = 0-25 connections). Harmonization and pooling improved sensitivity (k = 38), while also achieving higher precision when combined with thresholding (0.97).

CONCLUSION

We demonstrated that network consistency, precision and sensitivity to detect disease effects in SVD are improved by thresholding and harmonization. We recommend introducing these techniques to leverage large existing multicentre datasets to better understand the impact of disease on brain networks.

Association between Motor Planning and the Frontoparietal Network in Children: An Exploratory Multimodal Study

OBJECTIVE

Evidence from adult literature shows the involvement of cortical grey matter areas of the frontoparietal lobe and the white matter bundle, the superior longitudinal fasciculus (SLF) in motor planning. This is yet to be confirmed in children.

METHOD

A multimodal study was designed to probe the neurostructural basis of childhood motor planning. Behavioural (motor planning), magnetic resonance imaging (MRI) and diffusion weighted imaging (DWI) data were acquired from 19 boys aged 8-11 years. Motor planning was assessed using the one and two colour sequences of the octagon task. The MRI data were preprocessed and analysed using FreeSurfer 6.0. Cortical thickness and cortical surface area were extracted from the caudal middle frontal gyrus (MFG), superior frontal gyrus (SFG), precentral gyrus (PcG), supramarginal gyrus (SMG), superior parietal lobe (SPL) and the inferior parietal lobe (IPL) using the Desikan-Killiany atlas. The DWI data were preprocessed and analysed using ExploreDTI 4.8.6 and the white matter tract, the SLF was reconstructed.

RESULTS

Motor planning of the two colour sequence was associated with cortical thickness of the bilateral MFG and left SFG, PcG, IPL and SPL. The right SLF was related to motor planning for the two colour sequence as well as with the left cortical thickness of the SFG.

CONCLUSION

Altogether, morphology within frontodorsal circuity, and the white matter bundles that support communication between them, may be associated with individual differences in childhood motor planning.

Multimodal tract-based MRI metrics outperform whole brain markers in determining cognitive impact of small vessel disease-related brain injury

In cerebral small vessel disease (cSVD), whole brain MRI markers of cSVD-related brain injury explain limited variance to support individualized prediction. Here, we investigate whether considering abnormalities in brain tracts by integrating multimodal metrics from diffusion MRI (dMRI) and structural MRI (sMRI), can better capture cognitive performance in cSVD patients than established approaches based on whole brain markers. We selected 102 patients (73.7 ± 10.2 years old, 59 males) with MRI-visible SVD lesions and both sMRI and dMRI. Conventional linear models using demographics and established whole brain markers were used as benchmark of predicting individual cognitive scores. Multi-modal metrics of 73 major brain tracts were derived from dMRI and sMRI, and used together with established markers as input of a feed-forward artificial neural network (ANN) to predict individual cognitive scores. A feature selection strategy was implemented to reduce the risk of overfitting. Prediction was performed with leave-one-out cross-validation and evaluated with the R² of the correlation between measured and predicted cognitive scores. Linear models predicted memory and processing speed with R² = 0.26 and R² = 0.38, respectively. With ANN, feature selection resulted in 13 tract-specific metrics and 5 whole brain markers for predicting processing speed, and 28 tract-specific metrics and 4 whole brain markers for predicting memory. Leave-one-out ANN prediction with the selected features achieved R² = 0.49 and R² = 0.40 for processing speed and memory, respectively. Our results show proof-of-concept that combining tract-specific multimodal MRI metrics can improve the prediction of cognitive performance in cSVD by leveraging tract-specific multi-modal metrics.

Does Loss of Integrity of the Cingulum Bundle Link Amyloid-β Accumulation and Neurodegeneration in Alzheimer’s Disease?

Background: Alzheimer’s disease is characterized by the accumulation of amyloid-β (Aβ) into plaques, aggregation of tau into neurofibrillary tangles, and neurodegenerative processes including atrophy. However, there is a poorly understood spatial discordance between initial Aβ deposition and local neurodegeneration. Objective: Here, we test the hypothesis that the cingulum bundle links Aβ deposition in the cingulate cortex to medial temporal lobe (MTL) atrophy. Methods: 21 participants with mild cognitive impairment (MCI) from the UMC Utrecht memory clinic (UMCU, discovery sample) and 37 participants with MCI from Alzheimer’s Disease Neuroimaging Initiative (ADNI, replication sample) with available Aβ-PET scan, T1-weighted and diffusion-weighted MRI were included. Aβ load of the cingulate cortex was measured by the standardized uptake value ratio (SUVR), white matter integrity of the cingulum bundle was assessed by mean diffusivity and atrophy of the MTL by normalized MTL volume. Relationships were tested with linear mixed models, to accommodate multiple measures for each participant. Results: We found at most a weak association between cingulate Aβ and MTL volume (added R2 <0.06), primarily for the posterior hippocampus. In neither sample, white matter integrity of the cingulum bundle was associated with cingulate Aβ or MTL volume (added R2 <0.01). Various sensitivity analyses (Aβ-positive individuals only, posterior cingulate SUVR, MTL sub region volume) provided similar results. Conclusion: These findings, consistent in two independent cohorts, do not support our hypothesis that loss of white matter integrity of the cingulum is a connecting factor between cingulate gyrus Aβ deposition and MTL atrophy.

Post-hemorrhagic ventricular dilatation affects white matter maturation in extremely preterm infants

BACKGROUND

Data on microstructural white matter integrity in preterm infants with post-hemorrhagic ventricular dilatation (PHVD) using diffusion tensor imaging (DTI) are limited. Also, to date, no study has focused on the DTI changes in extremely preterm (EP) infants with PHVD.

METHODS

A case-control study of EP infants <28 weeks' gestation with PHVD was conducted. Diffusivity and fractional anisotropy (FA) values of corticospinal tracts (CST) and corpus callosum (CC) were measured using DTI at term-equivalent age. Outcomes were assessed at 2-years-corrected age.

RESULTS

Twenty-one infants with PHVD and 21 matched-controls were assessed. FA values in the CC were lower in infants with PHVD compared with controls (mean difference, 0.05 [95% confidence interval (CI), 0.02-0.08], p < 0.001). In infants with periventricular hemorrhagic infarction, FA values in the CC were lower than in controls (mean difference, 0.05 [95% CI, 0.02-0.09], p = 0.005). The composite cognitive and motor scores were associated with the FA value of the CC (coefficient 114, p = 0.01 and coefficient 147, p = 0.004; respectively).

CONCLUSIONS

Extremely preterm infants with PHVD showed lower FA values in CC. A positive correlation was also shown between the composite cognitive and motor scores and FA value of the CC at 2-years-corrected age.

IMPACT

Extremely preterm infants with post-hemorrhagic ventricular dilatation showed lower fractional anisotropy values in their corpus callosum compared with controls reflecting the impaired microstructure of these commissural nerve fibers that are adjacent to the dilated ventricles. Impaired microstructure of the corpus callosum was shown to be associated with cognitive and motor scores at 2-years-corrected age.

Cross-site harmonization of multi-shell diffusion MRI measures based on rotational invariant spherical harmonics (RISH)

Quantification methods based on the acquisition of diffusion magnetic resonance imaging (dMRI) with multiple diffusion weightings (e.g., multi-shell) are becoming increasingly applied to study the in-vivo brain. Compared to single-shell data for diffusion tensor imaging (DTI), multi-shell data allows to apply more complex models such as diffusion kurtosis imaging (DKI), which attempts to capture both diffusion hindrance and restriction effects, or biophysical models such as NODDI, which attempt to increase specificity by separating biophysical components. Because of the strong dependence of the dMRI signal on the measurement hardware, DKI and NODDI metrics show scanner and site differences, much like other dMRI metrics. These effects limit the implementation of multi-shell approaches in multicenter studies, which are needed to collect large sample sizes for robust analyses. Recently, a post-processing technique based on rotation invariant spherical harmonics (RISH) features was introduced to mitigate cross-scanner differences in DTI metrics. Unlike statistical harmonization methods, which require repeated application to every dMRI metric of choice, RISH harmonization is applied once on the raw data, and can be followed by any analysis. RISH features harmonization has been tested on DTI features but not its generalizability to harmonize multi-shell dMRI. In this work, we investigated whether performing the RISH features harmonization of multi-shell dMRI data removes cross-site differences in DKI and NODDI metrics while retaining longitudinal effects. To this end, 46 subjects underwent a longitudinal (up to 3 time points) two-shell dMRI protocol at 3 imaging sites. DKI and NODDI metrics were derived before and after harmonization and compared both at the whole brain level and at the voxel level. Then, the harmonization effects on cross-sectional and on longitudinal group differences were evaluated. RISH features averaged for each of the 3 sites exhibited prominent between-site differences in the frontal and posterior part of the brain. Statistically significant differences in fractional anisotropy, mean diffusivity and mean kurtosis were observed both at the whole brain and voxel level between all the acquisition sites before harmonization, but not after. The RISH method also proved effective to harmonize NODDI metrics, particularly in white matter. The RISH based harmonization maintained the magnitude and variance of longitudinal changes as compared to the non-harmonized data of all considered metrics. In conclusion, the application of RISH feature based harmonization to multi-shell dMRI data can be used to remove cross-site differences in DKI metrics and NODDI analyses, while retaining inherent relations between longitudinal acquisitions.

A diffusion tensor imaging analysis of white matter microstructures in non-operated craniosynostosis patients

PURPOSE

In 7 to 15-year-old operated syndromic craniosynostosis patients, we have shown the presence of microstructural anomalies in brain white matter by using DTI. To learn more about the cause of these anomalies, the aim of the study is to determine diffusivity values in white matter tracts in non-operated syndromic craniosynostosis patients aged 0-2 years compared to healthy controls.

METHODS

DTI datasets of 51 non-operated patients with syndromic craniosynostosis with a median [IQR] age of 0.40 [0.25] years were compared with 17 control subjects with a median of 1.20 [0.85] years. Major white matter tract pathways were reconstructed with ExploreDTI from MRI brain datasets acquired on a 1.5 T MRI system. Eigenvalues of these tract data were examined, with subsequent assessment of the affected tracts. Having syndromic craniosynostosis (versus control), gender, age, frontal occipital horn ratio (FOHR), and tract volume were treated as independent variables.

RESULTS

ʎ₂ and ʎ₃ of the tracts genu of the corpus callosum and the hippocampal segment of the cingulum bundle show a ƞ² > 0.14 in the comparison of patients vs controls, which indicates a large effect on radial diffusivity. Subsequent linear regressions on radial diffusivity of these tracts show that age and FOHR are significantly associated interacting factors on radial diffusivity (p < 0.025).

CONCLUSION

Syndromic craniosynostosis shows not to be a significant factor influencing the major white matter tracts. Enlargement of the ventricles show to be a significant factor on radial diffusivity in the tracts corpus callosum genu and the hippocampal segment of the cingulate bundle.

CLINICAL TRIAL REGISTRATION

MEC-2014-461.

Diffusion Kurtosis Imaging of Neonatal Spinal Cord in Clinical Routine

Diffusion kurtosis imaging (DKI) has undisputed advantages over the more classical diffusion magnetic resonance imaging (dMRI) as witnessed by the fast-increasing number of clinical applications and software packages widely adopted in brain imaging. However, in the neonatal setting, DKI is still largely underutilized, in particular in spinal cord (SC) imaging, because of its inherently demanding technological requirements. Due to its extreme sensitivity to non-Gaussian diffusion, DKI proves particularly suitable for detecting complex, subtle, fast microstructural changes occurring in this area at this early and critical stage of development, which are not identifiable with only DTI. Given the multiplicity of congenital anomalies of the spinal canal, their crucial effect on later developmental outcome, and the close interconnection between the SC region and the brain above, managing to apply such a method to the neonatal cohort becomes of utmost importance. This study will (i) mention current methodological challenges associated with the application of advanced dMRI methods, like DKI, in early infancy, (ii) illustrate the first semi-automated pipeline built on Spinal Cord Toolbox for handling the DKI data of neonatal SC, from acquisition setting to estimation of diffusion measures, through accurate adjustment of processing algorithms customized for adult SC, and (iii) present results of its application in a pilot clinical case study. With the proposed pipeline, we preliminarily show that DKI is more sensitive than DTI-related measures to alterations caused by brain white matter injuries in the underlying cervical SC.

Impact of thresholding on the consistency and sensitivity of diffusion MRI-based brain networks in patients with cerebral small vessel disease

INTRODUCTION

Thresholding of low-weight connections of diffusion MRI-based brain networks has been proposed to remove false-positive connections. It has been previously established that this yields more reproducible scan-rescan network architecture in healthy subjects. In patients with brain disease, network measures are applied to assess inter-individual variation and changes over time. Our aim was to investigate whether thresholding also achieves improved consistency in network architecture in patients, while maintaining sensitivity to disease effects for these applications.

METHODS

We applied fixed-density and absolute thresholding on brain networks in patients with cerebral small vessel disease (SVD, n = 86; ≈24 months follow-up), as a clinically relevant exemplar condition. In parallel, we applied the same methods in healthy young subjects (n = 44; scan-rescan interval ≈4 months) as a frame of reference. Consistency of network architecture was assessed with dice similarity of edges and intraclass correlation coefficient (ICC) of edge-weights and hub-scores. Sensitivity to disease effects in patients was assessed by evaluating interindividual variation, changes over time, and differences between those with high and low white matter hyperintensity burden, using correlation analyses and mixed ANOVA.

RESULTS

Compared to unthresholded networks, both thresholding methods generated more consistent architecture over time in patients (unthresholded: dice = .70; ICC: .70-.78; thresholded: dice = .77; ICC: .73-.83). However, absolute thresholding created fragmented nodes. Similar observations were made in the reference group. Regarding sensitivity to disease effects in patients, fixed-density thresholds that were optimal in terms of consistency (densities: .10-.30) preserved interindividual variation in global efficiency and node strength as well as the sensitivity to detect effects of time and group. Absolute thresholding produced larger fluctuations of interindividual variation.

CONCLUSIONS

Our results indicate that thresholding of low-weight connections, particularly when using fixed-density thresholding, results in more consistent network architecture in patients with longer rescan intervals, while preserving sensitivity to disease effects.

The relationship between quantitative magnetic resonance imaging of the ankle plantar flexors, muscle function during walking and maximal strength in people with neuromuscular diseases

BACKGROUND

Progression of plantar flexor weakness in neuromuscular diseases is usually monitored by muscle strength measurements, although they poorly relate to muscle function during walking. Pathophysiological changes such as intramuscular adipose tissue affect dynamic muscle function independent from isometric strength. Diffusion tensor imaging and T2 imaging are quantitative MRI measures reflecting muscular pathophysiological changes, and are therefore potential biomarkers to monitor plantar flexor functioning during walking in people with neuromuscular diseases.

METHODS

In fourteen individuals with plantar flexor weakness diffusion tensor imaging and T2 scans of the plantar flexors were obtained, and the diffusion indices fractional anisotropy and mean diffusivity calculated. With a dynamometer, maximal isometric plantar flexor strength was measured. 3D gait analysis was used to assess maximal ankle moment and power during walking.

FINDINGS

Fractional anisotropy, mean diffusivity and T2 relaxation time all moderately correlated with maximal plantar flexor strength (r > 0.512). Fractional anisotropy and mean diffusivity were not related with ankle moment or power (r < 0.288). T2 relaxation time was strongly related to ankle moment (r = -0.789) and ankle power (r = -0.798), and moderately related to maximal plantar flexor strength (r < 0.600).

INTERPRETATION

In conclusion, T2 relaxation time, indicative of multiple pathophysiological changes, was strongly related to plantar flexor function during walking, while fractional anisotropy and mean diffusivity, indicative of fiber size, only related to maximal plantar flexor strength. This indicates that these measures may be suitable to monitor muscle function and gain insights into the pathophysiological changes underlying a poor plantar flexor functioning during gait in people with neuromuscular diseases.

Multimodal MRI cerebral correlates of verbal fluency switching and its impairment in women with depression

BACKGROUND

The search of biomarkers in the field of depression requires easy implementable tests that are biologically rooted. Qualitative analysis of verbal fluency tests (VFT) are good candidates, but its cerebral correlates are unknown.

METHODS

We collected qualitative semantic and phonemic VFT scores along with grey and white matter anatomical MRI of depressed (n = 26) and healthy controls (HC, n = 25) women. Qualitative VFT variables are the "clustering score" (i.e. the ability to produce words within subcategories) and the "switching score" (i.e. the ability to switch between clusters). The clustering and switching scores were automatically calculated using a data-driven approach. Brain measures were cortical thickness (CT) and fractional anisotropy (FA). We tested for associations between CT, FA and qualitative VFT variables within each group.

RESULTS

Patients had reduced switching VFT scores compared to HC. Thicker cortex was associated with better switching score in semantic VFT bilaterally in the frontal (superior, rostral middle and inferior gyri), parietal (inferior parietal lobule including the supramarginal gyri), temporal (transverse and fusiform gyri) and occipital (lingual gyri) lobes in the depressed group. Positive association between FA and the switching score in semantic VFT was retrieved in depressed patients within the corpus callosum, right inferior fronto-occipital fasciculus, right superior longitudinal fasciculus extending to the anterior thalamic radiation (all p < 0.05, corrected).

CONCLUSION

Together, these results suggest that automatic qualitative VFT scores are associated with brain anatomy and reinforce its potential use as a surrogate for depression cerebral bases.

Inhibitory Control and the Structural Parcelation of the Right Inferior Frontal Gyrus

The right inferior frontal gyrus (rIFG) has most strongly, although not exclusively, been associated with response inhibition, not least based on covariations of behavioral performance measures and local gray matter characteristics. However, the white matter microstructure of the rIFG as well as its connectivity has been less in focus, especially when it comes to the consideration of potential subdivisions within this area. The present study reconstructed the structural connections of the three main subregions of the rIFG (i.e., pars opercularis, pars triangularis, and pars orbitalis) using diffusion tensor imaging, and further assessed their associations with behavioral measures of inhibitory control. The results revealed a marked heterogeneity of the three subregions with respect to the pattern and extent of their connections, with the pars orbitalis showing the most widespread inter-regional connectivity, while the pars opercularis showed the lowest number of interconnected regions. When relating behavioral performance measures of a stop signal task to brain structure, the data indicated an association between the dorsal opercular connectivity and the go reaction time and the stopping accuracy.

Jointly estimating parametric maps of multiple diffusion models from undersampled q-space data: A comparison of three deep learning approaches

PURPOSE

While advanced diffusion techniques have been found valuable in many studies, their clinical availability has been hampered partly due to their long scan times. Moreover, each diffusion technique can only extract a few relevant microstructural features. Using multiple diffusion methods may help to better understand the brain microstructure, which requires multiple expensive model fittings. In this work, we compare deep learning (DL) approaches to jointly estimate parametric maps of multiple diffusion representations/models from highly undersampled q-space data.

METHODS

We implement three DL approaches to jointly estimate parametric maps of diffusion tensor imaging (DTI), diffusion kurtosis imaging (DKI), neurite orientation dispersion and density imaging (NODDI), and multi-compartment spherical mean technique (SMT). A per-voxel q-space deep learning (1D-qDL), a per-slice convolutional neural network (2D-CNN), and a 3D-patch-based microstructure estimation with sparse coding using a separable dictionary (MESC-SD) network are considered.

RESULTS

The accuracy of estimated diffusion maps depends on the q-space undersampling, the selected network architecture, and the region and the parameter of interest. The smallest errors are observed for the MESC-SD network architecture (less than 10 % normalized RMSE in most brain regions).

CONCLUSION

Our experiments show that DL methods are very efficient tools to simultaneously estimate several diffusion maps from undersampled q-space data. These methods can significantly reduce both the scan ( ∼ 6-fold) and processing times ( ∼ 25-fold) for estimating advanced parametric diffusion maps while achieving a reasonable accuracy.

Identification of Proprioceptive Thalamocortical Tracts in Children: Comparison of fMRI, MEG, and Manual Seeding of Probabilistic Tractography

Studying white matter connections with tractography is a promising approach to understand the development of different brain processes, such as proprioception. An emerging method is to use functional brain imaging to select the cortical seed points for tractography, which is considered to improve the functional relevance and validity of the studied connections. However, it is unknown whether different functional seeding methods affect the spatial and microstructural properties of the given white matter connection. Here, we compared functional magnetic resonance imaging, magnetoencephalography, and manual seeding of thalamocortical proprioceptive tracts for finger and ankle joints separately. We showed that all three seeding approaches resulted in robust thalamocortical tracts, even though there were significant differences in localization of the respective proprioceptive seed areas in the sensorimotor cortex, and in the microstructural properties of the obtained tracts. Our study shows that the selected functional or manual seeding approach might cause systematic biases to the studied thalamocortical tracts. This result may indicate that the obtained tracts represent different portions and features of the somatosensory system. Our findings highlight the challenges of studying proprioception in the developing brain and illustrate the need for using multimodal imaging to obtain a comprehensive view of the studied brain process.

What's New and What's Next in Diffusion MRI Preprocessing

Diffusion MRI (dMRI) provides invaluable information for the study of tissue microstructure and brain connectivity, but suffers from a range of imaging artifacts that greatly challenge the analysis of results and their interpretability if not appropriately accounted for. This review will cover dMRI artifacts and preprocessing steps, some of which have not typically been considered in existing pipelines or reviews, or have only gained attention in recent years: brain/skull extraction, B-matrix incompatibilities w.r.t the imaging data, signal drift, Gibbs ringing, noise distribution bias, denoising, between- and within-volumes motion, eddy currents, outliers, susceptibility distortions, EPI Nyquist ghosts, gradient deviations, B₁ bias fields, and spatial normalization. The focus will be on “what’s new” since the notable advances prior to and brought by the Human Connectome Project (HCP), as presented in the predecessing issue on “Mapping the Connectome” in 2013. In addition to the development of novel strategies for dMRI preprocessing, exciting progress has been made in the availability of open source tools and reproducible pipelines, databases and simulation tools for the evaluation of preprocessing steps, and automated quality control frameworks, amongst others. Finally, this review will consider practical considerations and our view on “what’s next” in dMRI preprocessing.

Incorporating outlier information into diffusion-weighted MRI modeling for robust microstructural imaging and structural brain connectivity analyses

The white matter structures of the human brain can be represented using diffusion-weighted MRI tractography. Unfortunately, tractography is prone to find false-positive streamlines causing a severe decline in its specificity and limiting its feasibility in accurate structural brain connectivity analyses. Filtering algorithms have been proposed to reduce the number of invalid streamlines but the currently available filtering algorithms are not suitable to process data that contains motion artefacts which are typical in clinical research. We augmented the Convex Optimization Modelling for Microstructure Informed Tractography (COMMIT) algorithm to adjust for these signals drop-out motion artefacts. We demonstrate with comprehensive Monte-Carlo whole brain simulations and in vivo infant data that our robust algorithm is capable of properly filtering tractography reconstructions despite these artefacts. We evaluated the results using parametric and non-parametric statistics and our results demonstrate that if not accounted for, motion artefacts can have severe adverse effects in human brain structural connectivity analyses as well as in microstructural property mappings. In conclusion, the usage of robust filtering methods to mitigate motion related errors in tractogram filtering is highly beneficial, especially in clinical studies with uncooperative patient groups such as infants. With our presented robust augmentation and open-source implementation, robust tractogram filtering is readily available.

Along-tract analysis of the white matter is more informative about brain ageing, compared to whole-tract analysis

Diffusion Tensor Imaging (DTI) enabled the investigation of brain White Matter (WM), both qualitatively to study the macrostructure, and quantitatively to study the microstructure. The quantitative analyses are mostly performed at the whole-tract level, i.e., providing one measure of interest per tract; however, along-tract approaches may provide finer details of the quality of the WM tracts. In this study, using the DWI data collected from 40 young and 40 old individuals, we compared the DTI measures of FA, MD, AD, and RD, estimated by both whole-tract and along-tract approaches in 18 WM bundles, between the two groups. The results of the whole-tract quantitative analysis showed a statistically significant (p-FWER < 0.05) difference between the old and young groups in 6 tracts for FA, 8 tracts for MD, 1 tract for AD, and 7 tracts for RD. On the contrary, the along-tract approach showed differences between the two groups in 10 tracts for FA, 14 tracts for MD, 8 tracts for AD, and 11 tracts for RD. All the differences between the along-tract measures of the two groups had a large effect size (Cohen'd > 0.80). This study showed that the along-tract approach for the analysis of brain WM reveals changes in some WM tracts which had not shown any changes in the whole-tract approach, and therefore this finding emphasizes the utilization of the along-tract approach along with the whole-tract method for a more accurate study of the brain WM.

Non-parametric combination of multimodal MRI for lesion detection in focal epilepsy

•

Multivariate voxel-based analysis useful for lesion detection in focal epilepsy.

•

Non-parametric combination algorithm used to combine data from various MR sequences.

•

Successful lesion detection demonstrated in MRI-positive and MRI-negative patients.

•

Multimodal analysis detected abnormalities from diverse epileptogenic pathologies.

•

Sensitivity of multivariate analysis notably higher than univariate analyses.

One third of patients with medically refractory focal epilepsy have normal-appearing MRI scans. This poses a problem as identification of the epileptogenic region is required for surgical treatment. This study performs a multimodal voxel-based analysis (VBA) to identify brain abnormalities in MRI-negative focal epilepsy. Data was collected from 69 focal epilepsy patients (42 with discrete lesions on MRI scans, 27 with no visible findings on scans), and 62 healthy controls. MR images comprised T1-weighted, fluid-attenuated inversion recovery (FLAIR), fractional anisotropy (FA) and mean diffusivity (MD) from diffusion tensor imaging, and neurite density index (NDI) from neurite orientation dispersion and density imaging. These multimodal images were coregistered to T1-weighted scans, normalized to a standard space, and smoothed with 8 mm FWHM. Initial analysis performed voxel-wise one-tailed t-tests separately on grey matter concentration (GMC), FLAIR, FA, MD, and NDI, comparing patients with epilepsy to controls. A multimodal non-parametric combination (NPC) analysis was also performed simultaneously on FLAIR, FA, MD, and NDI. Resulting p-maps were family-wise error rate corrected, threshold-free cluster enhanced, and thresholded at p < 0.05. Sensitivity was established through visual comparison of results to manually drawn lesion masks or seizure onset zone (SOZ) from stereoelectroencephalography. A leave-one-out cross-validation with the same analysis protocols was performed on controls to determine specificity. NDI was the best performing individual modality, detecting focal abnormalities in 38% of patients with normal MRI and conclusive SOZ. GMC demonstrated the lowest sensitivity at 19%. NPC provided superior performance to univariate analyses with 50% sensitivity. Specificity in controls ranged between 96 and 100% for all analyses. This study demonstrated the utility of a multimodal VBA utilizing NPC for detecting epileptogenic lesions in MRI-negative focal epilepsy. Future work will apply this approach to datasets from other centres and will experiment with different combinations of MR sequences.

Toward more robust and reproducible diffusion kurtosis imaging

PURPOSE

The general utility of diffusion kurtosis imaging (DKI) is challenged by its poor robustness to imaging artifacts and thermal noise that often lead to implausible kurtosis values.

THEORY AND METHODS

A robust scalar kurtosis index can be estimated from powder-averaged diffusion-weighted data. We introduce a novel DKI estimator that uses this scalar kurtosis index as a proxy for the mean kurtosis to regularize the fit.

RESULTS

The regularized DKI estimator improves the robustness and reproducibility of the kurtosis metrics and results in parameter maps with enhanced quality and contrast.

CONCLUSION

Our novel DKI estimator promotes the wider use of DKI in clinical research and potentially diagnostics by improving the reproducibility and precision of DKI fitting and, as such, enabling enhanced visual, quantitative, and statistical analyses of DKI parameters.

High Inter-Rater Reliability of Manual Segmentation and Volume-Based Tractography in Healthy and Dystrophic Human Calf Muscle

BACKGROUND

Muscle diffusion tensor imaging (mDTI) is a promising surrogate biomarker in the evaluation of muscular injuries and neuromuscular diseases. Since mDTI metrics are known to vary between different muscles, separation of different muscles is essential to achieve muscle-specific diffusion parameters. The commonly used technique to assess DTI metrics is parameter maps based on manual segmentation (MSB). Other techniques comprise tract-based approaches, which can be performed in a previously defined volume. This so-called volume-based tractography (VBT) may offer a more robust assessment of diffusion metrics and additional information about muscle architecture through tract properties. The purpose of this study was to assess DTI metrics of human calf muscles calculated with two segmentation techniques-MSB and VBT-regarding their inter-rater reliability in healthy and dystrophic calf muscles.

METHODS

20 healthy controls and 18 individuals with different neuromuscular diseases underwent an MRI examination in a 3T scanner using a 16-channel Torso XL coil. DTI metrics were assessed in seven calf muscles using MSB and VBT. Coefficients of variation (CV) were calculated for both techniques. MSB and VBT were performed by two independent raters to assess inter-rater reliability by ICC analysis and Bland-Altman plots. Next to analysis of DTI metrics, the same assessments were also performed for tract properties extracted with VBT.

RESULTS

For both techniques, low CV were found for healthy controls (≤13%) and neuromuscular diseases (≤17%). Significant differences between methods were found for all diffusion metrics except for λ₁. High inter-rater reliability was found for both MSB and VBT (ICC ≥ 0.972). Assessment of tract properties revealed high inter-rater reliability (ICC ≥ 0.974).

CONCLUSIONS

Both segmentation techniques can be used in the evaluation of DTI metrics in healthy controls and different NMD with low rater dependency and high precision but differ significantly from each other. Our findings underline that the same segmentation protocol must be used to ensure comparability of mDTI data.

Diffusional Kurtosis Imaging in the Diffusion Imaging in Python Project

Diffusion-weighted magnetic resonance imaging (dMRI) measurements and models provide information about brain connectivity and are sensitive to the physical properties of tissue microstructure. Diffusional Kurtosis Imaging (DKI) quantifies the degree of non-Gaussian diffusion in biological tissue from dMRI. These estimates are of interest because they were shown to be more sensitive to microstructural alterations in health and diseases than measures based on the total anisotropy of diffusion which are highly confounded by tissue dispersion and fiber crossings. In this work, we implemented DKI in the Diffusion in Python (DIPY) project-a large collaborative open-source project which aims to provide well-tested, well-documented and comprehensive implementation of different dMRI techniques. We demonstrate the functionality of our methods in numerical simulations with known ground truth parameters and in openly available datasets. A particular strength of our DKI implementations is that it pursues several extensions of the model that connect it explicitly with microstructural models and the reconstruction of 3D white matter fiber bundles (tractography). For instance, our implementations include DKI-based microstructural models that allow the estimation of biophysical parameters, such as axonal water fraction. Moreover, we illustrate how DKI provides more general characterization of non-Gaussian diffusion compatible with complex white matter fiber architectures and gray matter, and we include a novel mean kurtosis index that is invariant to the confounding effects due to tissue dispersion. In summary, DKI in DIPY provides a well-tested, well-documented and comprehensive reference implementation for DKI. It provides a platform for wider use of DKI in research on brain disorders and in cognitive neuroscience.

Q-space trajectory imaging with positivity constraints (QTI+)

Q-space trajectory imaging (QTI) enables the estimation of useful scalar measures indicative of the local tissue structure. This is accomplished by employing generalized gradient waveforms for diffusion sensitization alongside a diffusion tensor distribution (DTD) model. The first two moments of the underlying DTD are made available by acquisitions at low diffusion sensitivity (b-values). Here, we show that three independent conditions have to be fulfilled by the mean and covariance tensors associated with distributions of symmetric positive semidefinite tensors. We introduce an estimation framework utilizing semi-definite programming (SDP) to guarantee that these conditions are met. Applying the framework on simulated signal profiles for diffusion tensors distributed according to non-central Wishart distributions demonstrates the improved noise resilience of QTI+ over the commonly employed estimation methods. Our findings on a human brain data set also reveal pronounced improvements, especially so for acquisition protocols featuring few number of volumes. Our method’s robustness to noise is expected to not only improve the accuracy of the estimates, but also enable a meaningful interpretation of contrast in the derived scalar maps. The technique’s performance on shorter acquisitions could make it feasible in routine clinical practice.

Improved reproducibility of diffusion MRI of the human brain with a four-way blip-up and down phase-encoding acquisition approach

PURPOSE

To assess the effects of blip-up and -down echo planar imaging (EPI) acquisition designs, with different choices of phase-encoding directions (PEDs) on the reproducibility of diffusion MRI (dMRI)-derived metrics in the human brain.

METHODS

Diffusion MRI data in seven subjects were acquired five times, each with five different protocols. The base design included 64 diffusion directions acquired with anterior-posterior (AP) PED, the first and second protocols added reverse phase-encoded b = 0 s / mm 2 posterior-anterior (PA) PED images. The third one included 32 directions all with PED acquisitions with opposite polarity (AP and PA). The fourth protocol, also with 32 unique directions used four PEDs (AP, PA, right-left (RL), and left-right (LR)). The scan time was virtually identical for all protocols. The variability of diffusion MRI metrics for each subject and each protocol was computed across the different sessions.

RESULTS

The highest reproducibility for all dMRI metrics was obtained with protocol four (AP/PA-RL/LR, ie, four-way PED). Protocols that used only b = 0 s / mm 2 for distortion correction, which are the most widely used designs, had the lowest reproducibility.

CONCLUSIONS

An acquisition design with four PEDs, including all DWIs in addition to b = 0 s / mm 2 images should be used to achieve high reproducibility in diffusion MRI studies.

Multimodal brain imaging connectivity analyses of emotional and motivational deficits in depression among women

BACKGROUND

Major depressive disorder (MDD) is characterized by impaired cortical-subcortical functional connectivity. Apathy adds to functional impairment, but its cerebral basis in MDD remains unknown. Our objective was to describe impairments in functional connectivity during emotional processing in MDD (with varying levels of congruency and attention), and to determine their correlation with apathy.

METHODS

We used the Variable Attention Affective Task during functional MRI, followed by diffusion-weighted MRI, to assess 55 right-handed women (30 with MDD and 25 healthy controls) between September 2012 and February 2015. We estimated functional connectivity using generalized psychophysiologic interaction and anatomic connectivity with tract-based spatial statistics. We measured apathy using the Apathy Evaluation Scale.

RESULTS

We found decreased functional connectivity between the left amygdala and the left anterior cingulate cortex (ACC) during negative stimuli in participants with MDD (t54 = 4.2; p = 0.035, family-wise error [FWE]-corrected). During high-attention stimuli, participants with MDD showed reduced functional connectivity between the right dorsolateral prefrontal cortex (dlPFC) and the right ACC (t54 = 4.06, pFWE = 0.02), but greater functional connectivity between the right dlPFC and the right amygdala (t54 = 3.35, p = 0.048). Apathy was associated with increased functional connectivity between the right dlPFC and the right ACC during high-attention stimuli (t28 = 5.2, p = 0.01) and increased fractional anisotropy in the right posterior cerebellum, the anterior and posterior cingulum and the bilateral internal capsule (all pFWE < 0.05).

LIMITATIONS

Limitations included a moderate sample size, concomitant antidepressant therapy and no directed connectivity.

CONCLUSION

We found that MDD was associated with impairments in cortical-subcortical functional connectivity during negative stimuli that might alter the recruitment of networks engaged in attention. Apathy-related features suggested networks similar to those observed in degenerative disorders, but possible different mechanisms.

The dissociative subtype of posttraumatic stress disorder is associated with subcortical white matter network alterations

Posttraumatic stress disorder (PTSD) is characterized by intrusions, avoidance, and hyperarousal while patients of the dissociative subtype (PTSD-D) experience additional dissociative symptoms. A neurobiological model proposes hyper-inhibition of limbic structures mediated by prefrontal cortices to underlie dissociation in PTSD. Here, we tested whether functional alterations in fronto-limbic circuits are underpinned by white matter network abnormalities on a network level. 23 women with PTSD-D and 19 women with classic PTSD participated. We employed deterministic diffusion tractography and graph theoretical analyses. Mean fractional anisotropy (FA) was chosen as a network weight and group differences assessed using network-based statistics. No significant white matter network alterations comprising both frontal and limbic structures in PTSD-D relative to classic PTSD were found. A subsequent whole brain exploratory analysis revealed relative FA alterations in PTSD-D in two subcortical networks, comprising connections between the left amygdala, hippocampus, and thalamus as well as links between the left ventral diencephalon, putamen, and pallidum, respectively. Dissociative symptom severity in the PTSD-D group correlated with FA values within both networks. Our findings suggest fronto-limbic inhibition in PTSD-D may present a dynamic neural process, which is not hard-wired via white matter tracts. Our exploratory results point towards altered fiber tract communication in a limbic-thalamic circuit, which may underlie (a) an initial strong emotional reaction to trauma reminders before conscious regulatory processes are enabled and (b) deficits in early sensory processing. In addition, aberrant structural connectivity in low-level motor regions may present neural correlates for dissociation as a passive threat-response.

Diffusion Tensor Imaging of Skeletal Muscle Contraction Using Oscillating Gradient Spin Echo

Diffusion tensor imaging (DTI) measures water diffusion in skeletal muscle tissue and allows for muscle assessment in a broad range of neuromuscular diseases. However, current DTI measurements, typically performed using pulsed gradient spin echo (PGSE) diffusion encoding, are limited to the assessment of non-contracted musculature, therefore providing limited insight into muscle contraction mechanisms and contraction abnormalities. In this study, we propose the use of an oscillating gradient spin echo (OGSE) diffusion encoding strategy for DTI measurements to mitigate the effect of signal voids in contracted muscle and to obtain reliable diffusivity values. Two OGSE sequences with encoding frequencies of 25 and 50 Hz were tested in the lower leg of five healthy volunteers with relaxed musculature and during active dorsiflexion and plantarflexion, and compared with a conventional PGSE approach. A significant reduction of areas of signal voids using OGSE compared with PGSE was observed in the tibialis anterior for the scans obtained in active dorsiflexion and in the soleus during active plantarflexion. The use of PGSE sequences led to unrealistically elevated axial diffusivity values in the tibialis anterior during dorsiflexion and in the soleus during plantarflexion, while the corresponding values obtained using the OGSE sequences were significantly reduced. Similar findings were seen for radial diffusivity, with significantly higher diffusivity measured in plantarflexion in the soleus muscle using the PGSE sequence. Our preliminary results indicate that DTI with OGSE diffusion encoding is feasible in human musculature and allows to quantitatively assess diffusion properties in actively contracting skeletal muscle. OGSE holds great potential to assess microstructural changes occurring in the skeletal muscle during contraction, and for non-invasive assessment of contraction abnormalities in patients with muscle diseases.

Thalamic energy dysfunction is associated with thalamo-cortical tract damage in multiple sclerosis: A diffusion spectroscopy study

BACKGROUND

Diffusion-weighted ¹H magnetic resonance spectroscopy (DW-MRS) allows to quantify creatine-phosphocreatine brain diffusivity (ADC(tCr)), whose reduction in multiple sclerosis (MS) has been proposed as a proxy of energy dysfunction.

OBJECTIVE

To investigate whether thalamic ADC(tCr) changes are associated with thalamo-cortical tract damage in MS.

METHODS

Twenty patients with MS and 13 healthy controls (HC) were enrolled in a DW-MRS and DW imaging (DWI) study. From DW-MRS, ADC(tCr) and total N-acetyl-aspartate diffusivity (ADC(tNAA)) were extracted in the thalami. Three thalamo-cortical tracts and one non-thalamic control tract were reconstructed from DWI. Fractional anisotropy (FA), mean (MD), axial (AD), and radial diffusivity (RD), reflecting microstructural integrity, were extracted for each tract. Associations between thalamic ADC(tCr) and tract metrics were assessed using linear regression models adjusting for age, sex, thalamic volume, thalamic ADC(tNAA), and tract-specific lesion load.

RESULTS

Lower thalamic ADC(tCr) was associated with higher MD and RD of thalamo-cortical projections in MS (MD: p = 0.029; RD: p = 0.017), but not in HC (MD: p = 0.625, interaction term between thalamic ADC(tCr) and group = 0.019; RD: p = 0.320, interaction term = 0.05). Thalamic ADC(tCr) was not associated with microstructural changes of the control tract.

CONCLUSION

Reduced thalamic ADC(tCr) correlates with thalamo-cortical tract damage in MS, showing that pathologic changes in thalamic energy metabolism are associated with structural degeneration of connected fibers.

Evaluation of interrater reliability of different muscle segmentation techniques in diffusion tensor imaging

INTRODUCTION

Muscle diffusion tensor imaging (mDTI) is a quantitative MRI technique that can provide information about muscular microstructure and integrity. Ultrasound and DTI studies have shown intramuscular differences, and therefore separation of different muscles for analysis is essential. The commonly used methods to assess DTI metrics in muscles are manual segmentation and tract-based analysis. Recently methods such as volume-based tractography have been applied to optimize muscle architecture estimation, but can also be used to assess DTI metrics.

PURPOSE

To evaluate diffusion metrics obtained using three different methods-volume-based tractography, manual segmentation-based analysis and tract-based analysis-with respect to their interrater reliability and their ability to detect intramuscular variance.

MATERIALS AND METHODS

30 volunteers underwent an MRI examination in a 3 T scanner using a 16-channel Torso XL coil. Diffusion-weighted images were acquired to obtain DTI metrics. These metrics were evaluated in six thigh muscles using volume-based tractography, manual segmentation and standard tractography. All three methods were performed by two independent raters to assess interrater reliability by ICC analysis and Bland-Altman plots. Ability to assess intramuscular variance was compared using an ANOVA with muscle as a between-subjects factor.

RESULTS

Interrater reliability for all methods was found to be excellent. The highest interrater reliability was found for volume-based tractography (ICC ≥ 0.967). Significant differences for the factor muscle in all examined diffusion parameters were shown in muscles using all methods (main effect p < 0.001).

CONCLUSIONS

Diffusion data can be assessed by volume tractography, standard tractography and manual segmentation with high interrater reliability. Each method produces different results for the investigated DTI parameters. Volume-based tractography was superior to conventional manual segmentation and tractography regarding interrater reliability and detection of intramuscular variance, while tract-based analysis showed the lowest coefficients of variation.

Chemotherapy-Induced Brain Effects in Small-Cell Lung Cancer Patients: A Multimodal MRI Study

The golden standard of treating Small Cell Lung Cancer (SCLC) entails application of platinum-based chemotherapy, is often accompanied by Prophylactic Cranial Irradiation (PCI), which have been linked to neurotoxic side-effects in cognitive functions. The related existing neuroimaging research mainly focuses on the effect of PCI treatment in life quality and expectancy, while little is known regarding the distinct adverse effects of chemotherapy. In this context, a multimodal MRI analysis based on structural and functional brain data is proposed in order to evaluate chemotherapy-specific effects on SCLC patients. Data from 20 patients (after chemotherapy and before PCI) and 14 healthy controls who underwent structural MRI, DTI and resting state fMRI were selected in this study. From a structural aspect, the proposed analysis included volumetry and thickness measurements on structural MRI data for assessing gray matter dissimilarities, as well as deterministic tractography and Tract-Based Spatial Statistics (TBSS) on DTI data, aiming to investigate potential white matter abnormalities. Functional data were also processed on the basis of connectivity analysis, evaluating brain network parameters to identify potential manifestation of functional inconsistencies. By comparing patients to healthy controls, the obtained results revealed statistically significant differences, with the patients' brains presenting reduced volumetry/thickness and fractional anisotropy values, accompanied by prominent differences in functional connectivity measurements. All above mentioned findings were observed in patients that underwent chemotherapy.

CMC is more than a measure of corticospinal tract integrity in acute stroke patients

•

CMC is weaker and occurs at lower frequencies in acute stroke patients.

•

Both afferent and efferent input signals contribute to CMC.

•

CMC should not be used as a direct measure of corticospinal tract integrity.

In healthy subjects, motor cortex activity and electromyographic (EMG) signals from contracting contralateral muscle show coherence in the beta (15–30 Hz) range. Corticomuscular coherence (CMC) is considered a sign of functional coupling between muscle and brain. Based on prior studies, CMC is altered in stroke, but functional significance of this finding has remained unclear. Here, we examined CMC in acute stroke patients and correlated the results with clinical outcome measures and corticospinal tract (CST) integrity estimated with diffusion tensor imaging (DTI). During isometric contraction of the extensor carpi radialis muscle, EMG and magnetoencephalographic oscillatory signals were recorded from 29 patients with paresis of the upper extremity due to ischemic stroke and 22 control subjects. CMC amplitudes and peak frequencies at 13–30 Hz were compared between the two groups. In the patients, the peak frequency in both the affected and the unaffected hemisphere was significantly (p < 0.01) lower and the strength of CMC was significantly (p < 0.05) weaker in the affected hemisphere compared to the control subjects. The strength of CMC in the patients correlated with the level of tactile sensitivity and clinical test results of hand function. In contrast, no correlation between measures of CST integrity and CMC was found. The results confirm the earlier findings that CMC is altered in acute stroke and demonstrate that CMC is bidirectional and not solely a measure of integrity of the efferent corticospinal tract.

Quantitative MRI of skeletal muscle in a cross-sectional cohort of patients with spinal muscular atrophy types 2 and 3

The aim of this study was to document upper leg involvement in spinal muscular atrophy (SMA) with quantitative MRI (qMRI) in a cross-sectional cohort of patients of varying type, disease severity and age. Thirty-one patients with SMA types 2 and 3 (aged 29.6 [7.6-73.9] years) and 20 healthy controls (aged 37.9 [17.7-71.6] years) were evaluated in a 3 T MRI with a protocol consisting of DIXON, T2 mapping and diffusion tensor imaging (DTI). qMRI measures were compared with clinical scores of motor function (Hammersmith Functional Motor Scale Expanded [HFMSE]) and muscle strength. Patients exhibited an increased fat fraction and fractional anisotropy (FA), and decreased mean diffusivity (MD) and T2 compared with controls (all P < .001). DTI parameters FA and MD manifest stronger effects than can be accounted for the effect of fatty replacement. Fat fraction, FA and MD show moderate correlation with muscle strength and motor function: FA is negatively associated with HFMSE and Medical Research Council sum score (τ = -0.56 and -0.59; both P < .001) whereas for fat fraction values are τ = -0.50 and -0.58, respectively (both P < .001). This study shows that DTI parameters correlate with muscle strength and motor function. DTI findings indirectly indicate cell atrophy and act as a measure independently of fat fraction. Combined these data suggest the potential of muscle DTI in monitoring disease progression and to study SMA pathogenesis in muscle.

Systemic multipotent adult progenitor cells protect the cerebellum after asphyxia in fetal sheep

Involvement of the cerebellum in the pathophysiology of hypoxic‐ischemic encephalopathy (HIE) in preterm infants is increasingly recognized. We aimed to assess the neuroprotective potential of intravenously administered multipotent adult progenitor cells (MAPCs) in the preterm cerebellum. Instrumented preterm ovine fetuses were subjected to transient global hypoxia‐ischemia (HI) by 25 minutes of umbilical cord occlusion at 0.7 of gestation. After reperfusion, two doses of MAPCs were administered intravenously. MAPCs are a plastic adherent bone‐marrow‐derived population of adult progenitor cells with neuroprotective potency in experimental and clinical studies. Global HI caused marked cortical injury in the cerebellum, histologically indicated by disruption of cortical strata, impeded Purkinje cell development, and decreased dendritic arborization. Furthermore, global HI induced histopathological microgliosis, hypomyelination, and disruption of white matter organization. MAPC treatment significantly prevented cortical injury and region‐specifically attenuated white matter injury in the cerebellum following global HI. Diffusion tensor imaging (DTI) detected HI‐induced injury and MAPC neuroprotection in the preterm cerebellum. This study has demonstrated in a preclinical large animal model that early systemic MAPC therapy improved structural injury of the preterm cerebellum following global HI. Microstructural improvement was detectable with DTI. These findings support the potential of MAPC therapy for the treatment of HIE and the added clinical value of DTI for the detection of cerebellar injury and the evaluation of cell‐based therapy.

Structural Network Analysis Using Diffusion MRI Tractography in Parkinson's Disease and Correlations With Motor Impairment

Functional impairment of spatially distributed brain regions in Parkinson's disease (PD) suggests changes in integrative and segregative network characteristics, for which novel analysis methods are available. To assess underlying structural network differences between PD patients and controls, we employed MRI T1 gray matter segmentation and diffusion MRI tractography to construct connectivity matrices to compare patients and controls with data originating from two different centers. In the Dutch dataset (Data-NL), 14 PD patients, and 15 healthy controls were analyzed, while 19 patients and 18 controls were included in the Canadian dataset (Data-CA). All subjects underwent T1 and diffusion-weighted MRI. Patients were assessed with Part 3 of the Unified Parkinson's Disease Rating Scale (UPDRS). T1 images were segmented using FreeSurfer, while tractography was performed using ExploreDTI. The regions of interest from the FreeSurfer segmentation were combined with the white matter streamline sets resulting from the tractography, to construct connectivity matrices. From these matrices, both global and local efficiencies were calculated, which were compared between the PD and control groups and related to the UPDRS motor scores. The connectivity matrices showed consistent patterns among the four groups, without significant differences between PD patients and control subjects, either in Data-NL or in Data-CA. In Data-NL, however, global and local efficiencies correlated negatively with UPDRS scores at both the whole-brain and the nodal levels [false discovery rate (FDR) 0.05]. At the nodal level, particularly, the posterior parietal cortex showed a negative correlation between UPDRS and local efficiency, while global efficiency correlated negatively with the UPDRS in the sensorimotor cortex. The spatial patterns of negative correlations between UPDRS and parameters for network efficiency seen in Data-NL suggest subtle structural differences in PD that were below sensitivity thresholds in Data-CA. These correlations are in line with previously described functional differences. The methodological approaches to detect such differences are discussed.

Multimodal principal component analysis to identify major features of white matter structure and links to reading

The role of white matter in reading has been established by diffusion tensor imaging (DTI), but DTI cannot identify specific microstructural features driving these relationships. Neurite orientation dispersion and density imaging (NODDI), inhomogeneous magnetization transfer (ihMT) and multicomponent driven equilibrium single-pulse observation of T1/T2 (mcDESPOT) can be used to link more specific aspects of white matter microstructure and reading due to their sensitivity to axonal packing and fiber coherence (NODDI) and myelin (ihMT and mcDESPOT). We applied principal component analysis (PCA) to combine DTI, NODDI, ihMT and mcDESPOT measures (10 in total), identify major features of white matter structure, and link these features to both reading and age. Analysis was performed for nine reading-related tracts in 46 neurotypical 6–16 year olds. We identified three principal components (PCs) which explained 79.5% of variance in our dataset. PC1 probed tissue complexity, PC2 described myelin and axonal packing, while PC3 was related to axonal diameter. Mixed effects regression models did not identify any significant relationships between principal components and reading skill. Bayes factor analysis revealed that the absence of relationships was not due to low power. Increasing PC1 in the left arcuate fasciculus with age suggest increases in tissue complexity, while increases of PC2 in the bilateral arcuate, inferior longitudinal, inferior fronto-occipital fasciculi, and splenium suggest increases in myelin and axonal packing with age. Multimodal white matter imaging and PCA provide microstructurally informative, powerful principal components which can be used by future studies of development and cognition. Our findings suggest major features of white matter undergo development during childhood and adolescence, but changes are not linked to reading during this period in our typically-developing sample.

Indirect frontocingulate structural connectivity predicts clinical response to accelerated rTMS in major depressive disorder

BACKGROUND

Repetitive transcranial magnetic stimulation (rTMS) is an established treatment for major depressive disorder (MDD), but its clinical efficacy remains rather modest. One reason for this could be that the propagation of rTMS effects via structural connections from the stimulated area to deeper brain structures (such as the cingulate cortices) is suboptimal.

METHODS

We investigated whether structural connectivity — derived from diffusion MRI data — could serve as a biomarker to predict treatment response. We hypothesized that stronger structural connections between the patient-specific stimulation position in the left dorsolateral prefrontal cortex (dlPFC) and the cingulate cortices would predict better clinical outcomes. We applied accelerated intermittent theta burst stimulation (aiTBS) to the left dlPFC in 40 patients with MDD. We correlated baseline structural connectivity, quantified using various metrics (fractional anisotropy, mean diffusivity, tract density, tract volume and number of tracts), with changes in depression severity scores after aiTBS.

RESULTS

Exploratory results (p < 0.05) showed that structural connectivity between the patient-specific stimulation site and the caudal and posterior parts of the cingulate cortex had predictive potential for clinical response to aiTBS.

LIMITATIONS

We used the diffusion tensor to perform tractography. A main limitation was that multiple fibre directions within voxels could not be resolved, which might have led to missing connections in some patients.

CONCLUSION

Stronger structural frontocingular connections may be of essence to optimally benefit from left dlPFC rTMS treatment in MDD. Even though the results are promising, further investigation with larger numbers of patients, more advanced tractography algorithms and classic daily rTMS treatment paradigms is warranted.

CLINICAL TRIAL REGISTRATION

http://clinicaltrials.gov/show/NCT01832805

Serial Diffusion Kurtosis Magnetic Resonance Imaging Study during Acute, Subacute, and Recovery Periods after Sport-Related Concussion

Sport-related concussion (SRC) is common in contact sports, but there remains a lack of reliable, unbiased biomarkers of brain injury and recovery. Although the symptoms of SRC generally resolve over a period of days to weeks, the lack of a biomarker impairs detection and return-to-play decisions. To this date, the pathophysiological recovery profile and relationships between brain changes and symptoms remained unclear. In the current study, diffusion kurtosis imaging (DKI) was used to monitor the effects of SRC on the brain and the trajectory of recovery in concussed American football players (n = 96) at <48 h, and 8, 15, and 45 days post-injury, who were compared with a matched group of uninjured players (n = 82). The concussed group reported significantly higher symptoms within 48 h after injury than controls, which resolved by the 8-day follow-up. The concussed group also demonstrated poorer performance on balance testing at <48 h and 8 days than controls. There were no significant differences between the groups in the Standardized Assessment of Concussion (SAC), a cognitive screening measure. DKI data were acquired with 3 mm isotropic resolution, and analyzed using tract-based spatial statistics (TBSS). Additionally, voxel- and region of interest-based analyses were also conducted. At <48 h, the concussed group showed significantly higher axial kurtosis than the control group. These differences increased in extent and magnitude at 8 days, then receded at 15 days, and returned to the normal levels by 45 days. Kurtosis fractional anisotropy (FA) exhibited a delayed response, with a consistent increase by days 15 and 45. The results indicate that changes detected in the acute period appear to be prolonged compared with clinical recovery, but additional brain changes not observable acutely appear to progress. Although further studies are needed to understand the pathological features of DKI changes after SRC, these findings highlight a potential disparity between clinical symptoms and pathophysiological recovery after SRC.