Fast magnetization transfer saturation imaging of the brain using MP2RAGE T1 mapping

PURPOSE

Magnetization transfer saturation (MTsat) mapping is commonly used to examine the macromolecular content of brain tissue. This study compared variable flip angle (VFA) T1 mapping against compressed-sensing MP2RAGE (csMP2RAGE) T1 mapping for accelerating MTsat imaging.

METHODS

VFA, MP2RAGE, and csMP2RAGE were compared against inversion-recovery T1 in an aqueous phantom at 3 T. The same 1-mm VFA, MP2RAGE, and csMP2RAGE protocols were acquired in 4 healthy subjects to compare T1 and MTsat. Bloch-McConnell simulations were used to investigate differences between the phantom and in vivo T1 results. Ten healthy controls were imaged twice with the csMP2RAGE MTsat protocol to quantify repeatability.

RESULTS

The MP2RAGE and csMP2RAGE protocols were 13.7% and 32.4% faster than the VFA protocol, respectively. At these scan times, all approaches provided strong repeatability and accurate T1 times (< 5% difference) in the phantom, but T1 accuracy was more impacted by T2 for VFA than for MP2RAGE. In vivo, VFA estimated longer T1 times than MP2RAGE and csMP2RAGE. Simulations suggest that the differences in the T1 measured using VFA, MP2RAGE, and inversion recovery could be explained by the magnetization-transfer effects. In the test-retest experiment, we found that the csMP2RAGE has a minimum detectable change of 2.3% for T1 mapping and 7.8% for MTsat imaging.

CONCLUSIONS

We demonstrated that MP2RAGE can be used in place of VFA T1 mapping in an MTsat protocol. Furthermore, a shorter scan time and high repeatability can be achieved using the csMP2RAGE sequence.

Decreased long-range temporal correlations in the resting-state functional magnetic resonance imaging blood-oxygen-level-dependent signal reflect motor sequence learning up to 2 weeks following training

Decreased long-range temporal correlations (LRTC) in brain signals can be used to measure cognitive effort during task execution. Here, we examined how learning a motor sequence affects long-range temporal memory within resting-state functional magnetic resonance imaging signal. Using the Hurst exponent (HE), we estimated voxel-wise LRTC and assessed changes over 5 consecutive days of training, followed by a retention scan 12 days later. The experimental group learned a complex visuomotor sequence while a complementary control group performed tightly matched movements. An interaction analysis revealed that HE decreases were specific to the complex sequence and occurred in well-known motor sequence learning associated regions including left supplementary motor area, left premotor cortex, left M1, left pars opercularis, bilateral thalamus, and right striatum. Five regions exhibited moderate to strong negative correlations with overall behavioral performance improvements. Following learning, HE values returned to pretraining levels in some regions, whereas in others, they remained decreased even 2 weeks after training. Our study presents new evidence of HE's possible relevance for functional plasticity during the resting-state and suggests that a cortical subset of sequence-specific regions may continue to represent a functional signature of learning reflected in decreased long-range temporal dependence after a period of inactivity.

The human brain connectome weighted by the myelin content and total intra-axonal cross-sectional area of white matter tracts

A central goal in neuroscience is the development of a comprehensive mapping between structural and functional brain features, which facilitates mechanistic interpretation of brain function. However, the interpretability of structure-function brain models remains limited by a lack of biological detail. Here, we characterize human structural brain networks weighted by multiple white matter microstructural features including total intra-axonal cross-sectional area and myelin content. We report edge-weight-dependent spatial distributions, variance, small-worldness, rich club, hubs, as well as relationships with function, edge length, and myelin. Contrasting networks weighted by the total intra-axonal cross-sectional area and myelin content of white matter tracts, we find opposite relationships with functional connectivity, an edge-length-independent inverse relationship with each other, and the lack of a canonical rich club in myelin-weighted networks. When controlling for edge length, networks weighted by either fractional anisotropy, radial diffusivity, or neurite density show no relationship with whole-brain functional connectivity. We conclude that the co-utilization of structural networks weighted by total intra-axonal cross-sectional area and myelin content could improve our understanding of the mechanisms mediating the structure-function brain relationship.

Optimization of acquisition parameters for cortical inhomogeneous magnetization transfer (ihMT) imaging using a rapid gradient echo readout

PURPOSE

Imaging biomarkers with increased myelin specificity are needed to better understand the complex progression of neurological disorders. Inhomogeneous magnetization transfer (ihMT) imaging is an emergent technique that has a high degree of specificity for myelin content but suffers from low signal to-noise ratio (SNR). This study used simulations to determine optimal sequence parameters for ihMT imaging for use in high-resolution cortical mapping.

METHODS

MT-weighted cortical image intensity and ihMT SNR were simulated using modified Bloch equations for a range of sequence parameters. The acquisition time was limited to 4.5 min/volume. A custom MT-weighted RAGE sequence with center-out k-space encoding was used to enhance SNR at 3 T. Pulsed MT imaging was studied over a range of saturation parameters, and the impact of the turbo factor on the effective ihMT resolution was investigated. 1 mm isotropic ihMTsat maps were generated in 25 healthy adults.

