The microstructural correlates of T1 in white matter

T1 mapping in patients with cervical spinal canal stenosis with and without decompressive surgery: A longitudinal study

BACKGROUND AND PURPOSE

Cervical spinal canal stenosis (cSCS) is a common cause of spinal impairment in the elderly. With conventional magnetic resonance imaging (MRI) suffering from various limitations, high-resolution single-shot T1 mapping has been proposed as a novel MRI technique in cSCS diagnosis. In this study, we investigated the effect of conservative and surgical treatment on spinal cord T1 relaxation times in cSCS.

METHODS

T1-mapping was performed in 54 patients with cSCS at 3 Tesla MRI at the maximum-, above and below the stenosis. Subsequently, intraindividual T1-differences (ΔT1) intrastenosis were calculated. Twenty-four patients received follow-up scans after 6 months.

RESULTS

Surgically treated patients showed higher ΔT1 at baseline (154.9 ± 81.6 vs. 95.3 ± 60.7), while absolute T1-values within the stenosis were comparable between groups (863.7 ± 89.3 milliseconds vs. 855.1 ± 62.2 milliseconds). In surgically treated patients, ΔT1 decreased inverse to stenosis severity. After 6 months, ΔT1 significantly decreased in the surgical group (154.9 ± 81.6 milliseconds to 85.7 ± 108.9 milliseconds, p = .021) and remained unchanged in conservatively treated patients. Both groups showed clinical improvement at the 6-month follow-up.

CONCLUSIONS

Baseline difference of T1 relaxation time (ΔT1) might serve as a supporting marker for treatment decision and change of T1 relaxation time might reflect relief of spinal cord narrowing indicating regenerative processes. Quantitative T1-mapping represents a promising additional imaging method to indicate a surgical treatment plan and to validate treatment success.

Estimating axon radius using diffusion-relaxation MRI: calibrating a surface-based relaxation model with histology

Axon radius is a potential biomarker for brain diseases and a crucial tissue microstructure parameter that determines the speed of action potentials. Diffusion MRI (dMRI) allows non-invasive estimation of axon radius, but accurately estimating the radius of axons in the human brain is challenging. Most axons in the brain have a radius below one micrometer, which falls below the sensitivity limit of dMRI signals even when using the most advanced human MRI scanners. Therefore, new MRI methods that are sensitive to small axon radii are needed. In this proof-of-concept investigation, we examine whether a surface-based axonal relaxation process could mediate a relationship between intra-axonal T2 and T1 times and inner axon radius, as measured using postmortem histology. A unique in vivo human diffusion-T1-T2 relaxation dataset was acquired on a 3T MRI scanner with ultra-strong diffusion gradients, using a strong diffusion-weighting (i.e., b = 6,000 s/mm2) and multiple inversion and echo times. A second reduced diffusion-T2 dataset was collected at various echo times to evaluate the model further. The intra-axonal relaxation times were estimated by fitting a diffusion-relaxation model to the orientation-averaged spherical mean signals. Our analysis revealed that the proposed surface-based relaxation model effectively explains the relationship between the estimated relaxation times and the histological axon radius measured in various corpus callosum regions. Using these histological values, we developed a novel calibration approach to predict axon radius in other areas of the corpus callosum. Notably, the predicted radii and those determined from histological measurements were in close agreement.

Axon fiber orientation as the source of T1 relaxation anisotropy in white matter: A study on corpus callosum in vivo and ex vivo

PURPOSE

Recent studies indicate that T₁ in white matter (WM) is influenced by fiber orientation in B₀ . The purpose of the study was to investigate the interrelationships between axon fiber orientation in corpus callosum (CC) and T₁ relaxation time in humans in vivo as well as in rat brain ex vivo.

METHODS

Volunteers were scanned for relaxometric and diffusion MRI at 3 T and 7 T. Angular T₁ plots from WM were computed using fractional anisotropy and fiber-to-field-angle maps. T₁ and fiber-to-field angle were measured in five sections of CC to estimate the effects of inherently varying fiber orientations on T₁ within the same tracts in vivo. Ex vivo rat-brain preparation encompassing posterior CC was rotated in B₀ and T₁ , and diffusion MRI images acquired at 9.4 T. T₁ angular plots were determined at several rotation angles in B₀ .

RESULTS

Angular T₁ plots from global WM provided reference for estimated fiber orientation-linked T₁ changes within CC. In anterior midbody of CC in vivo, where small axons are dominantly present, a shift in axon orientation is accompanied by a change in T₁ , matching that estimated from WM T₁ data. In CC, where large and giant axons are numerous, the measured T₁ change is about 2-fold greater than the estimated one. Ex vivo rotation of the same midsagittal CC region of interest produced angular T₁ plots at 9.4 T, matching those observed at 7 T in vivo.

CONCLUSION

These data causally link axon fiber orientation in B₀ to the T₁ relaxation anisotropy in WM.

Nonaqueous magnetization following adiabatic and selective pulses in brain: T1 and cross-relaxation dynamics

Inversion pulses are commonly employed in MRI for T 1 $$ {T}_1 $$ -weighted contrast and relaxation measurements. In the brain, it is often assumed that adiabatic pulses saturate the nonaqueous magnetization. We investigated this assumption using solid-state NMR to monitor the nonaqueous signal directly following adiabatic inversion and compared this with signals following hard and soft inversion pulses. The effects of the different preparations on relaxation dynamics were explored. Inversion recovery experiments were performed on ex vivo bovine and porcine brains using 360-MHz (8.4 T) and 200-MHz (4.7 T) NMR spectrometers, respectively, using broadband rectangular, adiabatic, and sinc inversion pulses as well as a long rectangular saturation pulse. Analogous human brain MRI experiments were performed at 3 T using single-slice echo-planar imaging. Relaxation data were fitted by mono- and biexponential decay models. Further fitting analysis was performed using only two inversion delay times. Adiabatic and sinc inversion left much of the nonaqueous magnetization along B 0 $$ {B}_0 $$ and resulted in biexponential relaxation. Saturation of both aqueous and nonaqueous magnetization components led to effectively monoexponential T 1 $$ {T}_1 $$ relaxation. Typical adiabatic inversion pulses do not, as has been widely assumed, saturate the nonaqueous proton magnetization in white matter. Unequal magnetization states in aqueous and nonaqueous ¹ H reservoirs prepared by soft and adiabatic pulses result in biexponential T 1 $$ {T}_1 $$ relaxation. Both pools must be prepared in the same magnetization state (e.g., saturated or inverted) in order to observe consistent monoexponential relaxation.

White matter microstructure and longitudinal relaxation time anisotropy in human brain at 3 and 7 T

A high degree of structural order by white matter (WM) fibre tracts creates a physicochemical environment where water relaxations are rendered anisotropic. Recently, angularly dependent longitudinal relaxation has been reported in human WM. We have characterised interrelationships between T1 relaxation and diffusion MRI microstructural indices at 3 and 7 T. Eleven volunteers consented to participate in the study. Multishell diffusion MR images were acquired with b-values of 0/1500/3000 and 0/1000/2000 s/mm² at 1.5 and 1.05 mm³ isotropic resolutions at 3 and 7 T, respectively. DTIFIT was used to compute DTI indices; the fibre-to-field angle (θ_FB ) maps were obtained using the principal eigenvector images. The orientations and volume fractions of multiple fibre populations were estimated using BedpostX in FSL, and the orientation dispersion index (ODI) was estimated using the NODDI protocol. MP2RAGE was used to acquire images for T1 maps at 1.0 and 0.9 mm³ isotropic resolutions at 3 and 7 T, respectively. At 3 T, T1 as a function of θ_FB in WM with high fractional anisotropy and one-fibre orientation volume fraction or low ODI shows a broad peak centred at 50^o , but a flat baseline at 0^o and 90^o . The broad peak amounted up to 7% of the mean T1. At 7 T, the broad peak appeared at 40^o and T1 in fibres running parallel to B0 was longer by up to 75 ms (8.3% of the mean T1) than in those perpendicular to the field. The peak at 40^o was approximately 5% of mean T1 (i.e., proportionally smaller than that at 54^o at 3 T). The data demonstrate T1 anisotropy in WM with high microstructural order at both fields. The angular patterns are indicative of the B0-dependency of T1 anisotropy. Thus myelinated WM fibres influence T1 contrast both by acting as a T1 contrast agent and rendering T1 dependent on fibre orientation with B0.

Quantitative evaluation of the spinal cord compression in patients with cervical spondylotic myelopathy using synthetic MRI

