Analysis and correction of biases in cross-relaxation MRI due to biexponential longitudinal relaxation

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

PURPOSE

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

Accurate actual flip angle imaging (AFI) in the presence of fat

PURPOSE

To investigate the effect of incomplete fat spoiling on the accuracy of B1 mapping with actual flip angle imaging (AFI) and to propose a method to minimize the errors using the chemical shift properties of fat.

THEORY AND METHODS

Diffusion-based dephasing is the main spoiling mechanism exploited in AFI. However, a very low diffusion in fat may make the spoiling insufficient, leading to ghosts in the B1 maps. As the errors retain the chemical-shift signature of fat, their impact can be minimized using chemical-shift-based fat signal removal from AFI acquisition modified to include multi-echo readout. The source of the errors and the proposed correction were studied in simulations and phantom and in-vivo imaging experiments.

RESULTS

Our results support that AFI artifacts are caused by the incomplete fat spoiling present in clinically attractive short TR acquisition regimes. The correction eliminated the ghosting and significantly improved the B1 mapping accuracy as well as the accuracy of R1 mapping performed with AFI-derived B1 maps.

CONCLUSIONS

The incomplete fat signal spoiling may be a source of AFI B1 mapping errors, especially in subjects with high fat content. Achieving complete fat spoiling requires longer TR, which is undesirable in clinical applications. The proposed approach based on fat signal removal can reduce errors without significant prolongation of the AFI pulse sequence. We propose that, when attaining complete fat spoiling is not feasible, AFI mapping should be performed in a multi-echo regime with appropriate fat separation or suppression to minimize these errors.

B1 inhomogeneity-corrected T1 mapping and quantitative magnetization transfer imaging via simultaneously estimating Bloch-Siegert shift and magnetization transfer effects

PURPOSE

To introduce a method of inducing Bloch-Siegert shift and magnetization Transfer Simultaneously (BTS) and demonstrate its utilization for measuring binary spin-bath model parameters free pool spin-lattice relaxation (T 1 F $$ {T}_1^{\mathrm{F}} $$ ), macromolecular fraction (f $$ f $$ ), magnetization exchange rate (k F $$ {k}_{\mathrm{F}} $$ ) and local transmit field (B 1 + $$ {B}_1^{+} $$ ).

THEORY AND METHODS

Bloch-Siegert shift and magnetization transfer is simultaneously induced through the application of off-resonance irradiation in between excitation and acquisition of an RF-spoiled gradient-echo scheme. Applying the binary spin-bath model, an analytical signal equation is derived and verified through Bloch simulations. Monte Carlo simulations were performed to analyze the method's performance. The estimation of the binary spin-bath parameters withB 1 + $$ {B}_1^{+} $$ compensation was further investigated through experiments, both ex vivo and in vivo.

RESULTS

Comparing BTS with existing methods, simulations showed that existing methods can significantly biasT 1 $$ {T}_1 $$ estimation when not accounting for transmitB 1 $$ {B}_1 $$ heterogeneity and MT effects that are present. Phantom experiments further showed that the degree of this bias increases with increasing macromolecular proton fraction. Multi-parameter fit results from an in vivo brain study generated values in agreement with previous literature. Based on these studies, we confirmed that BTS is a robust method for estimating the binary spin-bath parameters in macromolecule-rich environments, even in the presence ofB 1 + $$ {B}_1^{+} $$ inhomogeneity.

CONCLUSION

A method of estimating Bloch-Siegert shift and magnetization transfer effect has been developed and validated. Both simulations and experiments confirmed that BTS can estimate spin-bath parameters (T 1 F $$ {T}_1^{\mathrm{F}} $$ ,f $$ f $$ ,k F $$ {k}_{\mathrm{F}} $$ ) that are free fromB 1 + $$ {B}_1^{+} $$ bias.

Inhomogeneous magnetization transfer imaging: Concepts and directions for further development

Off-resonance radio frequency irradiation can induce the ordering of proton spins in the dipolar fields of their neighbors, in molecules with restricted mobility. This dipolar order decays with a characteristic relaxation time, T_1D , that is very different from the T₁ and T₂ relaxation of the nuclear alignment with the main magnetic field. Inhomogeneous magnetization transfer (ihMT) imaging is a refinement of magnetization transfer (MT) imaging that isolates the MT signal dependence on dipolar order relaxation times within motion-constrained molecules. Because T_1D relaxation is a unique contrast mechanism, ihMT may enable improved characterization of tissue. Initial work has stressed the high correlation between ihMT signal and myelin density. Dipolar order relaxation appears to be much longer in membrane lipids than other molecules. Recent work has shown, however, that ihMT acquisitions may also be adjusted to emphasize different ranges of T_1D . These newer approaches may be sensitive to other microstructural components of tissue. Here, we review the concepts and history of ihMT and outline the requirements for further development to realize its full potential.

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 (T_1D ) 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], T_2,f , T_2,b , R, and free pool T₁ ). 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 T₁ .

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 T₁ (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 ΔT₁  = -19.0 ms). The single-point methodology yielded consistent MPF and T₁ 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.

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.

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.

