Longitudinal motor system changes from acute to chronic spinal cord injury.

BACKGROUND AND PURPOSE: In acute spinal cord injury (SCI), magnetic resonance imaging (MRI) reveals tissue bridges and neurodegeneration for 2 years. This 5-year study aims to track initial lesion changes, subsequent neurodegeneration, and their impact on recovery. METHODS: This prospective longitudinal study enrolled acute SCI patients and healthy controls who were assessed clinically-and by MRI-regularly from 3 days postinjury up to 60 months. We employed histologically cross-validated quantitative MRI sequences sensitive to volume, myelin, and iron changes, thereby reflecting indirectly processes of neurodegeneration and neuroinflammation. General linear models tracked lesion and remote changes in volume, myelin- and iron-sensitive magnetic resonance indices over 5 years. Associations between lesion, degeneration, and recovery (using the Spinal Cord Independence Measure [SCIM] questionnaire and the International Standards for Neurological Classification of Spinal Cord Injury total motor score) were assessed. RESULTS: Patients' motor scores improved by an average of 12.86 (95% confidence interval [CI] = 6.70-19.00) points, and SCIM by 26.08 (95% CI = 17.00-35.20) points. Within 3-28 days post-SCI, lesion size decreased by more than two-thirds (3 days: 302.52 ± 185.80 mm2 , 28 days: 76.77 ± 88.62 mm2 ), revealing tissue bridges. Cervical cord and corticospinal tract volumes transiently increased in SCI patients by 5% and 3%, respectively, accompanied by cervical myelin decreases and iron increases. Over time, progressive atrophy was observed in both regions, which was linked to early lesion dynamics. Tissue bridges, reduced swelling, and myelin content decreases were predictive of long-term motor score recovery and improved SCIM score. CONCLUSIONS: Studying acute changes and their impact on longer follow-up provides insights into SCI trajectory, highlighting the importance of acute intervention while indicating the potential to influence outcomes in the later stages.

Predicting mortality in acutely hospitalised older patients: the impact of model dimensionality.

BACKGROUND: The prediction of long-term mortality following acute illness can be unreliable for older patients, inhibiting the delivery of targeted clinical interventions. The difficulty plausibly arises from the complex, multifactorial nature of the underlying biology in this population, which flexible, multimodal models based on machine learning may overcome. Here, we test this hypothesis by quantifying the comparative predictive fidelity of such models in a large consecutive sample of older patients acutely admitted to hospital and characterise their biological support. METHODS: A set of 804 admission episodes involving 616 unique patients with a mean age of 84.5 years consecutively admitted to the Acute Geriatric service at University College Hospital were identified, in whom clinical diagnoses, blood tests, cognitive status, computed tomography of the head, and mortality within 600 days after admission were available. We trained and evaluated out-of-sample an array of extreme gradient boosted trees-based predictive models of incrementally greater numbers of investigational modalities and modelled features. Both linear and non-linear associations with investigational features were quantified. RESULTS: Predictive models of mortality showed progressively increasing fidelity with greater numbers of modelled modalities and dimensions. The area under the receiver operating characteristic curve rose from 0.67 (sd = 0.078) for age and sex to 0.874 (sd = 0.046) for the most comprehensive model. Extracranial bone and soft tissue features contributed more than intracranial features towards long-term mortality prediction. The anterior cingulate and angular gyri, and serum albumin, were the greatest intracranial and biochemical model contributors respectively. CONCLUSIONS: High-dimensional, multimodal predictive models of mortality based on routine clinical data offer higher predictive fidelity than simpler models, facilitating individual level prognostication and interventional targeting. The joint contributions of both extracranial and intracranial features highlight the potential importance of optimising somatic as well as neural functions in healthy ageing. Our findings suggest a promising path towards a high-fidelity, multimodal index of frailty.

A hierarchical Bayesian model to find brain-behaviour associations in incomplete data sets.

