Brain microstructure by multi-modal MRI: Is the whole greater than the sum of its parts?

Metabolic checkpoints in glioblastomas: targets for new therapies and non-invasive detection

Glioblastoma (GBM) is a highly malignant tumor of the central nervous system that remains intractable despite advancements in current tumor treatment modalities, including immunotherapy. In recent years, metabolic checkpoints (aberrant metabolic pathways underlying the immunosuppressive tumor microenvironment) have gained attention as promising therapeutic targets and sensitive biomarkers across various cancers. Here, we briefly review the existing understanding of tumor metabolic checkpoints and their implications in the biology and management of GBM. Additionally, we discuss techniques that could evaluate metabolic checkpoints of GBM non-invasively, thereby potentially facilitating neo-adjuvant treatment and dynamic surveillance.

A genetic optimisation and iterative reconstruction framework for sparse multi-dimensional diffusion-relaxation correlation MRI

Multi-dimensional diffusion-relaxation correlation (DRC) magnetic resonance imaging (MRI) techniques have recently been developed to investigate tissue microstructures. Sub-voxel tissue heterogeneity is resolved from the local correlation distributions of relaxation time and molecular diffusivity. However, the implementation of these techniques considerably increases the total acquisition time, and simply reducing the scan time may be at the expense of detailed structural resolution. To overcome these limitations, an optimised framework was proposed for acquiring microstructural maps of the human brain on a clinically feasible timescale. First, the acquisition parameters of the multi-dimensional DRC MRI method were sparsely optimised using a genetic algorithm with a fitness function according to the spectral resolution of the correlation map, hardware requirements, and total scan time. Next, the acquired DRC MRI data were processed using a proposed numerical algorithm based on the dynamic inverse Laplace transform (ILT). Prior knowledge from one-dimensional data was then utilised in the iterative procedure to improve the spectral resolution. Finally, the proposed framework was validated using Monte Carlo simulations and experimental data acquired from healthy participants on an MRI scanner. The results demonstrated that the suggested approach is feasible for offering high-resolution DRC maps that correspond to distinct microstructures with a limited amount of optimised acquisition data from two-dimensional DRC sampling space. By significantly reducing scan time while retaining structural resolution, this approach may enable multi-dimensional DRC MRI to be more widely used for quantitative evaluation in biological and medical settings.

Optimisation of quantitative brain diffusion-relaxation MRI acquisition protocols with physics-informed machine learning

Diffusion-relaxation MRI aims to extract quantitative measures that characterise microstructural tissue properties such as orientation, size, and shape, but long acquisition times are typically required. This work proposes a physics-informed learning framework to extract an optimal subset of diffusion-relaxation MRI measurements for enabling shorter acquisition times, predict non-measured signals, and estimate quantitative parameters. In vivo and synthetic brain 5D-Diffusion-T1-T2∗-weighted MRI data obtained from five healthy subjects were used for training and validation, and from a sixth participant for testing. One fully data-driven and two physics-informed machine learning methods were implemented and compared to two manual selection procedures and Cramér-Rao lower bound optimisation. The physics-informed approaches could identify measurement-subsets that yielded more consistently accurate parameter estimates in simulations than other approaches, with similar signal prediction error. Five-fold shorter protocols yielded error distributions of estimated quantitative parameters with very small effect sizes compared to estimates from the full protocol. Selected subsets commonly included a denser sampling of the shortest and longest inversion time, lowest echo time, and high b-value. The proposed framework combining machine learning and MRI physics offers a promising approach to develop shorter imaging protocols without compromising the quality of parameter estimates and signal predictions.

A rapid multi-parametric quantitative MR imaging method to assess Parkinson's disease: a feasibility study

BACKGROUND

MULTIPLEX is a single-scan three-dimensional multi-parametric MRI technique that provides 1 mm isotropic T1-, T2*-, proton density- and susceptibility-weighted images and the corresponding quantitative maps. This study aimed to investigate its feasibility of clinical application in Parkinson's disease (PD).

METHODS

27 PD patients and 23 healthy control (HC) were recruited and underwent a MULTIPLEX scanning. All image reconstruction and processing were automatically performed with in-house C + + programs on the Automatic Differentiation using Expression Template platform. According to the HybraPD atlas consisting of 12 human brain subcortical nuclei, the region-of-interest (ROI) based analysis was conducted to extract quantitative parameters, then identify PD-related abnormalities from the T1, T2* and proton density maps and quantitative susceptibility mapping (QSM), by comparing patients and HCs.

RESULTS

The ROI-based analysis revealed significantly decreased mean T1 values in substantia nigra pars compacta and habenular nuclei, mean T2* value in subthalamic nucleus and increased mean QSM value in subthalamic nucleus in PD patients, compared to HCs (all p values < 0.05 after FDR correction). The receiver operating characteristic analysis showed all these four quantitative parameters significantly contributed to PD diagnosis (all p values < 0.01 after FDR correction). Furthermore, the two quantitative parameters in subthalamic nucleus showed hemicerebral differences in regard to the clinically dominant side among PD patients.

CONCLUSIONS

MULTIPLEX might be feasible for clinical application to assist in PD diagnosis and provide possible pathological information of PD patients' subcortical nucleus and dopaminergic midbrain regions.

Analyzing the cortical fine structure as revealed by ex-vivo anatomical MRI

