Optimization and validation of methods for mapping of the radiofrequency transmit field at 3T

Enduring differential patterns of neuronal loss and myelination along 6-month pulsatile gonadotropin-releasing hormone therapy in individuals with Down syndrome

Despite major progress in understanding the impact of the triplicated chromosome 21 on the brain and behaviour in Down syndrome, our knowledge of the underlying neurobiology in humans is still limited. We sought to address some of the pertinent questions about the drivers of brain structure differences and their associations with cognitive function in Down syndrome. To this aim, in a pilot magnetic resonance imaging (MRI) study, we monitored brain anatomy in individuals with Down syndrome receiving pulsatile gonadotropin-releasing hormone (GnRH) therapy over 6 months in comparison with typically developed age- and sex-matched healthy controls. We analysed cross-sectional (Down syndrome/healthy controls n  = 11/27; Down syndrome-2 females/9 males, age 26.7 ± 5.0 years old; healthy controls-8 females/19 males, age 24.1 ± 2.5 years old) and longitudinal (Down syndrome/healthy controls n  = 8/13; Down syndrome-1 female/7 males, age 26.4 ± 5.3 years old; healthy controls-4 females/9 males, 24.7 ± 2.2 years old) relaxometry and diffusion-weighted MRI data alongside standard cognitive assessment. The statistical tests looked for cross-sectional baseline differences and for differential changes over time between Down syndrome and healthy controls. The post hoc analysis confined to the Down syndrome group, tested for potential time-dependent interactions between individuals' overall cognitive performance and associated brain anatomy changes. The brain MRI statistical analyses covered both grey and white matter regions across the whole brain allowing for investigation of regional volume, macromolecular/myelin and iron content, additionally to diffusion tensor and neurite orientation and dispersion density characterization across major white matter tracts. The cross-sectional analysis showed reduced frontal, temporal and cerebellar volumes in Down syndrome with only the cerebellar differences remaining significant after adjustment for the presence of microcephaly (P family-wise-corrected < 0.05). The volume reductions were paralleled by decreased cortical and subcortical macromolecular/myelin content confined to the cortical motor system, thalamus and basal ganglia (P family-wise-corrected < 0.05). All major white matter tracts showed a ubiquitous mean diffusivity and intracellular volume fraction reduction contrasted with no differences in magnetization transfer saturation metrics (P family-wise-corrected < 0.05). Compared with healthy controls over the same period, Down syndrome individuals under GnRH therapy showed cognitive improvement (Montreal Cognitive Assessment from 11.4 ± 5.5 to 15.1 ± 5.6; P < 0.01) on the background of stability of the observed differential neuroanatomical patterns. Despite the lack of adequate Down syndrome control group, we interpret the obtained cross-sectional and longitudinal findings in young adults as evidence for predominant neurodevelopmental neuronal loss due to dysfunctional neurogenesis without signs for short-term myelin loss.

Longitudinal Changes of Quantitative Brain Tissue Properties Induced by Balance Training

The human brain can show remarkable experience-induced plasticity under conditions such as aging and pathology. However, the mapping of changes provided by many imaging approaches often lacks specificity with respect to biological tissue properties, which is relevant for treatment optimization and the evaluation of health-promoting lifestyle factors. Training-induced structural changes in cortical and subcortical gray matter likely reflect a mixture of various microstructural processes. In order to non-invasively map these different microstructural contributions, we used quantitative magnetic resonance imaging (qMRI) to measure clinically-relevant brain tissue property changes (such as iron, myelin, and water) in response to 4 weeks of motor balance training in 26 healthy young adults. Training resulted in a regionally-specific decrease in myelin-related magnetization transfer saturation (MTsat) in the left frontal cortex. We also found performance-related changes in iron-sensitive transverse relaxation rate (R2*) in visual cortical (signal increase along with positive performance correlation) and limbic subcortical (signal decrease along with negative performance correlation) brain areas. Our study contributes to a growing body of literature investigating motor training-induced microstructural brain plasticity. Specifically, we provide new insights into microstructural brain changes using whole-body motor learning (balance practice) and longitudinal quantitative mapping of brain tissue properties.

Circadian rapid eye movement sleep expression is associated with brain microstructural integrity in older adults

Rapid eye movement sleep (REMS) is increasingly suggested as a discriminant sleep state for subtle signs of age-related neurodegeneration. While REMS expression is under strong circadian control and circadian dysregulation increases with age, the association between brain aging and circadian REMS regulation has not yet been assessed. Here, we measure the circadian amplitude of REMS through a 40-h in-lab multiple nap protocol in controlled laboratory conditions, and brain microstructural integrity with quantitative multi-parameter mapping (MPM) imaging in 86 older individuals. We show that reduced circadian REMS amplitude is related to lower magnetization transfer saturation (MTsat), longitudinal relaxation rate (R1) and effective transverse relaxation rate (R2*) values in several white matter regions mostly located around the lateral ventricles, and with lower R1 values in grey matter clusters encompassing the hippocampus, parahippocampus, thalamus and hypothalamus. Our results further highlight the importance of considering circadian regulation for understanding the association between sleep and brain structure in older individuals.

The mouse motor system contains multiple premotor areas and partially follows human organizational principles

While humans are known to have several premotor cortical areas, secondary motor cortex (M2) is often considered to be the only higher-order motor area of the mouse brain and is thought to combine properties of various human premotor cortices. Here, we show that axonal tracer, functional connectivity, myelin mapping, gene expression, and optogenetics data contradict this notion. Our analyses reveal three premotor areas in the mouse, anterior-lateral motor cortex (ALM), anterior-lateral M2 (aM2), and posterior-medial M2 (pM2), with distinct structural, functional, and behavioral properties. By using the same techniques across mice and humans, we show that ALM has strikingly similar functional and microstructural properties to human anterior ventral premotor areas and that aM2 and pM2 amalgamate properties of human pre-SMA and cingulate cortex. These results provide evidence for the existence of multiple premotor areas in the mouse and chart a comparative map between the motor systems of humans and mice.

The Neurobiology of Life Course Socioeconomic Conditions and Associated Cognitive Performance in Middle to Late Adulthood

Despite major advances, our understanding of the neurobiology of life course socioeconomic conditions is still scarce. This study aimed to provide insight into the pathways linking socioeconomic exposures-household income, last known occupational position, and life course socioeconomic trajectories-with brain microstructure and cognitive performance in middle to late adulthood. We assessed socioeconomic conditions alongside quantitative relaxometry and diffusion-weighted magnetic resonance imaging indicators of brain tissue microstructure and cognitive performance in a sample of community-dwelling men and women (N = 751, aged 50-91 years). We adjusted the applied regression analyses and structural equation models for the linear and nonlinear effects of age, sex, education, cardiovascular risk factors, and the presence of depression, anxiety, and substance use disorders. Individuals from lower-income households showed signs of advanced brain white matter (WM) aging with greater mean diffusivity (MD), lower neurite density, lower myelination, and lower iron content. The association between household income and MD was mediated by neurite density (B = 0.084, p = 0.003) and myelination (B = 0.019, p = 0.009); MD partially mediated the association between household income and cognitive performance (B = 0.017, p < 0.05). Household income moderated the relation between WM microstructure and cognitive performance, such that greater MD, lower myelination, or lower neurite density was only associated with poorer cognitive performance among individuals from lower-income households. Individuals from higher-income households showed preserved cognitive performance even with greater MD, lower myelination, or lower neurite density. These findings provide novel mechanistic insights into the associations between socioeconomic conditions, brain anatomy, and cognitive performance in middle to late adulthood.

The effect of and correction for through-slice dephasing on 2D gradient-echo double angle B 1 + mapping

PURPOSE

