- Please check and comment entries here.
Pre-clinical Cerebral Perfusion MRI Techniques
Alterations to the cerebral microcirculation have been recognized to play a crucial role in the development of neurodegenerative disorders. However, the exact role of the microvascular alterations in the pathophysiological mechanisms often remains poorly understood. The early detection of changes in microcirculation and cerebral blood flow (CBF) can be used to get a better understanding of underlying disease mechanisms. This could be an important step towards the development of new treatment approaches. Animal models allow for the study of the disease mechanism at several stages of development, before the onset of clinical symptoms, and the verification with invasive imaging techniques. Specifically, pre-clinical magnetic resonance imaging (MRI) is an important tool for the development and validation of MRI sequences under clinically relevant conditions. This article reviews MRI strategies providing indirect non-invasive measurements of microvascular changes in the rodent brain that can be used for early detection and characterization of neurodegenerative disorders. The perfusion MRI techniques: Dynamic Contrast Enhanced (DCE), Dynamic Susceptibility Contrast Enhanced (DSC) and Arterial Spin Labeling (ASL), will be discussed, followed by less established imaging strategies used to analyze the cerebral microcirculation: Intravoxel Incoherent Motion (IVIM), Vascular Space Occupancy (VASO), Steady-State Susceptibility Contrast (SSC), Vessel size imaging, SAGE-based DSC, Phase Contrast Flow (PC) Quantitative Susceptibility Mapping (QSM) and quantitative Blood-Oxygenation-Level-Dependent (qBOLD). We will emphasize the advantages and limitations of each strategy, in particular on applications for high-field MRI in the rodent’s brain.
1. Dynamic Contrast Enhanced MRI
DCE-MRI is an exogenous contrast-based technique, mainly used to provide an estimation of the perfusion and the permeability of the cerebral microvasculature . Using DCE-MRI, T1 changes of tissue over time, following the introduction of a paramagnetic or superparamagnetic contrast agent, can be measured. For paramagnetic contrast agents, the acquisition consists of a baseline image (T10 map), followed by a series of T1-weighted images after an intravenous bolus of contrast agent. After injection, the contrast agent spreads through the tissue resulting in a change in the MR signal intensity proportional to the concentration of the contrast agent (T1 shortening). These temporal changes provide a signal intensity curve from which physiological parameters of the microvascular system such as perfusion, vessel permeability and the extravascular-extracellular space of the tissue can be derived .
There are two main approaches to analyze DCE measurements quantitatively, parametric and nonparametric approaches. The nonparametric (model free) or semi-quantitative techniques measure empirical metrics directly from the signal intensity curve such as the Bolus Arrival Time, Time-To-Peak, the maximal signal intensity or the area under the signal Attenuation Curve. If the MR signal is converted to a signal concentration curve, the latter is referred to as the initial area under the gadolinium curve. These metrics provide information about the shape and the structure of the MR signal intensity and can often be correlated to the underlying physiology . Nonparametric approaches are suitable for fast and simple diagnostics and have been used in various (pre-)clinical studies. However, nonparametric techniques have some important limitations. They are sensitive to several parameters, such as the MRI acquisition protocol (sequence and parameters), type of scanner, injection protocol and type of contrast agent, making comparison of the results difficult . Furthermore, the contribution of the physiological parameters to the MR signal intensity remains unclear .
Parametric or quantitative studies estimate the physiological parameters that can be directly related to the physiological properties . The parametric evaluation of the signal intensity attenuation curve is done by converting the MR signal intensity to a time contrast concentration curve and fitting it to a model. In order to perform this kinetic modeling, the Arterial Input Function (AIF), describing the concentration changes in a blood vessel entering the tissue as a function over time, has to be identified. The accurate identification of the AIF is far from straightforward. A detailed description of the determination of the AIF can be found elsewhere .
Most of the studies use the two-compartmental pharmacokinetic (PK) model comprising the intravascular-extracellular space and the extravascular-extracellular space introduced by Tofts et al. (Figure 1) . The two-compartment model estimates the volume transfer constant between the blood plasma and the extravascular-extracellular space (Ktrans); the reflux exchange rate between the extravascular-extracellular space and blood plasma (kep) and the volume of the EES per unit volume of tissue (νe = Ktrans/kep) .
Figure 1. Schematic representation of the two-compartmental pharmacokinetic (PK) Tofts model comprising an intra-vascular compartment (left) and extravascular -extracellular space (EES) compartment (right). The volume transfer between the two compartments is depicted by the transfer constant between the blood plasma and the EES (Ktrans) and the transfer constant between the EES and blood plasma (kep), respectively. kep can be calculated using kep = Ktrans /νe, whith νe volume of the EES per unit volume of tissue.
The transfer constant Ktrans is the most widely used kinetic parameter. It is a complex function depending on the blood flow, the vascular surface area per unit mass of tissue (representing the size and number of blood vessels) and the microvascular permeability, making it dependent on the tissue and its physiology . The dependence of Ktrans on both perfusion and microvascular permeability complicates the interpretation of the results, which is a major limiting factor of DCE-MRI. More complex DCE models, to describe the underlying tissue physiology more accurately, have been proposed . However, the higher degree of complexity of these models, and the increased number of parameters to be estimated, often results in increased errors in the parameters estimation .
DCE-MRI is frequently used in both clinical and pre-clinical studies. Parametric DCE-MRI requires the fast and accurate determination of the T1 values. In practice, a trade-off must often be made between the temporal resolution, spatial resolution, the SNR and the field of view . The limited spatial resolution in combination with increased susceptibility artifacts at higher field strength and fast changes in the blood plasma contrast agent often result in difficulties in the identification of the AIF in rodents . In pre-clinical studies, it is especially important to optimize the temporal resolution, allowing accurate T1 determination, without reducing the spatial resolution. To decrease the relatively long acquisition times in MRI, several acceleration techniques have been proposed . The use and implementation of acceleration techniques with high-field small-animal MRI systems has remained limited so far. Better implementations and increased use of acceleration techniques on high-field small-animal MRI scanners are needed to increase the reliability and reproducibility of the quantitative kinetic parameters in future pre-clinical DCE-MRI studies.
Since it is especially challenging to obtain a good identification of the AIF in small animals, several alternative approaches and models have been proposed. The AIF can also be determined by sampling arterial blood. However, due to their small blood volume, repeated blood sampling can have a significant effect on the physiological status of the animal . Furthermore, the low sampling rate and the inability to collect blood samples close to the tissue of interest often confound the identification of the AIF . An average population-based AIF is often used to avoid the need for a subject-specific AIF , but this approach does not take inter-subject variability into account.
Completely non-invasive methods have been developed to determine the AIF on an individual basis from the DCE-MRI . These methods need high temporal resolution and require often the presence of a large vessel within the field of view. A potentially interesting alternative approach, which compares the obtained curve with that of a reference region in healthy tissue, is the reference region model . This approach is frequently used in rodents and allows for quantitative analysis without the need of an AIF. Since this technique does not require a high temporal resolution to determine the AIF, a higher spatial resolution and/or SNR can be achieved. This has mainly been used in oncology to obtain quantitative measurements of perfusion in brain tumors .
Additional problems can arise due to the difficulty and lack of standardization of the contrast agent injection in rodents . The administration of the contrast agent is typically performed through injection in the jugular vein or tail vein, which is often challenging. Changes in the bolus of the contrast agent can have a significant effect on the quantitative parameters . Therefore, an automated injection system is preferred to increase the reproducibility compared to a manual injection.
2. Dynamic Susceptibility Contrast Enhanced MRI
The DSC MRI technique, also known as bolus tracking, is the most widely used MRI technique to measure brain perfusion . Similar to DCE MRI, the acquisition consists of a pre-contrast baseline image followed by a series of MR images after an intravenous bolus of contrast agent. In DSC MRI, the tissue perfusion is assessed by evaluation of a series of rapidly repeated T2- or T2*-weighted MR images resulting from the first pass of a contrast agent bolus through the tissue . During the first pass, the contrast agent is mainly confined to the intravascular space allowing a good estimation of the perfusion. The susceptibility induced signal loss over time (T2/T2* shortening), proportional to the contrast agent concentration, provides a signal intensity attenuation curve. The DSC-MRI technique is straightforward, has short acquisition times and has a high contrast to noise ratio compared to other perfusion methods. Its main limitation is the need for the identification of the AIF for quantitative analysis. In rodents, the capillary blood flow is about five times higher than in humans. Therefore, first pass measurements in rodents require extremely short bolus injections and acquisition times, making the MR signals much more sensitive to T1 changes and susceptibility artifacts .
Similar to DCE, the signal intensity attenuation curve can be used to derive semi-quantitative parameters. The transformation into the contrast concentration curve allows quantitative evaluation of the following physiological parameters: CBV, defined as the fraction of tissue volume occupied by blood (ml/g); CBF (ml/g/min) and MTT defined as the average time the contrast agent travels through the vasculature of the brain tissue (s) . Figure 2 shows an example of a typical tissue contrast concentration curve and explains how it can be used to derive the physiological parameters. The absolute quantification of the physiological parameters is strongly dependent on the AIF, and most of the artifacts in the quantification can be linked directly to an incorrect identification of the AIF . Since absolute quantification is difficult, the relative parameters calculated from the contrast concentration curve without the identification of the AIF are often used instead.
Figure 2. Parameters derivation in Dynamic Susceptibility Contrast Enhanced (DSC) MRI: (a) Semi-quantitative interpretation can be derived from the signal intensity curve. The Arrival Time (AT) of the bolus can be determined from the time interval between the injection of the contrast agent and the time point where the contrast agent is first detected in the tissue. The Time-To-Peak (TTP) is determined by the time interval between the injection of the contrast agent and the peak of the contrast agent in the tissue. The Full Width at Half Maximum (FWHM) is dependent on the Mean Transit Time (MTT) of the tissue. (b) the contrast concentration curve can be calculated by deconvolution of the Arterial Input Function (AIF) and the signal intensity curve. The Cerebral Blood Flow (CBF) is determined by the maximum height of the contrast concentration curve (Tmax). The Cerebral Blood Volume (CBV) is determined by the area under the contrast concentration curve, while the MTT can be calculated using MTT = CBV/CBF. Reproduced with permission from Leif Østergaard, Journal of Magnetic Resonance Imaging, published by John Wiley and Sons, 2005 .
The DSC MRI technique relies on the assumption that the contrast agent remains intravascular. As the contribution of the permeability and the volume of the extravascular-extracellular space increases with the time after the bolus injection, DSC MRI requires a high temporal resolution. The disruption of the BBB can introduce errors in the measurements. In order to correct for the leakage of the contrast agent, a preload of contrast agent can be administered , or a leakage-correction algorithm can be used . Recently, double or multi-echo approaches to remove the T1 effects caused by leakage of the contrast agent have gained increasing interest . In addition to the leakage-correction, multi-echo DSC MRI was shown to improve identification of the AIF , estimation of the transfer constant Ktrans and the volume of the extravascular-extracellular space per unit volume of tissue νe . The multi-echo approach thus allows the extraction of the kinetic parameters and the conventional DSC parameters without the need for multiple contrast agent injections (Table 2) .
Both spin-echo (SE, T2) and gradient-echo (GE, T2*) sequences can be used for DSC MRI. The susceptibility contrast in GE images arises from the contributions of both macro-and microvasculature . Therefore, it is preferable to use SE images, where the signal loss is mainly sensitive to the microvasculature . Recently, combination of the two techniques referred to as Spin- and Gradient-Echo Echo-Planar Imaging (SAGE EPI) was proposed by Schmiedeskamp et al. .
3. Arterial Spin Labeling
ASL is the most frequently used, completely non-invasive, cerebral perfusion MRI technique (Table 1). It is based on the labeling, often referred to as ‘tagging’, of water spins in the arterial blood supply. The labeled blood water spins are used as an endogenous contrast agent. Several techniques combining different preparations and readout schemes exist for ASL. The general ASL scheme consists of two consecutive acquisitions. Before the first acquisition, the arterial blood water spins are tagged ‘upstream’ from the tissue of interest using an inversion pulse. After a delay time (Inversion Time, TI), during which the tagged spins travel to the tissue of interest and exchange with the stationary spins in the extravascular-extracellular space, an image is acquired (Figure 3A). The second acquisition, which serves as a control, is identical to the first acquisition without the tagging of the arterial water spins . The difference between both images provides a signal proportional to the exchanged water magnetization and therefore the arterial blood flow to the tissue after the delay time. This difference is therefore directly proportional to the capillary blood flow in the tissue . The signal intensities from the perfusion weighted images can then be converted to quantitative measurements of CBF in physiological units of flow (mL/g/min) . The ASL technique was originally developed to provide information about CBF, but recently several methods have been developed that enable the estimation of other physiological parameters such as the bolus arrival time and arterial CBV (Table 2) .
Figure 3. Arterial Spin Labeling (ASL): (A) The general ASL scheme consists of two consecutive acquisitions. The first acquisition consists of the tagging of the arterial blood spins (Left), followed by a post-labeling delay time, during which the tagged spins flow to the tissue of interest and exchange with the stationary spins in the extravascular-extracellular space (Middle), and finally, the image is acquired (Right). The second acquisition, which serves as a control, is identical to the first acquisition without the tagging of the arterial water spins; (B) Schematic representation of the most widely used ASL sequences. Continuous ASL (CASL), uses a continuous Radiofrequency (RF) pulse in combination with a constant imaging gradient in the direction of the arterial blood flow to invert the water spins of arterial blood (Top). In Pulsed ASL (PASL), the labeling is achieved using a pulse or a pulse train of short RF inversion pulses (Middle). pseudo-Continuous ASL (pCASL) uses a train of short RF pulses that mimic the long continuous RF pulse from CASL (Top).
Several ASL sequences with different labeling schemes have been developed. The most widely used ASL sequences, continuous ASL (CASL) , pulsed ASL (PASL)  and pseudo-continuous ASL (pCASL), are shown in Figure 3B .
The CASL sequence uses a continuous Radiofrequency (RF) pulse in combination with a constant imaging gradient in the direction of the arterial blood flow to invert the arterial blood water spins. The main drawback of CASL is the requirement of long labeling RF pulses causing signal loss due to strong magnetization transfer (MT) effects and high specific absorption rates , historically limiting the applications to the acquisition of a single slice. To correct for MT effects in multiple slices, a separate labeling coil can be used (two-coil ASL)  or additional RF pulses during the control acquisition can be applied . In rats, a dedicated neck labeling coil has been used to minimize MT effects and improve the SNR . Due to their smaller size, two-coil CASL has proven to be more challenging in mice. When the dedicated coil is placed on the neck of a mouse, the short distance between the two coils will result in partial saturation of the brain signal. To overcome this limitation, an alternative approach was developed where the dedicated labeling coil is positioned at the heart region of the animal . Even though CASL has a higher sensitivity than its pulsed counterpart, the use of CASL is limited due to the need for hardware providing continuous RF, a dedicated labeling coil and a high specific absorption rate.
In PASL, the labeling is achieved using a pulse or a pulse train of short RF inversion pulses. Depending on how the labeling is applied, the PASL sequences can be subdivided into two groups: the techniques that perform the labeling symmetrical with respect to the measured plane, called flow alternating inversion recovery (FAIR) , and the asymmetrical PASL sequences . The PASL techniques have a higher inversion efficiency, a lower specific absorption rate and smaller MT effects than CASL. This, in combination with its ease of implementation and robustness, makes PASL the most widely used ASL sequence . Due to the ease of implementation, FAIR PASL is the most frequently used PASL technique in pre-clinical studies.
More recently, the pCASL sequence was developed to combine the advantages of CASL and PASL . It uses a train of short RF pulses that mimic the long continuous RF pulse from CASL. This way, it combines the high SNR of PASL with the high labeling efficiency of CASL without the need for specific hardware to generate a long continuous labeling pulse. pCASL is less dependent on flow velocity and shows better reproducibility than PASL and CASL . Furthermore, it also shows better control of the MT effects and increased label efficiency compared to CASL . The implementation of multi-slice pCASL is straightforward and does not require a dedicated labeling coil. Disadvantages are the lower label efficiency and higher specific absorption rate compared to PASL and CASL. Furthermore, pCASL has an increased sensitivity to magnetic field inhomogeneities in the labeling plane . This becomes an increasing problem at the higher field strengths used in pre-clinical small-animal MRI. Even though pCASL is the recommended implementation of ASL for clinical applications , it is rarely used in pre-clinical research due to the lack of commercially available pCASL sequences on high field pre-clinical systems .
The major disadvantage of ASL in general is the low intrinsic SNR. Only a small part of the brain tissue volume consists of blood (2–4%), while the rest of the volume is filled with stationary tissue. Therefore, the signal of the labelled blood water spins represents only of a very small fraction of the overall water volume in the tissue. Additionally, the longitudinal T1 relaxation time of tagged blood water spins at clinical field strengths is similar to the arterial transit time . This causes the labeling of the blood water spins to decrease by the time the measurement is performed resulting in lower SNR. The blood longitudinal relaxation time increases linearly with the magnetic field strength . Therefore, an increase in the magnetic field results in a decreased labeling loss due to the decay of the blood water spins. On the one hand, the high magnetic field strength and the increased gradient strength on pre-clinical scanners improve the SNR and the resolution of the ASL images . On the other hand, the higher magnetic field strength results in an increased T1 for gray matter, limiting SNR improvement . The reduction in the T2* of blood, due to stronger susceptibility effects of paramagnetic deoxygenated hemoglobin at higher fields strengths, can lead to a decreased contribution of the capillaries to the ASL signal, which can lead to an underestimation of the CBF . Furthermore, the increased magnetic field strength leads to increased contributions of magnetic field inhomogeneities in the labeling plane. The combination of these effects and several artifacts, such as subtraction errors, motion artifacts and susceptibility artifacts, make ASL a low SNR technique. To improve the SNR, several advanced pulse sequence or readout techniques, such as background suppression techniques , B0-correction , MT effect reduction  and motion correction techniques , have been proposed.
Nevertheless, ASL is increasingly used in pre-clinical research, in particular in models of neurological and neurodegenerative diseases . In contrast to exogenous-perfusion-weighted imaging methods, ASL can be used to measure relatively rapid changes in CBF, for example, in response to changed physiological parameters such as pCO2 or pH. Exogenous-perfusion-weighted imaging methods would require first a clearance of the contrast agent from the blood pool. Such ASL experiments have been used to determine the cerebral vascular response to hypercapnia/hypoventilation (Figure 4), which was proposed as an earlier marker for the detection of vascular disfunction in neurodegenerative diseases .
Figure 4. ALS experiments in control mice (Ctrl) and an Alzheimer disease model (biAT) according to . (A) Illustration of the experimental paradigm. Animals were intubated, and ventilation was controlled. ALS measurements commenced during a period of normoventilation. In order to assess the cerebral vascular response (CVR) to a hypercapnic challenge, animals were hypoventilated. This was followed by a control measurement under normoventilation. Expired CO2 levels (expCO2) were recorded during the experiments. More details on the protocol can be found in . (B) CBF maps and calculation of the CVR based on the baseline CBF and the CBF acquired under hypoventilation. (C) Comparison of CVR in the cortex (Ctx), hippocampus (hipp) and thalamus (thal) of control mice and biAT mice. This figure is based on data presented in .
Recently, a novel labeling scheme called velocity-selective ASL (vsASL) has been proposed . In vsASL, the labeling is based on the velocity of the blood flow instead of the spatial location of the arterial spins. This method was specifically developed to address the cerebral blood flow in situations of slow flow, where the arterial transit time can be longer than the T1 relaxation time. However, to our knowledge, this method has not been implemented on a high field pre-clinical scanner.
4. Intravoxel Incoherent Motion
Intravoxel Incoherent Motion (IVIM) MRI, first described by Le Bihan et al. , is a diffusion MRI technique that can provide information about both molecular diffusion and perfusion at the same time . The term IVIM stands for the motion of water within a voxel during an MR experiment. In biological tissue, this motion consists of contributions from the molecular diffusion of water and the microcirculation of blood in the microvascular network . Originally, IVIM MRI was introduced based on the assumption that blood flow in randomly oriented capillaries mimics a random walk, similar to the Brownian motion of water molecules. This so-called pseudo-diffusion process, characterized by the pseudo diffusion coefficient (D*), results in an additional signal attenuation in diffusion-weighted imaging measurements, making diffusion MRI sensitive to both diffusion and perfusion. In the presence of diffusion gradients, both phenomena result in a decay of the signal intensity. The signal attenuation increases with the degree of diffusion weighting that is applied (b-value) . However, during the typically applied short diffusion encoding times, the IVIM signal cannot be attributed solely to the randomized motion in the microvasculature network. The signal attenuation of the perfusion component will also have a contribution of phase dispersion caused by incoherent blood flow velocities within the microvascular network.
The pseudo-diffusion coefficient, D*, related to the perfusion in the microvascular network, is much higher than the molecular diffusion of water in the parenchymal tissue. Therefore, D* results in a faster decay of the signal attenuation. The faster signal decay, in combination with the small fraction capillary blood flow in the tissue, limits the contribution of the pseudo-diffusion to low b-values. Hence, the signal attenuation at high b-values is solely a consequence of diffusion effects in the extravascular-extracellular space. In order to capture both phenomena, the acquisition of an IVIM experiment consists of a conventional diffusion-weighted sequence, where a wide range of b-values are used. Compared to ‘conventional’ diffusion-weighted imaging, a number of measurements at low b-values are added to visualize the IVIM effect resulting from the pseudo-diffusion in the microcirculation (Figure 5).
Figure 5. The general principle of the two-step fitting approach in Intravoxel Incoherent Motion (IVIM) MRI: (A) In a first step, the signal attenuation at high b-values is fitted (red). Since the signal attenuation at high b-values only contains contributions from the molecular diffusion of water in the extravascular-extracellular space (EES), this can be used to estimate the diffusion coefficient. In the second step, the data from all b-values are used to extract the other IVIM parameters. The contribution of the pseudo-diffusion in the microvascular network (blue) can be visualized by subtracting the diffusion contribution from the total IVIM signal (black-red); (B) Coronal B0 map in a rat and; (C) Corresponding diffusion coefficient map, obtained from the first step in the two-step fitting approach. The images were acquired from a mouse brain using a single shot gradient spin-echo Echo Planar Imaging (EPI) sequence at a magnetic field strength of 9.4 T (Echo Time 37 ms, Repetition Time 1000 ms, flip angle 90°, 3 orthogonal diffusion gradient directions, gradient separation (Δ) 1.6 ms, gradient duration (δ) 7.9 ms and 30 b-values 0–1000 s/mm2).
In general, the signal attenuation can be expressed as:
where S(b) represents the signal intensity at a certain b-value; S0 the signal intensity without diffusion sensitizing gradients (b = 0 s/mm2); Fdiff (b) and Fperf (b) correspond to the diffusion and the perfusion component; fIVIM the fraction of flowing blood in a voxel (%) and (1 − fIVIM) the extravascular-extracellular space where only diffusion effects are considered, respectively.
Fitting the overall signal attenuation at multiple b-values allows the determination of fIVIM and D* (Table 2). The IVIM parameters fIVIM and D* provide information about the contribution of perfusion to the MRI signal and can be related to classical perfusion parameters . The original IVIM model describes both the perfusion and the diffusion component by a mono-exponential model resulting in the classical bi-exponential model for the signal attenuation. Over the past few years, several fitting approaches and different models for both the diffusion and the perfusion component have been explored.
The diffusion component can be approximated using a mono-exponential model. However, this is not correct in the case of biological tissue. This non-Gaussian behavior becomes visible at high b-values. Due to the small scale of the IVIM effects, this can result in significant errors in the parameter estimation. Examples of frequently used models that address this non-Gaussian behavior are the Kurtosis Model , which is known to fail at very high b-values , and the bi-exponential diffusion model containing a fast and a slow diffusion component .
Originally, two perfusion IVIM models were proposed by Le Bihan et al. . The mono-exponential model assumes that the direction of the blood flow changes several times during the measurement. This only holds true when long diffusion times are used or under conditions of fast flow and short vessel segments . At short diffusion times or under conditions of slow flow and long vessel segments, the blood will remain in a single segment during the measurement. Under these conditions, referred to as the ballistic limit, the signal attenuation can be more accurately fitted using a sinc model . However, in reality, the situation lies somewhere in between these two models. Several other models have been proposed to fit the more realistic intermediate situation .
Multi-compartment models have been shown to fit signal attenuation for the perfusion component more accurately , leading to several multi-compartment models being proposed . Recently, an in vivo rat brain study showed that the signal attenuation of the perfusion component can be more accurately fitted using a bi-exponential IVIM model that splits the perfusion component into a slow vascular pool, resulting from blood flow in the capillaries and a fast vascular pool, resulting from blood flow in the medium-sized vessel .
The fitting of the signal attention can be performed using several curve-fitting approaches. The most common methods are the one-step and the two-step approach. The one-step approach fits signal attenuation at all acquired b-values and simultaneously extracts the IVIM parameters using a nonlinear least-squares algorithm. This approach is easy to use, but is very sensitive to noise and local minima and often results in an inaccurate estimation of the parameters. The more robust two-step approach relies on the magnitude difference between the diffusion coefficient (D) and D*. Since the pseudo-diffusion only has a significant effect at low b-values, the signal attenuation at high b-values can be fitted to extract D (Figure 5). In the second step, the data from all b-values are used to get fIVIM and D*. The main limitation of these two-step approaches is that a threshold separating the perfusion and the diffusion component has to be chosen. The optimal threshold can vary strongly between the organ measured and different pathologies. To overcome the need to select a cut-off threshold manually, an iterative algorithm to select the optimal threshold was proposed . Recently, more complex fitting methods such as the Bayesian-fitting approach have gained interest .
The improvement of MRI hardware and software, improving the SNR and the acquisition time, has quickly expanded the field of IVIM MRI in recent years. One of the key features of IVIM is that it can provide quantitative information without the use of contrast agents and the need for identification of an AIF. The IVIM technique is a promising tool providing complementary information to the classical perfusion measurements (Table 1). In particular, the addition of IVIM to ASL measurements could provide additional diagnostic information without the administration of a contrast agent. The perfusion in several organs, such as the liver , kidney  and pancreas , has been successfully studied using IVIM MRI. Due to the small cerebral perfusion fraction, the implementation and the correlation with classical perfusion measurements in the brain has been proven to be more difficult .
Several factors influence the accuracy and reproducibility of IVIM. Apart from the different processing approaches, there is great diversity in the acquisition schemes used for IVIM MRI. In practice, the b-values and their distribution are often chosen heuristically. In general, an increase in the number of b-values results in a decreased error in the parameter estimation. However, the amount of b-values is often limited due to long acquisition times. Furthermore, it is often difficult to obtain accurate measurements at low b-values due to instabilities in the gradient amplifiers and MRI hardware limitations . Currently, an increasing number of efforts are ongoing to optimize the number and distribution of b-values .
Applications of IVIM in rodents are limited. The higher SNR resulting from the ultra-high magnetic field used in small animal MRI systems improves the accuracy of the fit and the parameter estimation. This allows for the identification of small differences, which is crucial when measuring cerebral perfusion. However, the limitations mentioned above in combination with the wide variations in magnetic field strength used in rodents result in significant variance in the parameters across studies. There is therefore need for more uniform acquisition schemes and post-processing approaches. This can lead to increased reliability and reproducibility of the IVIM parameters.
5. Vascular Space Occupancy
Vascular Space Occupancy (VASO) is a relatively new endogenous MRI technique that allows assessment of changes in the CBV based on the intrinsic T1 differences of tissue and blood (Table 1) . A schematic representation of the VASO principle is shown in Figure 6. First, the spins of both blood and tissue are inverted using a spatially non-slice-selective global inversion pulse. After the global inversion pulse, the longitudinal magnetization relaxes back to its equilibrium at the spin-specific T1 relaxation rate . The difference in T1 relaxation between tissue and blood is expressed in a time difference for their magnetization curves to cross zero. The VASO technique exploits this difference using an optimal inversion time (TI) for blood, referred to as blood nulling. When an image is acquired at the TI, the only signal which is detectable arises from the extracellular tissue. Assuming a constant water volume per voxel, the signal arising from the extravascular tissue is proportional to (1-CBV) . The optimal TI can be determined from the repetition time (TR) and the T1 of blood.
Figure 6. Schematic representation of the basic principle of Vascular Space Occupancy (VASO) MRI: First, the spins of both blood and tissue are inverted by the use of a spatially non-slice-selective global 180° inversion pulse. After the global inversion pulse, the longitudinal magnetization of blood and tissue relax back at their spin-specific T1 relaxation rate. At the optimal Inversion Time (TI), when the longitudinal magnetization of the tissue crosses zero, an additional 90° excitation pulse is applied to null the blood signal. Then, an image is acquired that contains only signals that arise from the extravascular tissue.
The VASO signal mainly arises from the differences in CBV. However, the contrast can be affected by several factors such as the CBF, the extravascular Blood Oxygenation Level Dependent (BOLD) effect, the contribution of cerebrospinal fluid, inflow effects and magnetization transfer effects. Some of these confounding contributions can be reduced by optimizing the acquisition parameters. The contribution of the BOLD effect can be minimized by using the shortest possible TE, while the TR should be chosen to be long enough to avoid strong CBF contributions. Inflow caused by non-steady-state inflowing blood spins can be minimized by the use of RF coils that cover the entire body or the application of magnetization-reset and crushing-gradient techniques . Due to exchange of water between blood and tissue, the inverted blood spins will affect the signal from the extravascular tissue. This will result in a lower SNR. This signal contrast between blood and tissue can be enhanced by the application of a magnetization transfer pulse .
The global inversion pulse used in VASO results in a low SNR. In general, the SNR is improved when higher magnetic field strengths are used. However, in VASO, there are several confounding factors that counteract this improvement. First of all, at higher magnetic field strengths, the difference between the T1 of tissue and blood will become smaller, reducing the sensitivity. Secondly, the higher field strength results in increased extravascular BOLD effects, counteracting the negative VASO signal and further reducing the sensitivity. This could underestimate CBV values. Furthermore, the efficiency of the blood nulling is limited due to increased field inhomogeneities, and the acquisition will become more prone to geometrical distortion and inflow artifacts . To overcome these limitations, an alternative method using a slab-selective inversion instead of a global inversion was proposed . This was further extended to slice-saturation slab-inversion VASO by applying additional RF pulses before the slab-selective inversion, which further improved the CNR compared to the original VASO method .
The main limitation of VASO is that it does not provide an absolute measurement of CBV. Therefore, VASO has been mainly used to measure CBV changes in functional MRI studies. In order to provide quantitative information about the CBV values, an alternative approach based on the T1 shortening effect of T1 contrast agents was developed . Contrast-based VASO is a steady-state technique that relies on the assumption that the contrast agent remains confined to the intravascular space. Therefore, the T1 shortening effect is limited to the blood compartment. The acquisition consists of identical pre-and post-contrast experiments. Due to the T1 shortening, the blood signal from the post-contrast measurement will no longer be nulled at TI. Since the contrast agent remains intravascular, the T1 weighting of the tissue will not be affected in between the two experiments. The difference between the pre-and post-contrast images can be used to calculate the absolute CBV .
Compared to DSC MRI, VASO MRI has several advantages and limitations. Since it is a steady-state technique, contrast-enhanced VASO does not require extremely short bolus injections and acquisition times. This allows a higher spatial resolution compared to a dynamic imaging method. Another advantage is that lower SNR of contrast-based VASO can be improved by signal averaging over different acquisitions, which is not possible in DSC . Finally, contrast-enhanced VASO is less prone to geometric distortions and susceptibility artifacts and does not require the estimation of the AIF. The main disadvantage of contrast-based VASO is that it only provides information about the CBV, while DSC also provides information on the CBF and MTT.
Recently, a new approach called inflow-VASO (iVASO) was introduced . In iVASO, only the water spins flowing into the slice are nulled using spatially slice-selective inversion. The TI for the nulling of the blood water spins is similar to the transit time of arterial blood water spins to reach the capillaries. Therefore, the iVASO technique is no longer sensitive to the total CBV, but is mainly sensitive to the contribution of the arterial and arteriolar compartment. Some of the main limitations of the original VASO technique can be overcome by the modification of the VASO sequence. By using a spatially slice-selective inversion pulse, the tissue signal is not affected by the blood nulling, which significantly increases the SNR. The iVASO technique does not depend on the difference between the T1 of the tissue, making it more interesting for use at higher magnetic fields strengths. Furthermore, since no inversion is applied in the imaging slice, iVASO is less dependent on T1 changes in the tissue. This makes iVASO a more suitable technique to study CBV in pathological conditions. Finally, iVASO significantly reduces the contribution of cerebrospinal fluid.
The iVASO signal changes are highly dependent on the choice of the TI and TR. When short TR is used, iVASO will become mainly sensitive to arteriolar blood water spin effects. At longer TR, capillary blood water spin effects will play a role. Due to its dependence on the TI, iVASO requires prior knowledge about the arterial transit time to provide absolute measurements of the arterial CBV . To overcome this limitation, a quantitative iVASO method was introduced . This approach uses the consecutive acquisition of a control scan without blood nulling and an iVASO acquisition where the inflowing blood water spins are nulled. Since the tissue magnetization is the same in both acquisitions, the subtraction of these two images provides an absolute measurement of the arterial CBV. The pulse sequence and the subtraction procedure are similar to ASL. However, the fundamental difference is that in iVASO, the water spins are nulled, whereas they are tagged in ASL .
The quantitative iVASO methods allow the completely non-invasive absolute measurements of CBV with a higher SNR than the originally proposed VASO technique. Quantitative iVASO is a promising tool to study both normal tissue and tissue with altered perfusion. It provides an interesting addition or alternative to the conventional methods for measuring CBV (Table 2). To date, VASO MRI was mainly applied in clinical research. However, due to the advantages of the more recently developed iVASO, it may become a useful tool in future pre-clinical research.
6. Steady-State Susceptibility Contrast MRI
Steady-state contrast enhanced MRI (SSC MRI) is an exogenous MRI technique that uses a contrast agent with a long (blood) half-life in order to maintain a more stable contrast agent concentration during the experiment (Table 1). Most SSC MRI studies use ultra-small super-paramagnetic iron oxide nanoparticles (USPIOs) as blood pool contrast agents. USPIOs have a longer (blood) half-life than the frequently used gadolinium chelates, allowing measurements closer to steady-state conditions. Furthermore, they have a larger hydrodynamic diameter than gadolinium chelates and are therefore more likely to remain confined to the intravascular space .
In SSC MRI, the susceptibility-induced contrast difference between the intravascular and extravascular compartments can be used to assess the CBV. First, the transverse relaxation rate R2 (SE) or R2* (GE) is measured before the injection of a contrast agent and under the steady-state conditions after the administration of a contrast agent is measured. Then, the linear relation between the changes in the transverse relaxation rate and the intravascular contrast agent concentration are used to generate an index, the relative CBV (rCBV).
Unlike DSC MRI, there is no need for rapid acquisition during the first pass of the bolus. Therefore, in the steady-state approach, a higher SNR and thus a higher spatial resolution can be achieved. The long half-life of the contrast agent allows for dynamic CBV measurements in functional MRI studies . The superparamagnetic properties of the USPIOs result in significant enhancement transverse relaxation rates compared to gadolinium base contrast agents. The larger diameter of iron oxide particles might result in slower leakage out of the vessel, even in case of BBB disruptions, reducing the error in the rCBV measurements .
A limitation of SSC MRI is that no information about the CBF and MTT can be obtained (Table 2). SSC MRI assumes that the susceptibility-induced contrast difference is independent of its neighboring voxels. However, this could lead to an overestimation of rCBV in the vicinity of large vessels. The steady-state approach has been commonly used in pre-clinical studies to measure rCBV. The major limiting factor for translation to the clinic is the lack of clinically approved contrast agents with a long enough half-life for steady-state measurements.
7. Vessel Size Imaging
Vessel size imaging is an exogenous MRI technique that provides a measurement of the average microvessel size within a voxel by comparing changes in the ΔR2 (GE) and ΔR2* (SE) relaxation rates induced by the injection of a blood pool contrast agent. The quantitative information about the average vessel diameter within a voxel is represented by the Vessel Size Index (VSI) (Table 2). Vessel size imaging can provide valuable information about the microvascular density that was previously only available from non-invasive measurements such as biopsies .
Vessel size imaging exploits the difference in response of the ΔR2 and ΔR2* transversal relaxation rate, as the vessel diameter increases . The dimensionless ratio ΔR2*/ΔR2, which increases with the vessel diameter, can be used to provide information on the average vessel size within a voxel . However, this dimensionless ratio depends on the concentration of the contrast agent. To avoid this dependence, the mean vessel density (Q), defined as the ratio ΔR2/(ΔR2*)2/3, was introduced. Under conditions of appropriately high contrast agent concentrations and long echo time, Q has shown a good correlation with histological measurements of the vessel density . Expanding on this, the VSI (µm) was introduced to provide a quantitative measurement of the average vessel size within a voxel. In order to measure the VSI from the ΔR2* and ΔR2 quantitatively, other parameters such as the water diffusion rate, D, the contrast agent concentration and the absolute blood volume fraction (CBV) should be considered.
Vessel size imaging has two types of applications, steady-state vessel size imaging  and dynamic vessel size imaging . Similar to SSC MRI, the steady-state method uses a contrast agent with a long blood half-life. The absolute CBV cannot be measured in the steady-state approach. However, under the assumption that the contrast agent remains intravascular during the measurement, and the concentration of the contrast agent in the blood is known, the blood volume fraction can be determined from ΔR2*. When the contrast agent concentration, which is proportional to the susceptibility difference between blood and brain tissue Δχ, and the diffusion coefficient D are determined, the absolute VSI can be determined using the relation (ΔR2/ΔR2*)3/2 (D/Δχ)1/2. The steady-state approach involves an invasive measurement of the blood contrast agent concentration, eliminating the need to measure absolute CBV. The value of D can be determined by performing a separate diffusion measurement or chosen based on prior knowledge. An advantage of the steady-state approach is the higher achievable spatial resolution using longer acquisition times. Similar to SSC MRI, the steady-state approach is mainly used for pre-clinical research, due to the lack of clinically approved contrast agents with a long enough blood half-life.
The dynamic approach is an extension of DSC MRI that acquires both GE and SE measurements during the first pass of the bolus of a contrast agent. In the dynamic approach it is impossible to measure the concentration of the contrast agent in the blood. However, in contrast to the steady-state approach, the absolute blood volume fraction can be determined from the susceptibility difference between blood and brain tissue Δχ obtained from ΔR2* measurements, as in conventional DSC, eliminating the need to measure the concentration of the contrast agent in the blood. When the absolute CBV and the diffusion coefficient D are determined, the quantitative VSI can be obtained using ΔR2/(ΔR2*)3/2 (CB D)1/2.
The dynamic approach requires the acquisition of both GE and SE measurements during the first pass of the contrast agent bolus. This limits the spatial resolution that can be achieved. In order to improve the spatial resolution, a method using separate GE and SE acquisitions and a dual injection of contrast agent was proposed . However, this technique has several other disadvantages such as the need of a higher dose of contrast agent, longer acquisition times and the potential influence of the first injection of the contrast agent on the second measurement.
To date, vessel size imaging has mainly been used in tumor research. However, vessel size imaging can also be explored as a promising tool in other neurodegenerative and neurological diseases where changes in the vessel size or density are expected. The VSI provides additional information on the mechanism behind perfusion changes in the brain complementing other perfusion measurements.
8. SAGE-Based DSC MRI
Recently, a combined SAGE Echo-Planer Imaging (SAGE EPI) MRI pulse sequence that allows the simultaneous measurement of spin-echo and gradient-echo DSC MRI without the need for sequential acquisition of two experiments was developed . Using a multi-echo approach, this technique provides measures of the typical DSC MRI parameters (CBF, CBV, MTT), the kinetic permeability parameters (Ktrans, νe) and the relative vessel size (VSI) (Table 2). Unlike conventional DSC measurements, the SAGE-based approach is independent of T1 effects .
The SAGE-based DSC MRI technique combines the advantages of several methods into a single acquisition. It provides complementary perfusion measures obtained from first pass measurements using a clinically approved contrast agent, making SAGE-based DSC MRI a promising tool to study neurodegenerative diseases (Table 1). However, SAGE-based DSC MRI also has its limitations.
Some of the limitations, such as the need to determine an AIF and geometric distortions, the dependence on the injection protocol and the need for extremely short acquisition times, have already been mentioned in the description of the conventional perfusion techniques. In addition to these limitations, the use of DSC MRI has been limited by the need for short TEs to measure changes in R2 and R2* accurately . In particular, the first TE used in SAGE-based MRI should be short enough, since the T2 and T2* become very short during the bolus of the contrast agent. Parallel imaging methods can be used to shorten the TE. However, on high-field pre-clinical small-animal systems, often a limited number of channels is available, limiting the acceleration that can be obtained using parallel imaging. Therefore, alternative approaches using partial Fourier encoding have been proposed . The implementation of other acceleration methods such as the keyhole method  and multiband or simultaneous multi-slice excitation  is highly beneficial for future pre-clinical studies, especially since several acceleration methods can be combined with each other and with parallel imaging, further reducing the scan time and the TE .
9. Phase Contrast Flow MRI
Phase Contrast MRI (PC MRI), also referred to as velocity mapping, is an endogenous MRI technique that uses flow-encoding gradients to visualize and quantify the velocity of moving fluids . It relies on the dephasing of moving spins when they are subjected to a bipolar gradient. The net phase shift of the moving spins will be proportional to the velocity of the spins along the direction of the bipolar gradient. Spins in the same direction of the bipolar gradient will obtain a positive net phase shift, while spins moving in the opposite direction will acquire a negative net phase shift . In order to compensate for unwanted phase shifts induced by other sequence parameters and to remove the background signal, the phase images are subtracted by a flow compensated reference image that was acquired with the same acquisition parameters but with an inverted bipolar gradient . The PC MRI technique is only sensitive to flow along the direction of the bipolar gradient. In order to obtain measurements of flow in arbitrary directions, flow encoding gradients can be applied along all three orthogonal directions. However, the addition of measurements along multiple axes significantly increases the acquisition time for PC MRI measurements.
Since the velocity-encoded images contain flow measured both in the direction off and in the opposite direction of the bipolar gradient, a phase range spanning from -π to π is chosen. Velocities corresponding to shifts larger than |π| will induce phase wrapping or aliasing in the velocity encoding image. Therefore, the maximum velocity that can be measured in a PC MRI experiment is limited. The maximum velocity that can be measured without the occurrence of aliasing at a certain gradient strength is expressed by the velocity encoding parameter (VENC). Before a PC MRI experiment, the VENC, which is inversely proportional to the gradient strength, should be chosen so that the maximum velocity corresponds to a phase shift of 180° to avoid aliasing. The correct choice of the VENC parameter, which should be estimated based on the velocities of interest, is crucial for the accuracy of PC MRI measurements. When the VENC is too high, the flow in the velocity-encoded image will be compressed in a small range of phase shifts. The inability to distinguish between small velocity differences will decrease the SNR and the accuracy of the images . This might lead to inaccurate measurements, especially in regions with slow flow . A too-low VENC will result in aliasing, making quantitative measurements difficult. The VENC is inversely proportional to the strength of the bipolar gradient and can be adapted by changing the bipolar gradient strength . However, a decrease in the field strength will be accompanied by a decreased SNR.
The main applications of PC MRI in the brain are flow measurements of cerebral spinal fluid  and vascular imaging using phase contrast MR angiography (PC MRA) . When introduced, the PC MRA technique was mainly used to visualize large vessels with active blood flow non-invasively . This can be achieved with 3D cine PC MRA or 4D PC MRA, which is discussed in detail elsewhere .
Rather than just generating angiograms, PC MRA can be used to measure the velocity and the direction of blood flow within the vascular network. In order to quantify the total global CBF in the brain non-invasively, the major cerebral feeding vessels 2D PC MRA are used . This is done by placing a thin imaging slice perpendicular to the main feeding arteries of the brain. In 2D PC MRA, the total CBF can be measured using either cardiac-gated or non-cardiac-gated acquisition, which significantly reduces the acquisition time . The total CBF, containing the entire blood supply to the brain, obtained with PC MRA can be used to normalize other CBF mapping techniques such as DSC  and ASL  to provide absolute CBF maps. This diminishes some confounding factors, such as the need to identify an AIF or the labeling efficiency, typically hindering the absolute CBF quantification . When used in combination with measurements for quantification of blood oxygenation, it can be used to estimate the cerebral metabolic rate of oxygen (CMRO2) . The CMRO2 is a key measure of cerebral functioning. Alterations in CMRO2 have been suggested to play an important role in several neurological disorders (19). In contrast to Positron Emission Tomography (PET), which is considered as the golden standard for CMRO2 quantification, PC MRI provides a completely non-invasive way to measure CMRO2.
Several limiting factors should be considered when 2D PC MRA is used to measure flow velocities. In order to describe the velocity accurately, the acquisition slice should be positioned perpendicular to the vessel orientation . When the slice is not positioned perpendicularly to the vessel, saturation effects might occur due to in-plane flow. Furthermore, the acquisition slice should be placed in a straight vessel segment with laminar flow to avoid intravoxel dephasing. Another limiting factor is the need to determine the optimal VENC, which requires prior knowledge or a good approximation of the maximum velocities in the vessel of interest.
10. Quantitative Susceptibility Mapping
Susceptibility Weighted Imaging (SWI) exploits the magnetic susceptibility differences between tissues, which are often the cause of artifacts in MR images, to obtain images with increased contrast . It exploits the complementary information about the structure and function of tissues contained in phase images by combining magnitude and phase data (Figure 7). The acquisition usually consists of a flow compensated T2*- weighted 3D gradient echo sequence. SWI is commonly used for tissue characterization based on susceptibility differences caused by deoxygenated blood, iron deposition and calcification in neurological disorders . The SWI approach has several limitations. It is highly dependent on imaging parameters, suffers from blooming artifacts and has the orientation dependence of the phase signal. This limits SWI to mainly qualitative measurements . To overcome these limitations, Quantitative Susceptibility Mapping (QSM) was developed (Table 1). QSM is a post-processing technique that quantifies the underlying magnetic susceptibilities based on the phase images.
Figure 7. Processing steps to generate Susceptibility Weighted Images (SWI) from magnitude and phase data. The raw phase data are filtered to remove low frequency fluctuations. Then, a phase mask that scales the filtered phase images to a range from 0 to 1 is created. This phase mask is then multiplied several times (n) with the magnitude image to generate the SWI with enhanced contrast.
The computation of the magnetic susceptibility consists of a three-step process: First, the magnetic field is estimated from the raw phase data using a phase unwrapping algorithm . Next, in order to get a map solely consisting of the susceptibility sources inside the brain or region of interest, the contributions from outside the field of view are removed by background field removal . Finally, the inverse problem from field perturbation to magnetic susceptibility has to be solved.
The magnetic field in perturbation in each voxel, caused by the magnetization of tissue when placed into an external magnetic field, can be approximated as a magnetic dipole producing a dipole field that extends beyond the voxel itself. Therefore, the magnetic field perturbation in a certain voxel is a superposition of its own dipole field and that of its neighboring voxels . When the susceptibility distribution is known, the field perturbation can be obtained by the convolution of the susceptibility distribution and the field of a unit dipole (i.e., dipole kernel) . However, multiple susceptibility distribution could result in the same field perturbation, making this is an ill-posed inversion problem. A simple kernel division would cause errors that would be represented as streaking artifacts in the reconstructed susceptibility map . Several methods to overcome this ill-posed inversion problem have been proposed.
The ill-posed problem can be solved analytically by repeated acquisition of the magnetic field perturbation at different orientations relative to the magnetic field. This method is referred to as calculation of susceptibility through multiple orientation sampling, or COSMOS , and is regarded as the golden standard for QSM. This requires, however, long acquisition times and the assumption of isotropic magnetic susceptibility, thereby making this technique impractical for in vivo measurements . Therefore, regularization approaches using prior information to determine a unique solution from a single acquisition are used. Several regularization algorithms for solving the ill-posed inversion problem from a single acquisition and minimizing the streaking artifacts have been proposed .
The QSM post-processing technique has been used to provide quantitative measurements of the cerebral mixed venous oxygenation saturation (SvO2) (Table 2). The SvO2 can be calculated by exploiting the magnetic properties of hemoglobin. When fully oxygenated, diamagnetic hemoglobin results in negative susceptibility, leading to a large susceptibility decrease in arteries. In contrast, deoxygenated hemoglobin is paramagnetic. Therefore, the venous vessels will result in increased susceptibility . The ratio of blood oxygen defined as the Oxygen Extraction Fraction (OEF) can be obtained by calculating the difference between the susceptibility in the veins and tissue (Table 2). Additionally, when used in combination with other MRI sequences that provides a measure of the CBF, QSM can provide a completely non-invasive quantitative measurement of CMRO2 (Table 2) . QSM can provide important complementary information to other perfusion MRI techniques and has high potential to be used as an early biomarker for small changes in the cerebral physiology .
A major limitation of QSM is the assumption of the isotropic magnetic susceptibility within a voxel. However, this assumption does not hold for all types of brain tissue. In order to account for tissue containing anisotropic magnetic susceptibilities, tensor imaging was developed . Another important limitation of QSM is that it can only provide a relative quantification of the magnetic susceptibility . So far, this has limited the inter-subject comparison of QSM measurements. The accuracy of QSM is strongly dependent on several factors, which requires careful standardization and optimization of the post-processing and acquisition parameters to increase its reproducibility and accuracy further . To date, QSM has mainly been used in the human brain, but it might also become useful in other parts of the body.
Only a few studies have investigated QSM in rodents . Pre-clinical QSM could provide interesting insight in brain functioning under altered physiological conditions. Comparing data under normoxic, hyperoxic and hypercapnic conditions, and verification with blood gas analysis, could improve our understanding of neurological disease mechanisms and cerebral functioning. QSM especially benefits from a high magnetic field strength, as higher field strength increases both the SNR and the contrast in the phase images.
11. Quantitative Blood-Oxygenation-Level-Depedent MRI
The BOLD MRI technique is based on the different magnetic properties of oxygenated and deoxygenated blood. During brain activation, the neurons consume an increased amount of oxygen. This results in a local increase in the regional CBF towards the site of activation. When placed in an external magnetic field, the deoxyhemoglobin in the blood vessels will induce a magnetic susceptibility difference between blood and tissue. This susceptibility difference results into signal changes of the local contrast in both T2- and T2*-weighted images and can therefore be used as an endogenous contrast agent to measure brain activation indirectly . The temporal resolution of BOLD MRI is too poor to perform direct measurements of neuronal activity in the brain. However, due to its high spatial resolution and the ability to measure the much slower response in regional blood flow, BOLD fMRI has been used extensively as a tool to investigate temporal MR signal changes non-invasively during brain activity as a response to a cognitive tasks or stimuli under both normal and pathophysiologic conditions .
However, the BOLD effect can also be used to study the resting state or baseline of the brain and its impairment by neurodegenerative diseases. This allows quantitative evaluation of cerebral blood oxygenation and was therefore named quantitative BOLD (qBOLD) . The deoxyhemoglobin present in the blood vessels creates tissue specific mesoscopic field inhomogeneities . Using a Gradient Echo Sampling of Spin Echo (GESSE) sequence, the mesoscopic field inhomogeneities can be separated from microscopic and macroscopic inhomogeneities . The mesoscopic, tissue specific BOLD signal can then be fitted to a multi-compartment model based on prior knowledge about the brain tissue. This allows derivation of quantitative hemodynamic parameters, such as the deoxygenated blood volume on the venous site of the blood vessel network (DBV) and the brain tissue OEF  . The results from qBOLD are in good correspondence with direct measurements of blood oxygenation in a validation study in a rat model .
The original qBOLD method has some limitations. An additional field map needs to be acquired to correct for signal loss near air-tissue interfaces. This becomes increasingly important when imaging rodents due to the smaller brain size, the close location of sinuses to the brain and the associated decreased shim quality high magnetic field strength. Furthermore, the use of multi-compartment models, removing the R2-weighting and nulling the cerebrospinal fluid signal, requires a high SNR due to the large amount of modeling parameters. An alternative method using a FLAIR-GASE method was recently proposed . However, this method has not been validated in animal models yet. Furthermore, a dynamic approach using a parallel acceleration technique to provide a higher temporal resolution compared to the GESSE sequence, named multi-echo asymmetric spin echo (MASE), was proposed . This method can be used to measure dynamic changes in the brain oxygenation during brain activation.
Similar to QSM, the combination of qBOLD with another technique that provides measurements of the CBF allows quantitative measurement of CMRO2. The estimation of CMRO2 in QSM and qBOLD relies on a number of assumptions. Recently, a combined QSM+qBOLD (QQ) approach was proposed . The QQ approach allows the quantification of OEF and CMRO2 without the need for altered physiological conditions or empirical assumptions. This combined approach showed clearer gray and white matter compared to separate acquisitions of QSM and qBOLD, and a more uniform OEF and better agreement with independent methods to estimate CMRO2, compared to QSM. This makes QQ a highly promising method for improvement in the accuracy of blood oxygenation measurements in the brain.
The entry is from 10.3390/diagnostics11060926
- Tofts, P.S.; Kermode, A.G. Measurement of the blood-brain barrier permeability and leakage space using dynamic MR imaging. 1. Fundamental concepts. Magn. Reson. Med. 1991, 17, 357–367.
- Essig, M.; Shiroishi, M.S.; Nguyen, T.B.; Saake, M.; Provenzale, J.M.; Enterline, D.; Anzalone, N.; Dörfler, A.; Rovira, À.; Wintermark, M.; et al. Perfusion MRI: The Five Most Frequently Asked Technical Questions. Am. J. Roentgenol. 2013, 200, 24–34.
- Gordon, Y.; Partovi, S.; Müller-Eschner, M.; Amarteifio, E.; Bäuerle, T.; Weber, M.-A.; Kauczor, H.-U.; Rengier, F. Dynamic contrast-enhanced magnetic resonance imaging: Fundamentals and application to the evaluation of the peripheral perfusion. Cardiovasc. Diagn. Ther. 2014, 4, 147–164.
- Nielsen, T.; Wittenborn, T.; Horsman, M.R. Dynamic Contrast-Enhanced Magnetic Resonance Imaging (DCE-MRI) in Preclinical Studies of Antivascular Treatments. Pharmaceutics 2012, 4, 563–589.
- Khalifa, F.; Soliman, A.; El-Baz, A.; El-Ghar, M.A.; El-Diasty, T.; Gimel’Farb, G.; Ouseph, R.; Dwyer, A.C. Models and methods for analyzing DCE-MRI: A review. Med Phys. 2014, 41, 124301.
- Walker-Samuel, S.; Leach, M.O.; Collins, D.J. Evaluation of response to treatment using DCE-MRI: The relationship between initial area under the gadolinium curve (IAUGC) and quantitative pharmacokinetic analysis. Phys. Med. Biol. 2006, 51, 3593–3602.
- Calamante, F. Arterial input function in perfusion MRI: A comprehensive review. Prog. Nucl. Magn. Reson. Spectrosc. 2013, 74, 1–32.
- Tofts, P.S.; Brix, G.; Buckley, D.L.; Evelhoch, J.L.; Henderson, E.; Knopp, M.V.; Larsson, H.B.W.; Lee, T.-Y.; Mayr, N.A.; Parker, G.J.M.; et al. Estimating kinetic parameters from dynamic contrast-enhanced t1-weighted MRI of a diffusable tracer: Standardized quantities and symbols. J. Magn. Reson. Imaging 1999, 10, 223–232.
- Duan, C.; Kallehauge, J.F.; Bretthorst, G.L.; Tanderup, K.; Ackerman, J.J.H.; Garbow, J.R. Are complex DCE-MRI models supported by clinical data? Magn. Reson. Med. 2017, 77, 1329–1339.
- Barnes, S.L.; Whisenant, J.G.; Loveless, M.E.; Yankeelov, T.E. Practical Dynamic Contrast Enhanced MRI in Small Animal Models of Cancer: Data Acquisition, Data Analysis, and Interpretation. Pharmaceutics 2012, 4, 442–478.
- Fruytier, A.-C.; Magat, J.; Colliez, F.; Jordan, B.; Cron, G.; Gallez, B. Dynamic contrast-enhanced MRI in mice at high field: Estimation of the arterial input function can be achieved by phase imaging. Magn. Reson. Med. 2014, 71, 544–550.
- Tsao, J.; Kozerke, S. MRI temporal acceleration techniques. J. Magn. Reson. Imaging 2012, 36, 543–560.
- Pain, F.; Lanièce, P.; Mastrippolito, R.; Gervais, P.; Hantraye, P.; Besret, L. Arterial Input Function Measurement Without Blood Sampling Using a β-Microprobe in Rats. J. Nucl. Med. 2004, 45, 1577–1582.
- Zhou, R.; Pickup, S.; Yankeelov, T.E.; Springer, C.S., Jr.; Glickson, J.D. Simultaneous measurement of arterial input function and tumor pharmacokinetics in mice by dynamic contrast enhanced imaging: Effects of transcytolemmal water exchange. Magn. Reson. Med. 2004, 52, 248–257.
- Yankeelov, T.E.; Luci, J.J.; Lepage, M.; Li, R.; Debusk, L.; Lin, P.C.; Price, R.R.; Gore, J.C. Quantitative pharmacokinetic analysis of DCE-MRI data without an arterial input function: A reference region model. Magn. Reson. Imaging 2005, 23, 519–529.
- Van Osch, M.J.P.; Vonken, E.-J.P.A.; Viergever, M.A.; van der Grond, J.; Bakker, C.J.G. Measuring the arterial input function with gradient echo sequences. Magn. Reson. Med. 2003, 49, 1067–1076.
- McGrath, D.M.; Bradley, D.P.; Tessier, J.L.; Lacey, T.; Taylor, C.J.; Parker, G.J.M. Comparison of model-based arterial input functions for dynamic contrast-enhanced MRI in tumor bearing rats. Magn. Reson. Med. 2009, 61, 1173–1184.
- Yankeelov, T.E.; de Busk, L.M.; Billheimer, D.D.; Luci, J.J.; Lin, P.C.; Price, R.R.; Gore, J.C. Repeatability of a reference region model for analysis of murine DCE-MRI data at 7T. J. Magn. Reson. Imaging 2006, 24, 1140–1147.
- Ortuño, J.E.; Ledesma-Carbayo, M.J.; Simões, R.V.; Candiota, A.P.; Arús, C.; Santos, A. : A dynamic contrast-enhanced MRI pharmacokinetic analysis tool for preclinical data. BMC Bioinform. 2013, 14, 316.
- Sedlacik, J.; Myers, A.; Loeffler, R.B.; Williams, R.F.; Davidoff, A.M.; Hillenbrand, C.M. A dedicated automated injection system for dynamic contrast-enhanced MRI experiments in mice. J. Magn. Reson. Imaging 2013, 37, 746–751.
- Calamante, F.; Gadian, D.G.; Connelly, A. Quantification of Perfusion Using Bolus Tracking Magnetic Resonance Imaging in Stroke. Stroke 2002, 33, 1146–1151.
- Stadler, K.L.; Pease, A.P.; Ballegeer, E.A. Dynamic Susceptibility Contrast Magnetic Resonance Imaging Protocol of the Normal Canine Brain. Front. Veter. Sci. 2017, 4, 41.
- Østergaard, L. Principles of cerebral perfusion imaging by bolus tracking. J. Magn. Reson. Imaging 2005, 22, 710–717.
- Jin, S.; Kang, M.; Cho, H. Cerebral blood perfusion deficits using dynamic susceptibility contrast MRI with gadolinium chelates in rats with post-ischemic reperfusion without significant dynamic contrast-enhanced MRI-derived vessel permeabilities: A cautionary note. PLoS ONE 2018, 13, e0201076.
- Boxerman, J.L.; Prah, D.E.; Paulson, E.S.; Machan, J.T.; Bedekar, D.; Schmainda, K.M. The Role of Preload and Leakage Correction in Gadolinium-Based Cerebral Blood Volume Estimation Determined by Comparison with MION as a Criterion Standard. Am. J. Neuroradiol. 2012, 33, 1081–1087.
- Skinner, J.T.; Moots, P.L.; Ayers, G.D.; Quarles, C.C. On the Use of DSC-MRI for Measuring Vascular Permeability. Am. J. Neuroradiol. 2016, 37, 80–87.
- Leu, K.; Boxerman, J.L.; Cloughesy, T.F.; Lai, A.; Nghiemphu, P.L.; Liau, L.M.; Pope, W.B.; Ellingson, B.M. Improved Leakage Correction for Single-Echo Dynamic Susceptibility Contrast Perfusion MRI Estimates of Relative Cerebral Blood Volume in High-Grade Gliomas by Accounting for Bidirectional Contrast Agent Exchange. Am. J. Neuroradiol. 2016, 37, 1440–1446.
- Uematsu, H.; Maeda, M.; Sadato, N.; Matsuda, T.; Ishimori, Y.; Koshimoto, Y.; Kimura, H.; Yamada, H.; Kawamura, Y.; Yonekura, Y.; et al. Blood volume of gliomas determined by double-echo dynamic perfusion-weighted MR imaging: A preliminary study. Am. J. Neuroradiol. 2001, 22, 1915–1919.
- Newton, A.T.; Pruthi, S.; Stokes, A.M.; Skinner, J.T.; Quarles, C.C. Improving Perfusion Measurement in DSC-MR Imaging with Multiecho Information for Arterial Input Function Determination. Am. J. Neuroradiol. 2016, 37, 1237–1243.
- Quarles, C.C.; Gore, J.C.; Xu, L.; Yankeelov, T.E. Comparison of dual-echo DSC-MRI- and DCE-MRI-derived contrast agent kinetic parameters. J. Magn. Reson. Imaging 2012, 30, 944–953.
- Dennie, J.; Mandeville, J.B.; Boxerman, J.L.; Packard, S.D.; Rosen, B.R.; Weisskoff, R.M. NMR imaging of changes in vascular morphology due to tumor angiogenesis. Magn. Reson. Med. 1998, 40, 793–799.
- Stokes, A.M.; Skinner, J.T.; Quarles, C.C. Assessment of a combined spin- and gradient-echo (SAGE) DSC-MRI method for preclinical neuroimaging. J. Magn. Reson. Imaging 2014, 32, 1181–1190.
- Calamante, F.; Thomas, D.L.; Pell, G.S.; Wiersma, J.; Turner, R. Measuring Cerebral Blood Flow Using Magnetic Resonance Imaging Techniques. J. Cereb. Blood Flow Metab. 1999, 19, 701–735.
- Schmiedeskamp, H.; Straka, M.; Newbould, R.D.; Zaharchuk, G.; Andre, J.B.; Olivot, J.-M.; Moseley, M.E.; Albers, G.W.; Bammer, R. Combined spin- and gradient-echo perfusion-weighted imaging. Magn. Reson. Med. 2012, 68, 30–40.
- Ferré, J.-C.; Bannier, E.; Raoult, H.; Mineur, G.; Carsin-Nicol, B.; Gauvrit, J.-Y. Arterial spin labeling (ASL) perfusion: Techniques and clinical use. Diagn. Interv. Imaging 2013, 94, 1211–1223.
- Alsaedi, A.; Thomas, D.; Bisdas, S.; Golay, X. Overview and Critical Appraisal of Arterial Spin Labelling Technique in Brain Perfusion Imaging. Contrast Media Mol. Imaging 2018, 2018, 1–15.
- Petersen, E.T.; Lim, T.; Golay, X. Model-free arterial spin labeling quantification approach for perfusion MRI. Magn. Reson. Med. 2006, 55, 219–232.
- Brookes, M.J.; Morris, P.G.; Gowland, P.A.; Francis, S.T. Noninvasive measurement of arterial cerebral blood volume using look-locker EPI and arterial spin labeling. Magn. Reson. Med. 2007, 58, 41–54.
- Wang, J.; Alsop, D.C.; Song, H.K.; Maldjian, J.A.; Tang, K.; Salvucci, A.E.; Detre, J.A. Arterial transit time imaging with flow encoding arterial spin tagging (FEAST). Magn. Reson. Med. 2003, 50, 599–607.
- Kim, T.; Kim, S.-G. Quantification of cerebral arterial blood volume and cerebral blood flow using MRI with modulation of tissue and vessel (MOTIVE) signals. Magn. Reson. Med. 2005, 54, 333–342.
- Williams, D.S.; Detre, J.A.; Leigh, J.S.; Koretsky, A.P. Magnetic resonance imaging of perfusion using spin inversion of arterial water. Proc. Natl. Acad. Sci. USA 1992, 89, 212–216.
- Edelman, R.R.; Siewert, B.; Darby, D.G.; Thangaraj, V.; Nobre, A.C.; Mesulam, M.M.; Warach, S. Qualitative mapping of cerebral blood flow and functional localization with echo-planar MR imaging and signal targeting with alternating radio frequency. Radiology 1994, 192, 513–520.
- Silva, A.C.; Kim, S.G. Pseudo-continuous arterial spin labeling technique for measuring CBF dynamics with high temporal resolution. Magn. Reson. Med. 1999, 42, 425–429.
- Dai, W.; Garcia, D.; de Bazelaire, C.; Alsop, D.C. Continuous flow-driven inversion for arterial spin labeling using pulsed radio frequency and gradient fields. Magn. Reson. Med. 2008, 60, 1488–1497.
- Silva, A.C.; Zhang, W.; Williams, D.S.; Koretsky, A.P. Multi-Slice MRI of Rat Brain Perfusion During Amphetamine Stimulation Using Arterial Spin Labeling. Magn. Reson. Med. 1995, 33, 209–214.
- Alsop, D.C.; Detre, J.A. Multisection cerebral blood flow MR imaging with continuous arterial spin labeling. Radiology 1998, 208, 410–416.
- Zhang, W.; Silva, A.C.; Williams, D.S.; Koretsky, A.P. NMR Measurement of Perfusion Using Arterial Spin Labeling Without Saturation of Macromolecular Spins. Magn. Reson. Med. 1995, 33, 370–376.
- Muir, E.R.; Shen, Q.; Duong, T.Q. Cerebral blood flow MRI in mice using the cardiac-spin-labeling technique. Magn. Reson. Med. 2008, 60, 744–748.
- Kim, S.-G. Quantification of relative cerebral blood flow change by flow-sensitive alternating inversion recovery (FAIR) technique: Application to functional mapping. Magn. Reson. Med. 1995, 34, 293–301.
- Wong, E.C.; Buxton, R.B.; Frank, L.R. Implementation of quantitative perfusion imaging techniques for functional brain mapping using pulsed arterial spin labeling. NMR Biomed. 1997, 10, 237–249.
- Chen, Y.; Wang, D.J.J.; Detre, J.A. Test-retest reliability of arterial spin labeling with common labeling strategies. J. Magn. Reson. Imaging 2011, 33, 940–949.
- Jahanian, H.; Noll, D.C.; Hernandez-Garcia, L. B0field inhomogeneity considerations in pseudo-continuous arterial spin labeling (pCASL): Effects on tagging efficiency and correction strategy. NMR Biomed. 2011, 24, 1202–1209.
- Alsop, D.C.; Detre, J.A.; Golay, X.; Günther, M.; Hendrikse, J.; Hernandez-Garcia, L.; Lu, H.; MacIntosh, B.J.; Parkes, L.M.; Smits, M.; et al. Recommended implementation of arterial spin-labeled perfusion MRI for clinical applications: A consensus of the ISMRM perfusion study group and the European consortium for ASL in dementia. Magn. Reson. Med. 2015, 73, 102–116.
- Hirschler, L.; Collomb, N.; Voiron, J.; Köhler, S.; Barbier, E.L.; Warnking, J.M. SAR comparison between CASL and pCASL at high magnetic field and evaluation of the benefit of a dedicated labeling coil. Magn. Reson. Med. 2020, 83, 254–261.
- Dobre, M.C.; Uğurbil, K.; Marjanska, M. Determination of blood longitudinal relaxation time (T1) at high magnetic field strengths. J. Magn. Reson. Imaging 2007, 25, 733–735.
- Franke, C.; van Dorsten, F.A.; Oláh, L.; Schwindt, W.; Hoehn, M. Arterial spin tagging perfusion imaging of rat brain: Dependency on magnetic field strength. J. Magn. Reson. Imaging 2000, 18, 1109–1113.
- Golay, X.; Petersen, E.T. Arterial Spin Labeling: Benefits and Pitfalls of High Magnetic Field. Neuroimaging Clin. North Am. 2006, 16, 259–268.
- St. Lawrence, K.S.; Wang, J. Effects of the apparent transverse relaxation time on cerebral blood flow measurements obtained by arterial spin labeling. Magn. Reson. Med. 2005, 53, 425–433.
- Maleki, N.; Dai, W.; Alsop, D.C. Optimization of background suppression for arterial spin labeling perfusion imaging. Magn. Reson. Mater. Physics Biol. Med. 2012, 25, 127–133.
- Shen, Q.; Duong, T.Q. Background suppression in arterial spin labeling MRI with a separate neck labeling coil. NMR Biomed. 2011, 24, 1111–1118.
- Damen, F.C.; Tain, R.-W.; Thomas, R.; Li, W.; Tai, L.; Cai, K. Evaluation of B0-correction of relative CBF maps using tagging distance dependent Z-spectrum (TADDZ). J. Magn. Reson. Imaging 2020, 65, 83–89.
- Berry, E.S.K.; Jezzard, P.; Okell, T.W. Off-resonance correction for pseudo-continuous arterial spin labeling using the optimized encoding scheme. NeuroImage 2019, 199, 304–312.
- Jahng, G.-H.; Weiner, M.W.; Schuff, N. Improved arterial spin labeling method: Applications for measurements of cerebral blood flow in human brain at high magnetic field MRI. Med Phys. 2007, 34, 4519–4525.
- Tanenbaum, A.B.; Snyder, A.Z.; Brier, M.R.; Ances, B.M. A Method for Reducing the Effects of Motion Contamination in Arterial Spin Labeling Magnetic Resonance Imaging. J. Cereb. Blood Flow Metab. 2015, 35, 1697–1702.
- Hébert, F.; Grand’Maison, M.; Ho, M.-K.; Lerch, J.P.; Hamel, E.; Bedell, B.J. Cortical atrophy and hypoperfusion in a transgenic mouse model of Alzheimer’s disease. Neurobiol. Aging 2013, 34, 1644–1652.
- Govaerts, K.; Lechat, B.; Struys, T.; Kremer, A.; Borghgraef, P.; van Leuven, F.; Himmelreich, U.; Dresselaers, T. Longitudinal assessment of cerebral perfusion and vascular response to hypoventilation in a bigenic mouse model of Alzheimer’s disease with amyloid and tau pathology. NMR Biomed. 2019, 32, e4037.
- Struys, T.; Govaerts, K.; Oosterlinck, W.; Casteels, C.; Bronckaers, A.; Koole, M.; van Laere, K.; Herijgers, P.; Lambrichts, I.; Himmelreich, U.; et al. In vivo evidence for long-term vascular remodeling resulting from chronic cerebral hypoperfusion in mice. J. Cereb. Blood Flow Metab. 2016, 37, 726–739.
- Holmes, H.E.; Colgan, N.; Ismail, O.; Ma, D.; Powell, N.M.; O’Callaghan, J.M.; Harrison, I.F.; Johnson, R.A.; Murray, T.K.; Ahmed, Z.; et al. Imaging the accumulation and suppression of tau pathology using multiparametric MRI. Neurobiol. Aging 2016, 39, 184–194.
- Oosterlinck, W.W.; Dresselaers, T.; Geldhof, V.; Van Santvoort, A.; Robberecht, W.; Herijgers, P.; Himmelreich, U. Response of mouse brain perfusion to hypo- and hyperventilation measured by arterial spin labeling. Magn. Reson. Med. 2011, 66, 802–811.
- Wong, E.C.; Cronin, M.; Wu, W.-C.; Inglis, B.; Frank, L.R.; Liu, T.T. Velocity-selective arterial spin labeling. Magn. Reson. Med. 2006, 55, 1334–1341.
- Le Bihan, D.; Breton, E.; Lallemand, D.; Aubin, M.L.; Vignaud, J.; Laval-Jeantet, M. Separation of diffusion and perfusion in intravoxel incoherent motion MR imaging. Radiology 1988, 168, 497–505.
- Le Bihan, D. What can we see with IVIM MRI? NeuroImage 2019, 187, 56–67.
- Jensen, J.H.; Helpern, J.A.; Ramani, A.; Lu, H.; Kaczynski, K. Diffusional kurtosis imaging: The quantification of non-gaussian water diffusion by means of magnetic resonance imaging. Magn. Reson. Med. 2005, 53, 1432–1440.
- Jensen, J.H.; Helpern, J.A. MRI quantification of non-Gaussian water diffusion by kurtosis analysis. NMR Biomed. 2010, 23, 698–710.
- Mulkern, R.V.; Haker, S.J.; Maier, S.E. On high b diffusion imaging in the human brain: Ruminations and experimental insights. J. Magn. Reson. Imaging 2009, 27, 1151–1162.
- Fournet, G.; Li, J.-R.; Cerjanic, A.M.; Sutton, B.P.; Ciobanu, L.; le Bihan, D. A two-pool model to describe the IVIM cerebral perfusion. J. Cereb. Blood Flow Metab. 2016, 37, 2987–3000.
- Kennan, R.P.; Gao, J.-H.; Zhong, J.; Gore, J.C. A general model of microcirculatory blood flow effects in gradient sensitized MRI. Med Phys. 1994, 21, 539–545.
- Wetscherek, A.; Stieltjes, B.; Laun, F.B. Flow-compensated intravoxel incoherent motion diffusion imaging. Magn. Reson. Med. 2014, 74, 410–419.
- Neil, J.J.; Scherrer, L.A.; Ackerman, J.J.H. An approach to solving the dynamic range problem in measurement of the pseudodilusion coefficient in vivo with spin echoes. J. Magn. Reson. 1991, 95, 607–614.
- Henkelman, R.M.; Neil, J.J.; Xiang, Q.-S. A quantitative interpretation of IVIM measurements of vascular perfusion in the rat brain. Magn. Reson. Med. 1994, 32, 464–469.
- Duong, T.Q.; Kim, S.-G. In vivo MR measurements of regional arterial and venous blood volume fractions in intact rat brain. Magn. Reson. Med. 2000, 43, 393–402.
- Wurnig, M.C.; Donati, O.F.; Ulbrich, E.; Filli, L.; Kenkel, D.; Thoeny, H.C.; Boss, A. Systematic analysis of the intravoxel incoherent motion threshold separating perfusion and diffusion effects: Proposal of a standardized algorithm. Magn. Reson. Med. 2014, 74, 1414–1422.
- While, P.T. A comparative simulation study of bayesian fitting approaches to intravoxel incoherent motion modeling in diffusion-weighted MRI. Magn. Reson. Med. 2017, 78, 2373–2387.
- Li, Y.T.; Cercueil, J.-P.; Yuan, J.; Chen, W.; Loffroy, R.; Wáng, Y.X.J. Liver intravoxel incoherent motion (IVIM) magnetic resonance imaging: A comprehensive review of published data on normal values and applications for fibrosis and tumor evaluation. Quant. Imaging Med. Surg. 2017, 7, 59–78.
- Ichikawa, S.; Motosugi, U.; Ichikawa, T.; Sano, K.; Morisaka, H.; Araki, T. Intravoxel incoherent motion imaging of the kidney: Alterations in diffusion and perfusion in patients with renal dysfunction. J. Magn. Reson. Imaging 2013, 31, 414–417.
- Schneider, M.J.; Dietrich, O.; Ingrisch, M.; Helck, A.; Winter, K.S.; Reiser, M.F.; Staehler, M.; Casuscelli, J.; Notohamiprodjo, M. Intravoxel Incoherent Motion Magnetic Resonance Imaging in Partially Nephrectomized Kidneys. Investig. Radiol. 2016, 51, 323–330.
- Klauß, M.; Maier-Hein, K.; Tjaden, C.; Hackert, T.; Grenacher, L.; Stieltjes, B. IVIM DW-MRI of autoimmune pancreatitis: Therapy monitoring and differentiation from pancreatic cancer. Eur. Radiol. 2016, 26, 2099–2106.
- Kang, K.M.; Lee, J.M.; Yoon, J.H.; Kiefer, B.; Han, J.K.; Choi, B.I. Intravoxel Incoherent Motion Diffusion-weighted MR Imaging for Characterization of Focal Pancreatic Lesions. Radiology 2014, 270, 444–453.
- Zhang, X.; Ingo, C.; Teeuwisse, W.M.; Chen, Z.; van Osch, M.J.P. Comparison of perfusion signal acquired by arterial spin labeling-prepared intravoxel incoherent motion (IVIM) MRI and conventional IVIM MRI to unravel the origin of the IVIM signal. Magn. Reson. Med. 2018, 79, 723–729.
- Paschoal, A.M.; Leoni, R.F.; dos Santos, A.C.; Paiva, F.F. Intravoxel incoherent motion MRI in neurological and cerebrovascular diseases. NeuroImage Clin. 2018, 20, 705–714.
- Jalnefjord, O.; Montelius, M.; Starck, G.; Ljungberg, M. Optimization of b-value schemes for estimation of the diffusion coefficient and the perfusion fraction with segmented intravoxel incoherent motion model fitting. Magn. Reson. Med. 2019, 82, 1541–1552.
- Perucho, J.A.U.; Chang, H.C.C.; Vardhanabhuti, V.; Wang, M.; Becker, A.S.; Wurnig, M.C.; Lee, E.Y.P. B-Value Optimization in the Estimation of Intravoxel Incoherent Motion Parameters in Patients with Cervical Cancer. Korean J. Radiol. 2020, 21, 218–227.
- Lu, H.; Golay, X.; Pekar, J.J.; van Zijl, P.C. Functional magnetic resonance imaging based on changes in vascular space occupancy. Magn. Reson. Med. 2003, 50, 263–274.
- Lu, H.; Hua, J.; van Zijl, P.C.M. Noninvasive functional imaging of cerebral blood volume with vascular-space-occupancy (VASO) MRI. NMR Biomed. 2013, 26, 932–948.
- Hua, J.; Jones, C.K.; Qin, Q.; van Zijl, P.C.M. Implementation of vascular-space-occupancy MRI at 7T. Magn. Reson. Med. 2013, 69, 1003–1013.
- Hua, J.; Donahue, M.J.; Zhao, J.M.; Grgac, K.; Huang, A.J.; Zhou, J.; van Zijl, P.C.M. Magnetization transfer enhanced vascular-space-occupancy (MT-VASO) functional MRI. Magn. Reson. Med. 2009, 61, 944–951.
- Huber, L.; Ivanov, D.; Krieger, S.N.; Streicher, M.N.; Mildner, T.; Poser, B.A.; Möller, H.E.; Turner, R. Slab-selective, BOLD-corrected VASO at 7 Tesla provides measures of cerebral blood volume reactivity with high signal-to-noise ratio. Magn. Reson. Med. 2014, 72, 137–148.
- Jin, T.; Kim, S.-G. Improved cortical-layer specificity of vascular space occupancy fMRI with slab inversion relative to spin-echo BOLD at 9.4 T. NeuroImage 2008, 40, 59–67.
- Lu, H.; Law, M.; Johnson, G.; Ge, Y.; van Zijl, P.C.M.; Helpern, J.A. Novel approach to the measurement of absolute cerebral blood volume using vascular-space-occupancy magnetic resonance imaging. Magn. Reson. Med. 2005, 54, 1403–1411.
- Hua, J.; Qin, Q.; Donahue, M.J.; Zhou, J.; Pekar, J.J.; van Zijl, P.C.M. Inflow-based vascular-space-occupancy (iVASO) MRI. Magn. Reson. Med. 2011, 66, 40–56.
- Donahue, M.J.; Strother, M.K.; Hendrikse, J. Novel MRI Approaches for Assessing Cerebral Hemodynamics in Ischemic Cerebrovascular Disease. Stroke 2012, 43, 903–915.
- Donahue, M.J.; Sideso, E.; MacIntosh, B.J.; Kennedy, J.; Handa, A.; Jezzard, P. Absolute Arterial Cerebral Blood Volume Quantification Using Inflow Vascular-Space-Occupancy with Dynamic Subtraction Magnetic Resonance Imaging. J. Cereb. Blood Flow Metab. 2010, 30, 1329–1342.
- Hua, J.; Qin, Q.; Pekar, J.J.; van Zijl, P.C.M. Measurement of absolute arterial cerebral blood volume in human brain without using a contrast agent. NMR Biomed. 2011, 24, 1313–1325.
- Kim, S.-G.; Harel, N.; Jin, T.; Kim, T.; Lee, P.; Zhao, F. Cerebral blood volume MRI with intravascular superparamagnetic iron oxide nanoparticles. NMR Biomed. 2012, 26, 949–962.
- Mandeville, J.B.; Marota, J.J.A. Vascular filters of functional MRI: Spatial localization using BOLD and CBV contrast. Magn. Reson. Med. 1999, 42, 591–598.
- Le Duc, G.; Péoc’h, M.; Rémy, C.; Charpy, O.; Muller, R.N.; Le Bas, J.F.; Décorps, M. Use of T2-weighted susceptibility contrast MRI for mapping the blood volume in the glioma-bearing rat brain. Magn. Reson. Med. 1999, 42, 754–761.
- Weinstein, J.S.; Varallyay, C.G.; Dosa, E.; Gahramanov, S.; Hamilton, B.; Rooney, W.D.; Muldoon, L.L.; Neuwelt, E.A. Superparamagnetic Iron Oxide Nanoparticles: Diagnostic Magnetic Resonance Imaging and Potential Therapeutic Applications in Neurooncology and Central Nervous System Inflammatory Pathologies, a Review. J. Cereb. Blood Flow Metab. 2010, 30, 15–35.
- Troprès, I.; Pannetier, N.; Grand, S.; Lemasson, B.; Moisan, A.; Peoc’h, M.; Rémy, C.; Barbier, E.L. Imaging the microvessel caliber and density: Principles and applications of microvascular MRI. J. Cereb. Blood Flow Metab. 2015, 73, 325–341.
- Boxerman, J.L.; Bandettini, P.A.; Kwong, K.K.; Baker, J.R.; Davis, T.L.; Rosen, B.R.; Weisskoff, R.M. The intravascular contribution to fMRI signal change: Monte Carlo modeling and diffusion-weighted studies in vivo. Magn. Reson. Med. 1995, 34, 4–10.
- Jensen, J.H.; Chandra, R. MR imaging of microvasculature. Magn. Reson. Med. 2000, 44, 224–230.
- Troprès, I.; Grimault, S.; Vaeth, A.; Grillon, E.; Julien, C.; Payen, J.-F.; Lamalle, L.; Décorps, M. Vessel size imaging. Magn. Reson. Med. 2001, 45, 397–408.
- Kiselev, V.G.; Strecker, R.; Ziyeh, S.; Speck, O.; Hennig, J. Vessel size imaging in humans. Magn. Reson. Med. 2005, 53, 553–563.
- Hsu, Y.-Y.; Yang, W.-S.; Lim, K.-E.; Liu, H.-L. Vessel size imaging using dual contrast agent injections. J. Magn. Reson. Imaging 2009, 30, 1078–1084.
- Schmiedeskamp, H.; Straka, M.; Bammer, R. Compensation of slice profile mismatch in combined spin- and gradient-echo echo-planar imaging pulse sequences. Magn. Reson. Med. 2012, 67, 378–388.
- Skinner, J.T.; Robison, R.K.; Elder, C.P.; Newton, A.T.; Damon, B.M.; Quarles, C.C. Evaluation of a multiple spin- and gradient-echo (SAGE) EPI acquisition with SENSE acceleration: Applications for perfusion imaging in and outside the brain.J. Magn. Reson. Imaging 2014, 32, 1171–1180.
- Van Vaals, J.J.; Brummer, M.E.; Dixon, W.T.; Tuithof, H.H.; Engels, H.; Nelson, R.C.; Gerety, B.M.; Chezmar, J.L.; Boer, J.A.D. “Keyhole” method for accelerating imaging of contrast agent uptake. J. Magn. Reson. Imaging 1993, 3, 671–675.
- Han, M.; Yang, B.; Fernandez, B.; Lafontaine, M.; Alcaide-Leon, P.; Jakary, A.; Burns, B.L.; Morrison, M.A.; Villanueva-Meyer, J.E.; Chang, S.M.; et al. Simultaneous multi-slice spin- and gradient-echo dynamic susceptibility-contrast perfusion-weighted MRI of gliomas. NMR Biomed. 2020, 34, e4399.
- Wymer, D.T.; Patel, K.P.; Burke, W.F., III; Bhatia, V.K. Phase-Contrast MRI: Physics, Techniques, and Clinical Applications. Radiographics 2020, 40, 122–140.
- Korbecki, A.; Zimny, A.; Podgórski, P.; Sąsiadek, M.; Bladowska, J. Imaging of cerebrospinal fluid flow: Fundamentals, techniques, and clinical applications of phase-contrast magnetic resonance imaging. Pol. J. Radiol. 2019, 84, e240–e250.
- Srichai, M.B.; Lim, R.P.; Wong, S.; Lee, V.S. Cardiovascular Applications of Phase-Contrast MRI. Am. J. Roentgenol. 2009, 192, 662–675.
- Nett, E.J.; Johnson, K.M.; Frydrychowicz, A.; del Rio, A.M.; Schrauben, E.; Francois, C.J.; Wieben, O. Four-dimensional phase contrast MRI with accelerated dual velocity encoding. J. Magn. Reson. Imaging 2012, 35, 1462–1471.
- Battal, B.; Kocaoglu, M.; Bulakbasi, N.; Husmen, G.; Sanal, H.T.; Tayfun, C. Cerebrospinal fluid flow imaging by using phase-contrast MR technique. Br. J. Radiol. 2011, 84, 758–765.
- Guo, G. The Quantification of Cerebral Blood Flow by Phase Contrast MRA: Basics and Applications. Neuroradiol. J. 2008, 21, 11–21.
- Saloner, D. The AAPM/RSNA physics tutorial for residents. An introduction to MR angiography. Radiographics 1995, 15, 453–465.
- Frydrychowicz, A.; François, C.J.; Turski, P.A. Four-dimensional phase contrast magnetic resonance angiography: Potential clinical applications. Eur. J. Radiol. 2011, 80, 24–35.
- Markl, M.; Frydrychowicz, A.; Kozerke, S.; Hope, M.; Wieben, O. 4D flow MRI. J. Magn. Reson. Imaging 2012, 36, 1015–1036.
- Spilt, A.; Box, F.M.A.; van der Geest, R.J.; Reiber, J.H.C.; Kunz, P.; Kamper, A.M.; Blauw, G.J.; van Buchem, M.A. Reproducibility of total cerebral blood flow measurements using phase contrast magnetic resonance imaging. J. Magn. Reson. Imaging 2002, 16, 1–5.
- Wei, Z.; Chen, L.; Lin, Z.; Jiang, D.; Xu, J.; Liu, P.; van Zijl, P.C.M.; Lu, H. Optimization of phase-contrast MRI for the estimation of global cerebral blood flow of mice at 11.7T. Magn. Reson. Med. 2019, 81, 2566–2575.
- Chiu, S.-C.; Hsu, S.-T.; Huang, C.-W.; Shen, W.-C.; Peng, S.-L. Phase Contrast Magnetic Resonance Imaging in the Rat Common Carotid Artery. J. Vis. Exp. 2018, e57304.
- Peng, S.-L.; Su, P.; Wang, F.-N.; Cao, Y.; Zhang, R.; Lu, H.; Liu, P. Optimization of phase-contrast MRI for the quantification of whole-brain cerebral blood flow. J. Magn. Reson. Imaging 2015, 42, 1126–1133.
- Enzmann, D.R.; Marks, M.P.; Pelc, N.J. Comparison of cerebral artery blood flow measurements with gated cine and ungated phase-contrast techniques. J. Magn. Reson. Imaging 1993, 3, 705–712.
- Hofman, M.B.M.; Kouwenhoven, M.; Sprenger, M.; van Rossum, A.C.; Valk, J.; Westerhof, N. Nontriggered magnetic resonance velocity measurement of the time-average of pulsatile velocity. Magn. Reson. Med. 1993, 29, 648–655.
- Bakker, C.J.G.; Hartkamp, M.J.; Mali, W.P.T.M. Measuring blood flow by nontriggered 2D phase-contrast MR angiography. J. Magn. Reson. Imaging 1996, 14, 609–614.
- Bonekamp, D.; Degaonkar, M.; Barker, P.B. Quantitative cerebral blood flow in dynamic susceptibility contrast MRI using total cerebral flow from phase contrast magnetic resonance angiography. Magn. Reson. Med. 2011, 66, 57–66.
- Aslan, S.; Xu, F.; Wang, P.L.; Uh, J.; Yezhuvath, U.S.; van Osch, M.; Lu, H. Estimation of labeling efficiency in pseudocontinuous arterial spin labeling. Magn. Reson. Med. 2010, 63, 765–771.
- Liu, P.; Lu, H.; Filbey, F.M.; Pinkham, A.E.; McAdams, C.J.; Adinoff, B.; Daliparthi, V.; Cao, Y. Automatic and Reproducible Positioning of Phase-Contrast MRI for the Quantification of Global Cerebral Blood Flow. PLoS ONE 2014, 9, e95721.
- Liu, P.; Huang, H.; Rollins, N.; Chalak, L.F.; Jeon, T.; Halovanic, C.; Lu, H. Quantitative assessment of global cerebral metabolic rate of oxygen (CMRO2) in neonates using MRI. NMR Biomed. 2014, 27, 332–340.
- Xu, F.; Ge, Y.; Lu, H. Noninvasive quantification of whole-brain cerebral metabolic rate of oxygen (CMRO2) by MRI. Magn. Reson. Med. 2009, 62, 141–148.
- Haacke, E.M.; Mittal, S.; Wu, Z.; Neelavalli, J.; Cheng, Y.-C.N. Susceptibility-Weighted Imaging: Technical Aspects and Clinical Applications, Part 1. Am. J. Neuroradiol. 2009, 30, 19–30.
- Liu, S.; Buch, S.; Chen, Y.; Choi, H.-S.; Dai, Y.; Habib, C.; Hu, J.; Jung, J.-Y.; Luo, Y.; Utriainen, D.; et al. Susceptibility-weighted imaging: Current status and future directions. NMR Biomed. 2017, 30, e3552.
- Nissi, M.J.; Tóth, F.; Wang, L.; Carlson, C.S.; Ellermann, J.M. Improved Visualization of Cartilage Canals Using Quantitative Susceptibility Mapping. PLoS ONE 2015, 10, e0132167.
- Jenkinson, M. Fast, automated, N-dimensional phase-unwrapping algorithm. Magn. Reson. Med. 2002, 49, 193–197.
- Li, W.; Wu, B.; Liu, C. Quantitative susceptibility mapping of human brain reflects spatial variation in tissue composition. NeuroImage 2011, 55, 1645–1656.
- Liu, T.; Khalidov, I.; de Rochefort, L.; Spincemaille, P.; Liu, J.; Tsiouris, A.J.; Wang, Y. A novel background field removal method for MRI using projection onto dipole fields (PDF). NMR Biomed. 2011, 24, 1129–1136.
- Schweser, F.; Deistung, A.; Lehr, B.W.; Reichenbach, J.R. Quantitative imaging of intrinsic magnetic tissue properties using MRI signal phase: An approach to in vivo brain iron metabolism? NeuroImage 2011, 54, 2789–2807.
- Li, W.; Avram, A.V.; Wu, B.; Xiao, X.; Liu, C. Integrated Laplacian-based phase unwrapping and background phase removal for quantitative susceptibility mapping. NMR Biomed. 2014, 27, 219–227.
- Liu, C.; Li, W.; Tong, K.A.; Yeom, K.W.; Kuzminski, S. Susceptibility-weighted imaging and quantitative susceptibility mapping in the brain. J. Magn. Reson. Imaging 2015, 42, 23–41.
- Liu, C.; Wei, H.; Gong, N.-J.; Cronin, M.; Dibb, R.; Decker, K. Quantitative Susceptibility Mapping: Contrast Mechanisms and Clinical Applications. Tomography 2015, 1, 3–17.
- Bollmann, S.; Rasmussen, K.G.B.; Kristensen, M.; Blendal, R.G.; Østergaard, L.R.; Plocharski, M.; O’Brien, K.; Langkammer, C.; Janke, A.; Barth, M. DeepQSM—Using deep learning to solve the dipole inversion for quantitative susceptibility mapping. NeuroImage 2019, 195, 373–383.
- Wang, Y.; Liu, T. Quantitative susceptibility mapping (QSM): Decoding MRI data for a tissue magnetic biomarker. Magn. Reson. Med. 2015, 73, 82–101.
- Liu, T.; Spincemaille, P.; de Rochefort, L.; Kressler, B.; Wang, Y. Calculation of susceptibility through multiple orientation sampling (COSMOS): A method for conditioning the inverse problem from measured magnetic field map to susceptibility source image in MRI. Magn. Reson. Med. 2008, 61, 196–204.
- Vinayagamani, S.; Sheelakumari, R.; Sabarish, S.; Senthilvelan, S.; Ros, R.; Thomas, B.; Kesavadas, C. Quantitative Susceptibility Mapping: Technical Considerations and Clinical Applications in Neuroimaging. J. Magn. Reson. Imaging 2021, 53, 23–37.
- Li, W.; Wang, N.; Yu, F.; Han, H.; Cao, W.; Romero, R.; Tantiwongkosi, B.; Duong, T.Q.; Liu, C. A method for estimating and removing streaking artifacts in quantitative susceptibility mapping. NeuroImage 2015, 108, 111–122.
- Chatnuntawech, I.; McDaniel, P.; Cauley, S.F.; Gagoski, B.A.; Langkammer, C.; Martin, A.; Grant, P.E.; Wald, L.L.; Setsompop, K.; Adalsteinsson, E.; et al. Single-step quantitative susceptibility mapping with variational penalties. NMR Biomed. 2017, 30, e3570.
- Liu, J.; Liu, T.; de Rochefort, L.; Ledoux, J.; Khalidov, I.; Chen, W.; Tsiouris, A.J.; Wisnieff, C.; Spincemaille, P.; Prince, M.R.; et al. Morphology enabled dipole inversion for quantitative susceptibility mapping using structural consistency between the magnitude image and the susceptibility map. NeuroImage 2012, 59, 2560–2568.
- Wharton, S.; Schäfer, A.; Bowtell, R. Susceptibility mapping in the human brain using threshold-based k-space division. Magn. Reson. Med. 2010, 63, 1292–1304.
- Birkl, C.; Langkammer, C.; Sati, P.; Enzinger, C.; Fazekas, F.; Ropele, S. Quantitative Susceptibility Mapping to Assess Cerebral Vascular Compliance. Am. J. Neuroradiol. 2019, 40, 460–463.
- Liu, C. Susceptibility tensor imaging. Magn. Reson. Med. 2010, 63, 1471–1477.
- Li, W.; Liu, C.; Duong, T.Q.; van Zijl, P.C.M.; Li, X. Susceptibility tensor imaging (STI) of the brain. NMR Biomed. 2017, 30, e3540.
- Deistung, A.; Schweser, F.; Reichenbach, J.R. Overview of quantitative susceptibility mapping. NMR Biomed. 2017, 30, e3569.
- Reichenbach, J.R.; Schweser, F.; Serres, B.; Deistung, A. Quantitative Susceptibility Mapping: Concepts and Applications. Clin. Neuroradiol. 2015, 25, 225–230.
- Hsieh, M.-C.; Kuo, L.-W.; Huang, Y.-A.; Chen, J.-H. Investigating hyperoxic effects in the rat brain using quantitative susceptibility mapping based on MRI phase. Magn. Reson. Med. 2017, 77, 592–602.
- Hsieh, M.-C.; Tsai, C.-Y.; Liao, M.-C.; Yang, J.-L.; Su, C.-H.; Chen, J.-H. Quantitative Susceptibility Mapping-Based Microscopy of Magnetic Resonance Venography (QSM-mMRV) for In Vivo Morphologically and Functionally Assessing Cerebromicrovasculature in Rat Stroke Model. PLoS ONE 2016, 11, e0149602.
- Wei, H.; Xie, L.; Dibb, R.; Li, W.; Decker, K.; Zhang, Y.; Johnson, G.A.; Liu, C. Imaging whole-brain cytoarchitecture of mouse with MRI-based quantitative susceptibility mapping. NeuroImage 2016, 137, 107–115.
- Ogawa, S.; Lee, T.M.; Kay, A.R.; Tank, D.W. Brain magnetic resonance imaging with contrast dependent on blood oxygenation. Proc. Natl. Acad. Sci. USA 1990, 87, 9868–9872.
- Jonckers, E.; Shah, D.; Hamaide, J.; Verhoye, M.; van der Linden, A. The power of using functional fMRI on small rodents to study brain pharmacology and disease. Front. Pharmacol. 2015, 6, 231.
- Pan, W.-J.; Billings, J.C.W.; Grooms, J.K.; Shakil, S.; Keilholz, S.D. Considerations for resting state functional MRI and functional connectivity studies in rodents. Front. Neurosci. 2015, 9, 269.
- He, X.; Yablonskiy, D.A. Quantitative BOLD: Mapping of human cerebral deoxygenated blood volume and oxygen extraction fraction: Default state. Magn. Reson. Med. 2006, 57, 115–126.
- Yablonskiy, D.A.; Haacke, E.M. Theory of NMR signal behavior in magnetically inhomogeneous tissues: The static dephasing regime. Magn. Reson. Med. 1994, 32, 749–763.
- Yablonskiy, D.A. Quantitation of intrinsic magnetic susceptibility-related effects in a tissue matrix. Phantom study. Magn. Reson. Med. 1998, 39, 417–428.
- He, X.; Zhu, M.; Yablonskiy, D.A. Validation of oxygen extraction fraction measurement by qBOLD technique. Magn. Reson. Med. 2008, 60, 882–888.
- Stone, A.J.; Blockley, N.P. A streamlined acquisition for mapping baseline brain oxygenation using quantitative BOLD. NeuroImage 2017, 147, 79–88.
- Yin, Y.; Zhang, Y.; Gao, J.-H. Dynamic measurement of oxygen extraction fraction using a multiecho asymmetric spin echo (MASE) pulse sequence. Magn. Reson. Med. 2018, 80, 1118–1124.
- Cho, J.; Kee, Y.; Spincemaille, P.; Nguyen, T.D.; Zhang, J.; Gupta, A.; Zhang, S.; Wang, Y. Cerebral metabolic rate of oxygen (CMRO2) mapping by combining quantitative susceptibility mapping (QSM) and quantitative blood oxygenation level-dependent imaging (qBOLD). Magn. Reson. Med. 2018, 80, 1595–1604.