The human cortex has a rich fiber structure as revealed by myelin-staining of histological slices. Myelin also contributes to the image contrast in Magnetic Resonance Imaging (MRI). Recent advances in Magnetic Resonance (MR) scanner and imaging technology allowed the acquisition of an ex-vivo data set at an isotropic resolution of 100 µm. This study focused on a computational analysis of this data set with the aim of bridging between histological knowledge and MRI-based results. This work highlights: (1) the design and implementation of a processing chain that extracts intracortical features from a high-resolution MR image; (2) a demonstration of the correspondence between MRI-based cortical intensity profiles and the myelo-architectonic layering of the cortex; (3) the characterization and classification of four basic myelo-architectonic profile types; (4) the distinction of cortical regions based on myelo-architectonic features; and (5) the segmentation of cortical modules in the entorhinal cortex.

Multiparametric MRI: From Simultaneous Rapid Acquisition Methods and Analysis Techniques Using Scoring, Machine Learning, Radiomics, and Deep Learning to the Generation of Novel Metrics

ABSTRACT

With the recent advancements in rapid imaging methods, higher numbers of contrasts and quantitative parameters can be acquired in less and less time. Some acquisition models simultaneously obtain multiparametric images and quantitative maps to reduce scan times and avoid potential issues associated with the registration of different images. Multiparametric magnetic resonance imaging (MRI) has the potential to provide complementary information on a target lesion and thus overcome the limitations of individual techniques. In this review, we introduce methods to acquire multiparametric MRI data in a clinically feasible scan time with a particular focus on simultaneous acquisition techniques, and we discuss how multiparametric MRI data can be analyzed as a whole rather than each parameter separately. Such data analysis approaches include clinical scoring systems, machine learning, radiomics, and deep learning. Other techniques combine multiple images to create new quantitative maps associated with meaningful aspects of human biology. They include the magnetic resonance g-ratio, the inner to the outer diameter of a nerve fiber, and the aerobic glycolytic index, which captures the metabolic status of tumor tissues.

BTMF-GAN: A multi-modal MRI fusion generative adversarial network for brain tumors

Image fusion techniques have been widely used for multi-modal medical image fusion tasks. Most existing methods aim to improve the overall quality of the fused image and do not focus on the more important textural details and contrast between the tissues of the lesion in the regions of interest (ROIs). This can lead to the distortion of important tumor ROIs information and thus limits the applicability of the fused images in clinical practice. To improve the fusion quality of ROIs relevant to medical implications, we propose a multi-modal MRI fusion generative adversarial network (BTMF-GAN) for the task of multi-modal MRI fusion of brain tumors. Unlike existing deep learning approaches which focus on improving the global quality of the fused image, the proposed BTMF-GAN aims to achieve a balance between tissue details and structural contrasts in brain tumor, which is the region of interest crucial to many medical applications. Specifically, we employ a generator with a U-shaped nested structure and residual U-blocks (RSU) to enhance multi-scale feature extraction. To enhance and recalibrate features of the encoder, the multi-perceptual field adaptive transformer feature enhancement module (MRF-ATFE) is used between the encoder and the decoder instead of a skip connection. To increase contrast between tumor tissues of the fused image, a mask-part block is introduced to fragment the source image and the fused image, based on which, we propose a novel salient loss function. Qualitative and quantitative analysis of the results on the public and clinical datasets demonstrate the superiority of the proposed approach to many other commonly used fusion methods.

Towards automated in vivo parcellation of the human cerebral cortex using supervised classification of magnetic resonance fingerprinting residuals

BACKGROUND

Accurate parcellation of the cerebral cortex in an individual is a guide to its underlying organization. The most promising in vivo quantitative magnetic resonance (MR)-based microstructural cortical mapping methods are yet to achieve a level of parcellation accuracy comparable to quantitative histology.

METHODS

We scanned 6 participants using a 3D echo-planar imaging MR fingerprinting (EPI-MRF) sequence on a 7T Siemens scanner. After projecting MRF signals to the individual-specific inflated model of the cortical surface, normalized autocorrelations of MRF residuals of vertices of 8 microstructurally distinct areas (BA1, BA2, BA4a, BA6, BA44, BA45, BA17, and BA18) from 3 cortical regions were used as feature vector inputs into linear support vector machine (SVM), radial basis function SVM (RBF-SVM), random forest, and k-nearest neighbors supervised classification algorithms. The algorithms' prediction performance was compared using: (i) features from each vertex or (ii) features from neighboring vertices.

RESULTS

The neighborhood-based RBF-SVM classifier achieved the highest prediction score of 0.85 for classification of MRF residuals in the central region from a held-out participant.

CONCLUSIONS

We developed an automated method of cortical parcellation using a combination of MR fingerprinting residual analysis and machine learning classification. Our findings provide the basis for employing unsupervised learning algorithms for whole-cortex structural parcellation in individuals.

Neuroprotection and neuroregeneration: roles for the white matter

Efficient strategies for neuroprotection and repair are still an unmet medical need for neurodegenerative diseases and lesions of the central nervous system. Over the last few decades, a great deal of attention has been focused on white matter as a potential therapeutic target, mainly due to the discovery of the oligodendrocyte precursor cells in the adult central nervous system, a cell type able to fully repair myelin damage, and to the development of advanced imaging techniques to visualize and measure white matter lesions. The combination of these two events has greatly increased the body of research into white matter alterations in central nervous system lesions and neurodegenerative diseases and has identified the oligodendrocyte precursor cell as a putative target for white matter lesion repair, thus indirectly contributing to neuroprotection. This review aims to discuss the potential of white matter as a therapeutic target for neuroprotection in lesions and diseases of the central nervous system. Pivot conditions are discussed, specifically multiple sclerosis as a white matter disease; spinal cord injury, the acute lesion of a central nervous system component where white matter prevails over the gray matter, and Alzheimer's disease, where the white matter was considered an ancillary component until recently. We first describe oligodendrocyte precursor cell biology and developmental myelination, and its regulation by thyroid hormones, then briefly describe white matter imaging techniques, which are providing information on white matter involvement in central nervous system lesions and degenerative diseases. Finally, we discuss pathological mechanisms which interfere with myelin repair in adulthood.