To show thatB 0 $$ {\mathrm{B}}_0 $$ variations through slice and slice profile effects are two major confounders affecting 2D dual angleB 1 + $$ {\mathrm{B}}_1^{+} $$ maps using gradient-echo signals and thus need to be corrected to obtain accurateB 1 + $$ {\mathrm{B}}_1^{+} $$ maps.

METHODS

The 2D gradient-echo transverse complex signal was Bloch-simulated and integrated across the slice dimension including nonlinear variations inB 0 $$ {\mathrm{B}}_0 $$ inhomogeneities through slice. A nonlinear least squares fit was used to find theB 1 + $$ {\mathrm{B}}_1^{+} $$ factor corresponding to the best match between the two gradient-echo signals experimental ratio and the Bloch-simulated ratio. The correction was validated in phantom and in vivo at 3T.

RESULTS

For our RF excitation pulse, the error in theB 1 + $$ {\mathrm{B}}_1^{+} $$ factor scales by approximately 3.8% for every 10 Hz/cm variation inB 0 $$ {\mathrm{B}}_0 $$ along the slice direction. Higher accuracy phantomB 1 + $$ {\mathrm{B}}_1^{+} $$ maps were obtained after applying the proposed correction; the root mean squareB 1 + $$ {\mathrm{B}}_1^{+} $$ error relative to the gold standardB 1 + $$ {\mathrm{B}}_1^{+} $$ decreased from 6.4% to 2.6%. In vivo whole-liverT 1 $$ {\mathrm{T}}_1 $$ maps using the correctedB 1 + $$ {\mathrm{B}}_1^{+} $$ map registered a significant decrease inT 1 $$ {\mathrm{T}}_1 $$ gradient through slice.

CONCLUSION

B 0 $$ {\mathrm{B}}_0 $$ inhomogeneities varying through slice were seen to have an impact on the accuracy of 2D double angleB 1 + $$ {\mathrm{B}}_1^{+} $$ maps using gradient-echo sequences. Consideration of this confounder is crucial for research relying on accurate knowledge of the true excitation flip angles, as is the case ofT 1 $$ {\mathrm{T}}_1 $$ mapping using a spoiled gradient recalled echo sequence.

Phosphorus-31: A table-top method for 3D B1-field amplitude and phase measurements

A novel method of high-spatial-resolution, 3D B1-field distribution measurements is presented. The method is independent of the MR-scanner, and it allows for automated acquisitions of complete maps of all magnetic field vector components for both proton and heteronuclear MR coils of arbitrary geometrical shapes. The advantage of the method proposed here, compared with methods based on measurements with an MR-scanner, is that a complete image of both receive and transmit B1-fields, including the phase of the B1-field, can be acquired. The B1 field maps obtained in this manner can be used for absolute quantification of metabolites in MRS experiments, as well as for intensity compensations in imaging experiments, both of which are important concepts in biological and medical MR applications. Another use might be in coil development and testing. A comparison with B1 field magnitude maps obtained with an MR-scanner was included to validate the accuracy of the proposed method.

Brain structure and function: a multidisciplinary pipeline to study hominoid brain evolution

To decipher the evolution of the hominoid brain and its functions, it is essential to conduct comparative studies in primates, including our closest living relatives. However, strong ethical concerns preclude in vivo neuroimaging of great apes. We propose a responsible and multidisciplinary alternative approach that links behavior to brain anatomy in non-human primates from diverse ecological backgrounds. The brains of primates observed in the wild or in captivity are extracted and fixed shortly after natural death, and then studied using advanced MRI neuroimaging and histology to reveal macro- and microstructures. By linking detailed neuroanatomy with observed behavior within and across primate species, our approach provides new perspectives on brain evolution. Combined with endocranial brain imprints extracted from computed tomographic scans of the skulls these data provide a framework for decoding evolutionary changes in hominin fossils. This approach is poised to become a key resource for investigating the evolution and functional differentiation of hominoid brains.

Multimodal imaging of microstructural cerebral alterations and loss of synaptic density in Alzheimer's disease

Multiple neuropathological events are involved in Alzheimer's disease (AD). The current study investigated the concurrence of neurodegeneration, increased iron content, atrophy, and demyelination in AD. Quantitative multiparameter magnetic resonance imaging (MRI) maps providing neuroimaging biomarkers for myelination and iron content along with synaptic density measurements using [18F] UCB-H PET were acquired in 24 AD and 19 Healthy controls (19 males). The whole brain voxel-wise group comparison revealed demyelination in the right hippocampus, while no significant iron content difference was detected. Bilateral atrophy and synaptic density loss were observed in the hippocampus and amygdala. The multivariate GLM (mGLM) analysis shows a bilateral difference in the hippocampus and amygdala, right pallidum, left fusiform and temporal lobe suggesting that these regions are the most affected despite the diverse differences in brain tissue properties in AD. Demyelination was identified as the most affecting factor in the observed differences. Here, the mGLM is introduced as an alternative for multiple comparisons between different modalities, reducing the risk of false positives while informing about the co-occurrence of neuropathological processes in AD.

B1 mapping using pre-learned subspaces for quantitative brain imaging

PURPOSE

To develop a machine learning-based method for estimation of both transmitter and receiver B1 fields desired for correction of the B1 inhomogeneity effects in quantitative brain imaging.

THEORY AND METHODS

A subspace model-based machine learning method was proposed for estimation of B1t and B1r fields. Probabilistic subspace models were used to capture scan-dependent variations in the B1 fields; the subspace basis and coefficient distributions were learned from pre-scanned training data. Estimation of the B1 fields for new experimental data was achieved by solving a linear optimization problem with prior distribution constraints. We evaluated the performance of the proposed method for B1 inhomogeneity correction in quantitative brain imaging scenarios, including T1 and proton density (PD) mapping from variable-flip-angle spoiled gradient-echo (SPGR) data as well as neurometabolic mapping from MRSI data, using phantom, healthy subject and brain tumor patient data.

RESULTS

In both phantom and healthy subject data, the proposed method produced high-quality B1 maps. B1 correction on SPGR data using the estimated B1 maps produced significantly improved T1 and PD maps. In brain tumor patients, the proposed method produced more accurate and robust B1 estimation and correction results than conventional methods. The B1 maps were also applied to MRSI data from tumor patients and produced improved neurometabolite maps, with better separation between pathological and normal tissues.

CONCLUSION

This work presents a novel method to estimate B1 variations using probabilistic subspace models and machine learning. The proposed method may make correction of B1 inhomogeneity effects more robust in practical applications.

Assessment of white matter hyperintensity severity using multimodal magnetic resonance imaging

White matter hyperintensities are radiological abnormalities reflecting cerebrovascular dysfunction detectable using MRI. White matter hyperintensities are often present in individuals at the later stages of the lifespan and in prodromal stages in the Alzheimer's disease spectrum. Tissue alterations underlying white matter hyperintensities may include demyelination, inflammation and oedema, but these are highly variable by neuroanatomical location and between individuals. There is a crucial need to characterize these white matter hyperintensity tissue alterations in vivo to improve prognosis and, potentially, treatment outcomes. How different MRI measure(s) of tissue microstructure capture clinically-relevant white matter hyperintensity tissue damage is currently unknown. Here, we compared six MRI signal measures sampled within white matter hyperintensities and their associations with multiple clinically-relevant outcomes, consisting of global and cortical brain morphometry, cognitive function, diagnostic and demographic differences and cardiovascular risk factors. We used cross-sectional data from 118 participants: healthy controls (n = 30), individuals at high risk for Alzheimer's disease due to familial history (n = 47), mild cognitive impairment (n = 32) and clinical Alzheimer's disease dementia (n = 9). We sampled the median signal within white matter hyperintensities on weighted MRI images [T1-weighted (T1w), T2-weighted (T2w), T1w/T2w ratio, fluid-attenuated inversion recovery (FLAIR)] as well as the relaxation times from quantitative T1 (qT1) and T2* (qT2*) images. qT2* and fluid-attenuated inversion recovery signals within white matter hyperintensities displayed different age- and disease-related trends compared to normal-appearing white matter signals, suggesting sensitivity to white matter hyperintensity-specific tissue deterioration. Further, white matter hyperintensity qT2*, particularly in periventricular and occipital white matter regions, was consistently associated with all types of clinically-relevant outcomes in both univariate and multivariate analyses and across two parcellation schemes. qT1 and fluid-attenuated inversion recovery measures showed consistent clinical relationships in multivariate but not univariate analyses, while T1w, T2w and T1w/T2w ratio measures were not consistently associated with clinical variables. We observed that the qT2* signal was sensitive to clinically-relevant microstructural tissue alterations specific to white matter hyperintensities. Our results suggest that combining volumetric and signal measures of white matter hyperintensity should be considered to fully characterize the severity of white matter hyperintensities in vivo. These findings may have implications in determining the reversibility of white matter hyperintensities and the potential efficacy of cardio- and cerebrovascular treatments.