Macromolecular Proton Fraction as a Myelin Biomarker: Principles, Validation, and Applications

Macromolecular proton fraction (MPF) is a quantitative MRI parameter describing the magnetization transfer (MT) effect and defined as a relative amount of protons bound to biological macromolecules with restricted molecular motion, which participate in magnetic cross-relaxation with water protons. MPF attracted significant interest during past decade as a biomarker of myelin. The purpose of this mini review is to provide a brief but comprehensive summary of MPF mapping methods, histological validation studies, and MPF applications in neuroscience. Technically, MPF maps can be obtained using a variety of quantitative MT methods. Some of them enable clinically reasonable scan time and resolution. Recent studies demonstrated the feasibility of MPF mapping using standard clinical MRI pulse sequences, thus substantially enhancing the method availability. A number of studies in animal models demonstrated strong correlations between MPF and histological markers of myelin with a minor influence of potential confounders. Histological studies validated the capability of MPF to monitor both demyelination and re-myelination. Clinical applications of MPF have been mainly focused on multiple sclerosis where this method provided new insights into both white and gray matter pathology. Besides, several studies used MPF to investigate myelin role in other neurological and psychiatric conditions. Another promising area of MPF applications is the brain development studies. MPF demonstrated the capabilities to quantitatively characterize the earliest stage of myelination during prenatal brain maturation and protracted myelin development in adolescence. In summary, MPF mapping provides a technically mature and comprehensively validated myelin imaging technology for various preclinical and clinical neuroscience applications.

Data-Driven Retrospective Correction of B1 Field Inhomogeneity in Fast Macromolecular Proton Fraction and R1 Mapping

Correction of B₁ field non-uniformity is critical for many quantitative MRI methods including variable flip angle (VFA) T₁ mapping and single-point macromolecular proton fraction (MPF) mapping. The latter method showed promising results as a fast and robust quantitative myelin imaging approach and involves VFA-based R₁=1/T₁ map reconstruction as an intermediate processing step. The need for B₁ correction restricts applications of the above methods, since B₁ mapping sequences increase the examination time and are not commonly available in clinics. A new algorithm was developed to enable retrospective data-driven simultaneous B₁ correction in VFA R₁ and single-point MPF mapping. The principle of the algorithm is based on different mathematical dependences of B₁ -related errors in R₁ and MPF allowing extraction of a surrogate B₁ field map from uncorrected R₁ and MPF maps. To validate the method, whole-brain R₁ and MPF maps with isotropic 1.25 mm³ resolution were obtained on a 3 T MRI scanner from 11 volunteers. Mean parameter values in segmented brain tissues were compared between three reconstruction options including the absence of correction, actual B₁ correction, and surrogate B₁ correction. Surrogate B₁ maps closely reproduced actual patterns of B₁ inhomogeneity. Without correction, B₁ non-uniformity caused highly significant biases in R₁ and MPF ( ). Surrogate B₁ field correction reduced the biases in both R₁ and MPF to a non-significant level ( 0.1 ≤ P ≤ 0.8 ). The described algorithm obviates the use of dedicated B₁ mapping sequences in fast single-point MPF mapping and provides an alternative solution for correction of B₁ non-uniformities in VFA R₁ mapping.

Microstructure-Weighted Connectomics in Multiple Sclerosis

Introduction: Graph theory has been applied to study the pathophysiology of multiple sclerosis (MS) since it provides global and focal measures of brain network properties that are affected by MS. Typically, the connection strength and, consequently, the network properties are computed by counting the number of streamlines (NOS) connecting couples of gray matter regions. However, recent studies have shown that this method is not quantitative.

Methods: We evaluated diffusion-based microstructural measures extracted from three different models to assess the network properties in a group of 66 MS patients and 64 healthy subjects. Besides, we assessed their correlation with patients' disability and with a biological measure of neuroaxonal damage.

Results: Graph metrics extracted from connectomes weighted by intra-axonal microstructural components were the most sensitive to MS pathology and the most related to clinical disability. In contrast, measures of network segregation extracted from the connectomes weighted by maps describing extracellular diffusivity were the most related to serum concentration of neurofilament light chain. Network properties assessed with NOS were neither sensitive to MS pathology nor correlated with clinical and pathological measures of disease impact in MS patients.

Conclusion: Using tractometry-derived graph measures in MS patients, we identified a set of metrics based on microstructural components that are highly sensitive to the disease and that provide sensitive correlates of clinical and biological deterioration in MS patients.

