Analysis of diffusion tensor magnetic resonance imaging data using principal component analysis

An analysis method for diffusion tensor (DT) magnetic resonance imaging data is described, which, contrary to the standard method (multivariate fitting), does not require a specific functional model for diffusion-weighted (DW) signals. The method uses principal component analysis (PCA) under the assumption of a single fibre per pixel. PCA and the standard method were compared using simulations and human brain data. The two methods were equivalent in determining fibre orientation. PCA-derived fractional anisotropy and DT relative anisotropy had similar signal-to-noise ratio (SNR) and dependence on fibre shape. PCA-derived mean diffusivity had similar SNR to the respective DT scalar, and it depended on fibre anisotropy. Appropriate scaling of the PCA measures resulted in very good agreement between PCA and DT maps. In conclusion, the assumption of a specific functional model for DW signals is not necessary for characterization of anisotropic diffusion in a single fibre.

Analyzing diffusion tensor images with ghosting artifacts: the effects of direct and indirect normalization

The current study aims to assess the applicability of direct or indirect normalization for the analysis of fractional anisotropy (FA) maps in the context of diffusion-weighted images (DWIs) contaminated by ghosting artifacts. We found that FA maps acquired by direct normalization showed generally higher anisotropy than indirect normalization, and the disparities were aggravated by the presence of ghosting artifacts in DWIs. The voxel-wise statistical comparisons demonstrated that indirect normalization reduced the influence of artifacts and enhanced the sensitivity of detecting anisotropy differences between groups. This suggested that images contaminated with ghosting artifacts can be sensibly analyzed using indirect normalization.

The role of corticospinal tract damage in chronic motor recovery and neurorehabilitation: a pilot study

Background. With diffusion-tensor imaging (DTi) it is possible to estimate the structural characteristics of fiber bundles in vivo. This study used DTi to infer damage to the corticospinal tract (CST) and relates this parameter to (a) the level of residual motor ability at least 1 year poststroke and (b) the outcome of intensive motor rehabilitation with constraint-induced movement therapy (CIMT). Objective. To explore the role of CST damage in recovery and CIMT efficacy. Methods. Ten patients with low-functioning hemiparesis were scanned and tested at baseline, before and after CIMT. Lesion overlap with the CST was indexed as reduced anisotropy compared with a CST variability map derived from 26 controls. Residual motor ability was measured through the Wolf Motor Function Test (WMFT) and the Motor Activity Log (MAL) acquired at baseline. CIMT benefit was assessed through the pre—post treatment comparison of WMFT and MAL performance. Results. Lesion overlap with the CST correlated with residual motor ability at baseline, with greater deficits observed in patients with more extended CST damage. Infarct volume showed no systematic association with residual motor ability. CIMT led to significant improvements in motor function but outcome was not associated with the extent of CST damage or infarct volume. Conclusion. The study gives in vivo support for the proposition that structural CST damage, not infarct volume, is a major predictor for residual functional ability in the chronic state. The results provide initial evidence for positive effects of CIMT in patients with varying, including more severe, CST damage.

Cortical thickness changes in the non-lesioned hemisphere associated with non-paretic arm immobilization in modified CI therapy

Recent evidence suggests that immobilization of the upper limb for 2–3 weeks induces changes in cortical thickness as well as motor performance. In constraint induced (CI) therapy, one of the most effective interventions for hemiplegia, the non-paretic arm is constrained to enforce the use of the paretic arm in the home setting. With the present study we aimed to explore whether non-paretic arm immobilization in CI therapy induces structural changes in the non-lesioned hemisphere, and how these changes are related to treatment benefit. 31 patients with chronic hemiparesis participated in CI therapy with (N = 14) and without (N = 17) constraint. Motor ability scores were acquired before and after treatment. Diffusion tensor imaging (DTI) data was obtained prior to treatment. Cortical thickness was measured with the Freesurfer software. In both groups cortical thickness in the contralesional primary somatosensory cortex increased and motor function improved with the intervention. However the cortical thickness change was not associated with the magnitude of motor function improvement. Moreover, the treatment effect and the cortical thickness change were not significantly different between the constraint and the non-constraint groups. There was no correlation between fractional anisotropy changes in the non-lesioned hemisphere and treatment outcome. CI therapy induced cortical thickness changes in contralesional sensorimotor regions, but this effect does not appear to be driven by the immobilization of the non-paretic arm, as indicated by the absence of differences between the constraint and the non-constraint groups. Our data does not suggest that the arm immobilization used in CI therapy is associated with noticeable cortical thinning.

Corticospinal tract integrity and lesion volume play different roles in chronic hemiparesis and its improvement through motor practice