Canonical Correlation Analysis (CCA) and its regularised versions have been widely used in the neuroimaging community to uncover multivariate associations between two data modalities (e.g., brain imaging and behaviour). However, these methods have inherent limitations: (1) statistical inferences about the associations are often not robust; (2) the associations within each data modality are not modelled; (3) missing values need to be imputed or removed. Group Factor Analysis (GFA) is a hierarchical model that addresses the first two limitations by providing Bayesian inference and modelling modality-specific associations. Here, we propose an extension of GFA that handles missing data, and highlight that GFA can be used as a predictive model. We applied GFA to synthetic and real data consisting of brain connectivity and non-imaging measures from the Human Connectome Project (HCP). In synthetic data, GFA uncovered the underlying shared and specific factors and predicted correctly the non-observed data modalities in complete and incomplete data sets. In the HCP data, we identified four relevant shared factors, capturing associations between mood, alcohol and drug use, cognition, demographics and psychopathological measures and the default mode, frontoparietal control, dorsal and ventral networks and insula, as well as two factors describing associations within brain connectivity. In addition, GFA predicted a set of non-imaging measures from brain connectivity. These findings were consistent in complete and incomplete data sets, and replicated previous findings in the literature. GFA is a promising tool that can be used to uncover associations between and within multiple data modalities in benchmark datasets (such as, HCP), and easily extended to more complex models to solve more challenging tasks.

Restoring statistical validity in group analyses of motion-corrupted MRI data.

Motion during the acquisition of magnetic resonance imaging (MRI) data degrades image quality, hindering our capacity to characterise disease in patient populations. Quality control procedures allow the exclusion of the most affected images from analysis. However, the criterion for exclusion is difficult to determine objectively and exclusion can lead to a suboptimal compromise between image quality and sample size. We provide an alternative, data-driven solution that assigns weights to each image, computed from an index of image quality using restricted maximum likelihood. We illustrate this method through the analysis of quantitative MRI data. The proposed method restores the validity of statistical tests, and performs near optimally in all brain regions, despite local effects of head motion. This method is amenable to the analysis of a broad type of MRI data and can accommodate any measure of image quality.

Correcting inter-scan motion artifacts in quantitative R1 mapping at 7T.

PURPOSE: Inter-scan motion is a substantial source of error in R1 estimation methods based on multiple volumes, for example, variable flip angle (VFA), and can be expected to increase at 7T where B1 fields are more inhomogeneous. The established correction scheme does not translate to 7T since it requires a body coil reference. Here we introduce two alternatives that outperform the established method. Since they compute relative sensitivities they do not require body coil images. THEORY: The proposed methods use coil-combined magnitude images to obtain the relative coil sensitivities. The first method efficiently computes the relative sensitivities via a simple ratio; the second by fitting a more sophisticated generative model. METHODS: R1 maps were computed using the VFA approach. Multiple datasets were acquired at 3T and 7T, with and without motion between the acquisition of the VFA volumes. R1 maps were constructed without correction, with the proposed corrections, and (at 3T) with the previously established correction scheme. The effect of the greater inhomogeneity in the transmit field at 7T was also explored by acquiring B1+ maps at each position. RESULTS: At 3T, the proposed methods outperform the baseline method. Inter-scan motion artifacts were also reduced at 7T. However, at 7T reproducibility only converged on that of the no motion condition if position-specific transmit field effects were also incorporated. CONCLUSION: The proposed methods simplify inter-scan motion correction of R1 maps and are applicable at both 3T and 7T, where a body coil is typically not available. The open-source code for all methods is made publicly available.

Ventralis intermedius nucleus anatomical variability assessment by MRI structural connectivity.