Statistical analyses of motion-corrupted MRI relaxometry data computed from multiple scans

BACKGROUND

Consistent noise variance across data points (i.e. homoscedasticity) is required to ensure the validity of statistical analyses of MRI data conducted using linear regression methods. However, head motion leads to degradation of image quality, introducing noise heteroscedasticity into ordinary-least square analyses.

NEW METHOD

The recently introduced QUIQI method restores noise homoscedasticity by means of weighted least square analyses in which the weights, specific for each dataset of an analysis, are computed from an index of motion-induced image quality degradation. QUIQI was first demonstrated in the context of brain maps of the MRI parameter R2 * , which were computed from a single set of images with variable echo time. Here, we extend this framework to quantitative maps of the MRI parameters R1, R2 * , and MTsat, computed from multiple sets of images.

RESULTS

QUIQI restores homoscedasticity in analyses of quantitative MRI data computed from multiple scans. QUIQI allows for optimization of the noise model by using metrics quantifying heteroscedasticity and free energy.

COMPARISON WITH EXISTING METHODS

QUIQI restores homoscedasticity more effectively than insertion of an image quality index in the analysis design and yields higher sensitivity than simply removing the datasets most corrupted by head motion from the analysis.

CONCLUSION

QUIQI provides an optimal approach to group-wise analyses of a range of quantitative MRI parameter maps that is robust to inherent homoscedasticity.

Release of cognitive and multimodal MRI data including real-world tasks and hippocampal subfield segmentations

We share data from N = 217 healthy adults (mean age 29 years, range 20-41; 109 females, 108 males) who underwent extensive cognitive assessment and neuroimaging to examine the neural basis of individual differences, with a particular focus on a brain structure called the hippocampus. Cognitive data were collected using a wide array of questionnaires, naturalistic tests that examined imagination, autobiographical memory recall and spatial navigation, traditional laboratory-based tests such as recalling word pairs, and comprehensive characterisation of the strategies used to perform the cognitive tests. 3 Tesla MRI data were also acquired and include multi-parameter mapping to examine tissue microstructure, diffusion-weighted MRI, T2-weighted high-resolution partial volume structural MRI scans (with the masks of hippocampal subfields manually segmented from these scans), whole brain resting state functional MRI scans and partial volume high resolution resting state functional MRI scans. This rich dataset will be of value to cognitive and clinical neuroscientists researching individual differences, real-world cognition, brain-behaviour associations, hippocampal subfields and more. All data are freely available on Dryad.

Emerging Perspectives on MRI Application in Multiple Sclerosis: Moving from Pathophysiology to Clinical Practice

MRI plays a central role in the diagnosis of multiple sclerosis (MS) and in the monitoring of disease course and treatment response. Advanced MRI techniques have shed light on MS biology and facilitated the search for neuroimaging markers that may be applicable in clinical practice. MRI has led to improvements in the accuracy of MS diagnosis and a deeper understanding of disease progression. This has also resulted in a plethora of potential MRI markers, the importance and validity of which remain to be proven. Here, five recent emerging perspectives arising from the use of MRI in MS, from pathophysiology to clinical application, will be discussed. These are the feasibility of noninvasive MRI-based approaches to measure glymphatic function and its impairment; T1-weighted to T2-weighted intensity ratio to quantify myelin content; classification of MS phenotypes based on their MRI features rather than on their clinical features; clinical relevance of gray matter atrophy versus white matter atrophy; and time-varying versus static resting-state functional connectivity in evaluating brain functional organization. These topics are critically discussed, which may guide future applications in the field.

Parkinson's disease may disrupt overlapping subthalamic nucleus and pallidal motor networks

There is an ongoing debate about differential clinical outcome and associated adverse effects of deep brain stimulation (DBS) in Parkinson's disease (PD) targeting the subthalamic nucleus (STN) or the globus pallidus pars interna (GPi). Given that functional connectivity profiles suggest beneficial DBS effects within a common network, the empirical evidence about the underlying anatomical circuitry is still scarce. Therefore, we investigate the STN and GPi-associated structural covariance brain patterns in PD patients and healthy controls. We estimate GPi's and STN's whole-brain structural covariance from magnetic resonance imaging (MRI) in a normative mid- to old-age community-dwelling cohort (n = 1184) across maps of grey matter volume, magnetization transfer (MT) saturation, longitudinal relaxation rate (R1), effective transversal relaxation rate (R2*) and effective proton density (PD*). We compare these with the structural covariance estimates in patients with idiopathic PD (n = 32) followed by validation using a reduced size controls' cohort (n = 32). In the normative data set, we observed overlapping spatially distributed cortical and subcortical covariance patterns across maps confined to basal ganglia, thalamus, motor, and premotor cortical areas. Only the subcortical and midline motor cortical areas were confirmed in the reduced size cohort. These findings contrasted with the absence of structural covariance with cortical areas in the PD cohort. We interpret with caution the differential covariance maps of overlapping STN and GPi networks in patients with PD and healthy controls as correlates of motor network disruption. Our study provides face validity to the proposed extension of the currently existing structural covariance methods based on morphometry features to multiparameter MRI sensitive to brain tissue microstructure.

Topography of associations between cardiovascular risk factors and myelin loss in the ageing human brain

Our knowledge of the mechanisms underlying the vulnerability of the brain's white matter microstructure to cardiovascular risk factors (CVRFs) is still limited. We used a quantitative magnetic resonance imaging (MRI) protocol in a single centre setting to investigate the cross-sectional association between CVRFs and brain tissue properties of white matter tracts in a large community-dwelling cohort (n = 1104, age range 46-87 years). Arterial hypertension was associated with lower myelin and axonal density MRI indices, paralleled by higher extracellular water content. Obesity showed similar associations, though with myelin difference only in male participants. Associations between CVRFs and white matter microstructure were observed predominantly in limbic and prefrontal tracts. Additional genetic, lifestyle and psychiatric factors did not modulate these results, but moderate-to-vigorous physical activity was linked to higher myelin content independently of CVRFs. Our findings complement previously described CVRF-related changes in brain water diffusion properties pointing towards myelin loss and neuroinflammation rather than neurodegeneration.

B 1 + $$ {B}_1^{+} $$ -correction of magnetization transfer saturation maps optimized for 7T postmortem MRI of the brain

PURPOSE

Magnetization transfer saturation ( MTsat $$ \mathrm{MTsat} $$ ) is a useful marker to probe tissue macromolecular content and myelination in the brain. The increased B 1 + $$ {B}_1^{+} $$ -inhomogeneity at ≥ 7 $$ \ge 7 $$  T and significantly larger saturation pulse flip angles which are often used for postmortem studies exceed the limits where previous MTsat $$ \mathrm{MTsat} $$ B 1 + $$ {B}_1^{+} $$ correction methods are applicable. Here, we develop a calibration-based correction model and procedure, and validate and evaluate it in postmortem 7T data of whole chimpanzee brains.