Objectives: This work aimed to investigate the feasibility and diagnostic value of synthetic MRI, including T1, T2 and PD values in determining the severity of cervical spondylotic myelopathy (CSM). Methods: All subjects (51 CSM patients and 9 healthy controls) underwent synthetic MRI scan on a 3.0T GE MR scanner. The cervical canal stenosis degree of subjects was graded 0-III based on the method of a MRI grading system. Regions of interest (ROIs) were manually drawn at the maximal compression level (MCL) by covering the whole spinal cord to generate T1_MCL, T2_MCL, and PD_MCL values in grade I-III groups. Besides, anteroposterior (AP) and transverse (Trans) diameters of the spinal cord at MCL were measured in grade II and grade III groups, and relative values were calculated as follows: rAP = AP_MCL/AP_normal, rTrans = Trans_MCL/Trans_normal. rMIN = rAP/rTrans. Results: T1_MCL value showed a decreasing trend with severity of grades (from grade 0 to grade II, p < 0.05), while it increased dramatically at grade III. T2_MCL value showed no significant difference among grade groups (from grade 0 to grade II), while it increased dramatically at grade III compared to grade II (p < 0.05). PD_MCL value showed no statistical difference among all grade groups. rMIN of grade III was significantly lower than that of grade II (p < 0.05). T2_MCL value was negatively correlated with rMIN, whereas positively correlated with rTrans. Conclusion: Synthetic MRI can provide not only multiple contrast images but also quantitative mapping, which is showed promisingly to be a reliable and efficient method in the quantitative diagnosis of CSM.

Decreased basal ganglia and thalamic iron in early psychotic spectrum disorders are associated with increased psychotic and schizotypal symptoms

Iron deficits have been reported as a risk factor for psychotic spectrum disorders (PSD). However, examinations of brain iron in PSD remain limited. The current study employed quantitative MRI to examine iron content in several iron-rich subcortical structures in 49 young adult individuals with PSD (15 schizophrenia, 17 schizoaffective disorder, and 17 bipolar disorder with psychotic features) compared with 35 age-matched healthy controls (HC). A parametric approach based on a two-pool magnetization transfer model was applied to estimate longitudinal relaxation rate (R₁), which reflects both iron and myelin, and macromolecular proton fraction (MPF), which is specific to myelin. To describe iron content, a synthetic effective transverse relaxation rate (R₂*) was modeled using a linear fitting of R₁ and MPF. PSD patients compared to HC showed significantly reduced R₁ and synthetic R₂* across examined regions including the pallidum, ventral diencephalon, thalamus, and putamen areas. This finding was primarily driven by decreases in the subgroup with schizophrenia, followed by schizoaffective disorder. No significant group differences were noted for MPF between PSD and HC while for regional volume, significant reductions in patients were only observed in bilateral caudate, suggesting that R₁ and synthetic R₂* reductions in schizophrenia and schizoaffective patients likely reflect iron deficits that either occur independently or precede structural and myelin changes. Subcortical R₁ and synthetic R₂* were also found to be inversely related to positive symptoms within the PSD group and to schizotypal traits across the whole sample. These findings that decreased iron in subcortical regions are associated with PSD risk and symptomatology suggest that brain iron deficiencies may play a role in PSD pathology and warrant further study.

Cortical myelin profile variations in healthy aging brain: A T1w/T2w ratio study

Demyelination is observed in both healthy aging and age-related neurodegenerative disorders. While the significance of myelin within the cortex is well acknowledged, studies focused on intracortical demyelination and depth-specific structural alterations in normal aging are lacking. Using the recently available Human Connectome Project Aging dataset, we investigated intracortical myelin in a normal aging population using the T1w/T2w ratio. To capture the fine changes across cortical depths, we employed a surface-based approach by constructing cortical profiles traveling perpendicularly through the cortical ribbon and sampling T1w/T2w values. The curvatures of T1w/T2w cortical profiles may be influenced by differences in local myeloarchitecture and other tissue properties, which are known to vary across cortical regions. To quantify the shape of these profiles, we parametrized the level of curvature using a nonlinearity index (NLI) that measures the deviation of the profile from a straight line. We showed that NLI exhibited a steep decline in aging that was independent of local cortical thinning. Further examination of the profiles revealed that lower T1w/T2w near the gray-white matter boundary and superficial cortical depths were major contributors to the apparent NLI variations with age. These findings suggest that demyelination and changes in other T1w/T2w related tissue properties in normal aging may be depth-specific and highlight the potential of NLI as a unique marker of microstructural alterations within the cerebral cortex.

Empirical transmit field bias correction of T1w/T2w myelin maps

T1-weighted divided by T2-weighted (T1w/T2w) myelin maps were initially developed for neuroanatomical analyses such as identifying cortical areas, but they are increasingly used in statistical comparisons across individuals and groups with other variables of interest. Existing T1w/T2w myelin maps contain radiofrequency transmit field (B1+) biases, which may be correlated with these variables of interest, leading to potentially spurious results. Here we propose two empirical methods for correcting these transmit field biases using either explicit measures of the transmit field or alternatively a 'pseudo-transmit' approach that is highly correlated with the transmit field at 3T. We find that the resulting corrected T1w/T2w myelin maps are both better neuroanatomical measures (e.g., for use in cross-species comparisons), and more appropriate for statistical comparisons of relative T1w/T2w differences across individuals and groups (e.g., sex, age, or body-mass-index) within a consistently acquired study at 3T. We recommend that investigators who use the T1w/T2w approach for mapping cortical myelin use these B1+ transmit field corrected myelin maps going forward.

Individual variability in the nonlinear development of the corpus callosum during infancy and toddlerhood: a longitudinal MRI analysis

The human brain spends several years bootstrapping itself through intrinsic and extrinsic modulation, thus gradually developing both spatial organization and functions. Based on previous studies on developmental patterns and inter-individual variability of the corpus callosum (CC), we hypothesized that inherent variations of CC shape among infants emerge, depending on the position within the CC, along the developmental timeline. Here we used longitudinal magnetic resonance imaging data from infancy to toddlerhood and investigated the area, thickness, and shape of the midsagittal plane of the CC by applying multilevel modeling. The shape characteristics were extracted using the Procrustes method. We found nonlinearity, region-dependency, and inter-individual variability, as well as intra-individual consistencies, in CC development. Overall, the growth rate is faster in the first year than in the second year, and the trajectory differs between infants; the direction of CC formation in individual infants was determined within six months and maintained to two years. The anterior and posterior subregions increase in area and thickness faster than other subregions. Moreover, we clarified that the growth rate of the middle part of the CC is faster in the second year than in the first year in some individuals. Since the division of regions exhibiting different tendencies coincides with previously reported divisions based on the diameter of axons that make up the region, our results suggest that subregion-dependent individual variability occurs due to the increase in the diameter of the axon caliber, myelination partly due to experience and axon elimination during the early developmental period.

Spin Lattice (T1) and Magnetization Transfer Saturation (MTsat) Imaging to Monitor Age-Related Differences in Skeletal Muscle Tissue

Background: The aim was to compare spin-lattice relaxation (T1) mapping from sequences with no fat suppression and three fat suppression methods and Magnetization Transfer Saturation (MT_sat) mapping, to identify regional and age-related differences in calf muscle. These differences may be of clinical significance in age-related loss of muscle force. Methods: Ten young and seven senior subjects were imaged on a 3T MRI scanner using a 3D Fast Low Angle Shot sequence without and with different fat suppression and with MT saturation pulse. Bland–Altman plots were used to assess T1 maps using the fat unsuppressed sequence as the reference image. Age and regional differences in T1 and in MT_sat were assessed using two-way factorial analyses of variance (ANOVAs) with Bonferroni-adjusted independent sample t-tests for post hoc analyses. Results: A significant age-related increase in T1 and decrease in MT_sat was seen in the calf muscles. The largest size effect was observed in the T1 sequence with fat saturation. Conclusions: T1 increase with age may reflect increase in inflammatory processes while the decrease in MTsat may indicate that magnetization transfer may also be associated with muscle fiber macromolecules. T1 and MT_sat maps of calf muscle have the potential to detect regional and age-related compositional differences in calf muscle.

Test-Retest Reproducibility of In Vivo Magnetization Transfer Ratio and Saturation Index in Mice at 9.4 Tesla

BACKGROUND

Magnetization transfer saturation (MTsat) imaging was developed to reduce T1 dependence and improve specificity to myelin, compared to the widely used MT ratio (MTR) approach, while maintaining a feasible scan time. As MTsat imaging is an emerging technique, the reproducibility of MTsat compared to MTR must be evaluated.

PURPOSE

To assess the test-retest reproducibility of MTR and MTsat in the mouse brain at 9.4 T and calculate sample sizes potentially required to detect effect sizes ranging from 6% to 14%.

STUDY TYPE

Prospective.

SUBJECTS

Twelve healthy C57Bl/6 mice.

FIELD STRENGTH/SEQUENCE

9.4 T; magnetization transfer imaging using FLASH-3D Gradient Echo; T2-weighted TurboRARE spin echo.

ASSESSMENT

All mice were scanned at two timepoints (5 days apart). MTR and MTsat maps were analyzed using mean region-of-interest (ROIs: corpus callosum [CC], internal capsule [IC], hippocampus [HC], cortex [CX], and thalamus [TH]), and whole brain voxel-wise analysis.