RESULTS

Greater SNR was observed for larger number of bursts consisting of 6-8 saturation pulses each, combined with a high readout turbo factor. However, that protocol suffered from a point spread function that was more than twice the nominal resolution. For high-resolution cortical imaging, we selected a protocol with a higher effective resolution at the cost of a lower SNR. We present the first group-average ihMTsat whole-brain map at 1 mm isotropic resolution.

CONCLUSION

This study presents the impact of saturation and excitation parameters on ihMTsat SNR and resolution. We demonstrate the feasibility of high-resolution cortical myelin imaging using ihMTsat in less than 20 min.

High-resolution diffusion-weighted imaging at 7 Tesla: single-shot readout trajectories and their impact on signal-to-noise ratio, spatial resolution and accuracy

Diffusion MRI (dMRI) is a valuable imaging technique to study the connectivity and microstructure of the brain in vivo. However, the resolution of dMRI is limited by the low signal-to-noise ratio (SNR) of this technique. Various multi-shot acquisition strategies have been developed to achieve sub-millimeter resolution, but they require long scan times which can be restricting for patient scans. Alternatively, the SNR of single-shot acquisitions can be increased by using a spiral readout trajectory to minimize the sequence echo time. Imaging at ultra-high fields (UHF) could further increase the SNR of single-shot dMRI; however, the shorter T2* of brain tissue and the greater field non-uniformities at UHFs will degrade image quality, causing image blurring, distortions, and signal loss. In this study, we investigated the trade-off between the SNR and resolution of different k-space trajectories, including echo planar imaging (EPI), partial Fourier EPI, and spiral trajectories, over a range of dMRI resolutions at 7T. The effective resolution, spatial specificity and sharpening effect were measured from the point spread function (PSF) of the simulated diffusion sequences for a nominal resolution range of 0.6-1.8 mm. In-vivo partial brain scans at a nominal resolution of 1.5 mm isotropic were acquired using the three readout trajectories to validate the simulation results. Field probes were used to measure dynamic magnetic fields offline up to the 3^rd order of spherical harmonics. Image reconstruction was performed using static ΔB₀ field maps and the measured trajectories to correct image distortions and artifacts, leaving T2* effects as the primary source of blurring. The effective resolution was examined in fractional anisotropy (FA) maps calculated from a multi-shell dataset with b-values of 300, 1000, and 2000 s/mm² in 5, 16, and 48 directions, respectively. In-vivo scans at nominal resolutions of 1, 1.2, and 1.5 mm were acquired and the SNR of the different trajectories calculated using the multiple replica method to investigate the SNR. Finally, in-vivo whole brain scans with an effective resolution of 1.5 mm isotropic were acquired to explore the SNR and efficiency of different trajectories at a matching effective resolution. FA and intra-cellular volume fraction (ICVF) maps calculated using neurite orientation dispersion and density imaging (NODDI) were used for the comparison. The simulations and in vivo imaging results showed that for matching nominal resolutions, EPI trajectories had the highest specificity and effective resolution with maximum image sharpening effect. However, spirals have a significantly higher SNR, in particular at higher resolutions and even when the effective image resolutions are matched. Overall, this work shows that the higher SNR of single-shot spiral trajectories at 7T allows us to achieve higher effective resolutions compared to EPI and PF-EPI to map the microstructure and connectivity of small brain structures.

High spatial overlap but diverging age-related trajectories of cortical magnetic resonance imaging markers aiming to represent intracortical myelin and microstructure

Statistical effects of cortical metrics derived from standard T1- and T2-weighted magnetic resonance imaging (MRI) images, such as gray-white matter contrast (GWC), boundary sharpness coefficient (BSC), T1-weighted/T2-weighted ratio (T1w/T2w), and cortical thickness (CT), are often interpreted as representing or being influenced by intracortical myelin content with little empirical evidence to justify these interpretations. We first examined spatial correspondence with more biologically specific microstructural measures, and second compared between-marker age-related trends with the underlying hypothesis that different measures primarily driven by similar changes in myelo- and microstructural underpinnings should be highly related. Cortical MRI markers were derived from MRI images of 127 healthy subjects, aged 18-81, using cortical surfaces that were generated with the CIVET 2.1.0 pipeline. Their gross spatial distributions were compared with gene expression-derived cell-type densities, histology-derived cytoarchitecture, and quantitative R1 maps acquired on a subset of participants. We then compared between-marker age-related trends in their shape, direction, and spatial distribution of the linear age effect. The gross anatomical distributions of cortical MRI markers were, in general, more related to myelin and glial cells than neuronal indicators. Comparing MRI markers, our results revealed generally high overlap in spatial distribution (i.e., group means), but mostly divergent age trajectories in the shape, direction, and spatial distribution of the linear age effect. We conclude that the microstructural properties at the source of spatial distributions of MRI cortical markers can be different from microstructural changes that affect these markers in aging.