THEORY

The B 1 + $$ {B}_1^{+} $$ dependence of MTsat $$ \mathrm{MTsat} $$ was investigated by varying the off-resonance saturation pulse flip angle. For the range of saturation pulse flip angles applied in typical experiments on postmortem tissue, the dependence was close to linear. A linear model with a single calibration constant C $$ C $$ is proposed to correct bias in MTsat $$ \mathrm{MTsat} $$ by mapping it to the reference value of the saturation pulse flip angle.

METHODS

C $$ C $$ was estimated voxel-wise in five postmortem chimpanzee brains. "Individual-based global parameters" were obtained by calculating the mean C $$ C $$ within individual specimen brains and "group-based global parameters" by calculating the means of the individual-based global parameters across the five brains.

RESULTS

The linear calibration model described the data well, though C $$ C $$ was not entirely independent of the underlying tissue and B 1 + $$ {B}_1^{+} $$ . Individual-based correction parameters and a group-based global correction parameter ( C = 1 . 2 $$ C=1.2 $$ ) led to visible, quantifiable reductions of B 1 + $$ {B}_1^{+} $$ -biases in high-resolution MTsat $$ \mathrm{MTsat} $$ maps.

CONCLUSION

The presented model and calibration approach effectively corrects for B 1 + $$ {B}_1^{+} $$ inhomogeneities in postmortem 7T data.

B1 mapping using an EPI-based double angle approach: A practical guide for correcting slice profile and B0 distortion effects

PURPOSE

Aim of this study was to develop a reliable B₁ mapping method for brain imaging based on vendor MR sequences available on clinical scanners. Correction procedures for B₀ distortions and slice profile imperfections are proposed, together with a phantom experiment for deriving the approximate time-bandwidth-product (TBP) of the excitation pulse, which is usually not known for vendor sequences.

METHODS

The double angle method was used, acquiring two gradient echo echo-planar imaging data sets with different excitation angles. A correction factor C (B₁ , TBP, B₀ ) was derived from simulations for converting double angle method signal quotients into bias-free B₁ maps. In vitro and in vivo tests compare results with reference B₁ maps based on an established in-house sequence.

RESULTS

The simulation shows that C has a negligible B₁ dependence, allowing for a polynomial approximation of C (TBP, B₀ ). Signal quotients measured in a phantom experiment with known TBP reconfirm the simulation results. In vitro and in vivo B₁ maps based on the proposed method, assuming TBP = 5.8 as derived from a phantom experiment, match closely the reference B₁ maps. Analysis without B₀ correction shows marked deviations in areas of distorted B₀ , highlighting the importance of this correction.

CONCLUSION

Double angle method-based B₁ mapping was set up for vendor gradient echo-echo-planar imaging sequences, using a correction procedure for slice profile imperfections and B₀ distortions. This will help to set up quantitative MRI studies on clinical scanners with release sequences, as the method does not require knowledge of the exact RF-pulse profiles or the use of in-house sequences.

High-resolution quantitative and functional MRI indicate lower myelination of thin and thick stripes in human secondary visual cortex

The characterization of cortical myelination is essential for the study of structure-function relationships in the human brain. However, knowledge about cortical myelination is largely based on post-mortem histology, which generally renders direct comparison to function impossible. The repeating pattern of pale-thin-pale-thick stripes of cytochrome oxidase (CO) activity in the primate secondary visual cortex (V2) is a prominent columnar system, in which histology also indicates different myelination of thin/thick versus pale stripes. We used quantitative magnetic resonance imaging (qMRI) in conjunction with functional magnetic resonance imaging (fMRI) at ultra-high field strength (7 T) to localize and study myelination of stripes in four human participants at sub-millimeter resolution in vivo. Thin and thick stripes were functionally localized by exploiting their sensitivity to color and binocular disparity, respectively. Resulting functional activation maps showed robust stripe patterns in V2 which enabled further comparison of quantitative relaxation parameters between stripe types. Thereby, we found lower longitudinal relaxation rates (R1) of thin and thick stripes compared to surrounding gray matter in the order of 1-2%, indicating higher myelination of pale stripes. No consistent differences were found for effective transverse relaxation rates (R2*). The study demonstrates the feasibility to investigate structure-function relationships in living humans within one cortical area at the level of columnar systems using qMRI.

Microstructural and neurochemical plasticity mechanisms interact to enhance human perceptual decision-making

Experience and training are known to boost our skills and mold the brain's organization and function. Yet, structural plasticity and functional neurotransmission are typically studied at different scales (large-scale networks, local circuits), limiting our understanding of the adaptive interactions that support learning of complex cognitive skills in the adult brain. Here, we employ multimodal brain imaging to investigate the link between microstructural (myelination) and neurochemical (GABAergic) plasticity for decision-making. We test (in males, due to potential confounding menstrual cycle effects on GABA measurements in females) for changes in MRI-measured myelin, GABA, and functional connectivity before versus after training on a perceptual decision task that involves identifying targets in clutter. We demonstrate that training alters subcortical (pulvinar, hippocampus) myelination and its functional connectivity to visual cortex and relates to decreased visual cortex GABAergic inhibition. Modeling interactions between MRI measures of myelin, GABA, and functional connectivity indicates that pulvinar myelin plasticity interacts-through thalamocortical connectivity-with GABAergic inhibition in visual cortex to support learning. Our findings propose a dynamic interplay of adaptive microstructural and neurochemical plasticity in subcortico-cortical circuits that supports learning for optimized decision-making in the adult human brain.

Quantitative MRI maps of human neocortex explored using cell type-specific gene expression analysis

Quantitative magnetic resonance imaging (qMRI) allows extraction of reproducible and robust parameter maps. However, the connection to underlying biological substrates remains murky, especially in the complex, densely packed cortex. We investigated associations in human neocortex between qMRI parameters and neocortical cell types by comparing the spatial distribution of the qMRI parameters longitudinal relaxation rate (${R_{1}}$), effective transverse relaxation rate (${R_{2}}^{\ast }$), and magnetization transfer saturation (MTsat) to gene expression from the Allen Human Brain Atlas, then combining this with lists of genes enriched in specific cell types found in the human brain. As qMRI parameters are magnetic field strength-dependent, the analysis was performed on MRI data at 3T and 7T. All qMRI parameters significantly covaried with genes enriched in GABA- and glutamatergic neurons, i.e. they were associated with cytoarchitecture. The qMRI parameters also significantly covaried with the distribution of genes enriched in astrocytes (${R_{2}}^{\ast }$ at 3T, ${R_{1}}$ at 7T), endothelial cells (${R_{1}}$ and MTsat at 3T), microglia (${R_{1}}$ and MTsat at 3T, ${R_{1}}$ at 7T), and oligodendrocytes and oligodendrocyte precursor cells (${R_{1}}$ at 7T). These results advance the potential use of qMRI parameters as biomarkers for specific cell types.

Abnormal brain iron accumulation in obstructive sleep apnea: A quantitative MRI study in the HypnoLaus cohort