STATISTICAL TESTS

Bland-Altman plots were used to assess biases between test-retest measurements. Test-retest reproducibility was evaluated via between and within-subject coefficients of variation (bsCV and wsCV, respectively). Sample sizes required were calculated (significance level: 95%; power: 80%), given effect sizes ranging from 6% to 14%, using both between and within-subject approaches. Results were considered statistically significant at P ≤ 0.05.

RESULTS

Bland-Altman plots showed negligible biases between test-retest sessions (MTR: 0.0009; MTsat: 0). ROI-based and voxel-wise CVs revealed high reproducibility for both MTR (ROI-bsCV/wsCV: CC-4.5/2.8%; IC-6.1/5.2%; HC-5.7/4.6%; CX-5.1/2.3%; TH-7.4/4.9%) and MTsat (ROI-bsCV/wsCV: CC-6.3/4.8%; IC-7.3/5.1%; HC-9.5/6.4%; CX-6.7/6.5%; TH-7.2/5.3%). With a sample size of 6, changes on the order of 15% could be detected in MTR and MTsat, both between and within subjects, while smaller changes (6%-8%) required sample sizes of 10-15 for MTR, and 15-20 for MTsat.

DATA CONCLUSION

MTsat exhibited comparable reproducibility to MTR, while providing sensitivity to myelin with less T1 dependence than MTR. EVIDENCE LEVEL: 2 TECHNICAL EFFICACY: Stage 1.

Multi-Contrast Magnetic Resonance Imaging of Visual White Matter Pathways in Patients With Glaucoma

Purpose

Glaucoma is a disorder that involves visual field loss caused by retinal ganglion cell damage. Previous diffusion magnetic resonance imaging (dMRI) studies have demonstrated that retinal ganglion cell damage affects tissues in the optic tract (OT) and optic radiation (OR). However, because previous studies have used a simple diffusion tensor model to analyze dMRI data, the microstructural interpretation of white matter tissue changes remains uncertain. In this study, we used a multi-contrast MRI approach to further clarify the type of microstructural damage that occurs in patients with glaucoma.

Methods

We collected dMRI data from 17 patients with glaucoma and 30 controls using 3-tesla (3T) MRI. Using the dMRI data, we estimated three types of tissue property metrics: intracellular volume fraction (ICVF), orientation dispersion index (ODI), and isotropic volume fraction (IsoV). Quantitative T1 (qT1) data, which may be relatively specific to myelin, were collected from all subjects.

Results

In the OT, all four metrics showed significant differences between the glaucoma and control groups. In the OR, only the ICVF showed significant between-group differences. ICVF was significantly correlated with qT1 in the OR of the glaucoma group, although qT1 did not show any abnormality at the group level.

Conclusions

Our results suggest that, at the group level, tissue changes in OR caused by glaucoma might be explained by axonal damage, which is reflected in the intracellular diffusion signals, rather than myelin damage. The significant correlation between ICVF and qT1 suggests that myelin damage might also occur in a smaller number of severe cases.

High Periventricular T1 Relaxation Times Predict Gait Improvement After Spinal Tap in Patients with Idiopathic Normal Pressure Hydrocephalus

PURPOSE

The diagnosis of idiopathic normal pressure hydrocephalus (iNPH) can be challenging. Aim of this study was to use a novel T1 mapping method to enrich the diagnostic work-up of patients with suspected iNPH.

METHODS

Using 3T magnetic resonance imaging (MRI) we prospectively evaluated rapid high-resolution T1 mapping at 0.5 mm resolution and 4 s acquisition time in 15 patients with suspected iNPH and 8 age-matched, healthy controls. T1 mapping in axial sections of the cerebrum, clinical and neuropsychological testing were performed prior to and after cerebrospinal fluid tap test (CSF-TT). T1 relaxation times were measured in 5 predefined periventricular regions.

RESULTS

All 15 patients with suspected iNPH showed gait impairment, 13 (86.6%) showed signs of cognitive impairment and 8 (53.3%) patients had urinary incontinence. Gait improvement was noted in 12 patients (80%) after CSF-TT. T1 relaxation times in all periventricular regions were elevated in patients with iNPH compared to controls with the most pronounced differences in the anterior (1006 ± 93 ms vs. 911 ± 77 ms; p = 0.023) and posterior horns (983 ± 103 ms vs. 893 ± 68 ms; p = 0.037) of the lateral ventricles. Montreal cognitive assessment (MoCA) scores at baseline were negatively correlated with T1 relaxation times (r < -0.5, p < 0.02). Higher T1 relaxation times were significantly correlated with an improvement of the 3‑m timed up and go test (r > 0.6 and p < 0.03) after CSF-TT.

CONCLUSION

In iNPH-patients, periventricular T1 relaxation times are increased compared to age-matched controls and predict gait improvement after CSF-TT. T1 mapping might enrich iNPH work-up and might be useful to indicate permanent shunting.

Mapping the human lateral geniculate nucleus and its cytoarchitectonic subdivisions using quantitative MRI

The human lateral geniculate nucleus (LGN) of the visual thalamus is a key subcortical processing site for visual information analysis. Due to its small size and deep location within the brain, a non-invasive characterization of the LGN and its microstructurally distinct magnocellular (M) and parvocellular (P) subdivisions in humans is challenging. Here, we investigated whether structural quantitative MRI (qMRI) methods that are sensitive to underlying microstructural tissue features enable MR-based mapping of human LGN M and P subdivisions. We employed high-resolution 7 Tesla in-vivo qMRI in N = 27 participants and ultra-high resolution 7 Tesla qMRI of a post-mortem human LGN specimen. We found that a quantitative assessment of the LGN and its subdivisions is possible based on microstructure-informed qMRI contrast alone. In both the in-vivo and post-mortem qMRI data, we identified two components of shorter and longer longitudinal relaxation time (T₁) within the LGN that coincided with the known anatomical locations of a dorsal P and a ventral M subdivision, respectively. Through ground-truth histological validation, we further showed that the microstructural MRI contrast within the LGN pertains to cyto- and myeloarchitectonic tissue differences between its subdivisions. These differences were based on cell and myelin density, but not on iron content. Our qMRI-based mapping strategy paves the way for an in-depth understanding of LGN function and microstructure in humans. It further enables investigations into the selective contributions of LGN subdivisions to human behavior in health and disease.

Measuring compartmental T2-orientational dependence in human brain white matter using a tiltable RF coil and diffusion-T2 correlation MRI

The anisotropy of brain white matter microstructure manifests itself in orientational-dependence of various MRI contrasts, and can result in significant quantification biases if ignored. Understanding the origins of this orientation-dependence could enhance the interpretation of MRI signal changes in development, ageing and disease and ultimately improve clinical diagnosis. Using a novel experimental setup, this work studies the contributions of the intra- and extra-axonal water to the orientation-dependence of one of the most clinically-studied parameters, apparent transverse relaxation T2. Specifically, a tiltable receive coil is interfaced with an ultra-strong gradient MRI scanner to acquire multidimensional MRI data with an unprecedented range of acquisition parameters. Using this setup, compartmental T2 can be disentangled based on differences in diffusional-anisotropy, and its orientation-dependence further elucidated by re-orienting the head with respect to the main magnetic field B→0. A dependence of (compartmental) T2 on the fibre orientation w.r.t. B→0 was observed, and further quantified using characteristic representations for susceptibility- and magic angle effects. Across white matter, anisotropy effects were dominated by the extra-axonal water signal, while the intra-axonal water signal decay varied less with fibre-orientation. Moreover, the results suggest that the stronger extra-axonal T2 orientation-dependence is dominated by magnetic susceptibility effects (presumably from the myelin sheath) while the weaker intra-axonal T2 orientation-dependence may be driven by a combination of microstructural effects. Even though the current design of the tiltable coil only offers a modest range of angles, the results demonstrate an overall effect of tilt and serve as a proof-of-concept motivating further hardware development to facilitate experiments that explore orientational anisotropy. These observations have the potential to lead to white matter microstructural models with increased compartmental sensitivity to disease, and can have direct consequences for longitudinal and group-wise T2- and diffusion-MRI data analysis, where the effect of head-orientation in the scanner is commonly ignored.

T1 Mapping for Microstructural Assessment of the Cervical Spinal Cord in the Evaluation of Patients with Degenerative Cervical Myelopathy

BACKGROUND AND PURPOSE

Although current radiologic evaluation of degenerative cervical myelopathy by conventional MR imaging accurately demonstrates spondylosis or degenerative disc disease causing spinal cord dysfunction, conventional MR imaging still fails to provide satisfactory anatomic and clinical correlations. In this context, we assessed the potential value of quantitative cervical spinal cord T1 mapping regarding the evaluation of patients with degenerative cervical myelopathy.

MATERIALS AND METHODS

