Magnetization Transfer (MT) has been extensively applied to study macromolecular biological tissue composition. The imaging contrast resides in the magnetization transfer between free-water protons and macromolecular proton compartments, through chemical exchange and dipolar interactions. In the two-pool tissue model, highly mobile protons are associated with the free-water pool while protons found in semisolid macromolecular sites are defined as the restricted pool (Sled 2018). A simple method for visualizing MT effects includes acquiring two images with and without an off-resonance MT pulse to calculate the MT ratio (MTR), which is the normalized difference of these two images (Wolff and Balaban 1989). Despite its proven usefulness to study multiple sclerosis (Zheng et al. 2018), Alzheimer’s disease (Fornari et al. 2012) and psychiatric disorders (Chen et al. 2015), the MTR is a semi-quantitative metric that depends critically on the imaging sequence parameters (Wood and Malik 2020). Another semi-quantitative approach is the estimation of MT saturation (MTsat) by fitting the MT signal obtained from an MT-weighted (MTw), a proton density (PD) weighted and T_{1}-weighted (T_{1}w) contrast (Helms et al. 2008). Quantitative MT (qMT) consists of fitting multiple images to a mathematical model to extract tissue-specific parameters related to physical quantities, such as pool sizes, magnetization exchange rates between pools, and T_{1},T_{2} relaxation times of each pool. Compared to semi-quantitative approaches (MTR, MTsat), qMT has long acquisition protocols and sometimes needs additional measurements (eg. B_{0}, B_{1}, T_{1}), and the complex models required to fit the quantitative maps makes it a challenging imaging technique.

In 2015, our lab published qMTLab (Cabana et al. 2015), an open-source software project seeking to unify three qMT methods in the same interface: qMT using spoiled gradient echo (qMT-SPGR), qMT using balanced steady-state free precession (qMT-bSSFP), and qMT using selective inversion recovery with fast spin echo (qMT-SIRFSE). qMTLab allowed users to simulate, evaluate, fit, and visualize qMT data with the possibility to share qMT protocols between researchers, allowing them to compare the performance of their methods (Cabana et al. 2015). Since then, we have extended the project and renamed it to qMRLab, (Karakuzu et al. 2020), which in addition to qMT now provides over 20 quantitative techniques under one umbrella, such as relaxation and diffusion models, quantitative susceptibility mapping, B_{0} and B_{1} mapping. In addition, we created interactive tutorials (Boudreau 2018a; Boudreau 2018b; Boudreau 2019) and blog posts for several qMRI techniques that were published under a creative commons license, and one of these tutorials (Boudreau 2018a) made it into a quantitative MRI book published by Elsevier (Seiberlich et al. 2020). This blog post is a continuation of our qMRLab outreach initiative, where we will focus on qMT and the tools we provide in qMRLab for this class of techniques.

Below is an introduction to qMT (with a focus on qMT-SPGR), where we will cover signal modelling and data fitting.

The magnetization transfer in a two-pool model is modelled by a set of coupled differential equations (Sled and Pike 2000):

where the magnetization at time $t$ is given by M = [M_{x,f}, M_{y,f}, M_{z,f}, M_{z,r}] and M_{x,f} and M_{y,f} are the transverse magnetization for the free pool in the $x$ and $y$ direction, respectively. The longitudinal magnetization for the free and restricted pool is denoted by M_{z,f} and M_{z,r}. Due to the very short T_{2} of the restricted pool (on the order of microseconds), the transverse magnetization of this pool is not explicitly modelled.

The constants k_{f} and k_{r} represent the exchange rate of the longitudinal magnetization from the free pool to the restricted pool (k_{f}) and from the restricted to the free pool (k_{r}). The ratio of these quantities is constrained by the size of the restricted to the free-water pool, expressed as k_{f}/k_{r} = M_{0,r}/M_{0,f}. This ratio is called the pool size ratio $F$ defined as F = M_{0,r}/M_{0,f}. The precessional frequency ω_{1} is a measure of the power of the off-resonance radiofrequency pulse and ∆ is a frequency offset at which the magnetic B_{1}^{+} field is applied. As shown in Figure 1, an on-resonance (∆ = 0) radiofrequency pulse is also part of the pulse sequence and is applied after the MT pulse. The relaxation time constants for the free and restricted pools are denoted by T_{1,f} and T_{1,r} for the longitudinal magnetization, and T_{2,f} and T_{2,r} for the transverse magnetization. Finally, $W$ is the saturation rate of the restricted pool, which is a function of the absorption lineshape $G$ that depends on the frequency offset and the transverse relaxation of the restricted pool (Figure 2).

