Investigating apparent differences between standard DKI and axisymmetric DKI and its consequences for biophysical parameter estimates.

PURPOSE: The purpose of the study is to identify differences between axisymmetric diffusion kurtosis imaging (DKI) and standard DKI, their consequences for biophysical parameter estimates, and the protocol choice influence on parameter estimation. METHODS: Noise-free and noisy, synthetic diffusion MRI human brain data is simulated using standard DKI for a standard and the fast "199" acquisition protocol. First the noise-free "baseline" difference between both DKI models is estimated and the influence of fiber complexity is investigated. Noisy data is used to establish the signal-to-noise ratio at which the baseline difference exceeds noise variability. The influence of protocol choices and denoising is investigated. The five axisymmetric DKI tensor metrics (AxTM), the parallel and perpendicular diffusivity and kurtosis and mean of the kurtosis tensor are used to compare both DKI models. Additionally, the baseline difference is also estimated for the five parameters of the WMTI-Watson model. RESULTS: The parallel and perpendicular kurtosis and all of the WMTI-Watson parameters had large baseline differences. Using a Westin or FA mask reduced the number of voxels with large baseline difference, that is, by selecting voxels with less complex fibers. For the noisy data, precision was worsened by the fast "199" protocol but adaptive denoising can help counteract these effects. CONCLUSION: For the diffusivities and mean of the kurtosis tensor, axisymmetric DKI with a standard protocol delivers similar results as standard DKI. Fiber complexity is one main driver of the baseline differences. Using the "199" protocol worsens precision in noisy data but adaptive denoising mitigates these effects.

Insights and improvements in correspondence between axonal volume fraction measured with diffusion-weighted MRI and electron microscopy.

Biophysical diffusion-weighted imaging (DWI) models are increasingly used in neuroscience to estimate the axonal water fraction ( f AW ), which in turn is key for noninvasive estimation of the axonal volume fraction ( f A ). These models require thorough validation by comparison with a reference method, for example, electron microscopy (EM). While EM studies often neglect the unmyelinated axons and solely report the fraction of myelinated axons, in DWI both myelinated and unmyelinated axons contribute to the DWI signal. However, DWI models often include simplifications, for example, the neglect of differences in the compartmental relaxation times or fixed diffusivities, which in turn might affect the estimation of f AW . We investigate whether linear calibration parameters (scaling and offset) can improve the comparability between EM- and DWI-based metrics of f A . To this end, we (a) used six DWI models based on the so-called standard model of white matter (WM), including two models with fixed compartmental diffusivities (e.g., neurite orientation dispersion and density imaging, NODDI) and four models that fitted the compartmental diffusivities (e.g., white matter tract integrity, WMTI), and (b) used a multimodal data set including ex vivo diffusion DWI and EM data in mice with a broad dynamic range of fibre volume metrics. We demonstrated that the offset is associated with the volume fraction of unmyelinated axons and the scaling factor is associated with different compartmental T 2 and can substantially enhance the comparability between EM- and DWI-based metrics of f A . We found that DWI models that fitted compartmental diffusivities provided the most accurate estimates of the EM-based f A . Finally, we introduced a more efficient hybrid calibration approach, where only the offset is estimated but the scaling is fixed to a theoretically predicted value. Using this approach, a similar one-to-one correspondence to EM was achieved for WMTI. The method presented can pave the way for use of validated DWI-based models in clinical research and neuroscience.

Fiber-orientation independent component of R2* obtained from single-orientation MRI measurements in simulations and a post-mortem human optic chiasm.