Obstructive sleep apnea syndrome (OSA) may be a risk factor for Alzheimer's disease. One of the hallmarks of Alzheimer's disease is disturbed iron homeostasis leading to abnormal iron deposition in brain tissue. To date, there is no empirical evidence to support the hypothesis of altered brain iron homeostasis in patients with obstructive sleep apnea as well. Data were analysed from 773 participants in the HypnoLaus study (mean age 55.9 ± 10.3 years) who underwent polysomnography and brain MRI. Cross-sectional associations were tested between OSA parameters and the MRI effective transverse relaxation rate (R2*) - indicative of iron content - in 68 grey matter regions, after adjustment for confounders. The group with severe OSA (apnea-hypopnea index ≥30/h) had higher iron levels in the left superior frontal gyrus (F3,760  = 4.79, p = 0.003), left orbital gyri (F3,760  = 5.13, p = 0.002), right and left middle temporal gyrus (F3,760  = 4.41, p = 0.004 and F3,760  = 13.08, p < 0.001, respectively), left angular gyrus (F3,760  = 6.29, p = 0.001), left supramarginal gyrus (F3,760  = 4.98, p = 0.003), and right cuneus (F3,760  = 7.09, p < 0.001). The parameters of nocturnal hypoxaemia were all consistently associated with higher iron levels. Measures of sleep fragmentation had less consistent associations with iron content. This study provides the first evidence of increased brain iron levels in obstructive sleep apnea. The observed iron changes could reflect underlying neuropathological processes that appear to be driven primarily by hypoxaemic mechanisms.

Reliability of multi-parameter mapping (MPM) in the cervical cord: A multi-center multi-vendor quantitative MRI study

MRI based multicenter studies which target neurological pathologies affecting the spinal cord and brain - including spinal cord injury (SCI) - require standardized acquisition protocols and image processing methods. We have optimized and applied a multi-parameter mapping (MPM) protocol that simultaneously covers the brain and the cervical cord within a traveling heads study across six clinical centers (Leutritz et al., 2020). The MPM protocol includes quantitative maps (magnetization transfer saturation (MT), proton density (PD), longitudinal (R1), and effective transverse (R2*) relaxation rates) sensitive to myelination, water content, iron concentration, and morphometric measures, such as cross-sectional cord area. Previously, we assessed the repeatability and reproducibility of the brain MPM data acquired in the five healthy participants who underwent two scan-rescans (Leutritz et al., 2020). This study focuses on the cervical cord MPM data derived from the same acquisitions to determine its repeatability and reproducibility in the cervical cord. MPM matrices of the cervical cord were generated and processed using the hMRI and the spinal cord toolbox. To determine reliability of the cervical MPM data, the intra-site (i.e., scan-rescan) coefficient of variation (CoV), inter-site CoV, and bias within region of interests (C1, C2 and C3 levels) were determined. The range of the mean intra- and inter-site CoV of MT, R1 and PD was between 2.5% and 12%, and between 1.1% and 4.0% for the morphometric measures. In conclusion, the cervical MPM data showed a high repeatability and reproducibility for key imaging biomarkers and hence can be employed as a standardized tool in multi-center studies, including clinical trials.

Universal pulses for homogeneous excitation using single channel coils

PURPOSE

Universal Pulses (UPs) are excitation pulses that reduce the flip angle inhomogeneity in high field MRI systems without subject-specific optimization, originally developed for parallel transmit (PTX) systems at 7 T. We investigated the potential benefits of UPs for single channel (SC) transmit systems at 3 T, which are widely used for clinical and research imaging, and for which flip angle inhomogeneity can still be problematic.

METHODS

SC-UPs were designed using a spiral nonselective k-space trajectory for brain imaging at 3 T using transmit field maps (B1+) and off-resonance maps (B0) acquired on two different scanner types: a 'standard' single channel transmit system and a system with a PTX body coil. The effect of training group size was investigated using data (200 subjects) from the standard system. The PTX system was used to compare SC-UPs to PTX-UPs (15 subjects). In two additional subjects, prospective imaging using SC-UP was studied.

RESULTS

Average flip angle homogeneity error fell from 9.5 ± 0.5 % for 'default' excitation to 3.0 ± 0.6 % using SC-UPs trained over 50 subjects. Performance of the UPs was found to steadily improve as training group size increased, but stabilized after ~15 subjects. On the PTX-enabled system, SC-UPs again outperformed default excitation in simulations (4.8 ± 0.6 % error versus 10.6 ± 0.8 % respectively) though greater homogenization could be achieved with PTX-UPs (3.9 ± 0.6 %) and personalized pulses (SC-PP 3.6 ± 1.0 %, PTX-PP 2.9 ± 0.6 %). MP-RAGE imaging using SC-UP resulted in greater separation between grey and white matter signal intensities than default excitation.

CONCLUSIONS

SC-UPs can improve excitation homogeneity in standard 3 T systems without further calibration and could be used instead of a default excitation pulse for nonselective neuroimaging at 3 T.

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.

GnRH replacement rescues cognition in Down syndrome

At the present time, no viable treatment exists for cognitive and olfactory deficits in Down syndrome (DS). We show in a DS model (Ts65Dn mice) that these progressive nonreproductive neurological symptoms closely parallel a postpubertal decrease in hypothalamic as well as extrahypothalamic expression of a master molecule that controls reproduction-gonadotropin-releasing hormone (GnRH)-and appear related to an imbalance in a microRNA-gene network known to regulate GnRH neuron maturation together with altered hippocampal synaptic transmission. Epigenetic, cellular, chemogenetic, and pharmacological interventions that restore physiological GnRH levels abolish olfactory and cognitive defects in Ts65Dn mice, whereas pulsatile GnRH therapy improves cognition and brain connectivity in adult DS patients. GnRH thus plays a crucial role in olfaction and cognition, and pulsatile GnRH therapy holds promise to improve cognitive deficits in DS.

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.

Test-retest reliability of multi-parametric maps (MPM) of brain microstructure

Multiparameter mapping (MPM) is a quantitative MRI protocol that is promising for studying microstructural brain changes in vivo with high specificity. Reliability values are an important prior knowledge for efficient study design and facilitating replicable findings in development, aging and neuroplasticity research. To explore longitudinal reliability of MPM we acquired the protocol in 31 healthy young subjects twice over a rescan interval of 4 weeks. We assessed the within-subject coefficient of variation (WCV), the between-subject coefficient of variation (BCV), and the intraclass correlation coefficient (ICC). Using these metrics, we investigated the reliability of (semi-) quantitative magnetization transfer saturation (MT_sat), proton density (PD), transversal relaxation (R2*) and longitudinal relaxation (R1). To increase relevance for explorative studies in development and training-induced plasticity, we assess reliability both on local voxel- as well as ROI-level. Finally, we disentangle contributions and interplay of within- and between-subject variability to ICC and assess the optimal degree of spatial smoothing applied to the data. We reveal evidence that voxelwise ICC reliability of MPMs is moderate to good with median values in cortex (subcortical GM): MT: 0.789 (0.447) PD: 0.553 (0.264) R1: 0.555 (0.369) R2*: 0.624 (0.477). The Gaussian smoothing kernel of 2 to 4 mm FWHM resulted in optimal reproducibility. We discuss these findings in the context of longitudinal intervention studies and the application to research designs in neuroimaging field.

In vivo Estimation of Axonal Morphology From Magnetic Resonance Imaging and Electroencephalography Data

Purpose

We present a novel approach that allows the estimation of morphological features of axonal fibers from data acquired in vivo in humans. This approach allows the assessment of white matter microscopic properties non-invasively with improved specificity.

Theory

The proposed approach is based on a biophysical model of Magnetic Resonance Imaging (MRI) data and of axonal conduction velocity estimates obtained with Electroencephalography (EEG). In a white matter tract of interest, these data depend on (1) the distribution of axonal radius [P(r)] and (2) the g-ratio of the individual axons that compose this tract [g(r)]. P(r) is assumed to follow a Gamma distribution with mode and scale parameters, M and θ, and g(r) is described by a power law with parameters α and β.

Methods

MRI and EEG data were recorded from 14 healthy volunteers. MRI data were collected with a 3T scanner. MRI-measured g-ratio maps were computed and sampled along the visual transcallosal tract. EEG data were recorded using a 128-lead system with a visual Poffenberg paradigm. The interhemispheric transfer time and axonal conduction velocity were computed from the EEG current density at the group level. Using the MRI and EEG measures and the proposed model, we estimated morphological properties of axons in the visual transcallosal tract.