Loss of age-related laminar differentiation of intracortical myelin in bipolar disorder

Age-related changes of intracortical myelin in bipolar disorder (BD) have been observed to deviate from the quadratic age curve observed in healthy controls (HC), but it is unclear if this holds at varying cortical depths. From BD (n = 44; age range = 17.6-45.5 years) and HC (n = 60; age range = 17.1-45.8 years) participants, we collected 3T T1-weighted (T1w) images with strong intracortical contrast. Signal values were sampled from 3 equivolume cortical depths. Linear mixed models were used to compare age-related changes in the T1w signal between depths and between groups at each depth. In HC, the age-related changes were significantly different between the superficial one-fourth depth and the deeper depths in the right ventral somatosensory (t = -4.63; FDRp = 0.00025), left dorsomedial somatosensory (t = -3.16; FDRp = 0.028), left rostral ventral premotor (t = -3.16; FDRp = 0.028), and right ventral inferior parietal cortex (t = -3.29; FDRp = 0.028). BD participants exhibited no differences in the age-related T1w signal between depths. Illness duration was negatively correlated with the T1w signal at the one-fourth depth in the right anterior cingulate cortex (rACC; rho = -0.50; FDRp = 0.029). Physiological age-related and depth-specific variation in the T1w signal were not observed in BD. The T1w signal in the rACC may reflect lifetime disease burden in the disorder.

Co-registration of Imaging Modalities (MRI, CT and PET) to Perform Frameless Stereotaxic Robotic Injections in the Common Marmoset

The common marmoset has emerged as a popular model in neuroscience research, in part due to its reproductive efficiency, genetic and neuroanatomical similarities to humans and the successful generation of transgenic lines. Stereotaxic procedures in marmosets are guided by 2D stereotaxic atlases, which are constructed with a limited number of animals and fail to account for inter-individual variability in skull and brain size. Here, we developed a frameless imaging-guided stereotaxic system that improves upon traditional approaches by using subject-specific registration of computed tomography (CT), magnetic resonance imaging (MRI) and positron emission tomography (PET) data to identify a surgical target, namely the putamen, in two marmosets. The skull surface was laser-scanned to create a point cloud that was registered to the 3D reconstruction of the skull from CT. Reconstruction of the skull, as well as of the brain from MR images, was crucial for surgical planning. Localisation and injection into the putamen was done using a 6-axis robotic arm controlled by a surgical navigation software (Brainsight™). Integration of subject-specific registration and frameless stereotaxic navigation allowed target localisation specific to each animal. Injection of alpha-synuclein fibrils into the putamen triggered progressive neurodegeneration of the nigro-striatal system, a key feature of Parkinson's disease. Four months post-surgery, a PET scan found evidence of nigro-striatal denervation, supporting accurate targeting of the putamen during co-registration and subsequent surgery. Our results suggest that this approach, coupled with frameless stereotaxic neuronavigation, is accurate in localising surgical targets and can be used to assess endpoints for longitudinal studies.

Motor sequences; separating the sequence from the motor. A longitudinal rsfMRI study

In motor learning, sequence specificity, i.e. the learning of specific sequential associations, has predominantly been studied using task-based fMRI paradigms. However, offline changes in resting state functional connectivity after sequence-specific motor learning are less well understood. Previous research has established that plastic changes following motor learning can be divided into stages including fast learning, slow learning and retention. A description of how resting state functional connectivity after sequence-specific motor sequence learning (MSL) develops across these stages is missing. This study aimed to identify plastic alterations in whole-brain functional connectivity after learning a complex motor sequence by contrasting an active group who learned a complex sequence with a control group who performed a control task matched for motor execution. Resting state fMRI and behavioural performance were collected in both groups over the course of 5 consecutive training days and at follow-up after 12 days to encompass fast learning, slow learning, overall learning and retention. Between-group interaction analyses showed sequence-specific decreases in functional connectivity during overall learning in the right supplementary motor area (SMA). We found that connectivity changes in a key region of the motor network, the superior parietal cortex (SPC) were not a result of sequence-specific learning but were instead linked to motor execution. Our study confirms the sequence-specific role of SMA that has previously been identified in online task-based learning studies, and extends it to resting state network changes after sequence-specific MSL.

Open science datasets from PREVENT-AD, a longitudinal cohort of pre-symptomatic Alzheimer's disease

To move Alzheimer Disease (AD) research forward it is essential to collect data from large cohorts, but also make such data available to the global research community. We describe the creation of an open science dataset from the PREVENT-AD (PResymptomatic EValuation of Experimental or Novel Treatments for AD) cohort, composed of cognitively unimpaired older individuals with a parental or multiple-sibling history of AD. From 2011 to 2017, 386 participants were enrolled (mean age 63 years old ± 5) for sustained investigation among whom 349 have retrospectively agreed to share their data openly. Repositories are findable through the unified interface of the Canadian Open Neuroscience Platform and contain up to five years of longitudinal imaging data, cerebral fluid biochemistry, neurosensory capacities, cognitive, genetic, and medical information. Imaging data can be accessed openly at https://openpreventad.loris.ca while most of the other information, sensitive by nature, is accessible by qualified researchers at https://registeredpreventad.loris.ca. In addition to being a living resource for continued data acquisition, PREVENT-AD offers opportunities to facilitate understanding of AD pathogenesis.