The effective transverse relaxation rate (R2*) is sensitive to the microstructure of the human brain like the g-ratio which characterises the relative myelination of axons. However, the fibre-orientation dependence of R2* degrades its reproducibility and any microstructural derivative measure. To estimate its orientation-independent part (R2,iso*) from single multi-echo gradient-recalled-echo (meGRE) measurements at arbitrary orientations, a second-order polynomial in time model (hereafter M2) can be used. Its linear time-dependent parameter, ß1, can be biophysically related to R2,iso* when neglecting the myelin water (MW) signal in the hollow cylinder fibre model (HCFM). Here, we examined the performance of M2 using experimental and simulated data with variable g-ratio and fibre dispersion. We found that the fitted ß1 can estimate R2,iso* using meGRE with long maximum-echo time (TEmax ≈ 54 ms), but not accurately captures its microscopic dependence on the g-ratio (error 84%). We proposed a new heuristic expression for ß1 that reduced the error to 12% for ex vivo compartmental R2 values. Using the new expression, we could estimate an MW fraction of 0.14 for fibres with negligible dispersion in a fixed human optic chiasm for the ex vivo compartmental R2 values but not for the in vivo values. M2 and the HCFM-based simulations failed to explain the measured R2*-orientation-dependence around the magic angle for a typical in vivo meGRE protocol (with TEmax ≈ 18 ms). In conclusion, further validation and the development of movement-robust in vivo meGRE protocols with TEmax ≈ 54 ms are required before M2 can be used to estimate R2,iso* in subjects.

Mitigating the impact of flip angle and orientation dependence in single compartment R2* estimates via 2-pool modeling.

PURPOSE: The effective transverse relaxation rate ( R 2 * $$ {\mathrm{R}}_2^{\ast } $$ ) is influenced by biological features that make it a useful means of probing brain microstructure. However, confounding factors such as dependence on flip angle (α) and fiber orientation with respect to the main field ( θ $$ \uptheta $$ ) complicate interpretation. The α- and θ $$ \uptheta $$ -dependence stem from the existence of multiple sub-voxel micro-environments (e.g., myelin and non-myelin water compartments). Ordinarily, it is challenging to quantify these sub-compartments; therefore, neuroscientific studies commonly make the simplifying assumption of a mono-exponential decay obtaining a single R 2 * $$ {\mathrm{R}}_2^{\ast } $$ estimate per voxel. In this work, we investigated how the multi-compartment nature of tissue microstructure affects single compartment R 2 * $$ {\mathrm{R}}_2^{\ast } $$ estimates. METHODS: We used 2-pool (myelin and non-myelin water) simulations to characterize the bias in single compartment R 2 * $$ {\mathrm{R}}_2^{\ast } $$ estimates. Based on our numeric observations, we introduced a linear model that partitions R 2 * $$ {\mathrm{R}}_2^{\ast } $$ into α-dependent and α-independent components and validated this in vivo at 7T. We investigated the dependence of both components on the sub-compartment properties and assessed their robustness, orientation dependence, and reproducibility empirically. RESULTS: R 2 * $$ {\mathrm{R}}_2^{\ast } $$ increased with myelin water fraction and residency time leading to a linear dependence on α. We observed excellent agreement between our numeric and empirical results. Furthermore, the α-independent component of the proposed linear model was robust to the choice of α and reduced dependence on fiber orientation, although it suffered from marginally higher noise sensitivity. CONCLUSION: We have demonstrated and validated a simple approach that mitigates flip angle and orientation biases in single-compartment R 2 * $$ {\mathrm{R}}_2^{\ast } $$ estimates.

Axisymmetric diffusion kurtosis imaging with Rician bias correction: A simulation study.