Primary Multiparametric Quantitative Brain MRI: State-of-the-Art Relaxometric and Proton Density Mapping Techniques

This review on brain multiparametric quantitative MRI (MP-qMRI) focuses on the primary subset of quantitative MRI (qMRI) parameters that represent the mobile ("free") and bound ("motion-restricted") proton pools. Such primary parameters are the proton densities, relaxation times, and magnetization transfer parameters. Diffusion qMRI is also included because of its wide implementation in complete clinical MP-qMRI application. MP-qMRI advances were reviewed over the past 2 decades, with substantial progress observed toward accelerating image acquisition and increasing mapping accuracy. Areas that need further investigation and refinement are identified as follows: (a) the biologic underpinnings of qMRI parameter values and their changes with age and/or disease and (b) the theoretical limitations implicitly built into most qMRI mapping algorithms that do not distinguish between the different spatial scales of voxels versus spin packets, the central physical object of the Bloch theory. With rapidly improving image processing techniques and continuous advances in computer hardware, MP-qMRI has the potential for implementation in a wide range of clinical applications. Currently, three emerging MP-qMRI applications are synthetic MRI, macrostructural qMRI, and microstructural tissue modeling.

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 R₁ and R₂^*, proton density PD, and magnetization transfer saturation MT_sat) 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, R₁, and MT_sat 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 R₁ and decreased from PD to MT_sat 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.

Rapid processing and quantitative evaluation of structural brain scans for adaptive multimodal imaging

Current neuroimaging acquisition and processing approaches tend to be optimised for quality rather than speed. However, rapid acquisition and processing of neuroimaging data can lead to novel neuroimaging paradigms, such as adaptive acquisition, where rapidly processed data is used to inform subsequent image acquisition steps. Here we first evaluate the impact of several processing steps on the processing time and quality of registration of manually labelled T1 -weighted MRI scans. Subsequently, we apply the selected rapid processing pipeline both to rapidly acquired multicontrast EPImix scans of 95 participants (which include T1 -FLAIR, T2 , T2 *, T2 -FLAIR, DWI and ADC contrasts, acquired in ~1 min), as well as to slower, more standard single-contrast T1 -weighted scans of a subset of 66 participants. We quantify the correspondence between EPImix T1 -FLAIR and single-contrast T1 -weighted scans, using correlations between voxels and regions of interest across participants, measures of within- and between-participant identifiability as well as regional structural covariance networks. Furthermore, we explore the use of EPImix for the rapid construction of morphometric similarity networks. Finally, we quantify the reliability of EPImix-derived data using test-retest scans of 10 participants. Our results demonstrate that quantitative information can be derived from a neuroimaging scan acquired and processed within minutes, which could further be used to implement adaptive multimodal imaging and tailor neuroimaging examinations to individual patients.

Virtual mouse brain histology from multi-contrast MRI via deep learning

¹H MRI maps brain structure and function non-invasively through versatile contrasts that exploit inhomogeneity in tissue micro-environments. Inferring histopathological information from magnetic resonance imaging (MRI) findings, however, remains challenging due to absence of direct links between MRI signals and cellular structures. Here, we show that deep convolutional neural networks, developed using co-registered multi-contrast MRI and histological data of the mouse brain, can estimate histological staining intensity directly from MRI signals at each voxel. The results provide three-dimensional maps of axons and myelin with tissue contrasts that closely mimic target histology and enhanced sensitivity and specificity compared to conventional MRI markers. Furthermore, the relative contribution of each MRI contrast within the networks can be used to optimize multi-contrast MRI acquisition. We anticipate our method to be a starting point for translation of MRI results into easy-to-understand virtual histology for neurobiologists and provide resources for validating novel MRI techniques.

Magnetic resonance fingerprinting residual signals can disassociate human grey matter regions

The importance of accurate structural discrimination of the human grey matter regions has motivated the development of observer-independent reproducible methods that account for inter-individual architectonic variations. We introduce a non-invasive statistical residual analysis framework, employing unique tissue-specific magnetic resonance fingerprinting (MRF) signals after adjusting for the effect of T₁ and T₂^* MR relaxometry parameters (here termed MRF residuals). A 7 T Siemens MR scanner was used to acquire MRF signals, quantitative transmit magnetic field (B₁⁺) maps and T₁-weighted anatomical images of eleven cortical areas (5L, 5M, 5Ci, 7A, 7P, 7PC, hIP3, BA2, BA4a, BA4p and BA6) from six female participants. MRF residual signal for each voxel was calculated as the difference between the actual and best matching MRF signal evolutions from a precomputed MRF dictionary covering a range of T₁, T₂^* and B₁⁺ values. To compare MRF residuals between regions of interest, normalised autocorrelation was used as a shape-based statistical signal characterisation method and the Euclidean distance between autocorrelation profiles of residuals was used to measure the interareal dissimilarity. In the eleven cortical areas in both cerebral hemispheres of six participants, the proposed MRF residual analysis consistently showed interareal dissimilarity profiles that concorded with histological studies, indicating that MRF residuals potentially contain tissue microstructural information. MRF residual signals provide additional area-specific information that is complementary to the MR relaxometry-based (T₁, T₂^*) information used previously for distinguishing microstructural differences between human cerebral cortex regions in vivo. The proposed approach led to more accurate identification of structural variations across cortical areas of interest.