A model-based framework for correcting B 1 + inhomogeneity effects in magnetization transfer saturation and inhomogeneous magnetization transfer saturation maps

PURPOSE

In this work, we propose that Δ B 1 + -induced errors in magnetization transfer (MT) saturation (MT_sat ) maps can be corrected with use of an R₁ and B 1 + map and through numerical simulations of the sequence.

THEORY AND METHODS

One healthy subject was scanned at 3.0T using a partial quantitative MT protocol to estimate the relationship between observed R₁ (R_1,obs ) and apparent bound pool size ( M 0 , a p p B ) in the brain. MT_sat values were simulated for a range of B 1 + , R_1,obs , and M 0 , a p p B . An equation was fit to the simulated MT_sat , then a linear relationship between R_1,obs and M 0 , a p p B was generated. These results were used to generate correction factor maps for the MT_sat acquired from single-point data. The proposed correction was compared to an empirical correction factor with different MT-preparation schemes.

RESULTS

M 0 , a p p B was highly correlated with R_1,obs (r > 0.96), permitting the use of R_1,obs to estimate M 0 , a p p B for B 1 + correction. All B 1 + corrected MT_sat maps displayed a decreased correlation with B 1 + compared to uncorrected MT_sat and MT_sat corrected with an empirical factor in the corpus callosum. There was good agreement between the proposed approach and the empirical correction with radiofrequency saturation at 2 kHz, with larger deviations seen when using saturation pulses further off-resonance and in inhomogeneous (ih) MT_sat maps.

CONCLUSION

The proposed correction decreases the dependence of MT_sat on B 1 + inhomogeneities. Furthermore, this flexible framework permits the use of different saturation protocols, making it useful for correcting B 1 + inhomogeneities in ihMT.

Hippocampal shape across the healthy lifespan and its relationship with cognition

The study of the hippocampus across the healthy adult lifespan has rendered inconsistent findings. While volumetric measurements have often been a popular technique for analysis, more advanced morphometric techniques have demonstrated compelling results that highlight the importance and improved specificity of shape-based measures. Here, the MAGeT Brain algorithm was applied on 134 healthy individuals aged 18-81 years old to extract hippocampal subfield volumes and hippocampal shape measurements, namely: local surface area (SA) and displacement. We used linear-, second- or third-order natural splines to examine the relationships between hippocampal measures and age. In addition, partial least squares analyses were performed to relate volume and shape measurements with cognitive and demographic information. Volumetric results indicated a relative preservation of the right cornus ammonis 1 with age and a global volume reduction linked with older age, female sex, lower levels of education and cognitive performance. Vertex-wise analysis demonstrated an SA preservation in the anterior hippocampus with a peak during the sixth decade, while the posterior hippocampal SA gradually decreased across lifespan. Overall, SA decrease was linked to older age, female sex and, to a lesser extent lower levels of education and cognitive performance. Outward displacement in the lateral hippocampus and inward displacement in the medial hippocampus were enlarged with older age, lower levels of cognition and education, indicating an accentuation of the hippocampal "C" shape with age. Taken together, our findings suggest that vertex-wise analyses have higher spatial specifity and that sex, education, and cognition are implicated in the differential impact of age on hippocampal subregions throughout its anteroposterior and medial-lateral axes. This article is part of the Virtual Special Issue titled COGNITIVE NEU- ROSCIENCE OF HEALTHY AND PATHOLOGICAL AGING. The full issue can be found on ScienceDirect at https://www.sciencedirect.com/journal/neurobiology-of-aging/special-issue/105379XPWJP.

White matter microstructural changes in short-term learning of a continuous visuomotor sequence

Efficient neural transmission is crucial for optimal brain function, yet the plastic potential of white matter (WM) has long been overlooked. Growing evidence now shows that modifications to axons and myelin occur not only as a result of long-term learning, but also after short training periods. Motor sequence learning (MSL), a common paradigm used to study neuroplasticity, occurs in overlapping learning stages and different neural circuits are involved in each stage. However, most studies investigating short-term WM plasticity have used a pre-post design, in which the temporal dynamics of changes across learning stages cannot be assessed. In this study, we used multiple magnetic resonance imaging (MRI) scans at 7 T to investigate changes in WM in a group learning a complex visuomotor sequence (LRN) and in a control group (SMP) performing a simple sequence, for five consecutive days. Consistent with behavioral results, where most improvements occurred between the two first days, structural changes in WM were observed only in the early phase of learning (d1-d2), and in overall learning (d1-d5). In LRNs, WM microstructure was altered in the tracts underlying the primary motor and sensorimotor cortices. Moreover, our structural findings in WM were related to changes in functional connectivity, assessed with resting-state functional MRI data in the same cohort, through analyses in regions of interest (ROIs). Significant changes in WM microstructure were found in a ROI underlying the right supplementary motor area. Together, our findings provide evidence for highly dynamic WM plasticity in the sensorimotor network during short-term MSL.