PURPOSE: To compare the estimation accuracy of axisymmetric diffusion kurtosis imaging (DKI) and standard DKI in combination with Rician bias correction (RBC). METHODS: Axisymmetric DKI is more robust against noise-induced variation in the measured signal than standard DKI because of its reduced parameter space. However, its susceptibility to Rician noise bias at low signal-to-noise ratios (SNR) is unknown. Here, we investigate two main questions: first, does RBC improve estimation accuracy of axisymmetric DKI?; second, is estimation accuracy of axisymmetric DKI increased compared to standard DKI? Estimation accuracy was investigated on the five axisymmetric DKI tensor metrics (AxTM): the parallel and perpendicular diffusivity and kurtosis and mean of the kurtosis tensor, using a noise simulation study based on synthetic data of tissues with varying fiber alignment and in-vivo data focusing on white matter. RESULTS: RBC mainly increased accuracy for the parallel AxTM in tissues with highly to moderately aligned fibers. For the perpendicular AxTM, axisymmetric DKI without RBC performed slightly better than with RBC. However, the combination of axisymmetric DKI with RBC was the overall best performing algorithm across all five AxTM in white matter and axisymmetric DKI itself substantially improved accuracy in axisymmetric tissues with low fiber alignment. CONCLUSION: Combining axisymmetric DKI with RBC facilitates accurate DKI parameter estimation at unprecedented low SNRs ( ≈ 15 $$ \approx 15 $$ ) in white matter, possibly making it a valuable tool for neuroscience and clinical research studies where scan time is a limited resource. The tools used here are available in the open-source ACID toolbox for SPM.

Conduction velocity along a key white matter tract is associated with autobiographical memory recall ability.

Conduction velocity is the speed at which electrical signals travel along axons and is a crucial determinant of neural communication. Inferences about conduction velocity can now be made in vivo in humans using a measure called the magnetic resonance (MR) g-ratio. This is the ratio of the inner axon diameter relative to that of the axon plus the myelin sheath that encases it. Here, in the first application to cognition, we found that variations in MR g-ratio, and by inference conduction velocity, of the parahippocampal cingulum bundle were associated with autobiographical memory recall ability in 217 healthy adults. This tract connects the hippocampus with a range of other brain areas. We further observed that the association seemed to be with inner axon diameter rather than myelin content. The extent to which neurites were coherently organised within the parahippocampal cingulum bundle was also linked with autobiographical memory recall ability. Moreover, these findings were specific to autobiographical memory recall and were not apparent for laboratory-based memory tests. Our results offer a new perspective on individual differences in autobiographical memory recall ability, highlighting the possible influence of specific white matter microstructure features on conduction velocity when recalling detailed memories of real-life past experiences.

Error quantification in multi-parameter mapping facilitates robust estimation and enhanced group level sensitivity.

Multi-Parameter Mapping (MPM) is a comprehensive quantitative neuroimaging protocol that enables estimation of four physical parameters (longitudinal and effective transverse relaxation rates R1 and R2*, proton density PD, and magnetization transfer saturation MTsat) that are sensitive to microstructural tissue properties such as iron and myelin content. Their capability to reveal microstructural brain differences, however, is tightly bound to controlling random noise and artefacts (e.g. caused by head motion) in the signal. Here, we introduced a method to estimate the local error of PD, R1, and MTsat maps that captures both noise and artefacts on a routine basis without requiring additional data. To investigate the method's sensitivity to random noise, we calculated the model-based signal-to-noise ratio (mSNR) and showed in measurements and simulations that it correlated linearly with an experimental raw-image-based SNR map. We found that the mSNR varied with MPM protocols, magnetic field strength (3T vs. 7T) and MPM parameters: it halved from PD to R1 and decreased from PD to MTsat by a factor of 3-4. Exploring the artefact-sensitivity of the error maps, we generated robust MPM parameters using two successive acquisitions of each contrast and the acquisition-specific errors to down-weight erroneous regions. The resulting robust MPM parameters showed reduced variability at the group level as compared to their single-repeat or averaged counterparts. The error and mSNR maps may better inform power-calculations by accounting for local data quality variations across measurements. Code to compute the mSNR maps and robustly combined MPM maps is available in the open-source hMRI toolbox.

Towards a representative reference for MRI-based human axon radius assessment using light microscopy.