Twenty patients diagnosed with mild and moderate-to-severe degenerative cervical myelopathy and 10 healthy subjects were enrolled in a multiparametric MR imaging protocol. Cervical spinal cord T1 mapping was performed with the MP2RAGE sequence procedure. Retrieved data were processed and analyzed regarding the global spinal cord and white and anterior gray matter on the basis of the clinical severity and the spinal canal stenosis grading.

RESULTS

Noncompressed levels in healthy controls demonstrated significantly lower T1 values than noncompressed, mild, moderate, and severe stenotic levels in patients. Concerning the entire spinal cord T1 mapping, patients with moderate-to-severe degenerative cervical myelopathy had higher T1 values compared with healthy controls. Regarding the specific levels, patients with moderate-to-severe degenerative cervical myelopathy demonstrated a T1 value increase at C1, C7, and the level of maximal compression compared with healthy controls. Patients with mild degenerative cervical myelopathy had lower T1 values than those with moderate-to-severe degenerative cervical myelopathy at the level of maximal compression. Analyses of white and anterior gray matter confirmed similar results. Strong negative correlations between individual modified Japanese Orthopaedic Association scores and T1 values were also observed.

CONCLUSIONS

In this preliminary study, 3D-MP2RAGE T1 mapping demonstrated increased T1 values in the pathology tissue samples, with diffuse medullary alterations in all patients with degenerative cervical myelopathy, especially relevant at C1 (nonstenotic level) and at the maximal compression level. Encouraging correlations observed with the modified Japanese Orthopaedic Association score make this novel approach a potential quantitative biomarker related to clinical severity in degenerative cervical myelopathy. Nevertheless, patients with mild degenerative cervical myelopathy demonstrated nonsignificant results compared with healthy controls and should now be studied in multicenter studies with larger patient populations.

The R1-weighted connectome: complementing brain networks with a myelin-sensitive measure

Myelin plays a crucial role in how well information travels between brain regions. Complementing the structural connectome, obtained with diffusion MRI tractography, with a myelin-sensitive measure could result in a more complete model of structural brain connectivity and give better insight into white-matter myeloarchitecture. In this work we weight the connectome by the longitudinal relaxation rate (R1), a measure sensitive to myelin, and then we assess its added value by comparing it with connectomes weighted by the number of streamlines (NOS). Our analysis reveals differences between the two connectomes both in the distribution of their weights and the modular organization. Additionally, the rank-based analysis shows that R1 can be used to separate transmodal regions (responsible for higher-order functions) from unimodal regions (responsible for low-order functions). Overall, the R1-weighted connectome provides a different perspective on structural connectivity taking into account white matter myeloarchitecture.

In the present work, we show that by using a myelin-sensitive measure we can complement the diffusion MRI-based connectivity and provide a different picture of the brain organization. We show that the R1-weighted average distribution does not follow the same trend as the number of streamlines strength distribution, and the two connectomes exhibit different modular organization. We also show that unimodal cortical regions tend to be connected by more streamlines, but the connections exhibit a lower R1-weighted average, while the transmodal regions have higher R1-weighted average but fewer streamlines. This could imply that the unimodal regions require more connections with lower myelination, whereas the transmodal regions rely on connections with higher myelination.

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.

Age dependency and lateralization in the three branches of the human superior longitudinal fasciculus

The superior longitudinal fascicle/fasciculus (SLF) is a major white matter tract connecting the frontal and parietal cortices in humans. Although the SLF has often been analyzed as a single entity, several studies have reported that the SLF is segregated into three distinct branches (SLF I, II, and III). They have also reported the right lateralization of the SLF III volume and discussed its relationship with lateralized cortical functions in the fronto-parietal network. However, to date, the homogeneity or heterogeneity of the age dependency and lateralization properties of SLF branches have not been fully clarified. Through this study, we aimed to clarify the age dependency and lateralization of SLF I-III by analyzing diffusion-weighted MRI (dMRI) and quantitative R1 (qR1) map datasets collected from a wide range of age groups, mostly comprising right-handed children, adolescents, adults, and seniors (6 to 81 years old). The age dependency in dMRI measurement (fractional anisotropy, FA) was heterogeneous among the three SLF branches, suggesting that these branches are regulated by distinct developmental and aging processes. Lateralization analysis on SLF branches revealed that the right SLF III was larger than the left SLF III in adults, replicating previous reports. FA measurement also suggested that, in addition to SLF III, SLF II was lateralized to the right hemisphere in adolescents and adults. We further found a left lateralization of SLF I in qR1 data, a microstructural measurement sensitive to myelin levels, in adults. These findings suggest that the SLF sub-bundles are distinct entities in terms of age dependency and lateralization.

The traveling heads 2.0: Multicenter reproducibility of quantitative imaging methods at 7 Tesla

OBJECT

This study evaluates inter-site and intra-site reproducibility at ten different 7 T sites for quantitative brain imaging.

MATERIAL AND METHODS

Two subjects - termed the "traveling heads" - were imaged at ten different 7 T sites with a harmonized quantitative brain MR imaging protocol. In conjunction with the system calibration, MP2RAGE, QSM, CEST and multi-parametric mapping/relaxometry were examined.

RESULTS

Quantitative measurements with MP2RAGE showed very high reproducibility across sites and subjects, and errors were in concordance with previous results and other field strengths. QSM had high inter-site reproducibility for relevant subcortical volumes. CEST imaging revealed systematic differences between the sites, but reproducibility was comparable to results in the literature. Relaxometry had also very high agreement between sites, but due to the high sensitivity, differences caused by different applications of the B1 calibration of the two RF coil types used were observed.

CONCLUSION

Our results show that quantitative brain imaging can be performed with high reproducibility at 7 T and with similar reliability as found at 3 T for multicenter studies of the supratentorial brain.

Microcephaly with altered cortical layering in GIT1 deficiency revealed by quantitative neuroimaging

G Protein-Coupled Receptor Kinase-Interacting Protein-1 (GIT1) regulates neuronal functions, including cell and axon migration and synapse formation and maintenance, and GIT1 knockout (KO) mice exhibit learning and memory deficits. We noted that male and female GIT1-KO mice exhibit neuroimaging phenotypes including microcephaly, and altered cortical layering, with a decrease in neuron density in cortical layer V. Micro-CT and magnetic resonance microscopy (MRM) were used to identify morphometric phenotypes for the skulls and throughout the GIT1-KO brains. High field MRM of actively-stained mouse brains from GIT1-KO and wild type (WT) controls (n = 6 per group) allowed segmenting 37 regions, based on co-registration to the Waxholm Space atlas. Overall brain size in GIT1-KO mice was ~32% smaller compared to WT controls. After correcting for brain size, several regions were significantly different in GIT1-KO mice relative to WT, including the gray matter of the ventral thalamic nuclei and the rest of the thalamus, the inferior colliculus, and pontine nuclei. GIT1-KO mice had reduced volume of white matter tracts, most notably in the anterior commissure (~26% smaller), but also in the cerebral peduncle, fornix, and spinal trigeminal tract. On the other hand, the basal ganglia appeared enlarged in GIT1-KO mice, including the globus pallidus, caudate putamen, and particularly the accumbens - supporting a possible vulnerability to addiction. Volume based morphometry based on high-resolution MRM (21.5 μm isotropic voxels) was effective in detecting overall, and local differences in brain volumes in GIT1-KO mice, including in white matter tracts. The reduced relative volume of specific brain regions suggests a critical, but not uniform, role for GIT1 in brain development, conducive to brain microcephaly, and aberrant connectivity.

Understanding aqueous and non-aqueous proton T1 relaxation in brain

A full picture of longitudinal relaxation in complex heterogeneous environments like white matter brain tissue remains elusive. In tissue, successive approximations, from the solvation layer model to the two pool model, have highlighted how longitudinal magnetization evolution depends on both inter-compartmental exchange and spin-lattice relaxation. In white matter, however, these models fail to capture the behaviour of the two distinct aqueous pools, myelin water and intra/extra-cellular water. A challenge with testing more comprehensive multi-pool models lies in directly observing all pools, both aqueous and non-aqueous. In this work, we advance these efforts by integrating three main experimental and analytical elements: direct observation of the longitudinal relaxation of both the aqueous and the non-aqueous protons in white matter, a wide range of different initial conditions, and application of an analysis pipeline which includes lineshape, CPMG, and fitting of a four pool model. An eigenvector interpretation of the four pool model highlights how longitudinal relaxation in white matter depends on initial conditions. We find that a single set of model parameters is able to describe the entire range of relaxation behaviour observed in all the separable aqueous and non-aqueous pools in experiments involving six different initial conditions. Understanding of the nature and connectedness of the tissue components is crucial in the design and interpretation of many MRI measurements, especially those based on magnetization transfer and longitudinal relaxation. In particular, the dependency of relaxation behaviour on initial conditions is likely the basis for understanding method-dependent discrepancies in in vivo T₁.