Efficient whole-brain tract-specific T1 mapping at 3T with slice-shuffled inversion-recovery diffusion-weighted imaging

PURPOSE

Most voxels in white matter contain multiple fiber populations with different orientations and levels of myelination. Conventional T₁ mapping measures 1 T₁ value per voxel, representing a weighted average of the multiple tract T₁ times. Inversion-recovery diffusion-weighted imaging (IR-DWI) allows the T₁ times of multiple tracts in a voxel to be disentangled, but the scan time is prohibitively long. Recently, slice-shuffled IR-DWI implementations have been proposed to significantly reduce scan time. In this work, we demonstrate that we can measure tract-specific T₁ values in the whole brain using simultaneous multi-slice slice-shuffled IR-DWI at 3T.

METHODS

We perform simulations to evaluate the accuracy and precision of our crossing fiber IR-DWI signal model for various fiber parameters. The proposed sequence and signal model are tested in a phantom consisting of crossing asparagus pieces doped with gadolinium to vary T₁ , and in 2 human subjects.

RESULTS

Our simulations show that tract-specific T₁ times can be estimated within 5% of the nominal fiber T₁ values. Tract-specific T₁ values were resolved in subvoxel 2 fiber crossings in the asparagus phantom. Tract-specific T₁ times were resolved in 2 different tract crossings in the human brain where myelination differences have previously been reported; the crossing of the cingulum and genu of the corpus callosum and the crossing of the corticospinal tract and pontine fibers.

CONCLUSION

Whole-brain tract-specific T₁ mapping is feasible using slice-shuffled IR-DWI at 3T. This technique has the potential to improve the microstructural characterization of specific tracts implicated in neurodevelopment, aging, and demyelinating disorders.

A novel ex vivo, in situ method to study the human brain through MRI and histology

BACKGROUND

MRI-histology correlation studies of the ex vivo brain mostly employ fresh, extracted (ex situ) specimens, aldehyde fixed by immersion, which has several disadvantages for MRI scanning (e.g. deformation of the organ). A minority of studies are done ex vivo-in situ (unfixed brain), requiring an MRI scanner readily available within a few hours of the time of death.

NEW METHOD

We propose a new technique, exploited by anatomists, for scanning the ex vivo brain: fixation by whole body perfusion, which implies fixation of the brain in situ. This allows scanning the brain surrounded by fluids, meninges, and skull, preserving the structural relationships of the brain in vivo. To evaluate the proposed method, five heads perfused-fixed with a saturated sodium chloride solution were employed. Three sequences were acquired on a 1.5 T MRI scanner: T1weighted, T2weighted-FLAIR, and Gradient-echo. Histology analysis included immunofluorescence for myelin basic protein and neuronal nuclei.

RESULTS

All MRIs were successfully processed through a validated pipeline used with in vivo MRIs. All cases exhibited positive antigenicity for myelin and neuronal nuclei.

COMPARISON WITH EXISTING METHODS

All scans registered to a standard neuroanatomical template in pseudo-Talairach space more accurately than an ex vivo-ex situ scan. The time interval to scan the ex vivo brain in situ was increased to at least 10 months.

CONCLUSIONS

MRI and histology study of the ex vivo-in situ brain fixed by perfusion is an alternative approach that has important procedural and practical advantages over the two standard methods to study the ex vivo brain.

Altered hippocampal centrality and dynamic anatomical covariance of intracortical microstructure in first episode psychosis

Hippocampal circuitry has been posited to be fundamental to positive symptoms in psychosis, but its contributions to other factors important for outcome remains unclear. We hypothesized that longitudinal changes in the hippocampal circuit and concomitant changes of intracortical microstructure are altered in first episode psychosis (FEP) patients and that such changes are associated with negative symptoms and verbal memory. Longitudinal brain scans (2-4 visits over 3-15 months) were acquired for 27 FEP and 29 age- and sex-matched healthy controls. Quantitative T1 maps, sensitive to myelin content, were used to sample the microstructure of the hippocampal subfields and output circuitry (fimbria, alveus, fornix, mammillary bodies), and intracortical regions. Dynamic anatomical covariance in pair-wise regional trajectories were assessed for each subject, and graph theory was used to calculate a participation coefficient metric that quantifies the similarity/divergence between hippocampal and intracortical microstructure. The mean participation coefficient of the hippocampus was significantly reduced in FEP patients compared with controls, reflecting differences in output hippocampal regions. Importantly, lower participation coefficient of the hippocampal circuit was associated with worse negative symptoms, a relationship that was mediated by changes in verbal memory. This study provides evidence for reduced hippocampal centrality in FEP and concomitant changes in intracortical anatomy. Myelin-rich output regions of the hippocampus may be an important biological trigger in early psychosis, with cascading effects on broader cortical networks and resultant clinical profiles.