Neurobiological underpinnings of rapid white matter plasticity during intensive reading instruction

Diffusion MRI is a powerful tool for imaging brain structure, but it is challenging to discern the biological underpinnings of plasticity inferred from these and other non-invasive MR measurements. Biophysical modeling of the diffusion signal aims to render a more biologically rich image of tissue microstructure, but the application of these models comes with important caveats. A separate approach for gaining biological specificity has been to seek converging evidence from multi-modal datasets. Here we use metrics derived from diffusion kurtosis imaging (DKI) and the white matter tract integrity (WMTI) model along with quantitative MRI measurements of T1 relaxation to characterize changes throughout the white matter during an 8-week, intensive reading intervention (160 total hours of instruction). Behavioral measures, multi-shell diffusion MRI data, and quantitative T1 data were collected at regular intervals during the intervention in a group of 33 children with reading difficulties (7-12 years old), and over the same period in an age-matched non-intervention control group. Throughout the white matter, mean 'extra-axonal' diffusivity was inversely related to intervention time. In contrast, model estimated axonal water fraction (AWF), overall diffusion kurtosis, and T1 relaxation time showed no significant change over the intervention period. Both diffusion and quantitative T1 based metrics were correlated with pre-intervention reading performance, albeit with distinct anatomical distributions. These results are consistent with the view that rapid changes in diffusion properties reflect phenomena other than widespread changes in myelin density. We discuss this result in light of recent work highlighting non-axonal factors in experience-dependent plasticity and learning.

Fetal brain issues in congenital heart disease

Following the improvements in the clinical management of patients with congenital heart disease (CHD) and their increased survival, neurodevelopmental outcome has become an emerging priority in pediatric cardiology. Large-scale efforts have been made to protect the brain during the postnatal, surgical, and postoperative period; however, the presence of brain immaturity and injury at birth suggests in utero and peripartum disturbances. Over the past decade, there has been considerable interest and investigations on fetal brain growth in the setting of CHD. Advancements in fetal brain imaging have identified abnormal brain development in fetuses with CHD from the macrostructural (brain volumes and cortical folding) down to the microstructural (biochemistry and water diffusivity) scale, with more severe forms of CHD showing worse disturbances and brain abnormalities starting as early as the first trimester. Anomalies in common genetic developmental pathways and diminished cerebral substrate delivery secondary to altered cardiovascular physiology are the forefront hypotheses, but other factors such as impaired placental function and maternal psychological stress have surfaced as important contributors to fetal brain immaturity in CHD. The characterization and timing of fetal brain disturbances and their associated mechanisms are important steps for determining preventative prenatal interventions, which may provide a stronger foundation for the developing brain during childhood.

Maternal hyperoxygenation in congenital heart disease

The importance of prenatal diagnosis and fetal intervention has been increasing as a preventative strategy for improving the morbidity and mortality in congenital heart disease (CHD). The advancements in medical imaging technology have greatly enhanced our understanding of disease progression, assessment, and impact in those with CHD. In particular, there has been a growing focus on improving the morbidity and mortality of fetuses diagnosed with left-sided lesions. The disruption of fetal hemodynamics resulting from poor structural developmental of the left outflow tract during cardiogenesis is considered a major factor in the progressive lethal underdevelopment of the left ventricle (LV). This positive feedback cycle of inadequate flow and underdevelopment of the LV leads to a disrupted fetal circulation, which has been described to impact fetal brain growth where systemic outflow is poor and, in some cases, the fetal lungs in the setting of a restrictive interatrial communication. For the past decade, maternal hyperoxygenation (MH) has been investigated as a diagnostic tool to assess the pulmonary vasculature and a therapeutic agent to improve the development of the heart and brain in fetuses with CHD with a focus on left-sided cardiac defects. This review discusses the findings of these studies as well as the utility of acute and chronic administration of MH in CHD.

Gray Matter Vulnerabilities Predict Longitudinal Development of Apathy in Huntington's Disease

BACKGROUND

Apathy, a common neuropsychiatric disturbance in Huntington's disease (HD), is subserved by a complex neurobiological network. However, no study has yet employed a whole-brain approach to examine underlying regional vulnerabilities that may precipitate apathy changes over time.

OBJECTIVES

To identify whole-brain gray matter volume (GMV) vulnerabilities that may predict longitudinal apathy development in HD.

METHODS

Forty-five HD individuals (31 female) were scanned and evaluated for apathy and other neuropsychiatric features using the short-Problem Behavior Assessment for a maximum total of six longitudinal visits (including baseline). In order to identify regions where changes in GMV may describe changes in apathy, we performed longitudinal voxel-based morphometry (VBM) on those 33 participants with a magnetic resonance imaging (MRI) scan on their second visit at 18 ± 6 months follow-up (78 MRI datasets). We next employed a generalized linear mixed-effects model (N = 45) to elucidate whether initial and specific GMV may predict apathy development over time.

RESULTS

Utilizing longitudinal VBM, we revealed a relationship between increases in apathy and specific GMV atrophy in the right middle cingulate cortex (MCC). Furthermore, vulnerability in the right MCC volume at baseline successfully predicted the severity and progression of apathy over time.