Background. Initial evidence suggests that the integrity of the ipsilesional corticospinal tract (CST) after stroke is strongly related to motor function in the chronic state but not the treatment gain induced by motor rehabilitation. Objective. We examined the association of motor status and treatment benefit by testing patients with a wide range of severity of hemiparesis of the left and right upper extremity. Method. Diffusion tensor imaging was performed in 22 patients beyond 12 months after onset of stroke with severe to moderate hemiparesis. Motor function was tested before and after 2 weeks of modified constraint-induced movement therapy. Results. CST integrity, but not lesion volume, correlated with the motor ability measures of the Wolf Motor Function Test and the Motor Activity Log. No differences were found between left and right hemiparesis. Motor performance improved significantly with the treatment regime, and did so equally for patients with left and right arm paresis. However, treatment benefit was not associated with either CST integrity or lesion volume. Conclusion. CST integrity correlated best in this small trial with chronic long-term status but not treatment-induced improvements. The CST may play a different role in the mechanisms mediating long-term outcome compared to those underlying practice-induced gains after a chronic plateau in motor function.

K-Surfer: a KNIME extension for the management and analysis of Human brain MRI FreeSurfer/FSL data

Human brain imaging techniques, such as Magnetic Resonance Imaging (MRI) or Diffusion Tensor Imaging (DTI), have been established as scientific and diagnostic tools and their adoption is growing in popularity. Statistical methods, machine learning and data mining algorithms have successfully been adopted to extract predictive and descriptive models from neuroimage data. However, the knowledge discovery process typically requires also the adoption of pre-processing, post-processing and visualisation techniques in complex data workflows. Currently, a main problem for the integrated preprocessing and mining of MRI data is the lack of comprehensive platforms able to avoid the manual invocation of preprocessing and mining tools, that yields to an error-prone and inefficient process. In this work we present K-Surfer, a novel plug-in of the Konstanz Information Miner (KNIME) workbench, that automatizes the preprocessing of brain images and leverages the mining capabilities of KNIME in an integrated way. K-Surfer supports the importing, filtering, merging and pre-processing of neuroimage data from FreeSurfer, a tool for human brain MRI feature extraction and interpretation. K-Surfer automatizes the steps for importing FreeSurfer data, reducing time costs, eliminating human errors and enabling the design of complex analytics workflow for neuroimage data by leveraging the rich functionalities available in the KNIME workbench.

Connectivity-based parcellation of the human frontal polar cortex

The frontal pole corresponds to Brodmann area (BA) 10, the largest single architectonic area in the human frontal lobe. Generally, BA10 is thought to contain two or three subregions that subserve broad functions such as multitasking, social cognition, attention, and episodic memory. However, there is a substantial debate about the functional and structural heterogeneity of this large frontal region. Previous connectivity-based parcellation studies have identified two or three subregions in the human frontal pole. Here, we used diffusion tensor imaging to assess structural connectivity of BA10 in 35 healthy subjects and delineated subregions based on this connectivity. This allowed us to determine the correspondence of structurally based subregions with the scheme previously defined functionally. Three subregions could be defined in each subject. However, these three subregions were not spatially consistent between subjects. Therefore, we accepted a solution with two subregions that encompassed the lateral and medial frontal pole. We then examined resting-state functional connectivity of the two subregions and found significant differences between their connectivities. The medial cluster was connected to nodes of the default-mode network, which is implicated in internally focused, self-related thought, and social cognition. The lateral cluster was connected to nodes of the executive control network, associated with directed attention and working memory. These findings support the concept that there are two major anatomical subregions of the frontal pole related to differences in functional connectivity.

Frontal networks in adults with autism spectrum disorder

It has been postulated that autism spectrum disorder is underpinned by an ‘atypical connectivity’ involving higher-order association brain regions. To test this hypothesis in a large cohort of adults with autism spectrum disorder we compared the white matter networks of 61 adult males with autism spectrum disorder and 61 neurotypical controls, using two complementary approaches to diffusion tensor magnetic resonance imaging. First, we applied tract-based spatial statistics, a ‘whole brain’ non-hypothesis driven method, to identify differences in white matter networks in adults with autism spectrum disorder. Following this we used a tract-specific analysis, based on tractography, to carry out a more detailed analysis of individual tracts identified by tract-based spatial statistics. Finally, within the autism spectrum disorder group, we studied the relationship between diffusion measures and autistic symptom severity. Tract-based spatial statistics revealed that autism spectrum disorder was associated with significantly reduced fractional anisotropy in regions that included frontal lobe pathways. Tractography analysis of these specific pathways showed increased mean and perpendicular diffusivity, and reduced number of streamlines in the anterior and long segments of the arcuate fasciculus, cingulum and uncinate—predominantly in the left hemisphere. Abnormalities were also evident in the anterior portions of the corpus callosum connecting left and right frontal lobes. The degree of microstructural alteration of the arcuate and uncinate fasciculi was associated with severity of symptoms in language and social reciprocity in childhood. Our results indicated that autism spectrum disorder is a developmental condition associated with abnormal connectivity of the frontal lobes. Furthermore our findings showed that male adults with autism spectrum disorder have regional differences in brain anatomy, which correlate with specific aspects of autistic symptoms. Overall these results suggest that autism spectrum disorder is a condition linked to aberrant developmental trajectories of the frontal networks that persist in adult life.