Non-invasive assessment of axon radii via MRI bears great potential for clinical and neuroscience research as it is a main determinant of the neuronal conduction velocity. However, there is a lack of representative histological reference data at the scale of the cross-section of MRI voxels for validating the MRI-visible, effective radius (reff). Because the current gold standard stems from neuroanatomical studies designed to estimate the bulk-determined arithmetic mean radius (rarith) on small ensembles of axons, it is unsuited to estimate the tail-weighted reff. We propose CNN-based segmentation on high-resolution, large-scale light microscopy (lsLM) data to generate a representative reference for reff. In a human corpus callosum, we assessed estimation accuracy and bias of rarith and reff. Furthermore, we investigated whether mapping anatomy-related variation of rarith and reff is confounded by low-frequency variation of the image intensity, e.g., due to staining heterogeneity. Finally, we analyzed the error due to outstandingly large axons in reff. Compared to rarith, reff was estimated with higher accuracy (maximum normalized-root-mean-square-error of reff: 8.5 %; rarith: 19.5 %) and lower bias (maximum absolute normalized-mean-bias-error of reff: 4.8 %; rarith: 13.4 %). While rarith was confounded by variation of the image intensity, variation of reff seemed anatomy-related. The largest axons contributed between 0.8 % and 2.9 % to reff. In conclusion, the proposed method is a step towards representatively estimating reff at MRI voxel resolution. Further investigations are required to assess generalization to other brains and brain areas with different axon radii distributions.

Finding the best clearing approach - Towards 3D wide-scale multimodal imaging of aged human brain tissue.

The accessibility of new wide-scale multimodal imaging techniques led to numerous clearing techniques emerging over the last decade. However, clearing mesoscopic-sized blocks of aged human brain tissue remains an extremely challenging task. Homogenizing refractive indices and reducing light absorption and scattering are the foundation of tissue clearing. Due to its dense and highly myelinated nature, especially in white matter, the human brain poses particular challenges to clearing techniques. Here, we present a comparative study of seven tissue clearing approaches and their impact on aged human brain tissue blocks (> 5 mm). The goal was to identify the most practical and efficient method in regards to macroscopic transparency, brief clearing time, compatibility with immunohistochemical processing and wide-scale multimodal microscopic imaging. We successfully cleared 26 × 26 × 5 mm3-sized human brain samples with two hydrophilic and two hydrophobic clearing techniques. Optical properties as well as light and antibody penetration depths highly vary between these methods. In addition to finding the best clearing approach, we compared three microscopic imaging setups (the Zeiss Laser Scanning Microscope (LSM) 880 , the Miltenyi Biotec Ultramicroscope ll (UM ll) and the 3i Marianas LightSheet microscope) regarding optimal imaging of large-scale tissue samples. We demonstrate that combining the CLARITY technique (Clear Lipid-exchanged Acrylamide-hybridized Rigid Imaging compatible Tissue hYdrogel) with the Zeiss LSM 880 and combining the iDISCO technique (immunolabeling-enabled three-dimensional imaging of solvent-cleared organs) with the Miltenyi Biotec UM ll are the most practical and efficient approaches to sufficiently clear aged human brain tissue and generate 3D microscopic images. Our results point out challenges that arise from seven clearing and three imaging techniques applied to non-standardized tissue samples such as aged human brain tissue.

Corrigendum: The Influence of Radio-Frequency Transmit Field Inhomogeneities on the Accuracy of G-ratio Weighted Imaging.

[This corrects the article DOI: 10.3389/fnins.2021.674719.].

Longitudinal changes of spinal cord grey and white matter following spinal cord injury.