Hypomyelinating leukodystrophies - unravelling myelin biology

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

The Utility of a Convolutional Neural Network for Generating a Myelin Volume Index Map from Rapid Simultaneous Relaxometry Imaging

Purpose:

A current algorithm to obtain a synthetic myelin volume fraction map (SyMVF) from rapid simultaneous relaxometry imaging (RSRI) has a potential problem, that it does not incorporate information from surrounding pixels. The purpose of this study was to develop a method that utilizes a convolutional neural network (CNN) to overcome this problem.

Methods:

RSRI and magnetization transfer images from 20 healthy volunteers were included. A CNN was trained to reconstruct RSRI-related metric maps into a myelin volume-related index (generated myelin volume index: GenMVI) map using the MVI map calculated from magnetization transfer images (MTMVI) as reference. The SyMVF and GenMVI maps were statistically compared by testing how well they correlated with the MTMVI map. The correlations were evaluated based on: (i) averaged values obtained from 164 atlas-based ROIs, and (ii) pixel-based comparison for ROIs defined in four different tissue types (cortical and subcortical gray matter, white matter, and whole brain).

Results:

For atlas-based ROIs, the overall correlation with the MTMVI map was higher for the GenMVI map than for the SyMVF map. In the pixel-based comparison, correlation with the MTMVI map was stronger for the GenMVI map than for the SyMVF map, and the difference in the distribution for the volunteers was significant (Wilcoxon sign-rank test, P < 0.001) in all tissue types.

Conclusion:

The proposed method is useful, as it can incorporate more specific information about local tissue properties than the existing method. However, clinical validation is necessary.

T1 Mapping Quantifies Spinal Cord Compression in Patients With Various Degrees of Cervical Spinal Canal Stenosis

Age-related degeneration of the cervical spinal column is the most common cause of spinal cord lesions. T1 mapping has been shown to indicate the grade and site of spinal cord compression in low grade spinal canal stenosis (SCS). Aim of our study was to further investigate the diagnostic potential of a novel T1 mapping method at 0.75 mm resolution and 4 s acquisition time in 31 patients with various grades of degenerative cervical SCS. T1 mapping was performed in axial sections of the stenosis as well as above and below. Included subjects received standard T2-weighted MRI of the cervical spine (including SCS-grading 0-III), electrophysiological, and clinical examination. We found that patients with cervical SCS showed a significant difference in T1 relaxation times within the stenosis (727 ± 66 ms, mean ± standard deviation) in comparison to non-stenotic segments above (854 ± 104 ms, p < 0.001) and below (893 ± 137 ms, p < 0.001). There was no difference in mean T1 in non-stenotic segments in patients (p = 0.232) or between segments in controls (p = 0.272). Mean difference of the T1 relaxation times was significantly higher in grade III stenosis (234 ± 45) vs. in grade II stenosis (176 ± 45, p = 0.037) vs. in grade I stenosis (90 ± 87 ms, p = 0.010). A higher difference in T1 relaxation time was associated with a central efferent conduction deficit. In conclusion, T1 mapping may be useful as a tool for SCS quantification in all grades of SCS, including high-grade stenosis with myelopathy signal in conventional T2-weighted imaging.

Methodological Considerations for Neuroimaging in Deep Brain Stimulation of the Subthalamic Nucleus in Parkinson's Disease Patients

Deep brain stimulation (DBS) of the subthalamic nucleus is a neurosurgical intervention for Parkinson's disease patients who no longer appropriately respond to drug treatments. A small fraction of patients will fail to respond to DBS, develop psychiatric and cognitive side-effects, or incur surgery-related complications such as infections and hemorrhagic events. In these cases, DBS may require recalibration, reimplantation, or removal. These negative responses to treatment can partly be attributed to suboptimal pre-operative planning procedures via direct targeting through low-field and low-resolution magnetic resonance imaging (MRI). One solution for increasing the success and efficacy of DBS is to optimize preoperative planning procedures via sophisticated neuroimaging techniques such as high-resolution MRI and higher field strengths to improve visualization of DBS targets and vasculature. We discuss targeting approaches, MRI acquisition, parameters, and post-acquisition analyses. Additionally, we highlight a number of approaches including the use of ultra-high field (UHF) MRI to overcome limitations of standard settings. There is a trade-off between spatial resolution, motion artifacts, and acquisition time, which could potentially be dissolved through the use of UHF-MRI. Image registration, correction, and post-processing techniques may require combined expertise of traditional radiologists, clinicians, and fundamental researchers. The optimization of pre-operative planning with MRI can therefore be best achieved through direct collaboration between researchers and clinicians.

Cortical aging - new insights with multiparametric quantitative MRI

Understanding the microstructural changes related to physiological aging of the cerebral cortex is pivotal to differentiate healthy aging from neurodegenerative processes. The aim of this study was to investigate the age-related global changes of cortical microstructure and regional patterns using multiparametric quantitative MRI (qMRI) in healthy subjects with a wide age range. 40 healthy participants (age range: 2^nd to 8^th decade) underwent high-resolution qMRI including T1, PD as well as T2, T2* and T2′ mapping at 3 Tesla. Cortical reconstruction was performed with the FreeSurfer toolbox, followed by tests for correlations between qMRI parameters and age. Cortical T1 values were negatively correlated with age (p=0.007) and there was a widespread age-related decrease of cortical T1 involving the frontal and the parietotemporal cortex, while T2 was correlated positively with age, both in frontoparietal areas and globally (p=0.004). Cortical T2′ values showed the most widespread associations across the cortex and strongest correlation with age (r= -0.724, p=0.0001). PD and T2* did not correlate with age. Multiparametric qMRI allows to characterize cortical aging, unveiling parameter-specific patterns. Quantitative T2′ mapping seems to be a promising imaging biomarker of cortical age-related changes, suggesting that global cortical iron deposition is a prominent process in healthy aging.

Toward quantitative neuroimaging biomarkers for Friedreich's ataxia at 7 Tesla: Susceptibility mapping, diffusion imaging, R2 and R1 relaxometry

Friedreich's ataxia (FRDA) is a rare genetic disorder leading to degenerative processes. So far, no effective treatment has been found. Therefore, it is important to assist the development of medication with imaging biomarkers reflecting disease status and progress. Ten FRDA patients (mean age 37 ± 14 years; four female) and 10 age- and sex-matched controls were included. Acquisition of magnetic resonance imaging (MRI) data for quantitative susceptibility mapping, R₁ , R₂ relaxometry and diffusion imaging was performed at 7 Tesla. Results of volume of interest (VOI)-based analyses of the quantitative data were compared with a voxel-based morphometry (VBM) evaluation. Differences between patients and controls were assessed using the analysis of covariance (ANCOVA; p < 0.01) with age and sex as covariates, effect size of group differences, and correlations with disease characteristics with Spearman correlation coefficient. For the VBM analysis, a statistical threshold of 0.001 for uncorrected and 0.05 for corrected p-values was used. Statistically significant differences between FRDA patients and controls were found in five out of twelve investigated structures, and statistically significant correlations with disease characteristics were revealed. Moreover, VBM revealed significant white matter atrophy within regions of the brainstem, and the cerebellum. These regions overlapped partially with brain regions for which significant differences between healthy controls and patients were found in the VOI-based quantitative MRI evaluation. It was shown that two independent analyses provided overlapping results. Moreover, positive results on correlations with disease characteristics were found, indicating that these quantitative MRI parameters could provide more detailed information and assist the search for effective treatments.

B0-field dependence of MRI T1 relaxation in human brain

Tissue longitudinal relaxation characterized by recovery time T1 or rate R1 is a fundamental MRI contrast mechanism that is increasingly being used to study the brain's myelination patterns in both health and disease. Nevertheless, the quantitative relationship between T1 and myelination, and its dependence on B0 field strength, is still not well known. It has been theorized that in much of brain tissue, T1 field-dependence is driven by that of macromolecular protons (MP) through a mechanism called magnetization transfer (MT). Despite the explanatory power of this theory and substantial support from in-vitro experiments at low fields (<3 ​T), in-vivo evidence across clinically relevant field strengths is lacking. In this study, T1-weighted MRI was acquired in a group of eight healthy volunteers at four clinically relevant field strengths (0.55, 1.5, 3 and 7 ​T) using the same pulse sequence at a single site, and jointly analyzed based on the two-pool model of MT. MP fraction and free-water pool T1 were obtained in several brain structures at 3 and 7 ​T, which allowed distinguishing between contributions from macromolecular content and iron to tissue T1. Based on this, the T1 of MP in white matter, indirectly determined by assuming a field independent T1 of free water, was shown to increase approximately linearly with B0. This study advances our understanding of the T1 contrast mechanism and its relation to brain myelin content across the wide range of currently available MRI strengths, and it has the potential to inform design of T1 mapping methods for improved reproducibility in the human brain.

Waves of Maturation and Senescence in Micro-structural MRI Markers of Human Cortical Myelination over the Lifespan