Results

The estimated interhemispheric transfer time was 11.72 ± 2.87 ms, leading to an average conduction velocity across subjects of 13.22 ± 1.18 m/s. Out of the 4 free parameters of the proposed model, we estimated θ – the width of the right tail of the axonal radius distribution – and β – the scaling factor of the axonal g-ratio, a measure of fiber myelination. Across subjects, the parameter θ was 0.40 ± 0.07 μm and the parameter β was 0.67 ± 0.02 μm^−α.

Conclusion

The estimates of axonal radius and myelination are consistent with histological findings, illustrating the feasibility of this approach. The proposed method allows the measurement of the distribution of axonal radius and myelination within a white matter tract, opening new avenues for the combined study of brain structure and function, and for in vivo histological studies of the human brain.

Reliability of quantitative multiparameter maps is high for magnetization transfer and proton density but attenuated for R1 and R2 * in healthy young adults

We investigate the reliability of individual differences of four quantities measured by magnetic resonance imaging‐based multiparameter mapping (MPM): magnetization transfer saturation (MT), proton density (PD), longitudinal relaxation rate (R₁), and effective transverse relaxation rate (R₂*). Four MPM datasets, two on each of two consecutive days, were acquired in healthy young adults. On Day 1, no repositioning occurred and on Day 2, participants were repositioned between MPM datasets. Using intraclass correlation effect decomposition (ICED), we assessed the contributions of session‐specific, day‐specific, and residual sources of measurement error. For whole‐brain gray and white matter, all four MPM parameters showed high reproducibility and high reliability, as indexed by the coefficient of variation (CoV) and the intraclass correlation (ICC). However, MT, PD, R₁, and R₂* differed markedly in the extent to which reliability varied across brain regions. MT and PD showed high reliability in almost all regions. In contrast, R₁ and R₂* showed low reliability in some regions outside the basal ganglia, such that the sum of the measurement error estimates in our structural equation model was higher than estimates of between‐person differences. In addition, in this sample of healthy young adults, the four MPM parameters showed very little variability over four measurements but differed in how well they could assess between‐person differences. We conclude that R₁ and R₂* might carry only limited person‐specific information in some regions of the brain in healthy young adults, and, by implication, might be of restricted utility for studying associations to between‐person differences in behavior in those regions.

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

PURPOSE

Inter-scan motion is a substantial source of error in R 1 estimation methods based on multiple volumes, for example, variable flip angle (VFA), and can be expected to increase at 7T where B 1 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

R 1 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. R 1 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 B 1 + 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 R 1 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.

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.

Recommendations of Choice of Head Coil and Prescan Normalize Filter Depend on Region of Interest and Task

The performance of MRI head coils together with the influence of the prescan normalize filter in different brain regions was evaluated. Functional and structural data were recorded from 26 participants performing motor, auditory, and visual tasks in different conditions: with the 20- and 64-channel Siemens head/neck coil and the prescan normalize filter turned ON or OFF. Data were analyzed with the MRIQC tool to evaluate data quality differences. The functional data were statistically evaluated by comparison of the β estimates and the time-course signal-to-noise ratio (tSNR) in four regions of interest, i.e., the auditory, visual, and motor cortices and the thalamus. The MRIQC tool indicated a better data quality for both functional and structural data with the prescan normalize filter, with an advantage for the 20-channel head coil in functional data and an advantage for the 64-channel head coil in structural measurements. Nevertheless, recommendations for the functional data regarding choice of head coils and prescan normalize filter depend on the brain regions of interest. Higher β estimates and tSNR values occurred in the auditory cortex and thalamus with the prescan normalize filter, whereas the contrary was true for the visual and motor cortices. Due to higher β estimates in the visual cortex in the 64-channel head coil, this head coil is recommended for studies investigating the visual cortex. For most of the research questions, the 20-channel head coil is better suited for functional experiments, with the prescan normalize filter, especially when investigating deep brain areas. For anatomical studies, the 64-channel head coil seemed to be the better choice.

Laminar dynamics of high amplitude beta bursts in human motor cortex

Motor cortical activity in the beta frequency range is one of the strongest and most studied movement-related neural signals. At the single trial level, beta band activity is often characterized by transient, high amplitude, bursting events rather than slowly modulating oscillations. The timing of these bursting events is tightly linked to behavior, suggesting a more dynamic functional role for beta activity than previously believed. However, the neural mechanisms underlying beta bursts in sensorimotor circuits are poorly understood. To address this, we here leverage and extend recent developments in high precision MEG for temporally resolved laminar analysis of burst activity, combined with a neocortical circuit model that simulates the biophysical generators of the electrical currents which drive beta bursts. This approach pinpoints the generation of beta bursts in human motor cortex to distinct excitatory synaptic inputs to deep and superficial cortical layers, which drive current flow in opposite directions. These laminar dynamics of beta bursts in motor cortex align with prior invasive animal recordings within the somatosensory cortex, and suggest a conserved mechanism for somatosensory and motor cortical beta bursts. More generally, we demonstrate the ability for uncovering the laminar dynamics of event-related neural signals in human non-invasive recordings. This provides important constraints to theories about the functional role of burst activity for movement control in health and disease, and crucial links between macro-scale phenomena measured in humans and micro-circuit activity recorded from animal models.

Model-based multi-parameter mapping

Quantitative MR imaging is increasingly favoured for its richer information content and standardised measures. However, computing quantitative parameter maps, such as those encoding longitudinal relaxation rate (R1), apparent transverse relaxation rate (R2*) or magnetisation-transfer saturation (MTsat), involves inverting a highly non-linear function. Many methods for deriving parameter maps assume perfect measurements and do not consider how noise is propagated through the estimation procedure, resulting in needlessly noisy maps. Instead, we propose a probabilistic generative (forward) model of the entire dataset, which is formulated and inverted to jointly recover (log) parameter maps with a well-defined probabilistic interpretation (e.g., maximum likelihood or maximum a posteriori). The second order optimisation we propose for model fitting achieves rapid and stable convergence thanks to a novel approximate Hessian. We demonstrate the utility of our flexible framework in the context of recovering more accurate maps from data acquired using the popular multi-parameter mapping protocol. We also show how to incorporate a joint total variation prior to further decrease the noise in the maps, noting that the probabilistic formulation allows the uncertainty on the recovered parameter maps to be estimated. Our implementation uses a PyTorch backend and benefits from GPU acceleration. It is available at https://github.com/balbasty/nitorch.

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.

Parkinson's disease multimodal imaging: F-DOPA PET, neuromelanin-sensitive and quantitative iron-sensitive MRI

Parkinson's disease (PD) is a neurodegenerative synucleinopathy characterized by the degeneration of neuromelanin (NM)-containing dopaminergic neurons and deposition of iron in the substantia nigra (SN). How regional NM loss and iron accumulation within specific areas of SN relate to nigro-striatal dysfunction needs to be clarified. We measured dopaminergic function in pre- and postcommissural putamen by [¹⁸F]DOPA PET in 23 Parkinson's disease patients and 23 healthy control (HC) participants in whom NM content and iron load were assessed in medial and lateral SN, respectively, by NM-sensitive and quantitative R2* MRI. Data analysis consisted of voxelwise regressions testing the group effect and its interaction with NM or iron signals. In PD patients, R2* was selectively increased in left lateral SN as compared to healthy participants, suggesting a local accumulation of iron in Parkinson's disease. By contrast, NM signal differed between PD and HC, without specific regional specificity within SN. Dopaminergic function in posterior putamen decreased as R2* increased in lateral SN, indicating that dopaminergic function impairment progresses with iron accumulation in the SN. Dopaminergic function was also positively correlated with NM signal in lateral SN, indicating that dopaminergic function impairment progresses with depigmentation in the SN. A complex relationship was detected between R2* in the lateral SN and NM signal in the medial SN. In conclusion, multimodal imaging reveals regionally specific relationships between iron accumulation and depigmentation within the SN of Parkinson's disease and provides in vivo insights in its neuropathology.

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.