OBJECTIVES: Traumatic and non-traumatic spinal cord injury produce neurodegeneration across the entire neuraxis. However, the spatiotemporal dynamics of spinal cord grey and white matter neurodegeneration above and below the injury is understudied. METHODS: We acquired longitudinal data from 13 traumatic and 3 non-traumatic spinal cord injury patients (8-8 cervical and thoracic cord injuries) within 1.5 years after injury and 10 healthy controls over the same period. The protocol encompassed structural and diffusion-weighted MRI rostral (C2/C3) and caudal (lumbar enlargement) to the injury level to track tissue-specific neurodegeneration. Regression models assessed group differences in the temporal evolution of tissue-specific changes and associations with clinical outcomes. RESULTS: At 2 months post-injury, white matter area was decreased by 8.5% and grey matter by 15.9% in the lumbar enlargement, while at C2/C3 only white matter was decreased (-9.7%). Patients had decreased cervical fractional anisotropy (FA: -11.3%) and increased radial diffusivity (+20.5%) in the dorsal column, while FA was lower in the lateral (-10.3%) and ventral columns (-9.7%) of the lumbar enlargement. White matter decreased by 0.34% and 0.35% per month at C2/C3 and lumbar enlargement, respectively, and grey matter decreased at C2/C3 by 0.70% per month. CONCLUSIONS: This study describes the spatiotemporal dynamics of tissue-specific spinal cord neurodegeneration above and below a spinal cord injury. While above the injury, grey matter atrophy lagged initially behind white matter neurodegeneration, in the lumbar enlargement these processes progressed in parallel. Tracking trajectories of tissue-specific neurodegeneration provides valuable assessment tools for monitoring recovery and treatment effects.

Reducing Susceptibility Distortion Related Image Blurring in Diffusion MRI EPI Data.

Diffusion magnetic resonance imaging (MRI) is an increasingly popular technique in basic and clinical neuroscience. One promising application is to combine diffusion MRI with myelin maps from complementary MRI techniques such as multi-parameter mapping (MPM) to produce g-ratio maps that represent the relative myelination of axons and predict their conduction velocity. Statistical Parametric Mapping (SPM) can process both diffusion data and MPMs, making SPM the only widely accessible software that contains all the processing steps required to perform group analyses of g-ratio data in a common space. However, limitations have been identified in its method for reducing susceptibility-related distortion in diffusion data. More generally, susceptibility-related image distortion is often corrected by combining reverse phase-encoded images (blip-up and blip-down) using the arithmetic mean (AM), however, this can lead to blurred images. In this study we sought to (1) improve the susceptibility-related distortion correction for diffusion MRI data in SPM; (2) deploy an alternative approach to the AM to reduce image blurring in diffusion MRI data when combining blip-up and blip-down EPI data after susceptibility-related distortion correction; and (3) assess the benefits of these changes for g-ratio mapping. We found that the new processing pipeline, called consecutive Hyperelastic Susceptibility Artefact Correction (HySCO) improved distortion correction when compared to the standard approach in the ACID toolbox for SPM. Moreover, using a weighted average (WA) method to combine the distortion corrected data from each phase-encoding polarity achieved greater overlap of diffusion and more anatomically faithful structural white matter probability maps derived from minimally distorted multi-parameter maps as compared to the AM. Third, we showed that the consecutive HySCO WA performed better than the AM method when combined with multi-parameter maps to perform g-ratio mapping. These improvements mean that researchers can conveniently access a wide range of diffusion-related analysis methods within one framework because they are now available within the open-source ACID toolbox as part of SPM, which can be easily combined with other SPM toolboxes, such as the hMRI toolbox, to facilitate computation of myelin biomarkers that are necessary for g-ratio mapping.

The Influence of Radio-Frequency Transmit Field Inhomogeneities on the Accuracy of G-ratio Weighted Imaging.

G-ratio weighted imaging is a non-invasive, in-vivo MRI-based technique that aims at estimating an aggregated measure of relative myelination of axons across the entire brain white matter. The MR g-ratio and its constituents (axonal and myelin volume fraction) are more specific to the tissue microstructure than conventional MRI metrics targeting either the myelin or axonal compartment. To calculate the MR g-ratio, an MRI-based myelin-mapping technique is combined with an axon-sensitive MR technique (such as diffusion MRI). Correction for radio-frequency transmit (B1+) field inhomogeneities is crucial for myelin mapping techniques such as magnetization transfer saturation. Here we assessed the effect of B1+ correction on g-ratio weighted imaging. To this end, the B1+ field was measured and the B1+ corrected MR g-ratio was used as the reference in a Bland-Altman analysis. We found a substantial bias (≈-89%) and error (≈37%) relative to the dynamic range of g-ratio values in the white matter if the B1+ correction was not applied. Moreover, we tested the efficiency of a data-driven B1+ correction approach that was applied retrospectively without additional reference measurements. We found that it reduced the bias and error in the MR g-ratio by a factor of three. The data-driven correction is readily available in the open-source hMRI toolbox (www.hmri.info) which is embedded in the statistical parameter mapping (SPM) framework.