Seminal human brain histology work has demonstrated developmental waves of myelination. Here, using a micro-structural magnetic resonance imaging (MRI) marker linked to myelin, we studied fine-grained age differences to deduce waves of growth, stability, and decline of cortical myelination over the life-cycle. In 484 participants, aged 8–85 years, we fitted smooth growth curves to T1- to T2-weighted ratio in each of 360 regions from one of seven cytoarchitectonic classes. From the first derivatives of these generally inverted-U trajectories, we defined three milestones: the age at peak growth; the age at onset of a stable plateau; and the age at the onset of decline. Age at peak growth had a bimodal distribution comprising an early (pre-pubertal) wave of primary sensory and motor cortices and a later (post-pubertal) wave of association, insular and limbic cortices. Most regions reached stability in the 30-s but there was a second wave reaching stability in the 50-s. Age at onset of decline was also bimodal: in some right hemisphere regions, the curve declined from the 60-s, but in other left hemisphere regions, there was no significant decline from the stable plateau. These results are consistent with regionally heterogeneous waves of intracortical myelinogenesis and age-related demyelination.

Inter-individual Differences in Occipital Alpha Oscillations Correlate with White Matter Tissue Properties of the Optic Radiation

Neural oscillations at ∼10 Hz, called alpha oscillations, are one of the most prominent components of neural oscillations in the human brain. In recent years, characteristics (power/frequency/phase) of occipital alpha oscillations have been correlated with various perceptual phenomena. However, the relationship between inter-individual differences in alpha oscillatory characteristics and the properties of the underlying brain structures, such as white matter pathways, is unclear. A possibility is that intrinsic occipital alpha oscillations are mediated by thalamocortical interaction; we hypothesized that the most promising candidate for characterizing the intrinsic alpha oscillation is optic radiation (OR), which is the geniculo-cortical pathway carrying signals between the lateral geniculate nucleus (LGN) and primary visual cortex (V1). We used resting-state magnetoencephalography (MEG) and diffusion-weighted/quantitative magnetic resonance imaging (MRI) (dMRI/qMRI) to correlate the frequency and power of occipital alpha oscillations with the tissue properties of the OR by focusing on the different characteristics across individuals. We found that the peak alpha frequency (PAF) negatively correlated with intracellular volume fraction (ICVF), reflecting diffusion properties in intracellular (axonal) space, whereas the peak alpha power was not correlated with any tissue properties measurements. No significant correlation was found between OR and beta frequency/amplitude or between other white matter tract connecting parietal and inferotemporal cortex and alpha frequency/amplitude. These results support the hypothesis that an interaction between thalamic nuclei and early visual areas is essential for the occipital alpha oscillatory rhythm.

Diffusion time dependency along the human corpus callosum and exploration of age and sex differences as assessed by oscillating gradient spin-echo diffusion tensor imaging

Conventional diffusion imaging uses pulsed gradient spin echo (PGSE) waveforms with diffusion times of tens of milliseconds (ms) to infer differences of white matter microstructure. The combined use of these long diffusion times with short diffusion times (<10 ​ms) enabled by oscillating gradient spin echo (OGSE) waveforms can enable more sensitivity to changes of restrictive boundaries on the scale of white matter microstructure (e.g. membranes reflecting the axon diameters). Here, PGSE and OGSE images were acquired at 4.7 ​T from 20 healthy volunteers aged 20-73 years (10 males). Mean, radial, and axial diffusivity, as well as fractional anisotropy were calculated in the genu, body and splenium of the corpus callosum (CC). Monte Carlo simulations were also conducted to examine the relationship of intra- and extra-axonal radial diffusivity with diffusion time over a range of axon diameters and distributions. The results showed elevated diffusivities with OGSE relative to PGSE in the genu and splenium (but not the body) in both males and females, but the OGSE-PGSE difference was greater in the genu for males. Females showed positive correlations of OGSE-PGSE diffusivity difference with age across the CC, whereas there were no such age correlations in males. Simulations of radial diffusion demonstrated that for axon sizes in human brain both OGSE and PGSE diffusivities were dominated by extra-axonal water, but the OGSE-PGSE difference nonetheless increased with area-weighted outer-axon diameter. Therefore, the lack of OGSE-PGSE difference in the body is not entirely consistent with literature that suggests it is composed predominantly of axons with large diameter. The greater OGSE-PGSE difference in the genu of males could reflect larger axon diameters than females. The OGSE-PGSE difference correlation with age in females could reflect loss of smaller axons at older ages. The use of OGSE with short diffusion times to sample the microstructural scale of restriction implies regional differences of axon diameters along the corpus callosum with preliminary results suggesting a dependence on age and sex.

Cerebral Gray and White Matter Involvement in Anorexia Nervosa Evaluated by T1, T2, and T2* Mapping

BACKGROUND AND PURPOSE

Changes in the brain composition of anorexics could potentially be expected, opening the door to new imaging approaches where quantitative and qualitative MRI have a role. Our purpose was to investigate anorexia-related brain dehydration and myelin depletion by analyzing T1, T2, and T2* relaxation times of different brain structures in anorexics and controls.

METHODS

Thirty-eight anorexic female patients (mean age, 26.2 years; age range, 16.2-48.7 years; mean BMI, 14.5 kg/m² ; BMI range, 10.0-18.4 kg/m² ) underwent brain MRI between August 2014 and August 2018. Controls were 16 healthy females (mean age, 28.0 years; age range, 22.3-34.7 years; mean BMI, 20.9 kg/m² ; BMI range, 18.4-26.6 kg/m² ). T1, T2, and T2* relaxation times were obtained for different brain structures in anorexics and controls as part of this retrospective case-control study.

RESULTS

The T1 relaxation times of gray and white matter were significantly lower in anorexics (P = .009), whereas the T2 relaxation times of gray matter were higher (P < .001). There were no statistically significant differences in gray matter T2* relaxation times or in white matter T2 and T2* relaxation times between anorexics and controls. Occipital lobe gray matter showed the shortest T1, T2, and T2* relaxation times of all brain regions (P < .05).

CONCLUSIONS

T1 shortening in anorexics suggests both dehydration and myelin loss, whereas T2 prolongation points toward myelin loss (myelin water has lower T2), which seems to be less discernible in white matter. Shorter overall relaxation times in the most posterior regions of the brain suggest higher iron content.

Reduced diffusion in normal appearing white matter of glioma patients following radio(chemo)therapy

BACKGROUND AND PURPOSE

Standard treatment of high grade gliomas includes gross tumour resection followed by radio(chemo)therapy. Radiotherapy inevitably leads to irradiation of normal brain tissue. The goal of this prospective, longitudinal study was to use MRI to quantify normal appearing white and grey matter changes following radiation treatment as a function of dose and time after radiotherapy.

MATERIALS AND METHODS

Pre-radiotherapy (proton or photon therapy) MRI and follow-up MRIs collected in 3 monthly intervals thereafter were analysed for 22 glioma patients and included diffusion tensor imaging, quantitative T1, T2* and proton density mapping. Abnormal tissue was excluded from analysis. MR signal changes were quantified within different dose bin regions for grey and white matter and subsequently for whole brain white matter.

RESULTS

We found significant reductions in mean diffusivity, radial diffusivity, axial diffusivity and T2* in normal appearing white matter regions receiving a radiation dose as low as 10-20 Gy within the observational period of up to 18 months. The magnitude of these changes increased with the received radiation dose and progressed with time after radiotherapy. Whole brain white matter also showed a significant reduction in radial diffusivity as a function of radiation dose and time after radiotherapy. No significant changes were observed in grey matter.

CONCLUSION

Diffusion tensor imaging and T2* imaging revealed normal appearing white matter changes following radiation treatment. The changes were dose dependant and progressed over time. Further work is needed to understand the underlying tissue changes and to correlate the observed diffusion changes with late brain malfunctions.

Rapid myelin water imaging for the assessment of cervical spinal cord myelin damage

Background

Rapid myelin water imaging (MWI) using a combined gradient and spin echo (GRASE) sequence can produce myelin specific metrics for the human brain. Spinal cord MWI could be similarly useful, but technical challenges have hindered routine application. GRASE rapid MWI was recently successfully implemented for imaging of healthy cervical spinal cord and may complement other advanced imaging methods, such as diffusion tensor imaging (DTI) and quantitative T₁ (qT₁).

Objective

To demonstrate the feasibility of cervical cord GRASE rapid MWI in multiple sclerosis (MS), primary lateral sclerosis (PLS) and neuromyelitis optica spectrum disorder (NMO), with comparison to DTI and qT₁ metrics.

Methods

GRASE MWI, DTI and qT₁ data were acquired in 2 PLS, 1 relapsing-remitting MS (RRMS), 1 primary-progressive MS (PPMS) and 2 NMO subjects, as well as 6 age (±3 yrs) and sex matched healthy controls (HC). Internal cord structure guided template registrations, used for region of interest (ROI) analysis. Z score maps were calculated for the difference between disease subject and mean HC metric values.