CONCLUSIONS

This study highlights that individual differences in apathy in HD may be explained by variability in atrophy and initial vulnerabilities in the right MCC, a region implicated in action-initiation. These findings thus serve to facilitate the prediction of an apathetic profile, permitting targeted, time-sensitive interventions in neurodegenerative disease with potential implications in otherwise healthy populations. © 2021 International Parkinson and Movement Disorder Society.

Advanced Analysis of Diffusion Tensor Imaging Along With Machine Learning Provides New Sensitive Measures of Tissue Pathology and Intra-Lesion Activity in Multiple Sclerosis

Tissue pathology in multiple sclerosis (MS) is highly complex, requiring multi-dimensional analysis. In this study, our goal was to test the feasibility of obtaining high angular resolution diffusion imaging (HARDI) metrics through single-shell modeling of diffusion tensor imaging (DTI) data, and investigate how advanced measures from single-shell HARDI and DTI tractography perform relative to classical DTI metrics in assessing MS pathology. We examined 52 relapsing-remitting MS patients who had 3T anatomical brain MRI and DTI. Single-shell HARDI modeling yielded 5 sub-voxel-based metrics, totalling 11 diffusion measures including 4 DTI and 2 tractography metrics. Based on machine learning of 3-dimensional regions of interest, we evaluated the importance of the measures through several tissue classification tasks. These included two within-subject comparisons: lesion versus normal appearing white matter (NAWM); and lesion core versus shell. Further, by stratifying patients as having high (above 75%^ile) and low (below 25%^ile) number of MS lesions, we also performed 2 classifications between subjects for lesions and NAWM respectively. Results showed that in lesion-NAWM analysis, HARDI orientation distribution function (ODF) energy, DTI fractional anisotropy (FA), and HARDI orientation dispersion index were the top three metrics, which together achieved 65.2% accuracy and 0.71 area under the receiver operating characteristic curve (AUROC). In core-shell analysis, DTI mean diffusivity (MD), radial diffusivity, and FA were the top three metrics, and MD dominated the classification, which achieved 59.3% accuracy and 0.59 AUROC alone. Between patients, FA was the leading feature in lesion comparisons, while ODF energy was the best in NAWM separation. Collectively, single-shell modeling of common diffusion data can provide robust orientation measures of lesion and NAWM pathology, and DTI metrics are most sensitive to intra-lesion abnormality. Combined analysis of both advanced and classical diffusion measures may be critical for improved understanding of MS pathology.

Hybrid diffusion imaging reveals altered white matter tract integrity and associations with symptoms and cognitive dysfunction in chronic traumatic brain injury

•

Hybrid Diffusion Imaging (HYDI) detects white matter associations in patients with cTBI.

•

The advanced diffusion model NODDI was more sensitive in detecting between-group differences than classic DTI.

•

DTI appeared to be just as sensitive as NODDI for detecting white matter correlations with self-reported symptoms.

•

This study highlights the advantages of acquiring both DTI and NODDI to fully characterize white matter microstructure in cTBI.

The detection and association of in vivo biomarkers in white matter (WM) pathology after acute and chronic mild traumatic brain injury (mTBI) are needed to improve care and develop therapies. In this study, we used the diffusion MRI method of hybrid diffusion imaging (HYDI) to detect white matter alterations in patients with chronic TBI (cTBI). 40 patients with cTBI presenting symptoms at least three months post injury, and 17 healthy controls underwent magnetic resonance HYDI. cTBI patients were assessed with a battery of neuropsychological tests. A voxel-wise statistical analysis within the white matter skeleton was performed to study between group differences in the diffusion models. In addition, a partial correlation analysis controlling for age, sex, and time after injury was performed within the cTBI cohort, to test for associations between diffusion metrics and clinical outcomes. The advanced diffusion modeling technique of neurite orientation dispersion and density imaging (NODDI) showed large clusters of between-group differences resulting in lower values in the cTBI across the brain, where the single compartment diffusion tensor model failed to show any significant results. However, the diffusion tensor model appeared to be just as sensitive in detecting self-reported symptoms in the cTBI population using a within-group correlation. To the best of our knowledge this study provides the first application of HYDI in evaluation of cTBI using combined DTI and NODDI, significantly enhancing our understanding of the effects of concussion on white matter microstructure and emphasizing the utility of full characterization of complex diffusion to diagnose, monitor, and treat brain injury.

The Neurobiology of Zika Virus: New Models, New Challenges

The Zika virus (ZIKV) attracted attention due to one striking characteristic: the ability to cross the placental barrier and infect the fetus, possibly causing severe neurodevelopmental disruptions included in the Congenital Zika Syndrome (CZS). Few years after the epidemic, the CZS incidence has begun to decline. However, how ZIKV causes a diversity of outcomes is far from being understood. This is probably driven by a chain of complex events that relies on the interaction between ZIKV and environmental and physiological variables. In this review, we address open questions that might lead to an ill-defined diagnosis of CZS. This inaccuracy underestimates a large spectrum of apparent normocephalic cases that remain underdiagnosed, comprising several subtle brain abnormalities frequently masked by a normal head circumference. Therefore, new models using neuroimaging and artificial intelligence are needed to improve our understanding of the neurobiology of ZIKV and its true impact in neurodevelopment.

Lesion distribution and substrate of white matter damage in myotonic dystrophy type 1: Comparison with multiple sclerosis

•

The supratentorial distribution of lesions is similar in DM1 and MS.

•