Towards in vivo g-ratio mapping using MRI: Unifying myelin and diffusion imaging.

BACKGROUND: The g-ratio, quantifying the comparative thickness of the myelin sheath encasing an axon, is a geometrical invariant that has high functional relevance because of its importance in determining neuronal conduction velocity. Advances in MRI data acquisition and signal modelling have put in vivo mapping of the g-ratio, across the entire white matter, within our reach. This capacity would greatly increase our knowledge of the nervous system: how it functions, and how it is impacted by disease. NEW METHOD: This is the second review on the topic of g-ratio mapping using MRI. RESULTS: This review summarizes the most recent developments in the field, while also providing methodological background pertinent to aggregate g-ratio weighted mapping, and discussing pitfalls associated with these approaches. COMPARISON WITH EXISTING METHODS: Using simulations based on recently published data, this review reveals caveats to the state-of-the-art calibration methods that have been used for in vivo g-ratio mapping. It highlights the need to estimate both the slope and offset of the relationship between these MRI-based markers and the true myelin volume fraction if we are really to achieve the goal of precise, high sensitivity g-ratio mapping in vivo. Other challenges discussed in this review further evidence the need for gold standard measurements of human brain tissue from ex vivo histology. CONCLUSIONS: We conclude that the quest to find the most appropriate MRI biomarkers to enable in vivo g-ratio mapping is ongoing, with the full potential of many novel techniques yet to be investigated.

Traumatic and nontraumatic spinal cord injury: pathological insights from neuroimaging.

Pathophysiological changes in the spinal cord white and grey matter resulting from injury can be observed with MRI techniques. These techniques provide sensitive markers of macrostructural and microstructural tissue integrity, which correlate with histological findings. Spinal cord MRI findings in traumatic spinal cord injury (tSCI) and nontraumatic spinal cord injury - the most common form of which is degenerative cervical myelopathy (DCM) - have provided important insights into the pathophysiological processes taking place not just at the focal injury site but also rostral and caudal to the spinal injury. Although tSCI and DCM have different aetiologies, they show similar degrees of spinal cord pathology remote from the injury site, suggesting the involvement of similar secondary degenerative mechanisms. Advanced quantitative MRI protocols that are sensitive to spinal cord pathology have the potential to improve diagnosis and, more importantly, predict outcomes in patients with tSCI or nontraumatic spinal cord injury. This Review describes the insights into tSCI and DCM that have been revealed by neuroimaging and outlines current activities and future directions for the field.

Biophysically motivated efficient estimation of the spatially isotropic R2* component from a single gradient-recalled echo measurement.

PURPOSE: To propose and validate an efficient method, based on a biophysically motivated signal model, for removing the orientation-dependent part of R2* using a single gradient-recalled echo (GRE) measurement. METHODS: The proposed method utilized a temporal second-order approximation of the hollow-cylinder-fiber model, in which the parameter describing the linear signal decay corresponded to the orientation-independent part of R2* . The estimated parameters were compared to the classical, mono-exponential decay model for R2* in a sample of an ex vivo human optic chiasm (OC). The OC was measured at 16 distinct orientations relative to the external magnetic field using GRE at 7T. To show that the proposed signal model can remove the orientation dependence of R2* , it was compared to the established phenomenological method for separating R2* into orientation-dependent and -independent parts. RESULTS: Using the phenomenological method on the classical signal model, the well-known separation of R2* into orientation-dependent and -independent parts was verified. For the proposed model, no significant orientation dependence in the linear signal decay parameter was observed. CONCLUSIONS: Since the proposed second-order model features orientation-dependent and -independent components at distinct temporal orders, it can be used to remove the orientation dependence of R2* using only a single GRE measurement.