The ventralis intermedius nucleus (Vim) is centrally placed in the dentato-thalamo-cortical pathway (DTCp) and is a key surgical target in the treatment of severe medically refractory tremor. It is not visible on conventional MRI sequences; consequently, stereotactic targeting currently relies on atlas-based coordinates. This fails to capture individual anatomical variability, which may lead to poor long-term clinical efficacy. Probabilistic tractography, combined with known anatomical connectivity, enables localisation of thalamic nuclei at an individual subject level. There are, however, a number of confounds associated with this technique that may influence results. Here we focused on an established method, using probabilistic tractography to reconstruct the DTCp, to identify the connectivity-defined Vim (cd-Vim) in vivo. Using 100 healthy individuals from the Human Connectome Project, our aim was to quantify cd-Vim variability across this population, measure the discrepancy with atlas-defined Vim (ad-Vim), and assess the influence of potential methodological confounds. We found no significant effect of any of the confounds. The mean cd-Vim coordinate was located within 1.88 mm (left) and 2.12 mm (right) of the average midpoint and 3.98 mm (left) and 5.41 mm (right) from the ad-Vim coordinates. cd-Vim location was more variable on the right, which reflects hemispheric asymmetries in the probabilistic DTC reconstructed. The method was reproducible, with no significant cd-Vim location differences in a separate test-retest cohort. The superior cerebellar peduncle was identified as a potential source of artificial variance. This work demonstrates significant individual anatomical variability of the cd-Vim that atlas-based coordinate targeting fails to capture. This variability was not related to any methodological confound tested. Lateralisation of cerebellar functions, such as speech, may contribute to the observed asymmetry. Tractography-based methods seem sensitive to individual anatomical variability that is missed by conventional neurosurgical targeting; these findings may form the basis for translational tools to improve efficacy and reduce side-effects of thalamic surgery for tremor.

Uncertainty analysis of MR-PET image registration for precision neuro-PET imaging.

Accurate regional brain quantitative PET measurements, particularly when using partial volume correction, rely on robust image registration between PET and MR images. We argue here that the precision, and hence the uncertainty, of MR-PET image registration is mainly driven by the registration implementation and the quality of PET images due to their lower resolution and higher noise compared to the structural MR images. We propose a dedicated uncertainty analysis for quantifying the precision of MR-PET registration, centred around the bootstrap resampling of PET list-mode events to generate multiple PET image realisations with different noise (count) levels. The effects of PET image reconstruction parameters, such as the use of attenuation and scatter corrections and different number of iterations, on the precision and accuracy of MR-PET registration were investigated. In addition, the performance of four software packages with their default settings for rigid inter-modality image registration were considered: NiftyReg, Vinci, FSL and SPM. Four distinct PET image distributions made of two early time frames (similar to cortical FDG) and two late frames using two amyloid PET dynamic acquisitions of one amyloid positive and one amyloid negative participants were investigated. For the investigated four PET frames, the biggest impact on the uncertainty was observed between registration software packages (up to 10-fold difference in precision) followed by the reconstruction parameters. On average, the lowest uncertainty for different PET frames and brain regions was observed with SPM and two iterations of fully quantitative image reconstruction. The observed uncertainty for the varying PET count-level (from 5% to 60%) was slightly lower than for the reconstruction parameters. We also observed that the registration uncertainty in quantitative PET analysis depends on amyloid status of the considered PET frames, with increased uncertainty (up to three times) when using post-reconstruction partial volume correction. This analysis is applicable for PET data obtained from either PET/MR or PET/CT scanners.

Simultaneous voxel-wise analysis of brain and spinal cord morphometry and microstructure within the SPM framework.

To validate a simultaneous analysis tool for the brain and cervical cord embedded in the statistical parametric mapping (SPM) framework, we compared trauma-induced macro- and microstructural changes in spinal cord injury (SCI) patients to controls. The findings were compared with results obtained from existing processing tools that assess the brain and spinal cord separately. A probabilistic brain-spinal cord template (BSC) was generated using a generative semi-supervised modelling approach. The template was incorporated into the pre-processing pipeline of voxel-based morphometry and voxel-based quantification analyses in SPM. This approach was validated on T1-weighted scans and multiparameter maps, by assessing trauma-induced changes in SCI patients relative to controls and comparing the findings with the outcome from existing analytical tools. Consistency of the MRI measures was assessed using intraclass correlation coefficients (ICC). The SPM approach using the BSC template revealed trauma-induced changes across the sensorimotor system in the cord and brain in SCI patients. These changes were confirmed with established approaches covering brain or cord, separately. The ICC in the brain was high within regions of interest, such as the sensorimotor cortices, corticospinal tracts and thalamus. The simultaneous voxel-wise analysis of brain and cervical spinal cord was performed in a unique SPM-based framework incorporating pre-processing and statistical analysis in the same environment. Validation based on a SCI cohort demonstrated that the new processing approach based on the brain and cord is comparable to available processing tools, while offering the advantage of performing the analysis simultaneously across the neuraxis.

Microstructural plasticity in nociceptive pathways after spinal cord injury.