MR-based age-related effects on the striatum, globus pallidus, and thalamus in healthy individuals across the adult lifespan

While numerous studies have used magnetic resonance imaging (MRI) to elucidate normative age‐related trajectories in subcortical structures across the human lifespan, there exists substantial heterogeneity among different studies. Here, we investigated the normative relationships between age and morphology (i.e., volume and shape), and microstructure (using the T1‐weighted/T2‐weighted [T1w/T2w] signal ratio as a putative index of myelin and microstructure) of the striatum, globus pallidus, and thalamus across the adult lifespan using a dataset carefully quality controlled, yielding a final sample of 178 for the morphological analyses, and 162 for the T1w/T2w analyses from an initial dataset of 253 healthy subjects, aged 18–83. In accordance with previous cross‐sectional studies of adults, we observed age‐related volume decrease that followed a quadratic relationship between age and bilateral striatal and thalamic volumes, and a linear relationship in the globus pallidus. Our shape indices consistently demonstrated age‐related posterior and medial areal contraction bilaterally across all three structures. Beyond morphology, we observed a quadratic inverted U‐shaped relationship between T1w/T2w signal ratio and age, with a peak value occurring in middle age (at around 50 years old). After permutation testing, the Akaike information criterion determined age relationships remained significant for the bilateral globus pallidus and thalamus, for both the volumetric and T1w/T2w analyses. Our findings serve to strengthen and expand upon previous volumetric analyses by providing a normative baseline of morphology and microstructure of these structures to which future studies investigating patients with various disorders can be compared.

High resolution atlas of the venous brain vasculature from 7 T quantitative susceptibility maps

The vascular organization of the human brain can determine neurological and neurophysiological functions, yet thus far it has not been comprehensively mapped. Aging and diseases such as dementia are known to be associated with changes to the vasculature and normative data could help detect these vascular changes in neuroimaging studies. Furthermore, given the well-known impact of venous vessels on the blood oxygen level dependent (BOLD) signal, information about the common location of veins could help detect biases in existing datasets. In this work, a quantitative atlas of the venous vasculature using quantitative susceptibility maps (QSM) acquired with a 0.6-mm isotropic resolution is presented. The Venous Neuroanatomy (VENAT) atlas was created from 5 repeated 7 Tesla MRI measurements in young and healthy volunteers (n = 20, 10 females, mean age = 25.1 ± 2.5 years) using a two-step registration method on 3D segmentations of the venous vasculature. This cerebral vein atlas includes the average vessel location, diameter (mean: 0.84 ± 0.33 mm) and curvature (0.11 ± 0.05 mm^-1) from all participants and provides an in vivo measure of the angio-architectonic organization of the human brain and its variability. This atlas can be used as a basis to understand changes in the vasculature during aging and neurodegeneration, as well as vascular and physiological effects in neuroimaging.

Physical Exercise and Spatial Training: A Longitudinal Study of Effects on Cognition, Growth Factors, and Hippocampal Plasticity

Physical exercise has been suggested to improve cognitive performance through various neurobiological mechanisms, mediated by growth factors such as BDNF, IGF-I, and VEGF. Moreover, animal research has demonstrated that combined physical and cognitive stimulation leads to increased adult neurogenesis as compared to either experimental condition alone. In the present study, we therefore investigated whether a sequential combination of physical and spatial training in young, healthy adults elicits an additive effect on training and transfer gains. To this end, we compared the effects of (i) eight 20-minute sessions of cycling, (ii) sixteen 30-minute sessions of spatial training, (iii) a combination of both, and included (iv) a passive control cohort. We assessed longitudinal changes in cognitive performance, growth factor levels, and T₁ relaxation of hippocampal subfields (acquired with 7 T MRI). While substantial physical and spatial training gains were elicited in all trained groups, longitudinal transfer changes did not differ between these groups. Notably, we found no evidence for an additive effect of sequential physical and spatial training. These results challenge the extrapolation from the findings reported in animals to young, healthy adults.

Regionally specific changes in the hippocampal circuitry accompany progression of cerebrospinal fluid biomarkers in preclinical Alzheimer's disease

Neuropathological and in vivo brain imaging studies agree that the cornu ammonis 1 and subiculum subfields of the hippocampus are most vulnerable to atrophy in the prodromal phases of Alzheimer's disease (AD). However, there has been limited investigation of the structural integrity of the components of the hippocampal circuit, including subfields and extra-hippocampal white matter structure, in relation to the progression of well-accepted cerebrospinal fluid (CSF) biomarkers of AD, amyloid-β 1-42 (Aβ) and total-tau (tau). We investigated these relationships in 88 aging asymptomatic individuals with a parental or multiple-sibling familial history of AD. Apolipoprotein (APOE) ɛ4 risk allele carriers were identified, and all participants underwent cognitive testing, structural magnetic resonance imaging, and lumbar puncture for CSF assays of tau, phosphorylated-tau (p-tau) and Aβ. Individuals with a reduction in CSF Aβ levels (an indicator of amyloid accretion into neuritic plaques) as well as evident tau pathology (believed to be linked to neurodegeneration) exhibited lower subiculum volume, lower fornix microstructural integrity, and a trend towards lower cognitive score than individuals who showed only reduction in CSF Aβ. In contrast, persons with normal levels of tau showed an increase in structural MR markers in relation to declining levels of CSF Aβ. These results suggest that hippocampal subfield volume and extra-hippocampal white matter microstructure demonstrate a complex pattern where an initial volume increase is followed by decline among asymptomatic individuals who, in some instances, may be a decade or more away from onset of cognitive or functional impairment.