Example dataset for the hMRI toolbox.

The hMRI toolbox is an open-source toolbox for the calculation of quantitative MRI parameter maps from a series of weighted imaging data, and optionally additional calibration data. The multi-parameter mapping (MPM) protocol, incorporating calibration data to correct for spatial variation in the scanner's transmit and receive fields, is the most complete protocol that can be handled by the toolbox. Here we present a dataset acquired with such a full MPM protocol, which is made freely available to be used as a tutorial by following instructions provided on the associated toolbox wiki pages, which can be found at http://hMRI.info, and following the theory described in: hMRI - A toolbox for quantitative MRI in neuroscience and clinical research [1].

In vivo evidence of remote neural degeneration in the lumbar enlargement after cervical injury.

OBJECTIVE: To characterize remote secondary neurodegeneration of spinal tracts and neurons below a cervical spinal cord injury (SCI) and its relation to the severity of injury, the integrity of efferent and afferent pathways, and clinical impairment. METHODS: A comprehensive high-resolution MRI protocol was acquired in 17 traumatic cervical SCI patients and 14 controls at 3T. At the cervical lesion, a sagittal T2-weighted scan provided information on the width of preserved midsagittal tissue bridges. In the lumbar enlargement, high-resolution T2*-weighted and diffusion-weighted scans were used to calculate tissue-specific cross-sectional areas and diffusion indices, respectively. Regression analyses determined associations between MRI readouts and the electrophysiologic and clinical measures. RESULTS: At the cervical injury level, preserved midsagittal tissue bridges were present in the majority of patients. In the lumbar enlargement, neurodegeneration-in terms of macrostructural and microstructural MRI changes-was evident in the white matter and ventral and dorsal horns. Patients with thinner midsagittal tissue bridges had smaller ventral horn area, higher radial diffusivity in the gray matter, smaller motor evoked potential amplitude from the lower extremities, and lower motor score. In addition, smaller width of midsagittal tissue bridges was also associated with smaller tibialis sensory evoked potential amplitude and lower light-touch score. CONCLUSIONS: This study shows extensive tissue-specific cord pathology in infralesional spinal networks following cervical SCI, its magnitude relating to lesion severity, electrophysiologic integrity, and clinical impairment of the lower extremity. The clinical eloquence of remote neurodegenerative changes speaks to the application of neuroimaging biomarkers in diagnostic workup and planning of clinical trials.

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.

Dynamic changes in white matter microstructure in anorexia nervosa: findings from a longitudinal study.

BACKGROUND: Gray matter (GM) 'pseudoatrophy' is well-documented in patients with anorexia nervosa (AN), but changes in white matter (WM) are less well understood. Here we investigated the dynamics of microstructural WM brain changes in AN patients during short-term weight restoration in a combined longitudinal and cross-sectional study design. METHODS: Diffusion-weighted images were acquired in young AN patients before (acAN-Tp1, n = 56) and after (acAN-Tp2, n = 44) short-term weight restoration as well as in age-matched healthy controls (HC, n = 60). Images were processed using Tract-Based-Spatial-Statistics to compare fractional anisotropy (FA) across groups and timepoints. RESULTS: In the cross-sectional comparison, FA was significantly reduced in the callosal body in acAN-Tp1 compared with HC, while no differences were found between acAN-Tp2 and HC. In the longitudinal arm, FA increased with weight gain in acAN-Tp2 relative to acAN-Tp1 in large parts of the callosal body and the fornix, while it decreased in the right corticospinal tract. CONCLUSIONS: Our findings reveal that dynamic, bidirectional changes in WM microstructure in young underweight patients with AN can be reversed with brief weight restoration therapy. These results parallel those previously observed in GM and suggest that alterations in WM in non-chronic AN are also state-dependent and rapidly reversible with successful intervention.