The Motor Basis for Misophonia

Misophonia is a common disorder characterized by the experience of strong negative emotions of anger and anxiety in response to certain everyday sounds, such as those generated by other people eating, drinking, and breathing. The commonplace nature of these "trigger" sounds makes misophonia a devastating disorder for sufferers and their families. How such innocuous sounds trigger this response is unknown. Since most trigger sounds are generated by orofacial movements (e.g., chewing) in others, we hypothesized that the mirror neuron system related to orofacial movements could underlie misophonia. We analyzed resting state fMRI (rs-fMRI) connectivity (N = 33, 16 females) and sound-evoked fMRI responses (N = 42, 29 females) in misophonia sufferers and controls. We demonstrate that, compared with controls, the misophonia group show no difference in auditory cortex responses to trigger sounds, but do show: (1) stronger rs-fMRI connectivity between both auditory and visual cortex and the ventral premotor cortex responsible for orofacial movements; (2) stronger functional connectivity between the auditory cortex and orofacial motor area during sound perception in general; and (3) stronger activation of the orofacial motor area, specifically, in response to trigger sounds. Our results support a model of misophonia based on "hyper-mirroring" of the orofacial actions of others with sounds being the "medium" via which action of others is excessively mirrored. Misophonia is therefore not an abreaction to sounds, per se, but a manifestation of activity in parts of the motor system involved in producing those sounds. This new framework to understand misophonia can explain behavioral and emotional responses and has important consequences for devising effective therapies.SIGNIFICANCE STATEMENT Conventionally, misophonia, literally "hatred of sounds" has been considered as a disorder of sound emotion processing, in which "simple" eating and chewing sounds produced by others cause negative emotional responses. Our data provide an alternative but complementary perspective on misophonia that emphasizes the action of the trigger-person rather than the sounds which are a byproduct of that action. Sounds, in this new perspective, are only a "medium" via which action of the triggering-person is mirrored onto the listener. This change in perspective has important consequences for devising therapies and treatment methods for misophonia. It suggests that, instead of focusing on sounds, which many existing therapies do, effective therapies should target the brain representation of movement.

Temporal trajectory of brain tissue property changes induced by electroconvulsive therapy

BACKGROUND

After more than eight decades of electroconvulsive therapy (ECT) for pharmaco-resistant depression, the mechanisms governing its anti-depressant effects remain poorly understood. Computational anatomy studies using longitudinal T1-weighted magnetic resonance imaging (MRI) data have demonstrated ECT effects on hippocampus volume and cortical thickness, but they lack the interpretational specificity about underlying neurobiological processes.

METHODS

We sought to fill in the gap of knowledge by acquiring quantitative MRI indicative for brain's myelin, iron and tissue water content at multiple time-points before, during and after ECT treatment. We adapted established tools for longitudinal spatial registration of MRI data to the relaxometry-based multi-parameter maps aiming to preserve the initial total signal amount and introduced a dedicated multivariate analytical framework.

RESULTS

The whole-brain voxel-based analysis based on a multivariate general linear model showed that there is no brain tissue oedema contributing to the predicted ECT-induced hippocampus volume increase neither in the short, nor in the long-term observations. Improvements in depression symptom severity over time were associated with changes in both volume estimates and brain tissue properties expanding beyond mesial temporal lobe structures to anterior cingulate cortex, precuneus and striatum.

CONCLUSION

The obtained results stemming from multi-contrast MRI quantitative data provided a fingerprint of ECT-induced brain tissue changes over time that are contrasted against the background of established morphometry findings. The introduced data processing and statistical testing algorithms provided a reliable analytical framework for longitudinal multi-parameter brain maps. The results, particularly the evidence of lack of ECT impact on brain tissue water, should be considered preliminary considering the small sample size of the study.

A standard system phantom for magnetic resonance imaging

Purpose

A standard MRI system phantom has been designed and fabricated to assess scanner performance, stability, comparability and assess the accuracy of quantitative relaxation time imaging. The phantom is unique in having traceability to the International System of Units, a high level of precision, and monitoring by a national metrology institute. Here, we describe the phantom design, construction, imaging protocols, and measurement of geometric distortion, resolution, slice profile, signal‐to‐noise ratio (SNR), proton‐spin relaxation times, image uniformity and proton density.

Methods

The system phantom, designed by the International Society of Magnetic Resonance in Medicine ad hoc committee on Standards for Quantitative MR, is a 200 mm spherical structure that contains a 57‐element fiducial array; two relaxation time arrays; a proton density/SNR array; resolution and slice‐profile insets. Standard imaging protocols are presented, which provide rapid assessment of geometric distortion, image uniformity, T₁ and T₂ mapping, image resolution, slice profile, and SNR.

Results

Fiducial array analysis gives assessment of intrinsic geometric distortions, which can vary considerably between scanners and correction techniques. This analysis also measures scanner/coil image uniformity, spatial calibration accuracy, and local volume distortion. An advanced resolution analysis gives both scanner and protocol contributions. SNR analysis gives both temporal and spatial contributions.

Conclusions

A standard system phantom is useful for characterization of scanner performance, monitoring a scanner over time, and to compare different scanners. This type of calibration structure is useful for quality assurance, benchmarking quantitative MRI protocols, and to transition MRI from a qualitative imaging technique to a precise metrology with documented accuracy and uncertainty.

The Myelin Content of the Human Precentral Hand Knob Reflects Interindividual Differences in Manual Motor Control at the Physiological and Behavioral Level

The primary motor cortex hand area (M1_HAND) and adjacent dorsal premotor cortex (PMd) form the so-called motor hand knob in the precentral gyrus. M1_HAND and PMd are critical for dexterous hand use and are densely interconnected via corticocortical axons, lacking a sharp demarcating border. In 24 young right-handed volunteers, we performed multimodal mapping to delineate the relationship between structure and function in the right motor hand knob. Quantitative structural magnetic resonance imaging (MRI) at 3 tesla yielded regional R1 maps as a proxy of cortical myelin content. Participants also underwent functional MRI (fMRI). We mapped task-related activation and temporal precision, while they performed a visuomotor synchronization task requiring visually cued abduction movements with the left index or little finger. We also performed sulcus-aligned transcranial magnetic stimulation of the motor hand knob to localize the optimal site (hotspot) for evoking a motor evoked potential (MEP) in two intrinsic hand muscles. Individual motor hotspot locations varied along the rostrocaudal axis. The more rostral the motor hotspot location in the precentral crown, the longer were corticomotor MEP latencies. “Hotspot rostrality” was associated with the regional myelin content in the precentral hand knob. Cortical myelin content also correlated positively with task-related activation of the precentral crown and temporal precision during the visuomotor synchronization task. Together, our results suggest a link among cortical myelination, the spatial cortical representation, and temporal precision of finger movements. We hypothesize that the myelination of cortical axons facilitates neuronal integration in PMd and M1_HAND and, hereby, promotes the precise timing of movements.

SIGNIFICANCE STATEMENT Here we used magnetic resonance imaging and transcranial magnetic stimulation of the precentral motor hand knob to test for a link among cortical myelin content, functional corticomotor representations, and manual motor control. A higher myelin content of the precentral motor hand knob was associated with more rostral corticomotor presentations, with stronger task-related activation and a higher precision of movement timing during a visuomotor synchronization task. We propose that a high precentral myelin content enables fast and precise neuronal integration in M1 (primary motor cortex) and dorsal premotor cortex, resulting in higher temporal precision during dexterous hand use. Our results identify the degree of myelination as an important structural feature of the neocortex that is tightly linked to the function and behavior supported by the cortical area.