OBJECTIVE: To track the interplay between (micro-) structural changes along the trajectories of nociceptive pathways and its relation to the presence and intensity of neuropathic pain (NP) after spinal cord injury (SCI). METHODS: A quantitative neuroimaging approach employing a multiparametric mapping protocol was used, providing indirect measures of myelination (via contrasts such as magnetisation transfer (MT) saturation, longitudinal relaxation (R1)) and iron content (via effective transverse relaxation rate (R2*)) was used to track microstructural changes within nociceptive pathways. In order to characterise concurrent changes along the entire neuroaxis, a combined brain and spinal cord template embedded in the statistical parametric mapping framework was used. Multivariate source-based morphometry was performed to identify naturally grouped patterns of structural variation between individuals with and without NP after SCI. RESULTS: In individuals with NP, lower R1 and MT values are evident in the primary motor cortex and dorsolateral prefrontal cortex, while increases in R2* are evident in the cervical cord, periaqueductal grey (PAG), thalamus and anterior cingulate cortex when compared with pain-free individuals. Lower R1 values in the PAG and greater R2* values in the cervical cord are associated with NP intensity. CONCLUSIONS: The degree of microstructural changes across ascending and descending nociceptive pathways is critically implicated in the maintenance of NP. Tracking maladaptive plasticity unravels the intimate relationships between neurodegenerative and compensatory processes in NP states and may facilitate patient monitoring during therapeutic trials related to pain and neuroregeneration.

A Symmetric Prior for the Regularisation of Elastic Deformations: Improved anatomical plausibility in nonlinear image registration.

Nonlinear registration is critical to many aspects of Neuroimaging research. It facilitates averaging and comparisons across multiple subjects, as well as reporting of data in a common anatomical frame of reference. It is, however, a fundamentally ill-posed problem, with many possible solutions which minimise a given dissimilarity metric equally well. We present a regularisation method capable of selectively driving solutions towards those which would be considered anatomically plausible by penalising unlikely lineal, areal and volumetric deformations. This penalty is symmetric in the sense that geometric expansions and contractions are penalised equally, which encourages inverse-consistency. We demonstrate that this method is able to significantly reduce local volume changes and shape distortions compared to state-of-the-art elastic (FNIRT) and plastic (ANTs) registration frameworks. Crucially, this is achieved whilst simultaneously matching or exceeding the registration quality of these methods, as measured by overlap scores of labelled cortical regions. Extensive leveraging of GPU parallelisation has allowed us to solve this highly computationally intensive optimisation problem while maintaining reasonable run times of under half an hour.

The influence of microsatellite polymorphisms in sex steroid receptor genes ESR1, ESR2 and AR on sex differences in brain structure.

The androgen receptor (AR), oestrogen receptor alpha (ESR1) and oestrogen receptor beta (ESR2) play essential roles in mediating the effect of sex hormones on sex differences in the brain. Using Voxel-based morphometry (VBM) and gene sizing in two independent samples (discovery n â€‹= â€‹173, replication â€‹= â€‹61), we determine the common and unique influences on brain sex differences in grey (GM) and white matter (WM) volume between repeat lengths (n) of microsatellite polymorphisms AR(CAG)n, ESR1(TA)n and ESR2(CA)n. In the hypothalamus, temporal lobes, anterior cingulate cortex, posterior insula and prefrontal cortex, we find increased GM volume with increasing AR(CAG)n across sexes, decreasing ESR1(TA)n across sexes and decreasing ESR2(CA)n in females. Uniquely, AR(CAG)n was positively associated with dorsolateral prefrontal and orbitofrontal GM volume and the anterior corona radiata, left superior fronto-occipital fasciculus, thalamus and internal capsule WM volume. ESR1(TA)n was negatively associated with the left superior corona radiata, left cingulum and left inferior longitudinal fasciculus WM volume uniquely. ESR2(CA)n was negatively associated with right fusiform and posterior cingulate cortex uniquely. We thus describe the neuroanatomical correlates of three microsatellite polymorphisms of steroid hormone receptors and their relationship to sex differences.

Efficacy of spoken word comprehension therapy in patients with chronic aphasia: a cross-over randomised controlled trial with structural imaging.