Patients with DM1 do not show infratentorial lesions.

•

Quantitative magnetization transfer supports the presence of demyelination in DM1 lesions, but not in the NAWM.

•

Anterior temporal lobe lesions in DM1 might have a different substrate than periventricular ones.

Myotonic Dystrophy type 1 (DM1) is an autosomal dominant condition caused by expansion of the CTG triplet repeats within the myotonic dystrophy protein of the kinase (DMPK) gene. The central nervous system is involved in the disease, with multiple symptoms including cognitive impairment. A typical feature of DM1 is the presence of widespread white matter (WM) lesions, whose total volume is associated with CTG triplet expansion. The aim of this study was to characterize the distribution and pathological substrate of these lesions as well as the normal appearing WM (NAWM) using quantitative magnetization transfer (qMT) MRI, and comparing data from DM1 patients with those from patients with multiple sclerosis (MS). Twenty-eight patients with DM1, 29 patients with relapsing-remitting MS, and 15 healthy controls had an MRI scan, including conventional and qMT imaging. The average pool size ratio (F), a proxy of myelination, was computed within lesions and NAWM for every participant. The lesion masks were warped into MNI space and lesion probability maps were obtained for each patient group. The lesion distribution, total lesion load and the tissue-specific mean F were compared between groups. The supratentorial distribution of lesions was similar in the 2 patient groups, although mean lesion volume was higher in MS than DM1. DM1 presented higher prevalence of anterior temporal lobe lesions, but none in the cerebellum and brainstem. Significantly reduced F values were found within DM1 lesions, suggesting a loss of myelin density. While F was reduced in the NAWM of MS patients, it did not differ between DM1 and controls. Our results provide further evidence for a need to compare histology and imaging using new MRI techniques in DM1 patients, in order to further our understanding of the underlying disease process contributing to WM disease.

Hypomyelinating leukodystrophies - unravelling myelin biology

Hypomyelinating leukodystrophies constitute a subset of genetic white matter disorders characterized by a primary lack of myelin deposition. Most patients with severe hypomyelination present in infancy or early childhood and develop severe neurological deficits, but the clinical presentation can also be mild with onset of symptoms in adolescence or adulthood. MRI can be used to visualize the process of myelination in detail, and MRI pattern recognition can provide a clinical diagnosis in many patients. Next-generation sequencing provides a definitive diagnosis in 80-90% of patients. Genes associated with hypomyelination include those that encode structural myelin proteins but also many that encode proteins involved in RNA translation and some lysosomal proteins. The precise pathomechanisms remain to be elucidated. Improved understanding of the process of myelination, the metabolic axonal support functions of myelin and the proposed contribution of myelin to CNS plasticity provide possible explanations as to why almost all patients with hypomyelination experience slow clinical decline after a long phase of stability. In this Review, we provide an overview of the hypomyelinating leukodystrophies, the advances in our understanding of myelin biology and of the genes involved in these disorders, and the insights these advances have provided into their clinical presentations and evolution.

Variability and Standardization of Quantitative Imaging: Monoparametric to Multiparametric Quantification, Radiomics, and Artificial Intelligence

Radiological images have been assessed qualitatively in most clinical settings by the expert eyes of radiologists and other clinicians. On the other hand, quantification of radiological images has the potential to detect early disease that may be difficult to detect with human eyes, complement or replace biopsy, and provide clear differentiation of disease stage. Further, objective assessment by quantification is a prerequisite of personalized/precision medicine. This review article aims to summarize and discuss how the variability of quantitative values derived from radiological images are induced by a number of factors and how these variabilities are mitigated and standardization of the quantitative values are achieved. We discuss the variabilities of specific biomarkers derived from magnetic resonance imaging and computed tomography, and focus on diffusion-weighted imaging, relaxometry, lung density evaluation, and computer-aided computed tomography volumetry. We also review the sources of variability and current efforts of standardization of the rapidly evolving techniques, which include radiomics and artificial intelligence.

NODDI in clinical research

Diffusion MRI (dMRI) has proven to be a useful imaging approach for both clinical diagnosis and research investigating the microstructures of nervous tissues, and it has helped us to better understand the neurophysiological mechanisms of many diseases. Though diffusion tensor imaging (DTI) has long been the default tool to analyze dMRI data in clinical research, acquisition with stronger diffusion weightings beyond the DTI regimen is now possible with modern clinical scanners, potentially enabling even more detailed characterization of tissue microstructures. To take advantage of such data, neurite orientation dispersion and density imaging (NODDI) has been proposed as a way to relate the dMRI signal to tissue features via biophysically inspired modeling. The number of reports demonstrating the potential clinical utility of NODDI is rapidly increasing. At the same time, the pitfalls and limitations of NODDI, and general challenges in microstructure modeling, are becoming increasingly recognized by clinicians. dMRI microstructure modeling is a rapidly evolving field with great promise, where people from different scientific backgrounds, such as physics, medicine, biology, neuroscience, and statistics, are collaborating to build novel tools that contribute to improving human healthcare. Here, we review the applications of NODDI in clinical research and discuss future perspectives for investigations toward the implementation of dMRI microstructure imaging in clinical practice.

Multi-parametric quantitative in vivo spinal cord MRI with unified signal readout and image denoising