Quantitative Macromolecular Proton Fraction Mapping Reveals Altered Cortical Myelin Profile in Schizophrenia Spectrum Disorders

Myelin abnormalities have been reported in schizophrenia spectrum disorders (SSD) in white matter. However, in vivo examinations of cortical myeloarchitecture in SSD, especially those using quantitative measures, are limited. Here, we employed macromolecular proton fraction (MPF) obtained from quantitative magnetization transfer imaging to characterize intracortical myelin organization in 30 SSD patients versus 34 healthy control (HC) participants. We constructed cortical myelin profiles by extracting MPF values at various cortical depths and quantified their shape using a nonlinearity index (NLI). To delineate the association of illness duration with myelin changes, SSD patients were further divided into 3 duration groups. Between-group comparisons revealed reduced NLI in the SSD group with the longest illness duration (>5.5 years) compared with HC predominantly in bilateral prefrontal areas. Within the SSD group, cortical NLI decreased with disease duration and was positively associated with a measure of spatial working memory capacity as well as with cortical thickness (CT). Layer-specific analyses suggested that NLI decreases in the long-duration SSD group may arise in part from significantly increased MPF values in the midcortical layers. The current study reveals cortical myelin profile changes in SSD with illness progression, which may reflect an abnormal compensatory mechanism of the disorder.

Global hypomyelination of the brain white and gray matter in schizophrenia: quantitative imaging using macromolecular proton fraction