OBJECTIVE: The efficacy of spoken language comprehension therapies for persons with aphasia remains equivocal. We investigated the efficacy of a self-led therapy app, 'Listen-In', and examined the relation between brain structure and therapy response. METHODS: A cross-over randomised repeated measures trial with five testing time points (12-week intervals), conducted at the university or participants' homes, captured baseline (T1), therapy (T2-T4) and maintenance (T5) effects. Participants with chronic poststroke aphasia and spoken language comprehension impairments completed consecutive Listen-In and standard care blocks (both 12 weeks with order randomised). Repeated measures analyses of variance compared change in spoken language comprehension on two co-primary outcomes over therapy versus standard care. Three structural MRI scans (T2-T4) for each participant (subgroup, n=25) were analysed using cross-sectional and longitudinal voxel-based morphometry. RESULTS: Thirty-five participants completed, on average, 85 hours (IQR=70-100) of Listen-In (therapy first, n=18). The first study-specific co-primary outcome (Auditory Comprehension Test (ACT)) showed large and significant improvements for trained spoken words over therapy versus standard care (11%, Cohen's d=1.12). Gains were largely maintained at 12 and 24 weeks. There were no therapy effects on the second standardised co-primary outcome (Comprehensive Aphasia Test: Spoken Words and Sentences). Change on ACT trained words was associated with volume of pretherapy right hemisphere white matter and post-therapy grey matter tissue density changes in bilateral temporal lobes. CONCLUSIONS: Individuals with chronic aphasia can improve their spoken word comprehension many years after stroke. Results contribute to hemispheric debates implicating the right hemisphere in therapy-driven language recovery. Listen-In will soon be available on GooglePlay. TRIAL REGISTRATION NUMBER: NCT02540889.

hMRI - A toolbox for quantitative MRI in neuroscience and clinical research.

Neuroscience and clinical researchers are increasingly interested in quantitative magnetic resonance imaging (qMRI) due to its sensitivity to micro-structural properties of brain tissue such as axon, myelin, iron and water concentration. We introduce the hMRI-toolbox, an open-source, easy-to-use tool available on GitHub, for qMRI data handling and processing, presented together with a tutorial and example dataset. This toolbox allows the estimation of high-quality multi-parameter qMRI maps (longitudinal and effective transverse relaxation rates R1 and R2⋆, proton density PD and magnetisation transfer MT saturation) that can be used for quantitative parameter analysis and accurate delineation of subcortical brain structures. The qMRI maps generated by the toolbox are key input parameters for biophysical models designed to estimate tissue microstructure properties such as the MR g-ratio and to derive standard and novel MRI biomarkers. Thus, the current version of the toolbox is a first step towards in vivo histology using MRI (hMRI) and is being extended further in this direction. Embedded in the Statistical Parametric Mapping (SPM) framework, it benefits from the extensive range of established SPM tools for high-accuracy spatial registration and statistical inferences and can be readily combined with existing SPM toolboxes for estimating diffusion MRI parameter maps. From a user's perspective, the hMRI-toolbox is an efficient, robust and simple framework for investigating qMRI data in neuroscience and clinical research.

Generative diffeomorphic modelling of large MRI data sets for probabilistic template construction.

In this paper we present a hierarchical generative model of medical image data, which can capture simultaneously the variability of both signal intensity and anatomical shapes across large populations. Such a model has a direct application for learning average-shaped probabilistic tissue templates in a fully automated manner. While in principle the generality of the proposed Bayesian approach makes it suitable to address a wide range of medical image computing problems, our work focuses primarily on neuroimaging applications. In particular we validate the proposed method on both real and synthetic brain MR scans including the cervical cord and demonstrate that it yields accurate alignment of brain and spinal cord structures, as compared to state-of-the-art tools for medical image registration. At the same time we illustrate how the resulting tissue probability maps can readily be used to segment, bias correct and spatially normalise unseen data, which are all crucial pre-processing steps for MR imaging studies.

A comparison of various MRI feature types for characterizing whole brain anatomical differences using linear pattern recognition methods.