In qMT imaging, the biophysical model relates the parameters observed in the two-pool tissue model to physical quantities such as the fractional size of the pools, relaxation times and magnetization exchange rates of the free and restricted pool (Sled and Pike 2001; Sled 2018). However, qMT experiments usually consist of long acquisition imaging protocols accompanied by complex data fitting. To this end, some software solutions have been proposed (Karakuzu et al. 2020; Tabelow et al. 2019; Wood 2018). qMRLab is an open-source project for quantitative MR analysis that is an extension of qMTLab, a software for data simulation and analysis of three MT models: qMT-SPGR, qMT-bSSFP and qMT-SIRFSE. In addition to the quantitative MT methods, qMRLab also contains semi-quantitative MT models including the magnetization transfer ratio (MTR) and magnetization transfer saturation (MTsat).

The qMT-SPGR method in qMRLab contains three fitting models: Sled and Pike, Ramani, and Yarnykh and Yuan (Sled and Pike 2001; Ramani et al. 2002; Yarnykh 2002). For the Sled and Pike model, the saturation fraction effect of the MT pulse on the free pool is pre-computed to accelerate the processing times. The MT effect of the pulse is approximated as an instantaneous fractional saturation of the longitudinal magnetization of the free pool, assuming the absence of chemical exchange processes (Pike 1996; Sled and Pike 2001). To fit the model, additional parameters related to the pulse sequence are required, namely timing parameters, the absorption lineshape, and the characteristics of the MT pulse, such as the shape and the bandwidth or the time-bandwidth product. In Figure 5, the qMT-SPGR method is used to show a single voxel curve simulation for the same MT data fitted by three different models. The fitted parameters were the pool size ratio F, the magnetization transfer rate from the restricted to the free-water pool (k_{r}), and the transverse relaxation time of the free-water (T_{2,f}) and restricted (T_{2,r}) pool.

_{0} inhomogeneity affects the actual off-resonance frequency experienced by the tissue at each voxel. In MT, B_{0} maps are computed to correct for off-resonance frequency values of the MT pulse in the presence of magnetic field non-uniformity. Radiofrequency field B_{1} inhomogeneity is another source of inaccuracies that depend on the operating frequency of the scanner, the pulse sequence, and the shape and electrical properties of the sample (Sled and Pike 1998). Therefore, B_{1} maps are typically needed to correct the radiofrequency amplitude variations that affect the nominal values of the MT pulse flip angle and the excitation flip angle (Boudreau et al. 2018c). Longitudinal relaxation T_{1} values vary naturally in biological tissue, but the choice of the T_{1} mapping method, has also been shown to influence the variability of T_{1} measurements (Stikov et al. 2015). In the context of a qMT experiment, T_{1} maps are acquired with an independent measurement of the apparent relaxation time T_{1} (T_{1}^{meas}). The T_{1} map is related to the relaxation rate of the free pool (R_{1,f}) as described by equation 6, expressed in terms of $F$, k_{f}, R_{1}^{meas} and R_{1,r}, where the relaxation rate of the restricted pool is arbitrarily set to 1 s^{-1} because it is insensitive to this kind of MT experiments (Sled and Pike 2001). Multiple qMRI maps with a range of B_{0} and B_{1} inaccuracies, as well as T_{1} maps with a variety of relaxation times, have been simulated to show the effect of the quality of these input maps on the qMT fitted parameters as shown in Figure 6.

As described above in Figure 6, inaccurate MT pulse flip angles and excitation flip angles affect the fitted MT parameters, and there is an additional error source related to the T_{1} mapping measurement. As Boudreau et al. (2018c) have shown, using specific acquisitions protocols, T_{1} values can be affected by B_{1} field non-uniformities, such as the variable flip angle method, while the inversion recovery method is insensitive to these field inhomogeneities (Stikov et al. 2015).

Figure 7 displays an example human dataset with the input qMRI maps used to fit the qMT parameters F, k_{f}, T_{2,f}, T_{2,r}.