Inhomogeneous Magnetization Transfer (ihMT) imaging in the acute cuprizone mouse model of demyelination/remyelination.

BACKGROUND: To investigate the association of ihMT (inhom signals with the demyelination and remyelination phases of the acute cuprizone mouse model in comparison with histology, and to assess the extent of tissue damage and repair from MRI data. METHODS: Acute demyelination by feeding 0.2% cuprizone for five weeks, followed by a four-week remyelination period was applied on genetically modified plp-GFP mice. Animals were scanned at different time points of the demyelination and remyelination phases of the cuprizone model using a multimodal MRI protocol, including ihMT T1D-filters, MPF (Macromolecular Proton Fraction) and R1 (longitudinal relaxation rate). For histology, plp-GFP (proteolipid protein - Green Fluorescent Protein) microscopy and LFB (Luxol Fast Blue) staining were employed as references for the myelin content. Comparison of MRI with histology was performed in the medial corpus callosum (mCC) and cerebral cortex (CTX) at two brain levels whereas ROI-wise and voxel-based analyses of the MRI metrics allowed investigating in vivo the spatial extent of myelin alterations. RESULTS: IhMT high-pass T1D-filters, targeted toward long T1D components, showed significant temporal variations in the mCC consistent with the effects induced by the cuprizone toxin. In addition, the corresponding signals correlated strongly and significantly with the myelin content assessed by GFP fluorescence and LFB staining over the demyelination and the remyelination phases. The signal of the band-pass T1D-filter, which isolates short T1D components, showed changes over time that were poorly correlated with histology, hence suggesting a sensitivity to pathological processes possibly not related to myelin. Although MPF was also highly correlated to histology, ihMT high-pass T1D-filters showed better capability to characterize the spatial-temporal patterns during the demyelination and remyelination phases of the acute cuprizone model (e.g., rostro-caudal gradient of demyelination in the mCC previously described in the literature). CONCLUSIONS: IhMT sequences selective for long T1D components are specific and sensitive in vivo markers of demyelination and remyelination and have successfully captured the spatially heterogeneous pattern of the demyelination and remyelination mechanisms in the cuprizone model. Interestingly, differences in signal variations between the ihMT high-pass and band-pass T1D-filter, suggest a sensitivity of the ihMT sequences targeted to short T1Ds to alterations other than those of myelin. Future studies will need to further address these differences by examining more closely the origin of the short T1D components and the variation of each T1D component in pathology.

Quantitative magnetization transfer MRI unbiased by on-resonance saturation and dipolar order contributions.

PURPOSE: To demonstrate the bias in quantitative MT (qMT) measures introduced by the presence of dipolar order and on-resonance saturation (ONRS) effects using magnetization transfer (MT) spoiled gradient-recalled (SPGR) acquisitions, and propose changes to the acquisition and analysis strategies to remove these biases. METHODS: The proposed framework consists of SPGR sequences prepared with simultaneous dual-offset frequency-saturation pulses to cancel out dipolar order and associated relaxation (T1D ) effects in Z-spectrum acquisitions, and a matched quantitative MT (qMT) mathematical model that includes ONRS effects of readout pulses. Variable flip angle and MT data were fitted jointly to simultaneously estimate qMT parameters (macromolecular proton fraction [MPF], T2,f , T2,b , R, and free pool T1 ). This framework is compared with standard qMT and investigated in terms of reproducibility, and then further developed to follow a joint single-point qMT methodology for combined estimation of MPF and T1 . RESULTS: Bland-Altman analyses demonstrated a systematic underestimation of MPF (-2.5% and -1.3%, on average, in white and gray matter, respectively) and overestimation of T1 (47.1 ms and 38.6 ms, on average, in white and gray matter, respectively) if both ONRS and dipolar order effects are ignored. Reproducibility of the proposed framework is excellent (ΔMPF = -0.03% and ΔT1 = -19.0 ms). The single-point methodology yielded consistent MPF and T1 values with respective maximum relative average bias of -0.15% and -3.5 ms found in white matter. CONCLUSION: The influence of acquisition strategy and matched mathematical model with regard to ONRS and dipolar order effects in qMT-SPGR frameworks has been investigated. The proposed framework holds promise for improved accuracy with reproducibility.

A strategy to reduce the sensitivity of inhomogeneous magnetization transfer (ihMT) imaging to radiofrequency transmit field variations at 3 T.

PURPOSE: To minimize the sensitivity of inhomogeneous magnetization transfer gradient-echo (ihMT-GRE) imaging to radiofrequency (RF) transmit field ( B1+ ) inhomogeneities at 3 T. METHODS: The ihMT-GRE sequence was optimized by varying the concentration of the RF saturation energy over time, obtained by increasing the saturation pulse power while extending the sequence repetition time (TR). Different protocols were tested using numerical simulations and human in vivo experiments in the brain white matter (WM) of healthy subjects at 3 T. The sensitivity of the ihMT ratio (ihMTR) to B1+ variations was investigated by comparing measurements obtained at nominal transmitter adjustments and following a 20% global B1+ drop. The resulting relative variations (δihMTR ) were evaluated voxelwise as a function of the local B1+ distribution. The reproducibility of the protocol providing minimal B1+ bias was assessed in a test-retest experiment. RESULTS: In line with simulations, ihMT-GRE experiments conducted at high concentration of the RF energy over time demonstrated strong reduction of the B1+ inhomogeneity effects in the human WM. Under the optimal conditions of 350-ms TR and 3-µT root mean square (RMS) saturation power, 73% of all WM voxels presented δihMTR below 10%. Reproducibility analysis yielded a close-to-zero systematic bias (ΔihMTR = -0.081%) and a high correlation (ρ² = 0.977) between test and retest experiments. CONCLUSION: Concentrating RF saturation energy in ihMT-GRE sequences mitigates the sensitivity of the ihMTR to B1+ variations and allows for clinical-ready ihMT imaging at 3 T. This feature is of particular interest for high and ultra-high field applications.

T1D -weighted ihMT imaging - Part II. Investigating the long- and short-T1D components correlation with myelin content. Comparison with R1 and the macromolecular proton fraction.

PURPOSE: To investigate the long- and short-T1D components correlation with myelin content using inhomogeneous magnetization transfer (ihMT) high-pass and band-pass T1D -filters and to compare ihMT, R1 , and the macromolecular proton fraction (MPF) for myelin specific imaging. METHODS: The 3D ihMT rapid gradient echo (ihMTRAGE) sequences with increasing switching times (Δt) were used to derive ihMT high-pass T1D -filters with increasing T1D cutoff values and an ihMT band-pass T1D -filter for components in the 100 µs to 1 ms range. 3D spoiled gradient echo quantitative MT (SPGR-qMT) protocols were used to derive R1 and MPF maps. The specificity of R1 , MPF, and ihMT T1D -filters was evaluated by comparison with two histological reference techniques for myelin imaging. RESULTS: The higher contribution of long-T1D s as compared to the short components as Δt got longer led to an increase in the specificity to myelination. In contrast, focusing on the signal originating from a narrow range of short-T1D s (< 1 ms) as isolated by the band-pass T1D -filter led to lower specificity. In addition, the significantly lower r2 correlation coefficient of the band-pass T1D -filter suggests that the origin of short-T1D components is mostly associated with non-myelin protons. Also, the important contribution of short-T1D s to the estimated MPF, explains its low specificity to myelination as compared to the ihMT high-pass T1D -filters. CONCLUSION: Long-T1D components imaging by means of ihMT high-pass T1D -filters is proposed as an MRI biomarker for myelin content. Future studies should enable the investigation of the sensitivity of ihMT T1D -filters for demyelinating processes.

T1D -weighted ihMT imaging - Part I. Isolation of long- and short-T1D components by T1D -filtering.

PURPOSE: To identify T1D -filtering methods, which can specifically isolate various ranges of T1D components as they may be sensitive to different microstructural properties. METHODS: Modified Bloch-Provotorov equations describing a bi-T1D component biophysical model were used to simulate the inhomogeneous magnetization transfer (ihMT) signal from ihMTRAGE sequences at high RF power and low duty-cycle with different switching time values for the dual saturation experiment: Δt = 0.0, 0.8, 1.6, and 3.2 ms. Simulations were compared with experimental signals on the brain gray and white matter tissues of healthy mice at 7T. RESULTS: The lengthening of Δt created ihMT high-pass T1D -filters, which efficiently eliminated the signal from T1D components shorter than 1 ms, while partially attenuating that of longer components (≥ 1 ms). Subtraction of ihMTR images obtained with Δt = 0.0 ms and Δt = 0.8 ms generated a new ihMT band-pass T1D -filter isolating short-T1D components in the 100-µs to 1-ms range. Simulated ihMTR values in central nervous system tissues were confirmed experimentally. CONCLUSION: Long- and short-T1D components were successfully isolated with high RF power and low duty-cycle ihMT filters in the healthy mouse brain. Future studies should investigate the various T1D -range microstructural correlations in in vivo tissues.

Determination of optimal parameters for 3D single-point macromolecular proton fraction mapping at 7T in healthy and demyelinated mouse brain.

PURPOSE: To determine optimal constrained tissue parameters and off-resonance sequence parameters for single-point macromolecular proton fraction (SP-MPF) mapping based on a comprehensive quantitative magnetization transfer (qMT) protocol in healthy and demyelinated living mice at 7T. METHODS: Using 3D spoiled gradient echo-based sequences, a comprehensive qMT protocol is performed by sampling the Z-spectrum of mice brains, in vivo. Provided additional T1 , B1+ and B0 maps allow for the estimation of qMT tissue parameters, among which three will be constrained, namely the longitudinal and transverse relaxation characteristics of the free pool (R1,f T2,f ), the cross-relaxation rate (R) and the bound pool transverse relaxation time (T2,r ). Different sets of constrained parameters are investigated to reduce the bias between the SP-MPF and its reference based on the comprehensive protocol. RESULTS: Based on a whole-brain histogram analysis about the constrained parameters, the optimal experimental parameters that minimize the global bias between reference and SP-MPF maps consist of a 600° and 6 kHz off-resonance irradiation pulse. Following a Bland-Altman analysis over regions of interest, optimal constrained parameters were R1,f T2,f  = 0.0129, R = 26.5 s-1 , and T2,r  = 9.1 µs, yielding an overall MPF bias of 10-4 (limits of agreement [-0.0068;0.0070]) and a relative variation of 0.64% ± 5.95% between the reference and the optimal single-point method across all mice. CONCLUSION: The necessity of estimating animal model- and field-dependent constrained parameters was demonstrated. The single-point MPF method can be reliably applied at 7T, as part of routine preclinical in vivo imaging protocol in mice.

Correlations of quantitative MRI metrics with myelin basic protein (MBP) staining in a murine model of demyelination.

Myelin imaging in the central nervous system is essential for monitoring pathologies involving white matter alterations. Various quantitative MRI protocols relying on the modeling of the interactions of water protons with myelinated tissues have shown sensitivities in case of myelin disruption. Some extracted model parameters are more sensitive to demyelination, such as the bound pool fraction (f) in quantitative magnetization transfer imaging (qMTI), the radial diffusivity in diffusion tensor imaging (DTI), and the myelin water fraction (MWF) in myelin water imaging (MWI). A 3D ultrashort echo time (UTE) sequence within an appropriate water suppression condition (Diff-UTE) is also considered for the direct visualization of the myelin semi-solid matrix (Diff-UTE normalized signal; rSPF). In this paper, we aimed at assessing the sensitivities and correlations of the parameters mentioned above to an immuno-histological study of the myelin basic protein (MBP) in a murine model of demyelination at 7 T. We demonstrated a high sensitivity of the MRI metrics to demyelination, and strong Spearman correlations in the corpus callosum between histology, macromolecular proton fraction (ρ>0.87) and Diff-UTE signal (ρ>0.76), but moderate ones with radial diffusivity and MWF (|ρ|<0.70).

A diffusion-based method for long-T2 suppression in steady state sequences: Validation and application for 3D-UTE imaging.

PURPOSE: To introduce a novel method for long-T2 signal physical suppression in steady-state based on configuration states combination and modulation using diffusion weighting. Its efficiency in yielding a high contrast in short-T2 structures using an ultrashort echo time acquisition module (Diff-UTE) is compared to the adiabatically prepared Inversion-Recovery-UTE sequence (IR-UTE). THEORY AND METHODS: Using a rectangular-pulse prepared 3D-UTE sequence, the possibility of long-T2 component signal cancellation through diffusion effects is addressed, and the condition met for sets of sequence parameters. Simultaneously, the short-T2 component signal is maximized using a Bloch equation-based optimization process. The method is evaluated from simulations, and experiments are conducted on a phantom composed of short and long-T2 components, as well as on an ex vivo mouse head. RESULTS: Within equal scan times, the proposed method allowed for an efficient long-T2 signal suppression, and expectedly yielded a higher signal to noise ratio in short-T2 structures compared to the IR-UTE technique, although an intrinsic short-T2 signal loss is expected through the preparation module. CONCLUSION: The Diff-UTE method represents an interesting alternative to the IR-UTE technique. Diffusion weighting allowing for a long-T2 suppression results in a less penalizing method to generate a high and selective contrast in short-T2 components. Magn Reson Med 80:548-559, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

A comprehensive protocol for quantitative magnetic resonance imaging of the brain at 3 Tesla.

Quantitative MRI (qMRI) has been shown to be clinically useful for numerous applications in the brain and body. The development of rapid, accurate, and reproducible qMRI techniques offers access to new multiparametric data, which can provide a comprehensive view of tissue pathology. This work introduces a multiparametric qMRI protocol along with full postprocessing pipelines, optimized for brain imaging at 3 Tesla and using state-of-the-art qMRI tools. The total scan time is under 50 minutes and includes eight pulse-sequences, which produce range of quantitative maps including T1, T2, and T2* relaxation times, magnetic susceptibility, water and macromolecular tissue fractions, mean diffusivity and fractional anisotropy, magnetization transfer ratio (MTR), and inhomogeneous MTR. Practical tips and limitations of using the protocol are also provided and discussed. Application of the protocol is presented on a cohort of 28 healthy volunteers and 12 brain regions-of-interest (ROIs). Quantitative values agreed with previously reported values. Statistical analysis revealed low variability of qMRI parameters across subjects, which, compared to intra-ROI variability, was x4.1 ± 0.9 times higher on average. Significant and positive linear relationship was found between right and left hemispheres' values for all parameters and ROIs with Pearson correlation coefficients of r>0.89 (P<0.001), and mean slope of 0.95 ± 0.04. Finally, scan-rescan stability demonstrated high reproducibility of the measured parameters across ROIs and volunteers, with close-to-zero mean difference and without correlation between the mean and difference values (across map types, mean P value was 0.48 ± 0.27). The entire quantitative data and postprocessing scripts described in the manuscript are publicly available under dedicated GitHub and Figshare repositories. The quantitative maps produced by the presented protocol can promote longitudinal and multi-center studies, and improve the biological interpretability of qMRI by integrating multiple metrics that can reveal information, which is not apparent when examined using only a single contrast mechanism.

Spinal cord and brain tissue impairments as long-term effects of rugby practice? An exploratory study based on T1 and ihMTsat measures.

Rugby players are subject to multiple impacts to their head and neck that could have adverse neurological effects and put them at increased risk of neurodegeneration. Previous studies demonstrated altered default mode network and diffusion metrics on brain, as well as more foraminal stenosis, disc protrusion and neck pain among players of contact sports as compared to healthy controls. However, the long-term effects of practice and repetitive impacts on brain and cervical spinal cord (cSC) of the rugby players have never been systematically investigated. In this study, 15 retired professional and amateur rugby players (R) and 15 age-matched healthy controls (HC) (all males; mean age R: 46.8 ± 7.6; and HC: 48.6 ± 9.5) were recruited both to investigate cord impairments and further characterize brain structure damage. Medical questionnaires including modified Japanese Orthopedic Association scale (mJOA) and Neck Disability Index (NDI) were filled by all participants. A 3 T multi-parametric MR protocol including conventional qualitative techniques such as T1-, T2-, and T2*-weighted sequences, as well as state-of-the art quantitative techniques including MP2RAGE T1 mapping and 3D ihMTRAGE, was used on both brain and cSC. Normalized brain WM and GM volumes, spine Overall Stenosis Score, cord cross-sectional area and regional T1 and ihMT metrics were derived from these acquisitions. Rugby players showed significantly higher NDI scores, as well as a faster decline of normalized brain GM volume with age as compared to HC. Moreover, higher T1 values on cSC suggestive of structural degeneration, together with higher T1 and lower ihMTsat on brain WM suggestive of demyelination, were observed in retired rugby players as compared to age-matched controls, which may suggest cumulative effects of long-term impacts on the tissues. Metrics also suggest early aging and different aging processes on brain tissue in the players. These preliminary observations provide new insights in the domain, which should now be further investigated on larger cohorts and multicentric longitudinal studies, and further correlated to the likelihood of neurodegenerative diseases and risk factors.

Non-invasive assessment of skeletal muscle fibrosis in mice using nuclear magnetic resonance imaging and ultrasound shear wave elastography.

Fibrosis is a key pathological feature in muscle disorders, but its quantification mainly relies on histological and biochemical assays. Muscle fibrosis most frequently is entangled with other pathological processes, as cell membrane lesions, inflammation, necrosis, regeneration, or fatty infiltration, making in vivo assessment difficult. Here, we (1) describe a novel mouse model with variable levels of induced skeletal muscle fibrosis displaying minimal inflammation and no fat infiltration, and (2) report how fibrosis affects non-invasive metrics derived from nuclear magnetic resonance (NMR) and ultrasound shear-wave elastography (SWE) associated with a passive biomechanical assay. Our findings show that collagen fraction correlates with multiple non-invasive metrics. Among them, muscle stiffness as measured by SWE, T2, and extracellular volume (ECV) as measured by NMR have the strongest correlations with histology. We also report that combining metrics in a multi-modality index allowed better discrimination between fibrotic and normal skeletal muscles. This study demonstrates that skeletal muscle fibrosis leads to alterations that can be assessed in vivo with multiple imaging parameters. Furthermore, combining NMR and SWE passive biomechanical assay improves the non-invasive evaluation of skeletal muscle fibrosis and may allow disentangling it from co-occurring pathological alterations in more complex scenarios, such as muscular dystrophies.

MRI assessment of multiple dipolar relaxation time (T1D) components in biological tissues interpreted with a generalized inhomogeneous magnetization transfer (ihMT) model.

T1D, the relaxation time of dipolar order, is sensitive to slow motional processes. Thus T1D is a probe for membrane dynamics and organization that could be used to characterize myelin, the lipid-rich membrane of axonal fibers. A mono-component T1D model associated with a modified ihMT sequence was previously proposed for in vivo evaluation of T1D with MRI. However, experiments have suggested that myelinated tissues exhibit multiple T1D components probably due to a heterogeneous molecular mobility. A bi-component T1D model is proposed and implemented. ihMT images of ex-vivo, fixed rat spinal cord were acquired with multiple frequency alternation rate. Fits to data yielded two T1Ds of about 500 µs and 10 ms. The proposed model seems to further explore the complexity of myelin organization compared to the previously reported mono-component T1D model.