Age-related mapping of intracortical myelin from late adolescence to middle adulthood using T1 -weighted MRI

Magnetic resonance imaging (MRI) studies in humans have reported that the T₁ -weighted signal in the cerebral cortex follows an inverted "U" trajectory over the lifespan. Here, we investigated the T₁ -weighted signal trajectory from late adolescence to middle adulthood in humans to characterize the age range when mental illnesses tend to present, and efficacy of treatments are evaluated. We compared linear to quadratic predictors of age on signal in 67 healthy individuals, 17-45 years old. We investigated ¼, ½, and ¾ depths in the cortex representing intracortical myelin (ICM), in the superficial white matter (SWM), and in a reference deep white matter tract. We found that the quadratic fit was superior in all regions of the cortex, while signal in the SWM and deep white matter showed no global dependence on age over this range. The signal trajectory in any region followed a similar shape regardless of cortical depth. The quadratic fit was analyzed in 70 cortical regions to obtain the age of maximum signal intensity. We found that visual, cingulate, and left ventromedial prefrontal cortices peak first around 34 years old, whereas motor and premotor areas peak latest at ∼38 years. Our analysis suggests that ICM trajectories over this range can be modeled well in small cohorts of subjects using quadratic functions, which are amenable to statistical analysis, thus suitable for investigating regional changes in ICM with disease. This study highlights a novel approach to map ICM trajectories using an age range that coincides with the onset of many mental illnesses. Hum Brain Mapp 38:3691-3703, 2017. © 2017 Wiley Periodicals, Inc.

B1 mapping for bias-correction in quantitative T1 imaging of the brain at 3T using standard pulse sequences

PURPOSE

B₁ mapping is important for many quantitative imaging protocols, particularly those that include whole-brain T₁ mapping using the variable flip angle (VFA) technique. However, B₁ mapping sequences are not typically available on many magnetic resonance imaging (MRI) scanners. The aim of this work was to demonstrate that B₁ mapping implemented using standard scanner product pulse sequences can produce B₁ (and VFA T₁ ) maps comparable in quality and acquisition time to advanced techniques.

MATERIALS AND METHODS

Six healthy subjects were scanned at 3.0T. An interleaved multislice spin-echo echo planar imaging double-angle (EPI-DA) B₁ mapping protocol, using a standard product pulse sequence, was compared to two alternative methods (actual flip angle imaging, AFI, and Bloch-Siegert shift, BS). Single-slice spin-echo DA B₁ maps were used as a reference for comparison (Ref. DA). VFA flip angles were scaled using each B₁ map prior to fitting T₁ ; the nominal flip angle case was also compared.

RESULTS

The pooled-subject voxelwise correlation (ρ) for B₁ maps (BS/AFI/EPI-DA) relative to the reference B₁ scan (Ref. DA) were ρ = 0.92/0.95/0.98. VFA T₁ correlations using these maps were ρ = 0.86/0.88/0.96, much better than without B₁ correction (ρ = 0.53). The relative error for each B₁ map (BS/AFI/EPI-DA/Nominal) had 95^th percentiles of 5/4/3/13%.

CONCLUSION

Our findings show that B₁ mapping implemented using product pulse sequences can provide excellent quality B₁ (and VFA T₁ ) maps, comparable to other custom techniques. This fast whole-brain measurement (∼2 min) can serve as an excellent alternative for researchers without access to advanced B₁ pulse sequences.

LEVEL OF EVIDENCE

1 Technical Efficacy: Stage 1 J. Magn. Reson. Imaging 2017;46:1673-1682.

A Systematic Relationship Between Functional Connectivity and Intracortical Myelin in the Human Cerebral Cortex

Research in the macaque monkey suggests that cortical areas with similar microstructure are more likely to be connected. Here, we examine this link in the human cerebral cortex using 2 magnetic resonance imaging (MRI) measures: quantitative  T1 maps, which are sensitive to intracortical myelin content and provide an in vivo proxy for cortical microstructure, and resting-state functional connectivity. Using ultrahigh-resolution MRI at 7 T and dedicated image processing tools, we demonstrate a systematic relationship between T1-based intracortical myelin content and functional connectivity. This effect is independent of the proximity of areas. We employ nonlinear dimensionality reduction to characterize connectivity components and identify specific aspects of functional connectivity that are linked to myelin content. Our results reveal a consistent spatial pattern throughout different analytic approaches. While functional connectivity and myelin content are closely linked in unimodal areas, the correspondence is lower in transmodal areas, especially in posteromedial cortex and the angular gyrus. Our findings are in agreement with comprehensive reports linking histologically assessed microstructure and connectivity in different mammalian species and extend them to the human cerebral cortex in vivo.