Multi-parametric quantitative MRI (qMRI) of the spinal cord is a promising non-invasive tool to probe early microstructural damage in neurological disorders. It is usually performed in vivo by combining acquisitions with multiple signal readouts, which exhibit different thermal noise levels, geometrical distortions and susceptibility to physiological noise. This ultimately hinders joint multi-contrast modelling and makes the geometric correspondence of parametric maps challenging. We propose an approach to overcome these limitations, by implementing state-of-the-art microstructural MRI of the spinal cord with a unified signal readout in vivo (i.e. with matched spatial encoding parameters across a range of imaging contrasts). We base our acquisition on single-shot echo planar imaging with reduced field-of-view, and obtain data from two different vendors (vendor 1: Philips Achieva; vendor 2: Siemens Prisma). Importantly, the unified acquisition allows us to compare signal and noise across contrasts, thus enabling overall quality enhancement via multi-contrast image denoising methods. As a proof-of-concept, here we provide a demonstration with one such method, known as Marchenko-Pastur (MP) Principal Component Analysis (PCA) denoising. MP-PCA is a singular value (SV) decomposition truncation approach that relies on redundant acquisitions, i.e. such that the number of measurements is large compared to the number of components that are maintained in the truncated SV decomposition. Here we used in vivo and synthetic data to test whether a unified readout enables more efficient MP-PCA denoising of less redundant acquisitions, since these can be denoised jointly with more redundant ones. We demonstrate that a unified readout provides robust multi-parametric maps, including diffusion and kurtosis tensors from diffusion MRI, myelin metrics from two-pool magnetisation transfer, and T1 and T2 from relaxometry. Moreover, we show that MP-PCA improves the quality of our multi-contrast acquisitions, since it reduces the coefficient of variation (i.e. variability) by up to 17% for mean kurtosis, 8% for bound pool fraction (myelin-sensitive), and 13% for T1, while enabling more efficient denoising of modalities limited in redundancy (e.g. relaxometry). In conclusion, multi-parametric spinal cord qMRI with unified readout is feasible and provides robust microstructural metrics with matched resolution and distortions, whose quality benefits from multi-contrast denoising methods such as MP-PCA.

The evolution of white matter microstructural changes after mild traumatic brain injury: A longitudinal DTI and NODDI study

Neuroimaging biomarkers that can detect white matter (WM) pathology after mild traumatic brain injury (mTBI) and predict long-term outcome are needed to improve care and develop therapies. We used diffusion tensor imaging (DTI) and neurite orientation dispersion and density imaging (NODDI) to investigate WM microstructure cross-sectionally and longitudinally after mTBI and correlate these with neuropsychological performance. Cross-sectionally, early decreases of fractional anisotropy and increases of mean diffusivity corresponded to WM regions with elevated free water fraction on NODDI. This elevated free water was more extensive in the patient subgroup reporting more early postconcussive symptoms. The longer-term longitudinal WM changes consisted of declining neurite density on NODDI, suggesting axonal degeneration from diffuse axonal injury for which NODDI is more sensitive than DTI. Therefore, NODDI is a more sensitive and specific biomarker than DTI for WM microstructural changes due to mTBI that merits further study for mTBI diagnosis, prognosis, and treatment monitoring.

The effect of network thresholding and weighting on structural brain networks in the UK Biobank

Whole-brain structural networks can be constructed using diffusion MRI and probabilistic tractography. However, measurement noise and the probabilistic nature of the tracking procedure result in an unknown proportion of spurious white matter connections. Faithful disentanglement of spurious and genuine connections is hindered by a lack of comprehensive anatomical information at the network-level. Therefore, network thresholding methods are widely used to remove ostensibly false connections, but it is not yet clear how different thresholding strategies affect basic network properties and their associations with meaningful demographic variables, such as age. In a sample of 3153 generally healthy volunteers from the UK Biobank Imaging Study (aged 44-77 years), we constructed whole-brain structural networks and applied two principled network thresholding approaches (consistency and proportional thresholding). These were applied over a broad range of threshold levels across six alternative network weightings (streamline count, fractional anisotropy, mean diffusivity and three novel weightings from neurite orientation dispersion and density imaging) and for four common network measures (mean edge weight, characteristic path length, network efficiency and network clustering coefficient). We compared network measures against age associations and found that: 1) measures derived from unthresholded matrices yielded the weakest age-associations (0.033 ​≤ ​|β| ​≤ ​0.409); and 2) the most commonly-used level of proportional-thresholding from the literature (retaining 68.7% of all possible connections) yielded significantly weaker age-associations (0.070 ​≤ ​|β| ​≤ ​0.406) than the consistency-based approach which retained only 30% of connections (0.140 ​≤ ​|β| ​≤ ​0.409). However, we determined that the stringency of the threshold was a stronger determinant of the network-age association than the choice of threshold method and the two thresholding approaches identified a highly overlapping set of connections (ICC ​= ​0.84), when matched at 70% network sparsity. Generally, more stringent thresholding resulted in more age-sensitive network measures in five of the six network weightings, except at the highest levels of sparsity (>90%), where crucial connections were then removed. At two commonly-used threshold levels, the age-associations of the connections that were discarded (mean β ​≤ ​|0.068|) were significantly smaller in magnitude than the corresponding age-associations of the connections that were retained (mean β ​≤ ​|0.219|, p ​< ​0.001, uncorrected). Given histological evidence of widespread degeneration of structural brain connectivity with increasing age, these results indicate that stringent thresholding methods may be most accurate in identifying true white matter connections.

Histological validation of per-bundle water diffusion metrics within a region of fiber crossing following axonal degeneration