Results

PLS subjects had low myelin water fraction (MWF) in the lateral funiculi compared to HC. RRMS subject MWF was heterogeneous within the cord. The PPMS subject showed no trends in ROI results but had a region of low MWF Z score corresponding to a focal lesion. The NMO subject with a longitudinally extensive transverse myelitis lesion had low values for whole cord mean MWF of 12.8% compared to 24.3% (standard deviation 2.2%) for HC. The NMO subject without lesions also had low MWF compared to HC. DTI and qT₁ metrics showed similar trends, corroborating the MWF results and providing complementary information.

Conclusion

GRASE is sufficiently sensitive to detect decreased myelin within MS spinal cord plaques, NMO lesions, and PLS diffuse spinal cord injury. Decreased MWF in PLS is consistent with demyelination secondary to motor neuron degeneration. GRASE MWI is a feasible method for rapid assessment of myelin content in the cervical spinal cord and provides complementary information to that of DTI and qT₁ measures.

•

Downstream myelin changes in motor tracts of primary lateral sclerosis spinal cord.

•

Low myelin water fraction in multiple sclerosis and neuromyelitis optica cord lesions.

•

Diffuse demyelination evidence in neuromyelitis optica normal-appearing white matter.

•

Myelin water imaging provides complementary information to diffusion and T₁ metrics.

Diffusivity and quantitative T1 profile of human visual white matter tracts after retinal ganglion cell damage

In patients with retinal ganglion cell diseases, recent diffusion tensor imaging (DTI) studies have revealed structural abnormalities in visual white matter tracts such as the optic tract, and optic radiation. However, the microstructural origin of these diffusivity changes is unknown as DTI metrics involve multiple biological factors and do not correlate directly with specific microstructural properties. In contrast, recent quantitative T1 (qT1) mapping methods provide tissue property measurements relatively specific to myelin volume fractions in white matter. This study aims to improve our understanding of microstructural changes in visual white matter tracts following retinal ganglion cell damage in Leber's hereditary optic neuropathy (LHON) patients by combining DTI and qT1 measurements. We collected these measurements from seven LHON patients and twenty age-matched control subjects. For all individuals, we identified the optic tract and the optic radiation using probabilistic tractography, and evaluated diffusivity and qT1 profiles along them. Both diffusivity and qT1 measurements in the optic tract differed significantly between LHON patients and controls. In the optic radiation, these changes were observed in diffusivity but were not evident in qT1 measurements. This suggests that myelin loss may not explain trans-synaptic diffusivity changes in the optic radiation as a consequence of retinal ganglion cell disease.

•

Retinal ganglion cell damage affects diffusivity and T1 along visual pathways.

•

DTI metric identified white matter change in both optic tract and optic radiation.

•

T1 measurement in optic radiation did not exhibit abnormality, unlike DTI metric.

•

Myelin loss may not be a major cause of diffusivity change along optic radiation.

MP2RAGEME: T1 , T2 * , and QSM mapping in one sequence at 7 tesla

Quantitative magnetic resonance imaging generates images of meaningful physical or chemical variables measured in physical units that allow quantitative comparisons between tissue regions and among subjects scanned at the same or different sites. Here, we show that we can acquire quantitative T₁, T₂^*, and quantitative susceptibility mapping (QSM) information in a single acquisition, using a multi‐echo (ME) extension of the second gradient‐echo image of the MP2RAGE sequence. This combination is called MP2RAGE ME, or MP2RAGEME. The simultaneous acquisition results in large time savings, perfectly coregistered data, and minimal image quality differences compared to separately acquired data. Following a correction for residual transmit B₁⁺‐sensitivity, quantitative T₁, T₂^*, and QSM values were in excellent agreement with those obtained from separately acquired, also B₁⁺‐corrected, MP2RAGE data and ME gradient echo data. The quantitative values from reference regions of interests were also in very good correspondence with literature values. From the MP2RAGEME data, we further derived a multiparametric cortical parcellation, as well as a combined arterial and venous map. In sum, our MP2RAGEME sequence has the benefit in large time savings, perfectly coregistered data and minor image quality differences.

Quantification of spinal cord compression using T1 mapping in patients with cervical spinal canal stenosis - Preliminary experience

Background

Degenerative changes of the cervical spinal column are the most common cause of spinal cord lesions in the elderly. Conventional clinical, electrophysiological and radiological diagnostics of spinal cord compression are often inconsistent.

Materials and methods

The feasibility and diagnostic potential of a novel T1 mapping method at 0.5 mm resolution and 4 s acquisition time was evaluated in 14 patients with degenerative cervical spinal canal stenosis (SCS) and 6 healthy controls. T1 mapping was performed in axial sections of the stenosis as well as above and below. All subjects received standard T2-weighted MRI of the cervical spine (including SCS-grading 0-III), electrophysiological and clinical examinations.

Results

Patients revealed significantly decreased T1 relaxation times of the compressed spinal cord within the SCS (912 ± 53 ms, mean ± standard deviation) in comparison to unaffected segments above (1027 ± 39 ms, p < .001) and below (1056 ± 93 ms, p < .001). There was no difference in mean T1 in unaffected segments in patients (p = .712) or between segments in controls (p = .443). Moreover, T1 values were significantly lower in grade II (881 ± 46 ms, p = .005) than in grade I SCS (954 ± 29 ms). Patients with central conduction deficit tended to have lower T1 values within the SCS than patients without (909 ± 50 ms vs 968 ± 7 ms, p = .069).

Conclusion

Rapid high-resolution T1 mapping is a robust MRI method for quantifying spinal cord compression in patients with cervical SCS. It promises additional diagnostic insights and warrants more extended patient studies.

•

Rapid T1 mapping at 0.5 mm resolution was tested in cervical spinal canal stenosis (SCS).

•

T1 relaxation times significantly decreased within the SCS.

•

T1 relaxation times were significantly lower in grade II vs grade I SCS.

•

Central conduction deficits were inversely correlated with T1 relaxation time.

•

Rapid T1 mapping robustly and accurately quantifies spinal cord compression.

High-resolution imaging of distinct human corpus callosum microstructure and topography of structural connectivity to cortices at high field

Characterization of the microstructural properties and topography of the human corpus callosum (CC) is key to understanding interhemispheric neural communication and brain function. In this work, we tested the hypothesis that high-resolution T₁ relaxometry at high field has adequate sensitivity and specificity for characterizing microstructural properties of the human CC, and elucidating the structural connectivity of the callosal fibers to the cortices of origin. The high-resolution parametric T₁ images acquired from healthy subjects (N = 16) at 7 T clearly showed a consistent T₁ distribution among individuals with substantial variation along the human CC axis, which is highly similar to the spatial patterns of myelin density and myelinated axon size based on the histology study. Compared to the anterior part of the CC, the posterior midbody and splenium had significantly higher T₁ values. In conjunction with T₁-based classification method, the splenial T₁ values were decoded more reliably compared to a conventional partitioning method, showing a much higher T₁ value in the inferior splenium than in the middle/superior splenium. Moreover, the T₁ profile of the callosal subdivision represented the topology of the fiber connectivity to the projected cortical regions: the fibers in the posterior midbody and inferior splenium with a higher T₁ (inferring a larger axon size) were mainly connected to motor-sensory and visual cortical areas, respectively; in contrast, the fibers in the anterior/posterior CC with a lower T₁ (inferring a smaller axon size) were primarily connected to the frontal/parietal-temporal areas. These findings indicate that high-resolution T₁ relaxometry imaging could provide a complementary and robust neuroimaging tool, useful for exploring the complex tissue properties and topographic organization of the human corpus callosum.

White matter intercompartmental water exchange rates determined from detailed modeling of the myelin sheath

PURPOSE

Magnetization exchange (ME) between hydrogen protons of water and large molecules (semisolids [SS]) in lipid bilayers is an important factor in MRI signal generation and can be exploited to study white matter pathology. Current models used to quantify ME in white matter generally consider water to reside in 1 or 2 distinct compartments, ignoring the complexities of the myelin sheath's multicompartment structure of alternating myelin SS and myelin water (MW) layers. Here, we investigated the effect of this by fitting ME data obtained from human brain at 7 T with a multilayer model of myelin.

METHODS

A multi-echo acquisition for a T₂ ^* -based separation of MW from other water signals was combined with various preparation pulses to change the (relative) state of the SS and water pools and analyzed by fitting with a multilayer exchange model.

RESULTS

The estimated lifetime within a single MW layer was 260 µs, corresponding to a lipid bilayer permeability of 6.7 µm/s. The magnetization lifetime of the aggregate of all MW was estimated at 13 ms, shorter than previously reported values in the range of 40 to 140 ms.

CONCLUSION