Assessing intracortical myelin in the living human brain using myelinated cortical thickness

Alterations in the myelination of the cerebral cortex may underlie abnormal cortical function in a variety of brain diseases. Here, we describe a technique for investigating changes in intracortical myelin in clinical populations on the basis of cortical thickness measurements with magnetic resonance imaging (MRI) at 3 Tesla. For this, we separately compute the thickness of the shallower, lightly myelinated portion of the cortex and its deeper, heavily myelinated portion (referred to herein as unmyelinated and myelinated cortex, respectively). Our expectation is that the thickness of the myelinated cortex will be a specific biomarker for disruptions in myeloarchitecture. We show representative atlases of total cortical thickness, T, unmyelinated cortical thickness, G, and myelinated cortical thickness, M, for a healthy group of 20 female subjects. We further demonstrate myelinated cortical thickness measurements in a preliminary clinical study of 10 bipolar disorder type-I subjects and 10 healthy controls, and report significant decreases in the middle frontal gyrus in T, G, and M in the disorder, with the largest percentage change occurring in M. This study highlights the potential of myelinated cortical thickness measurements for investigating intracortical myelin involvement in brain disease at clinically relevant field strengths and resolutions.

On the accuracy of T1 mapping: searching for common ground

PURPOSE

There are many T1 mapping methods available, each of them validated in phantoms and reporting excellent agreement with literature. However, values in literature vary greatly, with T1 in white matter ranging from 690 to 1100 ms at 3 Tesla. This brings into question the accuracy of one of the most fundamental measurements in quantitative MRI. Our goal was to explain these variations and look into ways of mitigating them.

THEORY AND METHODS

We evaluated the three most common T1 mapping methods (inversion recovery, Look-Locker, and variable flip angle) through Bloch simulations, a white matter phantom and the brains of 10 healthy subjects (single-slice). We pooled the T1 histograms of the subjects to determine whether there is a sequence-dependent bias and whether it is reproducible across subjects.

RESULTS

We found good agreement between the three methods in phantoms, but poor agreement in vivo, with the white matter T1 histogram peak in healthy subjects varying by more than 30% depending on the method used. We also found that the pooled brain histograms displayed three distinct white matter peaks, with Look-Locker consistently underestimating, and variable flip angle overestimating the inversion recovery T1 values. The Bloch simulations indicated that incomplete spoiling and inaccurate B1 mapping could account for the observed differences.

CONCLUSION

We conclude that the three most common T1 mapping protocols produce stable T1 values in phantoms, but not in vivo. To improve the accuracy of T1 mapping, we recommend that sites perform in vivo validation of their T1 mapping method against the inversion recovery reference method, as the first step toward developing a robust calibration scheme.

Regional impact of field strength on voxel-based morphometry results

The objective of this study was to characterize the sensitivity of voxel-based morphometry (VBM) results to choice field strength. We chose to investigate the two most widespread acquisition sequences for VBM, FLASH and MP-RAGE, at 1.5 and 3 T. We first evaluated image quality of the four acquisition protocols in terms of SNR and image uniformity. We then performed a VBM study on eight subjects scanned twice using the four protocols to evaluate differences in grey matter (GM) density and corresponding scan-rescan variability, and a power analysis for each protocol in the context a longitudinal and cross-sectional VBM study. As expected, the SNR increased significantly at 3 T for both FLASH and MP-RAGE. Image non-uniformity increased as well, in particular for MP-RAGE. The differences in CNR and contrast non-uniformity cause regional biases between protocols in the VBM results, in particular between sequences at 3 T. The power analysis results show an overall decrease in the number of subjects required in a longitudinal study to detect a difference in GM density at 3 T for MP-RAGE, but an increase for FLASH. The number of subjects required in a cross-sectional VBM study is higher at 3 T for both sequences. Our results show that each protocol has a distinct regional sensitivity pattern to morphometric change, which goes against the classical view of VBM as an unbiased whole brain analysis technique, complicates the combination of data within a VBM study and the direct comparison of VBM studies based on different protocols.

Segmentation of cortical MS lesions on MRI using automated laminar profile shape analysis

Cortical multiple sclerosis lesions are difficult to detect in magnetic resonance images due to poor contrast with surrounding grey matter, spatial variation in healthy grey matter and partial volume effects. We propose using an observer-independent laminar profile-based parcellation method to detect cortical lesions. Following cortical surface extraction, profiles are extended from the white matter surface to the grey matter surface. The cortex is parcellated according to profile intensity and shape features using a k-means classifier. The method is applied to a high-resolution quantitative magnetic resonance data set from a fixed post mortem multiple sclerosis brain, and validated using histology.