Myelin deficiency is commonly recognized as an important pathological feature of brain tissues in schizophrenia (SZ). In this pilot study, global myelin content abnormalities in white matter (WM) and gray matter (GM) of SZ patients were non-invasively investigated using a novel clinically-targeted quantitative myelin imaging technique, fast macromolecular proton fraction (MPF) mapping. MPF maps were obtained from 23 healthy subjects and 31 SZ patients using a clinical 1.5T magnetic resonance imaging (MRI) scanner. Mean MPF in WM and GM was compared between the healthy control subjects and SZ patients with positive and negative leading symptoms using the multivariate analysis of covariance. The SZ patients had significantly reduced MPF in GM (p < 0.001) and WM (p = 0.02) with the corresponding relative decrease of 5% and 3%, respectively. The effect sizes for the myelin content loss in SZ relative to the control group were 1.0 and 1.5 for WM and GM, respectively. The SZ patients with leading negative symptoms had significantly lower MPF in GM (p < 0.001) and WM (p = 0.003) as compared to the controls and showed a significant MPF decrease in WM (p = 0.03) relative to the patients with leading positive symptoms. MPF in WM significantly negatively correlated with the disease duration in SZ patients (Pearson's r = -0.51; p = 0.004). This study demonstrates that chronic SZ is characterized by global microscopic brain hypomyelination of both WM and GM, which is associated with the disease duration and negative symptoms. Myelin deficiency in SZ can be detected and quantified by the fast MPF mapping method.

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₁.

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 T₁ , B 1 + and B₀ 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 (R_1,f T_2,f ), the cross-relaxation rate (R) and the bound pool transverse relaxation time (T_2,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 R_1,f T_2,f  = 0.0129, R = 26.5 s^-1 , and T_2,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.

Three-dimensional fast single-point macromolecular proton fraction mapping of the human brain at 0.5 Tesla

Fast single-point macromolecular proton fraction (MPF) mapping is a recent magnetic resonance imaging (MRI) method enabling quantitative assessment of myelin content in neural tissues. To date, the reported technical implementations of MPF mapping utilized high-field MRI equipment (1.5 T or higher), while low-field applications might pose challenges due to signal-to-noise ratio (SNR) limitations and short T₁ . This study aimed to evaluate the feasibility of MPF mapping of the human brain at 0.5 T. The three-dimensional MPF mapping protocol was implemented according to the single-point synthetic-reference method, which includes three spoiled gradient-echo sequences providing proton density, T₁ , and magnetization transfer contrast weightings. Whole-brain MPF maps were obtained from three healthy volunteers with spatial resolution of 1.5×1.5×2 mm³ and the total scan time of 19 minutes. MPF values were measured in a series of white and gray matter structures and compared with literature data for 3 T magnetic field. MPF maps enabled high contrast between white and gray matter with notable insensitivity to paramagnetic effects in iron-rich structures, such as globus pallidus, substantia nigra, and dentate nucleus. MPF values at 0.5 T appeared in close agreement with those at 3 T. This study demonstrates the feasibility of fast MPF mapping with low-field MRI equipment and the independence of brain MPF values of magnetic field. The presented results confirm the utility of MPF as an absolute scale for MRI-based myelin content measurements across a wide range of magnetic field strengths and extend the applicability of fast MPF mapping to inexpensive low-field MRI hardware.

Quantitative Imaging of White and Gray Matter Remyelination in the Cuprizone Demyelination Model Using the Macromolecular Proton Fraction

Macromolecular proton fraction (MPF) has been established as a quantitative clinically-targeted MRI myelin biomarker based on recent demyelination studies. This study aimed to assess the capability of MPF to quantify remyelination using the murine cuprizone-induced reversible demyelination model. MPF was measured in vivo using the fast single-point method in three animal groups (control, cuprizone-induced demyelination, and remyelination after cuprizone withdrawal) and compared to quantitative immunohistochemistry for myelin basic protein (MBP), myelinating oligodendrocytes (CNP-positive cells), and oligodendrocyte precursor cells (OPC, NG2-positive cells) in the corpus callosum, caudate putamen, hippocampus, and cortex. In the demyelination group, MPF, MBP-stained area, and oligodendrocyte count were significantly reduced, while OPC count was significantly increased as compared to both control and remyelination groups in all anatomic structures (p < 0.05). All variables were similar in the control and remyelination groups. MPF and MBP-stained area strongly correlated in each anatomic structure (Pearson’s correlation coefficients, r = 0.80–0.90, p < 0.001). MPF and MBP correlated positively with oligodendrocyte count (r = 0.70–0.84, p < 0.01 for MPF; r = 0.81–0.92, p < 0.001 for MBP) and negatively with OPC count (r = −0.69–−0.77, p < 0.01 for MPF; r = −0.72–−0.89, p < 0.01 for MBP). This study provides immunohistological validation of fast MPF mapping as a non-invasive tool for quantitative assessment of de- and remyelination in white and gray matter and indicates the feasibility of using MPF as a surrogate marker of reparative processes in demyelinating diseases.

Whole knee joint T1 values measured in vivo at 3T by combined 3D ultrashort echo time cones actual flip angle and variable flip angle methods

PURPOSE

To measure T1 relaxations for the major tissues in whole knee joints on a clinical 3T scanner.

METHODS

The 3D UTE-Cones actual flip angle imaging (AFI) method was used to map the transmission radiofrequency field (B1 ) in both short and long T2 tissues, which was then used to correct the 3D UTE-Cones variable flip angle (VFA) fitting to generate accurate T1 maps. Numerical simulation was carried out to investigate the accuracy of T1 measurement for a range of T2 values, excitation pulse durations, and B1 errors. Then, the 3D UTE-Cones AFI-VFA method was applied to healthy volunteers (N = 16) to quantify the T1 of knee tissues including cartilage, meniscus, quadriceps tendon, patellar tendon, anterior cruciate ligament (ACL), posterior cruciate ligament (PCL), marrow, and muscles at 3T.

RESULTS

Numerical simulation showed that the 3D UTE-Cones AFI-VFA technique can provide accurate T1 measurements (error <1%) when the tissue T2 is longer than 1 ms and a 150 μs excitation RF pulse is used and therefore is suitable for most knee joint tissues. The proposed 3D UTE-Cones AFI-VFA method showed an average T1 of 1098 ± 67 ms for cartilage, 833 ± 47 ms for meniscus, 800 ± 66 ms for quadriceps tendon, 656 ± 43 ms for patellar tendon, 873 ± 38 ms for ACL, 832 ± 49 ms for PCL, 379 ± 18 ms for marrow, and 1393 ± 46 ms for muscles.

CONCLUSION

The 3D UTE-Cones AFI-VFA method allows volumetric T1 measurement of the major tissues in whole knee joints on a clinical 3T scanner.

Conventional and quantitative MRI in a novel feline model of demyelination and endogenous remyelination

Background

The feeding of irradiated food to healthy adult cats results in widespread, noninflammatory demyelination of the central nervous system (CNS); a return to a normal diet results in endogenous remyelination with functional recovery. This recently discovered, reversible disease might provide a compelling clinical neuroimaging model system for the development and testing of myelin‐directed MRI methods as well as future remyelination therapies.

Purpose

Identify the noninvasive imaging characteristics of this new disease model and determine whether it features measurable changes on conventional and quantitative MRI.

Study Type

Pilot study.

Animal Model

Ten adult cats at various stages of demyelinating disease induced by an irradiated diet (35–55 kGy), and during recovery following a return to a normal diet.

Field Strength/Sequence

Conventional (T₂‐weighted) and quantitative (diffusion tensor, magnetization transfer) at 3T.

Assessment

MRI of the brain, optic nerves, and cervical spinal cord; a subset of diseased cats was euthanized for comparative histopathology.

Statistical Tests

Descriptive statistics.

Results

Disease produced T₂ prolongation, progressing from patchy to diffuse throughout most of the cerebral white matter (eventually involving U‐fibers) and spinal cord (primarily dorsal columns, reminiscent of subacute combined degeneration but without evidence of B12 deficiency). Magnetization transfer parameters decreased by 50–53% in cerebral white matter and by 25–30% in optic nerves and spinal cord dorsal columns. Fractional diffusion anisotropy decreased by up to 20% in pyramidal tracts, primarily driven by increased radial diffusivity consistent with axon preservation. Histopathology showed scattered myelin vacuolation of major white matter tracts as well as many thin myelin sheaths consistent with remyelination in the recovery phase, which was detectable on magnetization transfer imaging.

Data Conclusion

Feline irradiated diet‐induced demyelination features noninvasively imageable and quantifiable demyelination and remyelination of the CNS. It is therefore a compelling clinical neuroimaging model system.

Level of Evidence: 4

Technical Efficacy: Stage 2

J. Magn. Reson. Imaging 2019;49:1304–1311.

Phase-sensitive B1 mapping: Effects of relaxation and RF spoiling

PURPOSE

To develop a phase-based B₁ mapping technique accounting for the effects of imperfect RF spoiling and magnetization relaxation.

THEORY AND METHODS

The technique is based on a multi-gradient-echo sequence with 2 successive orthogonal radiofrequency (RF) excitation pulses followed by the train of gradient echoes measurements. We have derived a theoretical expression relating the MR signal phase produced by the 2 successive RF pulses to the B₁ field and B₀ -related frequency shift. The expression takes into account effects of imperfections of RF spoiling and T₁ and T2* relaxations.

RESULTS

Our computer simulations and experiments revealed that imperfections of RF spoiling cause significant errors in B₁ mapping if not accounted for. By accounting for these effects along with effects of magnetization relaxation and frequency shift, we demonstrated the high accuracy of our approach. The technique has been tested on spherical phantoms and a healthy volunteer.

CONCLUSION

In this paper, we have proposed, implemented, and demonstrated the accuracy of a new phase-based technique for fast and robust B₁ mapping based on the measured MR signal phase, frequency, and relaxation. Because imperfect RF spoiling effects are accounted for, this technique can be applied with short TRs and therefore substantially reduces the scan time. Magn Reson Med 80:101-111, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Optimization of inhomogeneous magnetization transfer (ihMT) MRI contrast for preclinical studies using dipolar relaxation time (T1D ) filtering

A pulsed inhomogeneous magnetization transfer (ihMT)-prepared fast imaging sequence was implemented at 11.75 T for preclinical studies on mouse central nervous system. A strategy based on filtering the ihMT signal originating from short dipolar relaxation time (T_1D ) components is proposed. It involves increasing the repetition time of consecutive radiofrequency (RF) pulses of the dual saturation and allows improved signal specificity for long T_1D myelinated structures. Furthermore, frequency offset, power and timing saturation parameters were adjusted to optimize the ihMT sensitivity. The optimization of the ihMT sensitivity, whilst preserving the strong specificity for the long T_1D component of myelinated tissues, allowed measurements of ihMT ratios on the order of 4-5% in white matter (WM), 2.5% in gray matter (GM) and 1-1.3% in muscle. This led to high relative ihMT contrasts between myelinated tissues and others (~3-4 between WM and muscle, and ≥2 between GM and muscle). Conversely, higher ihMT ratios (~6-7% in WM) could be obtained using minimal T_1D filtering achieved with short saturation pulse repetition time or cosine-modulated pulses for the dual-frequency saturation. This study represents a first stage in the process of validating ihMT as a myelin biomarker by providing optimized ihMT preclinical sequences, directly transposable and applicable to other preclinical magnetic fields and scanners. Finally, ihMT ratios measured in various central nervous system areas are provided for future reference.

Propagation of error from parameter constraints in quantitative MRI: Example application of multiple spin echo T2 mapping

PURPOSE

Quantitative MRI may require correcting for nuisance parameters which can or must be constrained to independently measured or assumed values. The noise and/or bias in these constraints propagate to fitted parameters. For example, the case of refocusing pulse flip angle constraint in multiple spin echo T₂ mapping is explored.

METHODS

An analytical expression for the mean-squared error of a parameter of interest was derived as a function of the accuracy and precision of an independent estimate of a nuisance parameter. The expression was validated by simulations and then used to evaluate the effects of flip angle (θ) constraint on the accuracy and precision of T⁁2 for a variety of multi-echo T₂ mapping protocols.

RESULTS

Constraining θ improved T⁁2 precision when the θ-map signal-to-noise ratio was greater than approximately one-half that of the first spin echo image. For many practical scenarios, constrained fitting was calculated to reduce not just the variance but the full mean-squared error of T⁁2, for bias in θ⁁≲6%.

CONCLUSION

The analytical expression derived in this work can be applied to inform experimental design in quantitative MRI. The example application to T₂ mapping provided specific cases, depending on θ⁁ accuracy and precision, in which θ⁁ measurement and constraint would be beneficial to T⁁2 variance or mean-squared error. Magn Reson Med 79:673-682, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Fast Realistic MRI Simulations Based on Generalized Multi-Pool Exchange Tissue Model

We present MRiLab, a new comprehensive simulator for large-scale realistic MRI simulations on a regular PC equipped with a modern graphical processing unit (GPU). MRiLab combines realistic tissue modeling with numerical virtualization of an MRI system and scanning experiment to enable assessment of a broad range of MRI approaches including advanced quantitative MRI methods inferring microstructure on a sub-voxel level. A flexible representation of tissue microstructure is achieved in MRiLab by employing the generalized tissue model with multiple exchanging water and macromolecular proton pools rather than a system of independent proton isochromats typically used in previous simulators. The computational power needed for simulation of the biologically relevant tissue models in large 3D objects is gained using parallelized execution on GPU. Three simulated and one actual MRI experiments were performed to demonstrate the ability of the new simulator to accommodate a wide variety of voxel composition scenarios and demonstrate detrimental effects of simplified treatment of tissue micro-organization adapted in previous simulators. GPU execution allowed  ∼ 200× improvement in computational speed over standard CPU. As a cross-platform, open-source, extensible environment for customizing virtual MRI experiments, MRiLab streamlines the development of new MRI methods, especially those aiming to infer quantitatively tissue composition and microstructure.

High-resolution three-dimensional macromolecular proton fraction mapping for quantitative neuroanatomical imaging of the rodent brain in ultra-high magnetic fields

A well-known problem in ultra-high-field MRI is generation of high-resolution three-dimensional images for detailed characterization of white and gray matter anatomical structures. T₁-weighted imaging traditionally used for this purpose suffers from the loss of contrast between white and gray matter with an increase of magnetic field strength. Macromolecular proton fraction (MPF) mapping is a new method potentially capable to mitigate this problem due to strong myelin-based contrast and independence of this parameter of field strength. MPF is a key parameter determining the magnetization transfer effect in tissues and defined within the two-pool model as a relative amount of macromolecular protons involved into magnetization exchange with water protons. The objectives of this study were to characterize the two-pool model parameters in brain tissues in ultra-high magnetic fields and introduce fast high-field 3D MPF mapping as both anatomical and quantitative neuroimaging modality for small animal applications. In vivo imaging data were obtained from four adult male rats using an 11.7T animal MRI scanner. Comprehensive comparison of brain tissue contrast was performed for standard R₁ and T₂ maps and reconstructed from Z-spectroscopic images two-pool model parameter maps including MPF, cross-relaxation rate constant, and T₂ of pools. Additionally, high-resolution whole-brain 3D MPF maps were obtained with isotropic 170µm voxel size using the single-point synthetic-reference method. MPF maps showed 3-6-fold increase in contrast between white and gray matter compared to other parameters. MPF measurements by the single-point synthetic reference method were in excellent agreement with the Z-spectroscopic method. MPF values in rat brain structures at 11.7T were similar to those at lower field strengths, thus confirming field independence of MPF. 3D MPF mapping provides a useful tool for neuroimaging in ultra-high magnetic fields enabling both quantitative tissue characterization based on the myelin content and high-resolution neuroanatomical visualization with high contrast between white and gray matter.

Rapid multicomponent relaxometry in steady state with correction of magnetization transfer effects

PURPOSE

To study the effects of magnetization transfer (MT) on multicomponent T2 parameters obtained using mcDESPOT in macromolecule-rich tissues and to propose a new method called mcRISE to correct MT-induced biases.

METHODS

The two-pool mcDESPOT model was modified by the addition of an exchanging macromolecule proton pool to model the MT effect in cartilage. The mcRISE acquisition scheme was developed to provide sensitivity to all pools. An incremental fitting was applied to estimate MT and relaxometry parameters with minimized coupling. The interaction between MT and relaxometry parameters, efficacy of MT correction, and feasibility of mcRISE in vivo were investigated in simulations and in healthy volunteers.

RESULTS

The MT effect caused significant errors in multicomponent T1/T2 values and in fast-relaxing water fraction fF , which is consistent with previous experimental observations. fF increased significantly with macromolecule content if MT was ignored. mcRISE resulted in a multifold reduction of MT biases and yielded decoupled multicomponent T1/T2 relaxometry and quantitative MT parameters.

CONCLUSION

mcRISE is an efficient approach for correcting MT biases in multicomponent relaxometry based on steady state sequences. Improved specificity of mcRISE may help to elucidate the sources of the previously described high sensitivity of noncorrected mcDESPOT parameters to disease-related changes in cartilage and the brain.

MR fingerprinting using the quick echo splitting NMR imaging technique

PURPOSE

The purpose of the study is to develop a quantitative method for the relaxation properties with a reduced radio frequency (RF) power deposition by combining magnetic resonance fingerprinting (MRF) technique with quick echo splitting NMR imaging technique (QUEST).

METHODS

A QUEST-based MRF sequence was implemented to acquire high-order echoes by increasing the gaps between RF pulses. Bloch simulations were used to calculate a dictionary containing the range of physically plausible signal evolutions using a range of T₁ and T₂ values based on the pulse sequence. MRF-QUEST was evaluated by comparing to the results of spin-echo methods. The specific absorption rate (SAR) of MRF-QUEST was compared with the clinically available methods.

RESULTS

MRF-QUEST quantifies the relaxation properties with good accuracy at the estimated head SAR of 0.03 W/kg. T₁ and T₂ values estimated by MRF-QUEST are in good agreement with the traditional methods.

CONCLUSIONS

The combination of the MRF and the QUEST provides an accurate quantification of T₁ and T₂ simultaneously with reduced RF power deposition. The resulting lower SAR may provide a new acquisition strategy for MRF when RF energy deposition is problematic. Magn Reson Med 77:979-988, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Minimizing the effects of magnetization transfer asymmetry on inhomogeneous magnetization transfer (ihMT) at ultra-high magnetic field (11.75 T)

OBJECTIVES

The recently reported inhomogeneous magnetization transfer technique (ihMT) has been proposed for specific imaging of inhomogeneously broadened lines, and has shown great promise for characterizing myelinated tissues. The ihMT contrast is obtained by subtracting magnetization transfer images obtained with simultaneous saturation at positive and negative frequency offsets (dual frequency saturation experiment, MT (+/-)) from those obtained with single frequency saturation (MT (+)) at the same total power. Hence, ihMT may be biased by MT-asymmetry, especially at ultra-high magnetic field. Use of the average of single positive and negative frequency offset saturation MT images, i.e., (MT (+)+MT (-)) has been proposed to correct the ihMT signal from MT-asymmetry signal.

MATERIALS AND METHODS

The efficiency of this correction method was experimentally assessed in this study, performed at 11.75 T on mice. Quantitative corrected ihMT and MT-asymmetry ratios (ihMTR and MTRasym) were measured in mouse brain structures for several MT-asymmetry magnitudes and different saturation parameter sets.

RESULTS

Our results indicated a "safe" range of magnitudes (/MTRasym/<4 %) for which MT-asymmetry signal did not bias the corrected ihMT signal. Moreover, experimental evidence of the different natures of both MT-asymmetry and inhomogeneous MT contrasts were provided. In particular, non-zero ihMT ratios were obtained at zero MTRasym values.

CONCLUSION

MTRasym is not a confounding factor for ihMT quantification, even at ultra-high field, as long as MTRasym is restricted to ±4 %.

Time-efficient, high-resolution, whole brain three-dimensional macromolecular proton fraction mapping

PURPOSE

Macromolecular proton fraction (MPF) mapping is a quantitative MRI method that reconstructs parametric maps of a relative amount of macromolecular protons causing the magnetization transfer (MT) effect and provides a biomarker of myelination in neural tissues. This study aimed to develop a high-resolution whole brain MPF mapping technique using a minimal number of source images for scan time reduction.

METHODS

The described technique was based on replacement of an actually acquired reference image without MT saturation by a synthetic one reconstructed from R1 and proton density maps, thus requiring only three source images. This approach enabled whole brain three-dimensional MPF mapping with isotropic 1.25 × 1.25 × 1.25 mm(3) voxel size and a scan time of 20 min. The synthetic reference method was validated against standard MPF mapping with acquired reference images based on data from eight healthy subjects.

RESULTS

Mean MPF values in segmented white and gray matter appeared in close agreement with no significant bias and small within-subject coefficients of variation (<2%). High-resolution MPF maps demonstrated sharp white-gray matter contrast and clear visualization of anatomical details, including gray matter structures with high iron content.

CONCLUSIONS

The proposed synthetic reference method improves resolution of MPF mapping and combines accurate MPF measurements with unique neuroanatomical contrast features.

MPnRAGE: A technique to simultaneously acquire hundreds of differently contrasted MPRAGE images with applications to quantitative T1 mapping

PURPOSE

To introduce a new technique called MPnRAGE, which produces hundreds of images with different T1 contrasts and a B1 corrected T1 map.

THEORY AND METHODS

An interleaved three-dimensional radial k-space trajectory with a sliding window reconstruction is used in conjunction with magnetization preparation pulses. This work modifies the SNAPSHOT-FLASH T1 fitting equations for radial imaging with view-sharing and develops a new rapid B1 correction procedure. MPnRAGE is demonstrated in phantoms and volunteers, including two volunteers with eight scans each and eight volunteers with two scans each. T1 values from MPnRAGE were compared with those from fast spin echo inversion recovery (FSE-IR) in phantoms and a healthy human brain at 3 Tesla (T).

RESULTS

The T1 fit for human white and gray matter was T1MPnRAGE = 1.00 · T1FSE-IR + 24 ms, r(2) = 0.990. Voxel-wise coefficient of variation in T1 measurements across eight time points was between 0.02 and 0.08. Region of interest-based T1 values were reproducible to within 2% and agree well with literature values.

CONCLUSION

In the same amount of time as a traditional MPRAGE exam (7.5 min), MPnRAGE was shown to produce hundreds of images with alternate T1 contrasts as well as an accurate and reproducible T1 map that is robust to B1 errors.

MR fingerprinting using fast imaging with steady state precession (FISP) with spiral readout

PURPOSE

This study explores the possibility of using gradient echo-based sequences other than balanced steady-state free precession (bSSFP) in the magnetic resonance fingerprinting (MRF) framework to quantify the relaxation parameters .

METHODS

An MRF method based on a fast imaging with steady-state precession (FISP) sequence structure is presented. A dictionary containing possible signal evolutions with physiological range of T1 and T2 was created using the extended phase graph formalism according to the acquisition parameters. The proposed method was evaluated in a phantom and a human brain. T1 , T2 , and proton density were quantified directly from the undersampled data by the pattern recognition algorithm.

RESULTS

T1 and T2 values from the phantom demonstrate that the results of MRF FISP are in good agreement with the traditional gold-standard methods. T1 and T2 values in brain are within the range of previously reported values.

CONCLUSION

MRF-FISP enables a fast and accurate quantification of the relaxation parameters. It is immune to the banding artifact of bSSFP due to B0 inhomogeneities, which could improve the ability to use MRF for applications beyond brain imaging.

Removal of cerebrospinal fluid partial volume effects in quantitative magnetization transfer imaging using a three-pool model with nonexchanging water component

PURPOSE

Parameters of the two-pool model describing magnetization transfer (MT) in macromolecule-rich tissues may be significantly biased in partial volume (PV) voxels containing cerebrospinal fluid (CSF). The purpose of this study was to develop a quantitative MT (qMT) method that provides indices insensitive to CSF PV averaging.

THEORY AND METHODS

We propose a three-pool MT model, in which PV macro-compartment is modeled as an additional nonexchanging water pool. We demonstrate the feasibility of model parameter estimation from several MT-weighted spoiled gradient echo datasets. We validated the three-pool model in numerical, phantom, and in vivo studies.

RESULTS

PV averaging with the free water compartment reduces all qMT parameters, most significantly affecting macromolecular proton fraction (MPF) and cross-relaxation rate. Monte-Carlo simulations confirmed stability of the three-pool model fit. Unlike the standard two-pool model, the three-pool model qMT parameters were not affected by PV averaging in simulations and phantom studies. The three-pool model fit allowed CSF PV correction in brain PV voxels and resulted in good correlation with standard two-pool model parameters in non-PV voxels.

CONCLUSION

Quantitative MT imaging based on a three-pool model with a non-exchanging water component yields a set of CSF-insensitive qMT parameters, which may improve MPF-based assessment of myelination in structures strongly affected by CSF PV averaging such as brain gray matter.

Cross-relaxation imaging of human patellar cartilage in vivo at 3.0T

OBJECTIVE

To compare quantitative magnetization transfer (qMT) parameters of patellar cartilage measured using cross-relaxation imaging (CRI) in asymptomatic volunteers and patients with osteoarthritis.

DESIGN

The study was performed with Institutional Review Board approval and with all subjects signing informed consent. CRI of the knee joint was performed at 3.0T on 20 asymptomatic volunteers and 11 patients with osteoarthritis. The fraction of macromolecular bound protons (f), the exchange rate constant between macromolecular bound protons and free water protons (k), and the T2 relaxation time of macromolecular bound protons (T2(B)) of patellar cartilage were measured. Mann-Whitney-Wilcoxon rank-sum tests were used to compare qMT parameters between asymptomatic volunteers and patients with osteoarthritis.

RESULTS

Average f, k, and T2(B) of patellar cartilage was 12.46%, 7.22 s(-1), and 6.49 μs respectively for asymptomatic volunteers and 12.80%, 6.13 s(-1), and 6.80 μs respectively for patients with osteoarthritis. There were statistically significant differences between groups of subjects for k (P < 0.01) and T2(B) (P < 0.0001) but not f (P = 0.38) of patellar cartilage.

CONCLUSION

Patients with osteoarthritis had significantly lower k and significantly higher T2(B) of patellar cartilage than asymptomatic volunteers which suggests that qMT parameters can detect changes in the macromolecular matrix of degenerative cartilage.