Brain tissue properties link cardio-vascular risk factors, mood and cognitive performance in the CoLaus|PsyCoLaus epidemiological cohort

Given the controversy about the impact of modifiable risk factors on mood and cognition in ageing, we sought to investigate the associations between cardio-vascular risk, mental health, cognitive performance and brain anatomy in mid- to old age. We analyzed a set of risk factors together with multi-parameter magnetic resonance imaging (MRI) in the CoLaus|PsyCoLaus cohort (n > 1200). Cardio-vascular risk was associated with differences in brain tissue properties - myelin, free tissue water, iron content - and regional brain volumes that we interpret in the context of micro-vascular hypoxic lesions and neurodegeneration. The interaction between clinical subtypes of major depressive disorder and cardio-vascular risk factors showed differential associations with brain structure depending on individuals' lifetime trajectory. There was a negative correlation between melancholic depression, anxiety and MRI markers of myelin and iron content in the hippocampus and anterior cingulate. Verbal memory and verbal fluency performance were positively correlated with left amygdala volumes. The concomitant analysis of brain morphometry and tissue properties allowed for a neuro-biological interpretation of the link between modifiable risk factors and brain health.

The relationship between hippocampal-dependent task performance and hippocampal grey matter myelination and iron content

Individual differences in scene imagination, autobiographical memory recall, future thinking and spatial navigation have long been linked with hippocampal structure in healthy people, although evidence for such relationships is, in fact, mixed. Extant studies have predominantly concentrated on hippocampal volume. However, it is now possible to use quantitative neuroimaging techniques to model different properties of tissue microstructure in vivo such as myelination and iron. Previous work has linked such measures with cognitive task performance, particularly in older adults. Here we investigated whether performance on scene imagination, autobiographical memory, future thinking and spatial navigation tasks was associated with hippocampal grey matter myelination or iron content in young, healthy adult participants. Magnetic resonance imaging data were collected using a multi-parameter mapping protocol (0.8 mm isotropic voxels) from a large sample of 217 people with widely-varying cognitive task scores. We found little evidence that hippocampal grey matter myelination or iron content were related to task performance. This was the case using different analysis methods (voxel-based quantification, partial correlations), when whole brain, hippocampal regions of interest, and posterior:anterior hippocampal ratios were examined, and across different participant sub-groups (divided by gender and task performance). Variations in hippocampal grey matter myelin and iron levels may not, therefore, help to explain individual differences in performance on hippocampal-dependent tasks, at least in young, healthy individuals.

A practical method for post-acquisition reduction of bias in fast, whole-brain B1-maps

Large consistent differences have been observed between maps of the flip angle correction factor (commonly called "B₁-maps") produced with different fast methods in the human brain. We present an empirical procedure for first-order multiplicative bias correction that can be applied when more than one B₁-mapping method is available. We use a B₁-map measurement in a calibration phantom as a reference and the voxel-wise histogram mode between ratios of B₁-maps produced from different methods to calculate determine the bias as a multiplicative correcting scale factor. Institutional implementations of four common methods of B₁-mapping were assessed: Method of Slopes, FSE and EPI double angle methods (DAM), and Bloch-Siegert. In human subjects, the multiplicative bias used to correct for each of the four methods was: Method of Slopes = 1.005, FSE-DAM = 0.956, EPI-DAM = 1.080, and Bloch-Siegert = 1.128. Scaling to remove this bias between methods produces more consistent B₁-maps which enable more consistent values for any computations requiring flip angle correction. In addition, we present evidence that the corrected B₁ maps, using our calibration method, are also more accurate.

Mapping Human Laryngeal Motor Cortex during Vocalization

The representations of the articulators involved in human speech production are organized somatotopically in primary motor cortex. The neural representation of the larynx, however, remains debated. Both a dorsal and a ventral larynx representation have been previously described. It is unknown, however, whether both representations are located in primary motor cortex. Here, we mapped the motor representations of the human larynx using functional magnetic resonance imaging and characterized the cortical microstructure underlying the activated regions. We isolated brain activity related to laryngeal activity during vocalization while controlling for breathing. We also mapped the articulators (the lips and tongue) and the hand area. We found two separate activations during vocalization-a dorsal and a ventral larynx representation. Structural and quantitative neuroimaging revealed that myelin content and cortical thickness underlying the dorsal, but not the ventral larynx representation, are similar to those of other primary motor representations. This finding confirms that the dorsal larynx representation is located in primary motor cortex and that the ventral one is not. We further speculate that the location of the ventral larynx representation is in premotor cortex, as seen in other primates. It remains unclear, however, whether and how these two representations differentially contribute to laryngeal motor control.

Does hippocampal volume explain performance differences on hippocampal-dependant tasks?

Marked disparities exist across healthy individuals in their ability to imagine scenes, recall autobiographical memories, think about the future and navigate in the world. The importance of the hippocampus in supporting these critical cognitive functions has prompted the question of whether differences in hippocampal grey matter volume could be one source of performance variability. Evidence to date has been somewhat mixed. In this study we sought to mitigate issues that commonly affect these types of studies. Data were collected from a large sample of 217 young, healthy adult participants, including whole brain structural MRI data (0.8 mm isotropic voxels) and widely-varying performance on scene imagination, autobiographical memory, future thinking and navigation tasks. We found little evidence that hippocampal grey matter volume was related to task performance in this healthy sample. This was the case using different analysis methods (voxel-based morphometry, partial correlations), when whole brain or hippocampal regions of interest were examined, when comparing different sub-groups (divided by gender, task performance, self-reported ability), and when using latent variables derived from across the cognitive tasks. Hippocampal grey matter volume may not, therefore, significantly influence performance on tasks known to require the hippocampus in healthy people. Perhaps only in extreme situations, as in the case of licensed London taxi drivers, are measurable ability-related hippocampus volume changes consistently exhibited.

Multiparameter mapping of relaxation (R1, R2*), proton density and magnetization transfer saturation at 3 T: A multicenter dual-vendor reproducibility and repeatability study

Multicenter clinical and quantitative magnetic resonance imaging (qMRI) studies require a high degree of reproducibility across different sites and scanner manufacturers, as well as time points. We therefore implemented a multiparameter mapping (MPM) protocol based on vendor's product sequences and demonstrate its repeatability and reproducibility for whole‐brain coverage. Within ~20 min, four MPM metrics (magnetization transfer saturation [MT], proton density [PD], longitudinal [R1], and effective transverse [R2*] relaxation rates) were measured using an optimized 1 mm isotropic resolution protocol on six 3 T MRI scanners from two different vendors. The same five healthy participants underwent two scanning sessions, on the same scanner, at each site. MPM metrics were calculated using the hMRI‐toolbox. To account for different MT pulses used by each vendor, we linearly scaled the MT values to harmonize them across vendors. To determine longitudinal repeatability and inter‐site comparability, the intra‐site (i.e., scan‐rescan experiment) coefficient of variation (CoV), inter‐site CoV, and bias across sites were estimated. For MT, R1, and PD, the intra‐ and inter‐site CoV was between 4 and 10% across sites and scan time points for intracranial gray and white matter. A higher intra‐site CoV (16%) was observed in R2* maps. The inter‐site bias was below 5% for all parameters. In conclusion, the MPM protocol yielded reliable quantitative maps at high resolution with a short acquisition time. The high reproducibility of MPM metrics across sites and scan time points combined with its tissue microstructure sensitivity facilitates longitudinal multicenter imaging studies targeting microstructural changes, for example, as a quantitative MRI biomarker for interventional clinical trials.