Disrupted white matter in language and motor tracts in developmental stuttering

White matter tractsc onnecting areas involved in speech and motor control were examined using diffusion-tensor imagingingin a sample of peoplewhostutter (n=29) who were heterogeneous with respect to age, sex, handedness and stuttering severity. The goals were to replicate previous findings in developmental stuttering and to extend ourknowledge by evaluating the relationship between white matter differences in people who stutter and factors such as age, sex, handedness and stuttering severity. We replicated previous findings that showed reduced integrity in white matter underlying ventral premotorcortex, cerebral peduncles and posteriorcorpus callosum in people who stutter, relative to controls. Tractography analysis additionally revealed significantly reduced white matter integrity in the arcuate fasciculus bilaterally and the left corticospinal tract and significantly reduced connectivity within theleft corticobulbar tract in people who stutter. Region-of-interest analyses revealed reduced white matter integrity in people whostutter in the three pairs ocerebellar peduncles thatcarry the afferent and efferent fibers of the cerebellum. Within thegroup of people who stutter, the higher the stuttering severity index, the lower the white matter integrity in the leftangular gyrus but the greater the white matter connectivity in theleft corticobulbartract. Also,in people who stutter, handedness and age predicted the integrity of the corticospinal tract and peduncles, respectively. Further studies are needed to determine which of these white matter differences relate to the neural basis of stuttering and which reflect experience-dependent plasticity.

Centennial changes in solar activity and the response of galactic cosmic rays

There is a growing consensus that the eleven year modulation of galactic cosmic rays (GCRs) resulting from solar activity is related to interplanetary propagating diffusive barriers (PDBs). The source of these PDBs is not well understood and numerical models describing GCR modulation simulate their effect by scaling the diffusion tensor to the interplanetary magnetic field strength (IMF). The implications of a century-scale change in solar wind speed and open solar flux, for numerical modelling of GCR modulation and the reconstruction of GCR variations over the last hundred years are discussed. The dominant role of the solar non-axisymmetric magnetic field in both forcing longitudinal solar wind speed fluctuations at solar maximum and in increasing the IMF is discussed in the context of a long-term rise in the open solar magnetic flux.

Effects of Thomson-scattering geometry on white-light imaging of an interplanetary shock: synthetic observations from forward magnetohydrodynamic modelling

Stereoscopic white-light imaging of a large portion of the inner heliosphere has been used to track interplanetary coronal mass ejections. At large elongations from the Sun, the white-light brightness depends on both the local electron density and the efficiency of the Thomson-scattering process. To quantify the effects of the Thomson-scattering geometry, we study an interplanetary shock using forward magnetohydrodynamic simulation and synthetic white-light imaging. Identifiable as an inclined streak of enhanced brightness in a time–elongation map, the travelling shock can be readily imaged by an observer located within a wide range of longitudes in the ecliptic. Different parts of the shock front contribute to the imaged brightness pattern viewed by observers at different longitudes. Moreover, even for an observer located at a fixed longitude, a different part of the shock front will contribute to the imaged brightness at any given time. The observed brightness within each imaging pixel results from a weighted integral along its corresponding ray-path. It is possible to infer the longitudinal location of the shock from the brightness pattern in an optical sky map, based on the east–west asymmetry in its brightness and degree of polarisation. Therefore, measurement of the interplanetary polarised brightness could significantly reduce the ambiguity in performing three-dimensional reconstruction of local electron density from white-light imaging.

Influence of mass diffusion in sea ice dynamical models

The mixing of floes of different thickness caused by repeated deformation of the ice cover is modeled as diffusion, and the mass balance equation for sea ice accounting for mass diffusion is developed. The effect of deformational diffusion on the ice thickness balance is shown to reach 1% of the divergence effect, which describes ridging and lead formation. This means that with the same accuracy the mass balance equation can be written in terms of mean velocity rather than mean mass-weighted velocity, which one should correctly use for a multicomponent fluid such as sea ice with components identified by floe thickness. Mixing (diffusion) of sea ice also occurs because of turbulent variations in wind and ocean drags that are unresolved in models. Estimates of the importance of turbulent mass diffusion on the dynamic redistribution of ice thickness are determined using empirical data for the turbulent diffusivity. For long-time-scale prediction (≫5 days), where unresolved atmospheric motion may have a length scale on the order of the Arctic basin and the time scale is larger than the synoptic time scale of atmospheric events, turbulent mass diffusion can exceed 10% of the divergence effect. However, for short-time-scale prediction, for example, 5 days, the unresolved scales are on the order of 100 km, and turbulent diffusion is about 0.1% of the divergence effect. Because inertial effects are small in the dynamics of the sea ice pack, diffusive momentum transfer can be disregarded.