There is a widespread interest in applying pattern recognition methods to anatomical neuroimaging data, but so far, there has been relatively little investigation into how best to derive image features in order to make the most accurate predictions. In this work, a Gaussian Process machine learning approach was used for predicting age, gender and body mass index (BMI) of subjects in the IXI dataset, as well as age, gender and diagnostic status using the ABIDE and COBRE datasets. MRI data were segmented and aligned using SPM12, and a variety of feature representations were derived from this preprocessing. We compared classification and regression accuracy using the different sorts of features, and with various degrees of spatial smoothing. Results suggested that feature sets that did not ignore the implicit background tissue class, tended to result in better overall performance, whereas some of the most commonly used feature sets performed relatively poorly.

Multivariate dynamical modelling of structural change during development.

Here we introduce a multivariate framework for characterising longitudinal changes in structural MRI using dynamical systems. The general approach enables modelling changes of states in multiple imaging biomarkers typically observed during brain development, plasticity, ageing and degeneration, e.g. regional gray matter volume of multiple regions of interest (ROIs). Structural brain states follow intrinsic dynamics according to a linear system with additional inputs accounting for potential driving forces of brain development. In particular, the inputs to the system are specified to account for known or latent developmental growth/decline factors, e.g. due to effects of growth hormones, puberty, or sudden behavioural changes etc. Because effects of developmental factors might be region-specific, the sensitivity of each ROI to contributions of each factor is explicitly modelled. In addition to the external effects of developmental factors on regional change, the framework enables modelling and inference about directed (potentially reciprocal) interactions between brain regions, due to competition for space, or structural connectivity, and suchlike. This approach accounts for repeated measures in typical MRI studies of development and aging. Model inversion and posterior distributions are obtained using earlier established variational methods enabling Bayesian evidence-based comparisons between various models of structural change. Using this approach we demonstrate dynamic cortical changes during brain maturation between 6 and 22 years of age using a large openly available longitudinal paediatric dataset with 637 scans from 289 individuals. In particular, we model volumetric changes in 26 bilateral ROIs, which cover large portions of cortical and subcortical gray matter. We account for (1) puberty-related effects on gray matter regions; (2) effects of an early transient growth process with additional time-lag parameter; (3) sexual dimorphism by modelling parameter differences between boys and girls. There is evidence that the regional pattern of sensitivity to dynamic hidden growth factors in late childhood is similar across genders and shows a consistent anterior-posterior gradient with strongest impact to prefrontal cortex (PFC) brain changes. Finally, we demonstrate the potential of the framework to explore the coupling of structural changes across a priori defined subnetworks using an example of previously established resting state functional connectivity.

OpenNFT: An open-source Python/Matlab framework for real-time fMRI neurofeedback training based on activity, connectivity and multivariate pattern analysis.

Neurofeedback based on real-time functional magnetic resonance imaging (rt-fMRI) is a novel and rapidly developing research field. It allows for training of voluntary control over localized brain activity and connectivity and has demonstrated promising clinical applications. Because of the rapid technical developments of MRI techniques and the availability of high-performance computing, new methodological advances in rt-fMRI neurofeedback become possible. Here we outline the core components of a novel open-source neurofeedback framework, termed Open NeuroFeedback Training (OpenNFT), which efficiently integrates these new developments. This framework is implemented using Python and Matlab source code to allow for diverse functionality, high modularity, and rapid extendibility of the software depending on the user's needs. In addition, it provides an easy interface to the functionality of Statistical Parametric Mapping (SPM) that is also open-source and one of the most widely used fMRI data analysis software. We demonstrate the functionality of our new framework by describing case studies that include neurofeedback protocols based on brain activity levels, effective connectivity models, and pattern classification approaches. This open-source initiative provides a suitable framework to actively engage in the development of novel neurofeedback approaches, so that local methodological developments can be easily made accessible to a wider range of users.

Subthalamic deep brain stimulation sweet spots and hyperdirect cortical connectivity in Parkinson's disease.