Contrary to expectations and previous reports, ME between protons in myelin SS and water is not limited by the myelin sheath but rather by the exchange between SS and water protons. The analysis of ME contrast should account for the relatively short MW lifetime and affects the interpretation of tissue compartmentalization from MRI contrasts such as T₁ - and diffusion-weighting.

Myelin Measurement: Comparison Between Simultaneous Tissue Relaxometry, Magnetization Transfer Saturation Index, and T1w/T2w Ratio Methods

Magnetization transfer (MT) imaging has been widely used for estimating myelin content in the brain. Recently, two other approaches, namely simultaneous tissue relaxometry of R1 and R2 relaxation rates and proton density (SyMRI) and the ratio of T1-weighted to T2-weighted images (T1w/T2w ratio), were also proposed as methods for measuring myelin. SyMRI and MT imaging have been reported to correlate well with actual myelin by histology. However, for T1w/T2w ratio, such evidence is limited. In 20 healthy adults, we examined the correlation between these three methods, using MT saturation index (MTsat) for MT imaging. After calibration, white matter (WM) to gray matter (GM) contrast was the highest for SyMRI among these three metrics. Even though SyMRI and MTsat showed strong correlation in the WM (r = 0.72), only weak correlation was found between T1w/T2w and SyMRI (r = 0.45) or MTsat (r = 0.38) (correlation coefficients significantly different from each other, with p values < 0.001). In subcortical and cortical GM, these measurements showed moderate to strong correlations to each other (r = 0.54 to 0.78). In conclusion, the high correlation between SyMRI and MTsat indicates that both methods are similarly suited to measure myelin in the WM, whereas T1w/T2w ratio may be less optimal.

Bihemispheric alterations in myelination in children following unilateral perinatal stroke

Background

Stroke is a leading cause of perinatal brain injury with variable outcomes including cerebral palsy and epilepsy. The biological processes that underlie these heterogeneous outcomes are poorly understood. Alterations in developmental myelination are recognized as a major determinant of outcome in preterm brain injury but have not been explored in perinatal stroke. We aimed to characterize myelination in hemiparetic children after arterial perinatal stroke, hypothesizing that ipsilesional myelination would be impaired, the degree of which would correlate with poor outcome.

Methods

Retrospective, controlled cohort study. Participants were identified through the Alberta Perinatal Stroke Project (APSP), a population-based research cohort (n > 400). Inclusion criteria were: 1) MRI-confirmed, unilateral arterial perinatal stroke, 2) T1-weighted MRI after 6 months of age, 3) absence of other neurological disorders, 4) neurological outcome that included at least one of the following tests - Pediatric Stroke Outcome Measure (PSOM), Assisting Hand Assessment (AHA), Melbourne Assessment (MA), neuropsychological evaluation (NPE), and robotic sensorimotor measurements. FreeSurfer software measured hemispheric asymmetry in myelination intensity (primary outcome). A second method using ImageJ software validated the detection of myelination asymmetry. A repeated measures ANOVA was used to compare perilesional, ipsilesional remote, and contralesional homologous region myelination between stroke cases and typically developing controls. Myelination metrics were compared to clinical outcome measures (t-test, Pearson's correlation).

Results

Twenty youth with arterial stroke (mean age: 13.4 ± 4.2yo) and 27 typically developing controls (mean age: 12.5 ± 3.7yo) were studied in FreeSurfer. Participants with stroke demonstrated lower myelination in the ipsilesional hemisphere (p < 0.0001). Myelination in perilesional regions had lower intensity compared to ipsilesional remote areas (p < .00001) and contralesional homologous areas (p < 0.00001). Ipsilesional remote regions had decreased myelination compared to homologous regions on the contralesional hemisphere (p = 0.016). Contralesional myelination was decreased compared to controls (p < 0.00001). Myelination metrics were not strongly associated with clinical motor, robotic sensorimotor, or neuropsychological outcomes though some complex tests requiring speeded responses had moderate effect sizes.

Conclusion

Myelination of apparently uninjured brain in both the ipsilesional and contralesional hemispheres is decreased after perinatal stroke. Differences appear to radiate outward from the lesion. Further study is needed to determine clinical significance.

•

Myelination is altered in the lesioned hemisphere after perinatal stroke.

•

The uninjured, contralesional hemisphere also demonstrates differences in myelination.

•

Simple software can estimate MRI myelination abnormalities in children with perinatal brain injury.

Scan-rescan of axcaliber, macromolecular tissue volume, and g-ratio in the spinal cord

PURPOSE

Recent MRI techniques have been introduced that can extract microstructural information in the white matter, such as the density or macromolecular content. Translating quantitative MRI to the clinic raises many challenges in terms of acquisition strategy, modeling of the MRI signal, artifact corrections, and metric extraction (template registration and partial volume effects). In this work, we investigated the scan-rescan repeatability of several quantitative MRI techniques in the human spinal cord.

METHODS

AxCaliber metrics, macromolecular tissue volume, and the fiber g-ratio were estimated in the spinal cord of eight healthy subjects, scanned and rescanned the same day in two different sessions.

RESULTS

Scan-rescan repeatability deviation was 3% for all metrics, in average in the white matter of all subjects. Intraclass correlation coefficient was up to 0.9. A three-way analysis of variance showed significant effects of white matter pathway, laterality, and subject.

CONCLUSION

The present study suggests that quantitative MRI gives stable measurements of white matter microstructure in the spinal cord of healthy subjects. Our findings remain to be evaluated in diseased populations. Magn Reson Med 79:2759-2765, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Microstructural imaging in the spinal cord and validation strategies

In vivo histology using magnetic resonance imaging (MRI) is a newly emerging research field that aims to non-invasively characterize tissue microstructure. The implications of in vivo histology are many, from discovering novel biomarkers to studying human development, to providing tools for disease diagnosis and monitoring the effects of novel treatments on tissue. This review focuses on quantitative MRI (qMRI) techniques that are used to map spinal cord microstructure. Opening with a rationale for non-invasive imaging of the spinal cord, this article continues with a brief overview of the existing MRI techniques for axon and myelin imaging, followed by the specific challenges and potential solutions for acquiring and processing such data. The final part of this review focuses on histological validation, with suggested tissue preparation, acquisition and processing protocols for large-scale microscopy.

Studying brain microstructure with magnetic susceptibility contrast at high-field

A rapidly developing application of high field MRI is the study of brain anatomy and function with contrast based on the magnetic susceptibility of tissues. To study the subtle variations in susceptibility contrast between tissues and with changes in brain activity, dedicated scan techniques such as susceptibility-weighted MRI and blood-oxygen level dependent functional MRI have been developed. Particularly strong susceptibility contrast has been observed with systems that operate at 7T and above, and their recent widespread use has led to an improved understanding of contributing sources and mechanisms. To interpret magnetic susceptibility contrast, analysis approaches have been developed with the goal of extracting measures that report on local tissue magnetic susceptibility, a physical quantity that, under certain conditions, allows estimation of blood oxygenation, local tissue iron content, and quantification of its changes with disease. Interestingly, high field studies have also brought to light that not only the makeup of tissues affects MRI susceptibility contrast, but that also a tissue's sub-voxel structure at scales all the way down to the molecular level plays an important role as well. In this review, various ways will be discussed by which sub-voxel structure can affect the MRI signal in general, and magnetic susceptibility in particular, sometimes in a complex fashion. In the light of this complexity, it appears likely that accurate, brain-wide quantification of iron will require the combination of multiple contrasts that may include diffusion and magnetization transfer information with susceptibility-weighted contrast. On the other hand, this complexity also offers opportunities to use magnetic susceptibility contrast to inform about specific microstructural aspects of brain tissue. Details and several examples will be presented in this review.

Inferring brain tissue composition and microstructure via MR relaxometry

MRI relaxometry is sensitive to a variety of tissue characteristics in a complex manner, which makes it both attractive and challenging for characterizing tissue. This article reviews the most common water proton relaxometry measures, T₁, T₂, and T₂^*, and reports on their development and current potential to probe the composition and microstructure of brain tissue. The development of these relaxometry measures is challenged by the need for suitably accurate tissue models, as well as robust acquisition and analysis methodologies. MRI relaxometry has been established as a tool for characterizing neural tissue, particular with respect to myelination, and the potential for further development exists.

Promise and pitfalls of g-ratio estimation with MRI

The fiber g-ratio is the ratio of the inner to the outer diameter of the myelin sheath of a myelinated axon. It has a limited dynamic range in healthy white matter, as it is optimized for speed of signal conduction, cellular energetics, and spatial constraints. In vivo imaging of the g-ratio in health and disease would greatly increase our knowledge of the nervous system and our ability to diagnose, monitor, and treat disease. MRI based g-ratio imaging was first conceived in 2011, and expanded to be feasible in full brain white matter with preliminary results in 2013. This manuscript reviews the growing g-ratio imaging literature and speculates on future applications. It details the methodology for imaging the g-ratio with MRI, and describes the known pitfalls and challenges in doing so.