Micro-architectural characteristics of white matter can be inferred through analysis of diffusion-weighted magnetic resonance imaging (dMRI). The diffusion-dependent signal can be analyzed through several methods, with the tensor model being the most frequently used due to its straightforward interpretation and low requirements for acquisition parameters. While valuable information can be gained from the tensor-derived metrics in regions of homogeneous tissue organization, this model does not provide reliable microstructural information at crossing fiber regions, which are pervasive throughout human white matter. Several multiple fiber models have been proposed that seem to overcome the limitations of the tensor, with few providing per-bundle dMRI-derived metrics. However, biological interpretations of such metrics are limited by the lack of histological confirmation. To this end, we developed a straightforward biological validation framework. Unilateral retinal ischemia was induced in ten rats, which resulted in axonal (Wallerian) degeneration of the corresponding optic nerve, while the contralateral was left intact; the intact and injured axonal populations meet at the optic chiasm as they cross the midline, generating a fiber crossing region in which each population has different diffusion properties. Five rats served as controls. High-resolution ex vivo dMRI was acquired five weeks after experimental procedures. We correlated and compared histology to per-bundle descriptors derived from three methodologies for dMRI analysis (constrained spherical deconvolution and two multi-tensor representations). We found a tight correlation between axonal density (as evaluated through automatic segmentation of histological sections) with per-bundle apparent fiber density and fractional anisotropy (derived from dMRI). The multi-fiber methods explored were able to correctly identify the damaged fiber populations in a region of fiber crossings (chiasm). Our results provide validation of metrics that bring substantial and clinically useful information about white-matter tissue at crossing fiber regions. Our proposed framework is useful to validate other current and future dMRI methods.

Dimensionality reduction of diffusion MRI measures for improved tractometry of the human brain

Various diffusion MRI (dMRI) measures have been proposed for characterising tissue microstructure over the last 15 years. Despite the growing number of experiments using different dMRI measures in assessments of white matter, there has been limited work on: 1) examining their covariance along specific pathways; and on 2) combining these different measures to study tissue microstructure. Indeed, it quickly becomes intractable for existing analysis pipelines to process multiple measurements at each voxel and at each vertex forming a streamline, highlighting the need for new ways to visualise or analyse such high-dimensional data. In a sample of 36 typically developing children aged 8–18 years, we profiled various commonly used dMRI measures across 22 brain pathways. Using a data-reduction approach, we identified two biologically-interpretable components that capture 80% of the variance in these dMRI measures. The first derived component captures properties related to hindrance and restriction in tissue microstructure, while the second component reflects characteristics related to tissue complexity and orientational dispersion. We then demonstrate that the components generated by this approach preserve the biological relevance of the original measurements by showing age-related effects across developmentally sensitive pathways. In summary, our findings demonstrate that dMRI analyses can benefit from dimensionality reduction techniques, to help disentangling the neurobiological underpinnings of white matter organisation.

Iron, Myelin, and the Brain: Neuroimaging Meets Neurobiology

Although iron is crucial for neuronal functioning, many aspects of cerebral iron biology await clarification. The ability to quantify specific iron forms in the living brain would open new avenues for diagnosis, therapeutic monitoring, and understanding pathogenesis of diseases. A modality that allows assessment of brain tissue composition in vivo, in particular of iron deposits or myelin content on a submillimeter spatial scale, is magnetic resonance imaging (MRI). Multimodal strategies combining MRI with complementary analytical techniques ex vivo have emerged, which may lead to improved specificity. Interdisciplinary collaborations will be key to advance beyond simple correlative analyses in the biological interpretation of MRI data and to gain deeper insights into key factors leading to iron accumulation and/or redistribution associated with neurodegeneration.

NODDI-DTI: Estimating Neurite Orientation and Dispersion Parameters from a Diffusion Tensor in Healthy White Matter

The NODDI-DTI signal model is a modification of the NODDI signal model that formally allows interpretation of standard single-shell DTI data in terms of biophysical parameters in healthy human white matter (WM). The NODDI-DTI signal model contains no CSF compartment, restricting application to voxels without CSF partial-volume contamination. This modification allowed derivation of analytical relations between parameters representing axon density and dispersion, and DTI invariants (MD and FA) from the NODDI-DTI signal model. These relations formally allow extraction of biophysical parameters from DTI data. NODDI-DTI parameters were estimated by applying the proposed analytical relations to DTI parameters estimated from the first shell of data, and compared to parameters estimated by fitting the NODDI-DTI model to both shells of data (reference dataset) in the WM of 14 in vivo diffusion datasets recorded with two different protocols, and in simulated data. The first two datasets were also fit to the NODDI-DTI model using only the first shell (as for DTI) of data. NODDI-DTI parameters estimated from DTI, and NODDI-DTI parameters estimated by fitting the model to the first shell of data gave similar errors compared to two-shell NODDI-DTI estimates. The simulations showed the NODDI-DTI method to be more noise-robust than the two-shell fitting procedure. The NODDI-DTI method gave unphysical parameter estimates in a small percentage of voxels, reflecting voxelwise DTI estimation error or NODDI-DTI model invalidity. In the course of evaluating the NODDI-DTI model, it was found that diffusional kurtosis strongly biased DTI-based MD values, and so, making assumptions based on healthy WM, a novel heuristic correction requiring only DTI data was derived and used to mitigate this bias. Since validations were only performed on healthy WM, application to grey matter or pathological WM would require further validation. Our results demonstrate NODDI-DTI to be a promising model and technique to interpret restricted datasets acquired for DTI analysis in healthy white matter with greater biophysical specificity, though its limitations must be borne in mind.