OBJECTIVES: Firstly, to identify subthalamic region stimulation clusters that predict maximum improvement in rigidity, bradykinesia and tremor, or emergence of side-effects; and secondly, to map-out the cortical fingerprint, mediated by the hyperdirect pathways which predict maximum efficacy. METHODS: High angular resolution diffusion imaging in twenty patients with advanced Parkinson's disease was acquired prior to bilateral subthalamic nucleus deep brain stimulation. All contacts were screened one-year from surgery for efficacy and side-effects at different amplitudes. Voxel-based statistical analysis of volumes of tissue activated models was used to identify significant treatment clusters. Probabilistic tractography was employed to identify cortical connectivity patterns associated with treatment efficacy. RESULTS: All patients responded well to treatment (46% mean improvement off medication UPDRS-III [p < 0.0001]) without significant adverse events. Cluster corresponding to maximum improvement in tremor was in the posterior, superior and lateral portion of the nucleus. Clusters corresponding to improvement in bradykinesia and rigidity were nearer the superior border in a further medial and posterior location. The rigidity cluster extended beyond the superior border to the area of the zona incerta and Forel-H2 field. When the clusters where averaged, the coordinates of the area with maximum overall efficacy was X = -10(-9.5), Y = -13(-1) and Z = -7(-3) in MNI(AC-PC) space. Cortical connectivity to primary motor area was predictive of higher improvement in tremor; whilst that to supplementary motor area was predictive of improvement in bradykinesia and rigidity; and connectivity to prefrontal cortex was predictive of improvement in rigidity. INTERPRETATION: These findings support the presence of overlapping stimulation sites within the subthalamic nucleus and its superior border, with different cortical connectivity patterns, associated with maximum improvement in tremor, rigidity and bradykinesia.

Spinal cord grey matter segmentation challenge.

An important image processing step in spinal cord magnetic resonance imaging is the ability to reliably and accurately segment grey and white matter for tissue specific analysis. There are several semi- or fully-automated segmentation methods for cervical cord cross-sectional area measurement with an excellent performance close or equal to the manual segmentation. However, grey matter segmentation is still challenging due to small cross-sectional size and shape, and active research is being conducted by several groups around the world in this field. Therefore a grey matter spinal cord segmentation challenge was organised to test different capabilities of various methods using the same multi-centre and multi-vendor dataset acquired with distinct 3D gradient-echo sequences. This challenge aimed to characterize the state-of-the-art in the field as well as identifying new opportunities for future improvements. Six different spinal cord grey matter segmentation methods developed independently by various research groups across the world and their performance were compared to manual segmentation outcomes, the present gold-standard. All algorithms provided good overall results for detecting the grey matter butterfly, albeit with variable performance in certain quality-of-segmentation metrics. The data have been made publicly available and the challenge web site remains open to new submissions. No modifications were introduced to any of the presented methods as a result of this challenge for the purposes of this publication.

l-Dopa responsiveness is associated with distinctive connectivity patterns in advanced Parkinson's disease.

BACKGROUND: Neuronal loss and dopamine depletion alter motor signal processing between cortical motor areas, basal ganglia, and the thalamus, resulting in the motor manifestations of Parkinson's disease. Dopamine replacement therapy can reverse these manifestations with varying degrees of improvement. METHODS: To evaluate functional connectivity in patients with advanced Parkinson's disease and changes in functional connectivity in relation to the degree of response to l-dopa, 19 patients with advanced Parkinson's disease underwent resting-state functional magnetic resonance imaging in the on-medication state. Scans were obtained on a 3-Tesla scanner in 3 × 3 × 2.5 mm3 voxels. Seed-based bivariate regression analyses were carried out with atlas-defined basal ganglia regions as seeds, to explore relationships between functional connectivity and improvement in the motor section of the UPDRS-III following an l-dopa challenge. False discovery rate-corrected P was set at < 0.05 for a 2-tailed t test. RESULTS: A greater improvement in UPDRS-III scores following l-dopa administration was characterized by higher resting-state functional connectivity between the prefrontal cortex and the striatum (P = 0.001) and lower resting-state functional connectivity between the pallidum (P = 0.001), subthalamic nucleus (P = 0.003), and the paracentral lobule (supplementary motor area, mesial primary motor, and primary sensory areas). CONCLUSIONS: Our findings show characteristic basal ganglia resting-state functional connectivity patterns associated with different degrees of l-dopa responsiveness in patients with advanced Parkinson's disease. l-Dopa exerts a graduated influence on remapping connectivity in distinct motor control networks, potentially explaining some of the variance in treatment response. © 2017 International Parkinson and Movement Disorder Society.