Robust MR assessment of cerebral blood volume

NeuroImage 112 (2015) 382–389
Contents lists available at ScienceDirect
NeuroImage
journal homepage: www.elsevier.com/locate/ynimg
Technical Note
Robust MR assessment of cerebral blood volume and mean vessel size
using SPION-enhanced ultrashort echo acquisition
S.H. Han a, J.H. Cho a, H.S. Jung a, J.Y. Suh b, J.K. Kim c, Y.R. Kim d, G. Cho b,⁎, H. Cho a,⁎
a
Ulsan National Institute of Science and Technology, Ulsan, South Korea
Korea Basic Science Institute, Ochang, South Korea
Asan Medical Center, Seoul, South Korea
d
Martinos Center for Biomedical Imaging, Charlestown, MA, USA
b
c
a r t i c l e
i n f o
Article history:
Received 8 December 2014
Accepted 16 March 2015
Available online 24 March 2015
Keywords:
Cerebral blood volume
Mean Vessel Size Index
SPION
UTE
a b s t r a c t
Intravascular superparamagnetic iron oxide nanoparticles (SPION)-enhanced MR transverse relaxation rates
(ΔR2⁎ and ΔR2) are widely used to investigate in vivo vascular parameters, such as the cerebral blood volume
(CBV), microvascular volume (MVV), and mean vessel size index (mVSI, ΔR2⁎/ΔR2). Although highly efficient, regional comparison of vascular parameters acquired using gradient-echo based ΔR2⁎ is hampered by its high sensitivity to magnetic field perturbations arising from air–tissue interfaces and large vessels. To minimize such
demerits, we took advantage of the dual contrast property of SPION and both theoretically and experimentally
verified the direct benefit of replacing gradient-echo based ΔR2⁎ measurement with ultra-short echo time
(UTE)-based ΔR1 contrast to generate the robust CBV and mVSI maps. The UTE acquisition minimized the local
measurement errors from susceptibility perturbations and enabled dose-independent CBV measurement using
the vessel/tissue ΔR1 ratio, while independent spin-echo acquisition enabled simultaneous ΔR2 measurement
and mVSI calculation of the cortex, cerebellum, and olfactory bulb, which are animal brain regions typified by significant susceptibility-associated measurement errors.
© 2015 Elsevier Inc. All rights reserved.
Introduction
Alternations in microvascular parameters, such as cerebral blood
volume (CBV) and microvascular size, are often caused by cerebrovascular disease (Gibbs et al., 1984; Hamberg et al., 1993; Kucharczyk et al.,
1993; Warach et al., 1992), vascular malfunction (Spetzler et al., 1977),
and neovascularization within tumors (Aronen et al., 1994; Aronen
et al., 1995). Detection of these microvascular variations plays an important role in assessing brain diseases and monitoring the therapeutic response in both preclinical and clinical settings (Dennie et al., 1998).
Many non-invasive estimation methods of microvascular anatomy
have been investigated, such as iodinated contrast radiography (Grubb
et al., 1973; Risberg et al., 1969), positron emission tomography studies
using oxygen-15-labeled carbon monoxide (Eichling et al., 1975;
Leggett and Williams, 1991; Phelps et al., 1979; Powers and Raichle,
1985), and single photon emission computed tomography studies
using 99 mTc-pertechnetate labeled red blood cells (Toyama et al.,
1990). However, exposure to ionizing radiation and the lack of spatial
resolution in radioactive tracer approaches pose critical limitations for
the accurate quantification of CBV and vessel information. Recently,
⁎ Corresponding author at: Unist-gil 50 (100 Banyeon-ri), Eonyang-eup, Uljugun, Ulsan
Metropolitan City, 689-798, Republic of Korea. Fax: +82 52 217 2509.
E-mail address: [email protected] (H. Cho).
http://dx.doi.org/10.1016/j.neuroimage.2015.03.042
1053-8119/© 2015 Elsevier Inc. All rights reserved.
several methods to obtain microvascular information using steadystate MRI and various contrast agents have been proposed to overcome
the limitations of radioactive tracer approaches (Kuppusamy et al.,
1996; Moseley et al., 1992; Schwarzbauer et al., 1993). One of these
methods is the use of superparamagnetic iron oxide nanoparticles
(SPION), which have a long blood half-life (Wu et al., 2004a). SPION induces sustained intravascular susceptibility perturbations in the surrounding tissue, as reflected by an increase in the transverse magnetic
resonance (MR) relaxation rates, i.e., ΔR2⁎ and ΔR2 from gradient and
spin echoes, respectively (Loubeyre et al., 1999; Weissleder et al.,
1990). The ΔR2⁎ and ΔR2 values are often used to calculate the CBV and
microvascular volume (MVV) and have been widely used for microvascular characterizations (Yablonskiy and Haacke, 1994). However, these
MR transverse relaxation measurements alone may not provide a representative regional comparison of vascular parameters. For example, the
capillary direction with respect to main B0 magnetic field changes the
magnitude of susceptibility-induced field inhomogeneities around the
vessel and confounds the quantification of CBV from ΔR2⁎ and ΔR2 measurements (Boxerman et al., 1995; Tropres et al., 2001). Moreover, the
CBV of the local tissue from gradient-echo based ΔR2⁎ measurements
can be further biased in tissue neighboring large vessels or air–tissue interfaces (Kim et al., 2013). These geometrical and environmental variations often complicate the interpretation of MR transverse relaxation
rates based on CBV measurements in regions such as the cortex, where
S.H. Han et al. / NeuroImage 112 (2015) 382–389
vessels are fairly well aligned, and in the olfactory bulb and cerebellum,
which are close to air–tissue interfaces. In addition, an accurate knowledge of the in vivo SPION concentration is required to compensate the
dose dependency of transverse relaxation rates. Accordingly, when regional and dose-independent CBV measurement is necessary, the ratio
between tissue and vessel signal enhancements (i.e., fractional blood
volume) from the ΔR1-contrast is preferable. A positive ΔR1 ratio between the tissue and vessel is not influenced by local susceptibility
field variations from the capillary structure, and reflects the doseindependent concentration-ratio of the intravascular contrast agent.
Furthermore, when obtaining size information on the microvasculature,
the ΔR2⁎/ΔR2 ratio is often used to estimate the mean vessel size index
(mVSI), by taking advantage of the sustained strong transverse relaxation rates of SPION and molecular diffusion modeling in the presence
of a susceptibility-induced field (Dennie et al., 1998; Tropres et al.,
2001; Prinster et al., 1997). However, when biased susceptibilityinduced field inhomogeneities from variable capillary orientation or
air–tissue interface impair local gradient-based ΔR2⁎ measurements,
the mVSI estimations may also become unreliable.
Previously, for robust simultaneous measurements of the regional
CBV from ΔR1 contrast and the mVSI from ΔR2⁎/ΔR2 contrast, in vivo
administration of Gd-based intravascular T1 and SPION-based T2/T2*
contrast agents was commonly employed in order to extract both
microvascular parameters (Beaumont et al., 2009; Huang et al., 2013;
Kim et al., 2008). The recently introduced ultra-short echo time (UTE)
pulse sequence minimizes susceptibility perturbations (Glover et al.,
1992; Holmes and Bydder, 2005; Pauly et al., 1989; Robson et al.,
2003), and provides an unique opportunity for dual-contrast microvascular acquisitions using a single intravascular agent, namely SPION,
when combined with conventional spin-echo acquisition (Jung et al.,
2014). However, despite the potential importance, direct and systematic comparisons between UTE-based and ΔR2⁎-based CBV measurements
via UTE acquisitions have not been reported as of yet.
In this study, Monte-Carlo spin dynamics simulations were performed to show the dependence of transverse relaxation rates on vessel
orientation and the unwanted bias in ΔR2⁎ of tissues near air-tissue interfaces. Additionally, the feasibility of alternative mVSI measurements
using UTE-based regional CBV values instead of ΔR2⁎ measurements was
proposed and validated. Experimentally, direct comparisons were performed between in vivo UTE-based and ΔR2⁎-based CBV maps in the cortex, cerebellum, and olfactory bulb regions in nude mice, which clearly
showed inconsistent ΔR2⁎ measurements at the vicinity of the air-tissue
interface and large vessels. Correspondingly, for the validation of mVSI
mapping with UTE acquisition, in vivo ΔR2⁎/ΔR2 versus UTE-based CBV/
ΔR2 maps were compared in the same brain regions. Next, the regional
relationships between the UTE-based and the ΔR2⁎-based CBV were
investigated in the olfactory bulb, cortex, caudate putamen, thalamus,
superior colliculus, and cerebellum in a large cohort. The regional relationships between the ΔR2⁎/ΔR2 versus UTE-based CBV/ΔR2 were investigated as well. The UTE-based CBV was compared with previously
reported regional CBVs in multiple studies. The dose dependency of
the UTE-based and the ΔR2⁎-based regional CBV were compared.
Method
Monte Carlo (MC) simulations
General methodology
All MC simulations in this study were performed using an in-house
Matlab script, which is similar to the method described by Pathak
et al. (2008). First, 3 dimensional-binary matrices were generated to
emulate the vascular structure using randomly oriented cylinders. The
number of cylindrical vessels was set to match the blood volume fraction (Vb), e.g., approximately 100 cylinders were used for a 2% Vb, and
the number of cylinders was increased proportionally for 5% and 10%
Vb. The vessel length was set to be 40 times greater than the vessel
383
radius (Rv). The binary matrix for the vascular structure was generated
three times to obtain the mean and standard deviation for each data
point in the relaxation simulations. The, finite perturber method
(FPM) was then used to calculate the magnetic field shift caused by
the intravascular contrast agent filling the cylinder using the Fourier
convolution theorem (Pathak et al., 2008). In FPM, the magnetic field
shift by an arbitrary structure is calculated as the sum of the magnetic
field shifts (ΔB) by uniformly packed small perturbers within the structure. The FPM inputs were the binary matrix of the vascular structure,
the magnetic susceptibility difference between the intravascular contrast agent and surrounding tissue (Δχ), and the main magnetic field
(B0). Finally, the Monte Carlo simulation was performed to calculate
the MR transverse relaxation rates (ΔR2 and ΔR2⁎) caused by proton diffusion in the presence of a susceptibility induced inhomogeneous magnetic field. Initially, 200,000 protons were positioned randomly within
the diffusion space that was filled by the cylinders. Each proton
moved to an adjacent position with each increment of the simulation
time (Δt) in a random trajectory and at a standard deviation of displacepffiffiffiffiffiffiffiffiffiffiffiffi
ment following the diffusion length ( 2DΔt , D is the diffusion coefficient). Subsequently each spin experienced different magnetic field
shifts at corresponding positions. This process induced the phase accumulation of each spins at each diffusion step. The total accumulated
phase of the nth proton during time t is given by the following equation:
Xt=Δt
γΔBðP n ðiΔt ÞÞΔt
ð1Þ
ϕn ðt Þ ¼
i¼1
where Δt is the diffusion time step, γ is the proton gyromagnetic ratio
(2.675 × 108 rad/T/s), Pn (t) is the position of the nth proton at time t,
and ΔB (P) is the magnetic field shift at position P. The MR signal was
then calculated from the mean accumulated phase of N total diffusing
protons:
SðtÞ ¼
1 XN iΦn ðt Þ
e
:
n¼1
N
ð2Þ
The change in the associated relaxation rates were calculated as
follows:
ΔR2 ; ΔR2 ¼ −
ln ½SðTEÞ
:
TE
ð3Þ
Several parameters were kept consistent in all the simulations. Each
diffusion process used 1000 randomly placed protons, and was repeated
200 times to obtain the mean MR signal. Hence, the number of total
spins (Nt) was 2 × 105. The diffusion coefficient (D) was set to
10−5 cm2/s, the diffusion time step (Δt) to 0.2 ms, and the vessels
were assumed impermeable to molecular diffusion. The Vb, Rv, B0, Δχ,
TE (SE), and TE (GE) parameters were set as variables for the particular
simulation.
Validation of the MC simulations
The simulation of the ΔR2 and ΔR2⁎ was performed as a function of
vessel radius (Supplementary Fig. 1A) (Boxerman et al., 1995; Pathak
et al., 2008). The parameters were identical to the previously reported
values (Boxerman et al., 1995; Pathak et al., 2008): Vb was set to 2%,
B0 to 1.5 T, Δχ to 10−7, TE (SE) to 100 ms, and TE (GE) to 60 ms. The
Rv of each vascular structure varied from 1 μm to 23 μm at a fixed Vb.
To determine representative experimental conditions, the dependency
of ΔR2 and ΔR2⁎ as a function of TE was investigated at two different
magnetic susceptibilities of the intravascular contrast agent (Δχ),
2 × 10−7 and 4 × 10−7 (Supplementary Fig. 1B). The main magnetic
field (B0), blood volume (Vb), and vessel diameter (VD) were set to
3 T, 2%, and 2 μm, respectively.
Vulnerability of transverse relaxation rates
The dependency of ΔR2 and ΔR2⁎ values on the vessel orientation and
air–tissue interface is shown in Fig. 1. In both simulations, the Vb was set
384
S.H. Han et al. / NeuroImage 112 (2015) 382–389
(A) Effect of vessel orientation
(A-1)
(A-2)
20
40
along x
along y
along z
B0, z-axis
14
35
ΔR2* [1/s ]
16
12
10
140
along x
along y
along z
120
ΔR2* , ΔR 2 [1/s]
18
ΔR2 [1/s ]
(B) Effect of air-tissue interface
30
25
20
15
8
6
−20
0
20
40
60
80
Relative portion of randomly
oriented vessels [%]
100
120
100
80
60
40
20
10
−20
0
20
40
60
80
100
0
−0.05
120
0
0.05
0.1
0.15
0.2
0.25
0.2
Bubble size [Radius/Space size]
Relative portion of randomly
oriented vessels [%]
Fig. 1. ΔR2 and ΔR⁎2 curves due to vessel orientations as a function of randomness (Figs. 1A-1 and A-2) and air–tissue interface within the tissue as a function of bubble size (Fig. 1B). The
total blood volume and vessel radius (Rv) were set to 2% and 2 μm, respectively. Magnetic susceptibility differences between pre- and post-SPION injection (ΔχSPION), and between the
captured air-bubble and tissue (Δχair-bubble) were set to 4 × 10−7 and 7 × 10−7 ,respectively. The blue, red, and black lines indicate the vessel alignment along the x, y, and zdirections, respectively. SPION, superparamagnetic iron oxide nanoparticles.
to 2%, Rv to 2 μm, B0 to 3 T along the z-direction, magnetic susceptibility
difference of the intravascular contrast agent to 4 × 10−7, the TE (SE) to
100 ms, and the TE (GE) to 15 ms. In the simulation of vessel orientation
effect for ΔR2 (Fig. 1A-1) and ΔR2⁎ (Fig. 1A-2), all the vessels were initially set to be aligned along a certain direction, e.g. x-axis, as shown in the
vessel structure of Fig. 1A-1. Then, the relative portion of randomly oriented vessels was increased from 0% (all aligned) to 100% (all randomly
oriented), while keeping the total Vb. This simulation was performed for
three cases with the vessels aligned in the x, y, and z-directions. Next, in
the simulation investigating the effect of the air–tissue interface on the
transverse relaxation (Fig. 1B), the air-bubble was set as a sphere and
was located in the center of the vascular space, with all the vessels
(the total Vb of 2%) set to random orientation. The bubble radius was
then increased as shown in vascular structure of Fig. 1B. More specifically, the ratio of the bubble radius to the one directional length of total
vascular matrix was varied from 0 to 1/5. The magnetic susceptibility
difference between the captured air-bubble and tissue (Δχair-bubble)
was set to 7 × 10−7, which was estimated based on a tissue magnetic
susceptibility of ~ 7 × 10−7 (Schenck, 1996), and an air (80% Nitrogen
and 20% Oxygen gas) magnetic susceptibility of ~ 5.6 × 10− 9 (Glick,
1961). Additionally, each proton was not allowed to exist inside the
air-bubble initially and the air-bubble was assumed impermeable.
Mean vessel size simulations
The vessel size dependency was examined in four MR vascular
parameters, ΔR2, ΔR2⁎, ΔR2⁎/ΔR2, and the UTE-based CBV/ΔR2. The
strength of the magnetic field (B0) was 3 T and the magnetic susceptibility difference between pre- and post-SPION injection (ΔχSPION) was set
to 2 × 10−7 and 4 × 10−7 for the simulation. These ΔχSPION values were
estimated based on the administered SPION dose in the animal experiment (Kim et al., 2013). The fractional blood volume (Vb) was set to 2%,
5%, and 10%, and these values were used for the corresponding UTEbased CBV in the simulations. As shown in Fig. 2, the vessel radius (Rv)
varied from 1 μm to 23 μm for each given Vb (2%, 5%, and 10%). For the
calculation of ΔR2 and ΔR2⁎, the TE (SE) and TE (GE) were set to
100 ms and 15 ms, respectively.
In vivo experiments
All studies were performed on a 3T MRI scanner (Philips Achieva,
Best, Netherlands). The Institutional Animal Care and Use Committee
of the Korea Basic Science Institute approved the in vivo animal experiments conducted in this study.
Animal and contrast agent protocols
A total 29 male nude mouse (20g ~ 28 g) were used in these studies.
The contrast agent used in this study was SPION, whose particle size,
half-life and MR relaxivity values were characterized in a previous study
(Jung et al., 2014). In most of the studies, UTE, MEGE and MSME MR images were acquired before and after SPION was intravenously administered at 70 μmol/kg (n = 21). For the dose dependency
characterization, higher SPION doses were administered at 120 μmol/kg
(n = 4) and 140 μmol/kg (n = 4).
Regional CBV measurements
An MRI anatomical brain atlas template was registered to the T2⁎ images from the MEGE sequence in order to obtain regional CBV information. The MRI anatomical brain template was accessed at http://
scalablebrainatlas.incf.org/main/coronal3d.php?template=WHS10.
Slice averaging of the template was performed due to the differing slice
thickness between the template and T2⁎ images. Representative six brain
regional masks were shown in Supplementary Fig. 2.
CBV measurements from UTE acquisitions
The UTE images were acquired under the following parameters in
order to maximize the intravascular T1 vessel/tissue contrast while
minimizing the T2⁎ effect of SPION: FA = 90°, TR = 50 ms, TE = 90 μs,
FOV = 5 × 5 × 8 cm3, and resolution = 195 × 195 × 1000 μm3. The
signal intensity from UTE is presented by before (Stissue _ pre = (1 −
Vb)Sextra + VbSintra) and after (Stissue _ post = (1 − Vb)Sextra + VbSintra′)
SPION administration based on a two-compartment tissue model with
intravascular and extravascular spaces. The Vb is the blood volume
fraction. The Sextra and Sintra are the signal intensities of the extravascular
and intravascular compartments before SPION administration, respectively. After administering the intravascular contrast agent, the signal in
the intravascular compartment is only altered due to the T1-shortening
effect of the intravascular contrast agent. (FA = 90° was used to minimize
the extravascular signal alteration from trans-vascular water exchange in
the presence of T1-shortening intravascular contrast agent.) The Sintra′
represents the intravascular space signal after injection. Since Sblood_pre
and Sblood_post can be replaced by Sintra and Sintra′ respectively, and they
are measurable on images, Vb was quantified from the following
equation:
Vb ¼
Stissue
Sblood
post −Stissue pre
post −Sblood pre
:
ð4Þ
S.H. Han et al. / NeuroImage 112 (2015) 382–389
-7
= 2 x 10 , B0 = 3T
(A-2)
20
15
50
40
30
10
20
5
10
0
0
0
100
10
10
0
10
0.1
−0.1
10
−0.3
10
−0.5
10
100
0
Radius [ µm]
10
100
Radius [ µm]
= 4 x 10 , B0 = 3T
(B-1)
(B-3)
(B-2)
50
40
30
20
−0.1
10
140
120
R2*/ R2
BV 2%
BV 5%
BV 10%
60
R2* [1/s]
160
100
80
60
40
10
20
0
0
0
10
-7
(B)
R2 [1/s]
0
100
Radius [ µm]
Radius [µm]
70
10
UTE-based CBV/ R2
R2 [1/s]
25
(A-4)
(A-3)
80
70
60
UTE-based CBV/ R2
BV 2%
BV 5%
BV 10%
30
R2* [1/s]
(A-1)
35
R2*/ R2
(A)
0
385
10
Radius [µm]
100
0
10
100
0
Radius [ µm]
10
100
(B-4)
10
−0.3
10
−0.5
10
−0.7
10
−0.9
10
0
Radius [ µm]
10
100
Radius [ µm]
Fig. 2. Vessel size dependence of four MR parameters, ΔR2 (Figs. 2A-1 and B-1), ΔR⁎2 (Figs. 2A-2 and B-2), ΔR⁎2 /ΔR2 (Figs. 2A-3 and B-3), and UTE-based CBV/ΔR2 (Figs. 2A-4 and B-4).
Figs. 2A and B show the effect under two different magnetic susceptibility values, 2 × 10−7 and 4 × 10−7, respectively. The red, blue, and black lines correspond to simulation performed
at a blood volume of 2%, 5%, and 10%, respectively. MR, magnetic resonance.
CBV measurements from MEGE acquisitions
The MEGE images were acquired using the following parameters to
determine ΔR2⁎: TR = 3000 ms, FOV = 5 × 5 × 2.4 cm3, resolution =
195 × 195 × 1000 μm3, and TE = 7, 13, 19, 25, 31, and 37 ms. To extract
ΔR2⁎ values in accordance with simulation, data points whose signal amplitudes were 5 times above the noise level were linearly fitted on a
semi-log plot and the signals before and after SPION administration
were obtained at TE = 15 ms. Then, the signal differences before and
after SPION administration were converted to ΔR2⁎ values at TE =
15 ms, using the single exponential model (ΔR2⁎ = ln(Spre/Spost)/TE)
(Tropres et al., 2001). Assuming that the static dephasing was
dominant, the fractional CBV (%) was estimated from the linear relationship relative to the ΔR2⁎ as follows (Tropres et al., 2001):
ΔR2 ¼
4
π ð1−Hct Þ V b Δχ γB0
3
ð5Þ
where Vb is the blood volume fraction, γ is the gyromagnetic ratio
(2.675 × 108 rad/T/s), Hct is the hematocrit set at 0.4, and Δχ is the
susceptibility difference of the contrast agent (0.36 ppm) at 3 T and
15 mg Fe/kg (Kim et al., 2013).
ΔR2⁎-based and UTE-based mVSI measurements
The MSME images were acquired using the following parameters to
determine the ΔR2 : TR = 3800 ms, FOV = 5 × 5 × 2.4 cm3, resolution =
195 × 195 × 2000 μm3, and TE = 60 and 100 ms. The signal differences
before and after SPION administration at TE = 60 ms were converted to
ΔR2 values, using the single exponential model (ΔR2 = ln(Spre/Spost)/TE)
(Dennie et al., 1998; Tropres et al., 2001; Prinster et al., 1997). The signals at TE = 100 ms were not used due to lowered signal to noise ratio.
The ΔR2⁎-based mVSI was calculated from the ratio of ΔR2⁎ and ΔR2, obtained with MEGE and MSME acquisitions, respectively (Dennie et al.,
1998; Tropres et al., 2001; Prinster et al., 1997). The UTE-based mVSI
was calculated from the ratio of the UTE-based CBV and ΔR2, obtained
with UTE and MSME acquisitions, respectively. As MSME measurements were performed with slice thickness of 2 mm, slice averaging
of the UTE and MEGE images was performed for the following mVSI
analyses.
Results
Validation of the MC simulations
The dependency of ΔR2 and ΔR2⁎ on vessel size was evaluated to validate the computational accuracy of our simulation, compared to previous results (Boxerman et al., 1995; Pathak et al., 2008). As shown in
Supplementary Fig. 1A, the ΔR2 and ΔR2⁎ in the current simulation
agree with the previous values obtained under the same conditions
(B0 = 1.5 T and Δχ = 1 × 10−7). The dependency of ΔR2 and ΔR2⁎ on
TE was also evaluated to optimize experimental TE values for the
MSME and MEGE acquisitions at 3 T. As shown in Supplementary
Fig. 1B, both ΔR2 and ΔR2⁎ drastically increased as the TE increased,
and rate of increase decreased at TE values of ~ 30 ms for ΔR2
and ~ 5 ms for ΔR2⁎.
Vulnerability of transverse relaxation rates
The effects of vessel orientation (Fig. 1A) and air–tissue interface
(Fig. 1B) on ΔR2 and ΔR2⁎ were examined. Fig. 1A showed the behavior
of ΔR2 and ΔR2⁎ as the relative portion of randomly oriented vessels
increased, while maintaining the total Vb at 2%. The blue, red, and black
lines indicate the initial vessel alignment along the x-, y-, and zdirections, respectively. Similar trends were observed between ΔR2 and
ΔR2⁎ as the relative portion of randomly oriented vessels increased in all
386
S.H. Han et al. / NeuroImage 112 (2015) 382–389
directions (Fig. 1A). The ΔR2 and ΔR2⁎ decreased as the relative portion of
randomly oriented vessels increased that were initially aligned along the
x- and y-directions, but the opposite relationship was observed shown in
vessels that were initially aligned along the main magnetic field direction
(z-direction). Fig. 1B shows the effect of the air–tissue interface on the ΔR2
and ΔR2⁎. Generally, both ΔR2 and ΔR2⁎ (blue and red lines, respectively)
increased as the radius of the air-bubble increased from zero to onefifth of the total space size. However, the increasing rate of ΔR2⁎ was
much higher than that of ΔR2, which shows that the gradient-echo
based ΔR2⁎ value is more sensitive to influence by the air–tissue interface
than the spin-echo based ΔR2 value is.
Mean vessel size index (mVSI) simulations
The vessel size dependency of the four MR parameters, ΔR2, ΔR2⁎,
ΔR2⁎/ΔR2, and UTE-based CBV/ΔR2 at 3T was shown in Fig. 2. Figs. 2A
and B showed the corresponding results at both magnetic susceptibility
values of 2 × 10−7 and 4 × 10−7, respectively. The red, blue, and black
lines correspond to the cases performed at a total Vb of 2%, 5%, and 10%,
respectively. At Δχ = 2 × 10−7, the ΔR2 peaked at an approximately
2 μm-radius and decreased as the radius increased (Fig. 2A-1). The
ΔR2⁎ increased and plateaued at an approximately 3 μm-radius
(Fig. 2A-2). This pattern was observed at all three Vb conditions (2%,
5%, and 10%). At the increased value of Δχ = 4 × 10−7, the ΔR2 steadily
decreased as the radius increased (Fig. 2B-1). The ΔR2⁎ quickly plateaued
to a constant value over the entire radius range for all Vb (Fig. 2B-2).
Next, the ΔR2⁎/ΔR2 and UTE-based CBV/ΔR2 were plotted as a
function of vessel size for all Vb. At both Δχ values (2 × 10− 7 and
4 × 10−7), neither the ΔR2⁎/ΔR2 nor the UTE-based CBV/ΔR2 showed
any significant dependency on the blood volume (Figs. 2A-3, A-4, B-3,
and B-4). For Δχ = 2 × 10−7, the UTE-based CBV/ΔR2 pattern was in
a good agreement with the ΔR2⁎/ΔR2 pattern, except at a radius of
1 μm; the UTE-based CBV/ΔR2 decreased from 1 μm to 2 μm and became
proportional to the increase in radius (Fig. 2A-3). The ΔR2⁎/ΔR2 steadily
increased as the radius increased (Fig. 2A-4). At Δχ = 4 × 10− 7, the
UTE-based CBV/ΔR2 showed a steady increase, which was consistent
with the ΔR2⁎/ΔR2 pattern (Figs. 2B-3 and B-4). This steady increase as
a function of vessel size and its independency from blood volume
make the UTE-based CBV/ΔR2 suitable for robust mVSI estimation at
our experimental conditions.
(A-1)
(B-1)
(C-1)
UTE-based CBV(%)
Direct comparison of CBV and mVSI maps from UTE, MEGE, and MSME
Representative voxel-wise maps of the regional CBV and mVSI calculated from the UTE, MEGE and MSME sequences were shown in Fig. 3.
The T2⁎-weighted images of the cortex, cerebellum, and olfactory bulb
regions were shown in Fig. 3A. The UTE-based CBV maps from Eq. (4)
are shown in Fig. 3B. The MEGE- and MSME-derived ΔR2⁎ and ΔR2 values
were shown in Figs. 3C and D, respectively. Conventional mVSI maps
using ΔR2⁎/ΔR2 were shown in Fig. 3E. The alternative mVSI maps
constructed using UTE-based CBV/ΔR2 were shown in Fig. 3F.
There were several notable differences between the ΔR2⁎-based and
UTE-based CBV maps in Figs. 3B and C. As indicated in Figs. 3B-1 and C-1
with black arrows, the UTE-based CBV showed positive values as the
blood volume increased, but the ΔR2⁎ showed negative values in the corresponding region. For the region indicated by black arrows in Fig. 3B-2,
no signal enhancement was observed in either map of the air–tissue interface region near both ears; this area was also negatively enhanced on
the gradient-echo acquisition shown in Fig. 3C-2. The UTE-based CBV
was able to clearly distinguish the high blood volume region from the
air–tissue interface, while the gradient-echo-based ΔR2⁎ failed to identify the corresponding regions. In addition, a spatial shift was observed in
the image on the ΔR2⁎ maps due to the field inhomogeneity compared
with the UTE-based CBV and ΔR2 maps, which would require an additional registration process for the direct voxel comparison and analysis.
Consequently, the resulting conventional mVSI (ΔR2⁎/ΔR2) showed an
unstable voxel map.
As indicated in Fig. 3E-1 with a white arrow, the mVSI was decreased
in cortical regions containing large vessels. However, the corresponding
mVSI showed a higher mVSI value than in other regions (Fig. 3F-1,
white arrow). The conventional mVSI map in the cerebellum was severely inconsistent between the hemispheres due to ΔR2⁎ artifacts near
the air-tissue interface (Fig. 3E-2, white arrow) compared to the UTEbased mVSI (Fig. 3F-2, white arrow). The olfactory bulb region showed
similar artifacts due to the large vessel penetrating perpendicular to
the main field (Fig. 3E-3, white arrow) compared to the UTE-based
mVSI (Fig. 3F-3, white arrow).
Comparison of regional CBV and mVSI values
Direct comparisons of the ΔR2⁎-based CBV (%) and UTE-based CBV (%)
values for the olfactory bulb, cortex, caudate putamen, thalamus, superior
(D-1)
10
0
−10
(E-1)
R2 [1/s]
R2* [1/s ]
50
0
−50
R2*/ R2 [1/s]
5
0
−5
(F-1)
UTE-based CBV(%)/ R 2
20
0
−20
5
0
−5
Cortex
(A-2)
(B-2)
(C-2)
10
0
−10
(D-2)
50
0
−50
(E-2)
5
0
−5
(F-2)
20
0
−20
5
0
−5
Cerebellum
(A-3)
(B-3)
(C-3)
10
0
−10
(D-3)
50
0
−50
(E-3)
5
0
−5
(F-3)
20
0
−20
5
0
−5
Olfactory bulb
Fig. 3. MR imaging of nude mice before and after SPION administration in the cortex, cerebellum, and olfactory bulb. (A) T⁎2 MR image before SPION administration. (B) UTE-based CBV
map. (C) ΔR⁎2 -based CBV map. (D) ΔR2 map. (E) ΔR⁎2 -based mVSI (ΔR⁎2 /ΔR2). (F) UTE-CBV based mVSI (UTE-CBV/ΔR2). MR, magnetic resonance; SPION, superparamagnetic iron oxide
nanoparticles; CBV, cerebral blood volume; mVSI, mean vessel size index.
S.H. Han et al. / NeuroImage 112 (2015) 382–389
colliculus, and cerebellum were shown in Fig. 4A. The standard deviation
of the measured CBV values was marked on the plot. The gradient-echobased ΔR2⁎-based CBV (%) showed a significantly larger standard deviation than that of the UTE-based CBV (%). By contrast, ΔR2⁎-based and the
UTE-based CBVs across the different brain regions were linearly correlated when a large sample size was included. To further validate its feasibility, the UTE-based CBV in various brain regions were compared with
previously reported values shown in Table 1. Each CBV ratio was calculated by dividing the regional CBV by the cortical CBV in order to identify the
correlation between regional CBV values in the different studies. The
cortex, thalamus, caudate putamen, and cerebellum generally showed
CBV values consistent with earlier reports; our CBV was close to or lied
between the minimum and maximum values in prior studies. For instance, the measured CBV in the cortex was 2.03, which ranged from 1
to 4 in earlier reports, indicating agreement between our cortical CBV
the earlier values. The same trend was observed in the CBV of the thalamus, caudate putamen, and cerebellum. However, in the olfactory bulb
and superior colliculus, our CBV was greater than the earlier values. The
CBV ratio was consistent with previous results in every brain region except the superior colliculus. Accordingly, regional ΔR2⁎/ΔR2 and UTEbased CBV/ΔR2 values were compared as shown in Fig. 4B. Numerical
values of respective experimental mVSI estimations consistently fell within the range of simulational calculations as shown in Figs. 2A-3 and A-4.
Considering apparent good coefficient of determination (r2 = 0.90) between UTE-based CBV and ΔR2⁎-based CBV across the different brain regions from Fig. 4A, UTE-based mVSI showed reduced linearity (r2 =
0.04, r2 w/o olfactory = 0.49, and r2 w/o olfactory & thalamus = 0.85)
with respect to ΔR2⁎-based mVSI. The reduced linearity between both
mVSI values is likely from the relative spatial shift of MEGE images with
respect to MSME and UTE images, which is conspicuous in olfactory
bulb region as shown in Fig. 3C-3 due to susceptibility mismatch. Significantly larger standard deviations of ΔR2⁎-based mVSI were also observed
for olfactory bulb and thalamus regions, where air–tissue interfaces and
macro-vessels exist.
387
Table 1
Regional CBV compared between previous and the current study for the cortex, thalamus,
caudate putamen, olfactory bulb, superior colliculus, and cerebellum. The CBV ratio was
calculated by dividing the regional CBV with the cortex CBV in order to identify the correlation of the regional CBV to those of the previous investigations. CBV, cerebral blood
volume.
Cortex
Thalamus
Caudate
putamen
Olfactory
bulb
Superior
colliculus
Cerebellum
CBV
(ml∙100g−1)
CBV
ratio
Method
animal
4.07
1.63
1.13
1.03
0.934
3.26
2
1
1
1
1
1
1
1
Rat
Rat
Rat
Mouse
Rat
Rat
Rat
2.44
2.03
3.03
1.11
1.1
0.785
3.56
2.87
1
1
1.86
0.98
1.07
0.84
1.75
0.71
MRI (Tropres et al., 2001)
MRI (Schwarzbauer et al., 1997)
Autoradiography (Lear et al., 1990)
MRI (Wu et al., 2004b)
CT (Cremer and Seville, 1983)
MRI (Chugh et al., 2009)
Autoradiography
(Julien‐Dolbec et al., 2002)
CT (Everett et al., 1956)
MRI (our study)
MRI (Schwarzbauer et al., 1997)
Autoradiography (Lear et al., 1990)
MRI (Wu et al., 2004b)
CT (Cremer and Seville, 1983)
MRI (our study)
MRI (Tropres et al., 2001)
0.88
0.629
2
1.3
1.07
0.78
0.673
0.82
0.64
1.04
Autoradiography (Lear et al., 1990)
CT (Cremer and Seville, 1983)
CT (Everett et al., 1956)
MRI (our study)
MRI (Wu et al., 2004b)
Rat
Rat
Rat
Mouse
Mouse
2.3
1.31
1.13
1.27
MRI (our study)
MRI (Wu et al., 2004b)
Mouse
Mouse
1.03
4.82
1.17
1.09
4.28
3.78
1.10
2.37
1.14
1.17
1.31
1.86
CT (Cremer and Seville, 1983)
MRI (our study)
MRI (Wu et al., 2004b)
CT (Cremer and Seville, 1983)
MRI (Chugh et al., 2009)
MRI (our study)
Rat
Mouse
Mouse
Rat
Rat
Mouse
Rat
Mouse
Rat
Rat
Mouse
Rat
Mouse
Rat
Comparison of dose dependency in CBV from UTE and MEGE
The dose dependent behaviors of the ΔR2⁎-(1/s) and UTE-based
CBV (%) values in the different regions of interest (ROI) of the mouse
brain are shown in Figs. 4C and D, respectively. An increased contrast
dose elevated the Δχ (Kim et al., 2013), thus increasing the ΔR2⁎-(1/s),
as shown in Fig. 4C. Dose correction was required to calculate the
(B)
rte
x
pu
tam
en
tha
su
lam
pe
rio
us
rc
oll
icu
lus
ce
reb
ell
um
co
ca
ud
ate
ac
tor
yb
ulb
0
MEGE−CBV(1/s)
30
20
8
6
4
2
10
0
rte
x
pu
tam
en
tha
su
lam
pe
rio
us
rc
oll
icu
lus
ce
reb
ell
um
2
5
10
80
100
120
concentration (µmol/kg)
140
0
80
100
120
140
concentration (µmol/kg)
ud
ate
4
ulb
6
10
40
co
8
15
olfactory bulb
cortex
caudate putamen
thalamus
superior colliculus
cerebellum
12
2
ca
10
mean Vessel Size Index (mVSI)
12
3 x UTE-CBV/ΔR2(TE = 60ms)
ΔR2*(TE = 15ms)/ΔR2(TE = 60ms)
r (all) = 0.04, r (w/o olfactory) = 0.49
2
r (w/o olfactory & thalamus) = 0.85
tor
yb
r2 (all) = 0.90
14
(D)
50
2
ac
16
olf
Blood volume fraction (%)
UTE-CBV
MEGE-CBV(TE = 15ms)
18
olf
(C)
20
20
UTE−CBV(%)
(A)
ΔR2⁎-based CBV (%) shown in Fig. 4A, which was generally larger than
the value calculated from the UTE-based CBV (%). By contrast, the
UTE-based CBV showed constant behavior at different doses. This
dose independency proves the stability of measuring the CBV using
the UTE sequence.
Fig. 4. (A) UTE-based CBV and MEGE (ΔR2⁎)-based CBV in the olfactory bulb, cortex, caudate putamen, thalamus, superior colliculus, and cerebellum regions (slice thickness = 1 mm). The
coefficient of determination (r2) was 0.90. (B) Corresponding (3 × UTE-based CBV)/ΔR2 and ΔR2⁎/ΔR2 derived mVSI values in the olfactory bulb, cortex, caudate putamen, thalamus, superior colliculus, and cerebellum regions (slice thickness = 2 mm). The coefficient of determination (r2) were 0.04, 0.49, and 0.85 for all regions, regions without olfactory, and regions
without olfactory & thalamus, respectively. (C) The MEGE (ΔR2⁎)-based and (D) UTE-based CBV at SPION doses of 70 μmol/kg (n = 21), 120 μmol/kg (n = 4), and 140 μmol/kg (n = 4) in
the different regions. CBV, cerebral blood volume; mVSI, mean vessel size index; SPION, superparamagnetic iron oxide nanoparticles.
388
S.H. Han et al. / NeuroImage 112 (2015) 382–389
Discussion
Microvascular parameters such as the CBV and mVSI are typically
measured using transverse relaxation changes, i.e., the ΔR2⁎ and ΔR2 in
the presence of intravascular SPION, usually known as a T2 contrast
agent. However, transverse relaxation measurements on a gradient
echo sequence are easily biased by geometrical and environmental conditions, and may cause several potential artifacts when estimating microvascular parameters. In this study, by taking advantage of the dual
contrast ability of SPION, we theoretically and experimentally verified
the direct benefit of replacing the gradient-echo-based ΔR2⁎ measurement with an ultrashort echo time (UTE)-based ΔR1 contrast to generate robust CBV and mVSI maps.
As suggested by the Monte Carlo simulations, the experimental ΔR2⁎based CBV showed larger standard deviations within the cohort than
those from the UTE-CBV. Different geometrical alignments relative to
the main magnetic field, vascular environments, and shimming value
for each subject may account for the ΔR2⁎ variations. This larger variation
in the ΔR2⁎-based CBV should be considered when longitudinal studies
of small cohorts are performed. A ΔR1-based CBV should be less
influenced by these external parameters in a large cohort, because R1
measurement does not include spin signal dephasing due to molecular
diffusion in the presence of field inhomogeneity. This was verified by
the stable CBV measurements in the animal experiment (n = 29)
with UTE acquisition. Furthermore, the conventional mVSI (ΔR2⁎/ΔR2)
map showed significant artifacts near the air–tissue interfaces and large
vessels due to unaccounted variations in the ΔR2⁎ measurements. In comparison, the newly proposed mVSI method (UTE-based CBV/ΔR2)
showed consistent values in both hemispheres and distinguished the
high blood pool area from the air–tissue interface region, truly reflecting
the in vivo concentration variation in the SPION.
Another possible advantage of dual contrast, SPION-assisted MR lies
in the unambiguous CBV fraction (UTE-based CBV) that is obtained independent of the size of underlying vessels. When traditional vessel
size index (ΔR2⁎/ΔR2) information is necessary for the specific purpose,
the UTE-based CBV and ΔR2⁎/ΔR2 can be separately measured using a
single agent without the confounding effect resulting from multiple
contrast agents. In this way, accurate ΔR1 and sensitive ΔR2⁎ contrasts
may be synergistically combined to enhance in vivo SPION contrast for
tailored applications (Jung et al.,2014 ; Girard et al., 2011).
It should be cautioned that the UTE-based CBV mapping with SPION,
which is independent of the hematocrit value and susceptibility effect
still provides fractional CBV, while dynamic susceptibility contrast
(DSC)-MRI is known to provide absolute CBV with suitable pharmacokinetic modeling with arterial input function information (Calamante
et al., 1999; Cha et al., 2002). Direct comparisons of these values will
shed light on the potential opportunities and limitations of UTE-based
CBV. Additionally, even though minimized from UTE acquisition with
90° flip angle excitation, the underlying assumption of R1-insensitivity
of extravascular signal from intravascular SPION or perfusion effect is
not fully experimentally proven for the determination of UTE-based
CBV, and future investigations may be focused on either compensating
or utilizing these potential errors of apparent CBV for the quantification
of trans-vascular water exchange using UTE acquisition. Technically,
variations of achievable shortest UTE echo time and worsening of
temporal resolution is still challenging for the wide-spread use of UTE
acquisition and future technical advances in UTE-related hardware
and acquisitions/reconstructions may further facilitate UTE applications
to microvascular characterizations.
A potential translational limitation of our study is that the SPION we
used is not yet clinically approved. However, a chemically equivalent
iron oxide nanoparticle (Feraheme, AMAG Pharmaceuticals, Lexington,
MA, USA) has been clinically approved for the treatment of irondeficiency anemia and is under consideration for use in MRI (Christen
et al., 2013; D'Arceuil et al., 2013; Li et al., 2005). Additional characterization of the MR contrast properties of this nanoparticle may validate
the clinical applications of dual-contrast microvascular quantification
described in this study.
In summary, by taking advantage of the dual contrast property of
SPION, UTE acquisition minimizes susceptibility perturbations and enables regional CBV measurements from the vessel/tissue ΔR1 ratio,
while consecutive spin-echo acquisition provides robust vessel size
information through simultaneous ΔR2 measurements in the cortex,
cerebellum, and olfactory bulb areas of the animal brain, which is subject to significant susceptibility perturbation. Application of this method
in longitudinal follow-up investigations of microvasculature brain disease may enable clear quantification of microvascular parameters.
Supplementary data to this article can be found online at http://dx.
doi.org/10.1016/j.neuroimage.2015.03.042.
Acknowledgments
This work was supported by the National Research Foundation of
Korea Grants funded by the Korean Government (no. 2010-0028684
and no. 2014 R1A1A1 008255).
References
Aronen, H.J., Gazit, I.E., Louis, D.N., Buchbinder, B.R., Pardo, F.S., Weisskoff, R.M., Harsh,
G.R., Cosgrove, G., Halpern, E.F., Hochberg, F.H., 1994. Cerebral blood volume maps
of gliomas: comparison with tumor grade and histologic findings. Radiology 191,
41–51.
Aronen, H.J., Glass, J., Pardo, F., Belliveau, J., Gruber, M., Buchbinder, B., Gazit, I., Linggood,
R., Fischman, A., Rosen, B., 1995. Echo-planar MR cerebral blood volume mapping of
gliomas: clinical utility. Acta Radiol. 36, 520–528.
Beaumont, M., Lemasson, B., Farion, R., Segebarth, C., Rémy, C., Barbier, E.L., 2009. Characterization of tumor angiogenesis in rat brain using iron-based vessel size index MRI in
combination with gadolinium-based dynamic contrast-enhanced MRI. J. Cereb. Blood
Flow Metab. 29, 1714–1726.
Boxerman, J.L., Hamberg, L.M., Rosen, B.R., Weisskoff, R.M., 1995. MR contrast due to intravascular magnetic susceptibility perturbations. Magn. Reson. Med. 34, 555–566.
Calamante, F., Thomas, D.L., Pell, G.S., Wiersma, J., Turner, R., 1999. Measuring cerebral
blood flow using magnetic resonance imaging techniques. J. Cereb. Blood Flow
Metab. 19 (7), 701–735.
Cha, S., Knopp, E.A., Johnson, G., Wetzel, S.G., Litt, A.W., Zagzag, D., 2002. Intracranial mass
lesions: dynamic contrast-enhanced susceptibility-weighted echo-planar perfusion
MR imaging. Radiology 223, 11–29.
Christen, T., Ni, W., Qiu, D., Schmiedeskamp, H., Bammer, R., Moseley, M., Zaharchuk, G.,
2013. High‐resolution cerebral blood volume imaging in humans using the blood
pool contrast agent ferumoxytol. Magn. Reson. Med. 70, 705–710.
Chugh, B.P., Lerch, J.P., Yu, L.X., Pienkowski, M., Harrison, R.V., Henkelman, R.M., Sled, J.G.,
2009. Measurement of cerebral blood volume in mouse brain regions using microcomputed tomography. Neuroimage 47, 1312–1318.
Cremer, J.E., Seville, M.P., 1983. Regional brain blood flow, blood volume, and haematocrit
values in the adult rat. J. Cereb. Blood Flow Metab. 3, 254–256.
D'Arceuil, H., Coimbra, A., Triano, P., Dougherty, M., Mello, J., Moseley, M., Glover, G.,
Lansberg, M., Blankenberg, F., 2013. Ferumoxytol enhanced resting state fMRI and
relative cerebral blood volume mapping in normal human brain. Neuroimage 83,
200–209.
Dennie, J., Mandeville, J.B., Boxerman, J.L., Packard, S.D., Rosen, B.R., Weisskoff, R.M., 1998.
NMR imaging of changes in vascular morphology due to tumor angiogenesis. Magn.
Reson. Med. 40, 793–799.
Eichling, J.O., Raichle, M.E., Grubb, R., Larson, K.B., Ter-Pogossian, M., 1975. In vivo determination of cerebral blood volume with radioactive oxygen-15 in the monkey. Circ.
Res. 37, 707–714.
Everett, N.B., Simmons, B., Lasher, E.P., 1956. Distribution of blood (Fe59) and plasma
(I131) volumes of rats determined by liquid nitrogen freezing. Circ. Res. 4, 419–424.
Gibbs, J., Leenders, K., Wise, R., Jones, T., 1984. Evaluation of cerebral perfusion reserve in
patients with carotid-artery occlusion. Lancet 323, 310–314.
Girard, O.M., Du, J., Agemy, L., Sugahara, K.N., Kotamraju, V.R., Ruoslahti, E., Bydder, G.M.,
Mattrey, R.F., 2011. Optimization of iron oxide nanoparticle detection using
ultraschort echo time pulse sequences: comparisons of T1, T2*, and synergistic T1T2* contrast mechanisms. Magn. Reson. Med. 65, 1649–1660.
Glick, R., 1961. On the diamagnetic suceptibility of gases. J. Phys. Chem. 65, 1552–1555.
Glover, G.H., Pauly, J.M., Bradshaw, K.M., 1992. Boron‐11 imaging with a three‐dimensional reconstruction method. J. Magn. Reson. Imaging 2, 47–52.
Grubb Jr., R.L., Phelps, M.E., Ter-Pogossian, M.M., 1973. Regional cerebral blood volume in
humans: X-ray fluorescence studies. Arch. Neurol. 28, 38.
Hamberg, L., Macfarlane, R., Tasdemiroglu, E., Boccalini, P., Hunter, G., Belliveau, J.,
Moskowitz, M., Rosen, B., 1993. Measurement of cerebrovascular changes in cats
after transient ischemia using dynamic magnetic resonance imaging. Stroke 24,
444–450.
Holmes, J.E., Bydder, G.M., 2005. MR imaging with ultrashort TE (UTE) pulse sequences:
basic principles. Radiography 11, 163–174.
S.H. Han et al. / NeuroImage 112 (2015) 382–389
Huang, S., Kim, J.K., Atochin, D.N., Farrar, C.T., Huang, P.L., Suh, J.Y., Kwon, S.J., Shim, W.H.,
Cho, H., Cho, G., 2013. Cerebral blood volume affects blood–brain barrier integrity in
an acute transient stroke model. J. Cereb. Blood Flow Metab. 33, 898–905.
Julien‐Dolbec, C., Tropres, I., Montigon, O., Reutenauer, H., Ziegler, A., Decorps, M., Payen,
J.F., 2002. Regional response of cerebral blood volume to graded hypoxic hypoxia in
rat brain. Br. J. Anaesth. 89, 287–293.
Jung, H., Park, B., Lee, C., Cho, J., Suh, J., Park, J., Kim, Y., Kim, J., Cho, G., Cho, H., 2014. Dual
MRI T1 and T2(⁎) contrast with size-controlled iron oxide nanoparticles.
Nanomedicine: nanotechnology, Biology and Medicine 10, 1679–1689.
Kim, Y., Tejima, E., Huang, S., Atochin, D., Dai, G., Lo, E., Huang, P., Bogdanov, A., Rosen, B.,
2008. In vivo quantification of transvascular water exchange during the acute phase
of permanent stroke. Magn. Reson. Med. 60, 813–821.
Kim, S.G., Harel, N., Jin, T., Kim, T., Lee, P., Zhao, F., 2013. Cerebral blood volume MRI with
intravascular superparamagnetic iron oxide nanoparticles. NMR Biomed. 26,
949–962.
Kucharczyk, J., Roberts, T., Moseley, M.E., Watson, A., 1993. Contrast‐enhanced perfusion‐
sensitive MR imaging in the diagnosis of cerebrovascular disorders. J. Magn. Reson.
Imaging 3, 241–245.
Kuppusamy, K., Lin, W., Cizek, G.R., Haacke, E.M., 1996. In vivo regional cerebral blood volume: quantitative assessment with 3D T1-weighted pre-and postcontrast MR imaging. Radiology 201, 106–112.
Lear, J.L., Kasliwal, R., Feyerabend, A., 1990. Mapping regional cerebral vascular transit
time by simultaneous determination of local cerebral blood flow and local cerebral
blood volume. Metab. Brain Dis. 5, 155–165.
Leggett, R., Williams, L., 1991. Suggested reference values for regional blood volumes in
humans. Health Phys. 60, 139–154.
Li, W., Tutton, S., Vu, A.T., Pierchala, L., Li, B.S., Lewis, J.M., Prasad, P.V., Edelman, R.R., 2005.
First‐pass contrast‐enhanced magnetic resonance angiography in humans using
ferumoxytol, a novel ultrasmall superparamagnetic iron oxide (USPIO)‐based blood
pool agent. J. Magn. Reson. Imaging 21, 46–52.
Loubeyre, P., De Jaegere, T., Bosmans, H., Miao, Y., Ni, Y., Landuyt, W., Marchal, G., 1999.
Comparison of iron oxide particles (AMI 227) with a gadolinium complex (Gd‐
DOTA) in dynamic susceptibility contrast MR imagings (FLASH and EPI) for both
phantom and rat brain at 1.5 tesla. J. Magn. Reson. Imaging 9, 447–453.
Moseley, M.E., Chew, W.M., White, D.L., Kucharczyk, J., Litt, L., Derugin, N., Dupon, J.,
Brasch, R.C., Norman, D., 1992. Hypercarbia‐induced changes in cerebral blood volume in the cat: a 1H MRI and intravascular contrast agent study. Magn. Reson.
Med. 23, 21–30.
Pathak, A.P., Ward, B.D., Schmainda, K.M., 2008. A novel technique for modeling
susceptibility-based contrast mechanisms for arbitrary microvascular geometries:
the finite perturber method. Neuroimage 40, 1130–1143.
Pauly, J., Conolly, S., Nishimura, D., Macovski, A., 1989. Slice-Selective Excitation for Very
Short T2 Species. Proceedings of SMRM, p. 28.
389
Phelps, M., Huang, S., Hoffman, E., Kuhl, D., 1979. Validation of tomographic measurement
of cerebral blood volume with C-11-labeled carboxyhemoglobin. J. Nucl. Med. (U. S.
A.) 20.
Powers, W.J., Raichle, M.E., 1985. Positron emission tomography and its application to the
study of cerebrovascular disease in man. Stroke 16, 361–376.
Prinster, A., Pierpaoli, C., Turner, R., Jezzard, P., 1997. Simultaneous measurement of ΔR2
and ΔR2* in cat brain during hypoxia and hypercapnia. Neuroimage 6, 191–200.
Risberg, J., Lundberg, N., Ingvar, D.H., 1969. Regional cerebral blood volume during acute
transient rises of the intracranial pressure (plateau waves)*. J. Neurosurg. 31,
303–310.
Robson, M.D., Gatehouse, P.D., Bydder, M., Bydder, G.M., 2003. Magnetic resonance: an introduction to ultrashort TE (UTE) imaging. J. Comput. Assist. Tomogr. 27, 825–846.
Schenck, J.F., 1996. The role of magnetic susceptibility in magnetic resonance imaging:
MRI magnetic compatibility of the first and second kinds. Med. Phys. 23, 815–850.
Schwarzbauer, C., Syha, J., Haase, A., 1993. Quantification of regional blood volumes by
rapid T1 mapping. Magn. Reson. Med. 29, 709–712.
Schwarzbauer, C., Morrissey, S.P., Deichmann, R., Hillenbrand, C., Syha, J., Adolf, H., Nöth,
U., Haase, A., 1997. Quantitative magnetic resonance imaging of capillary water permeability and regional blood volume with an intravascular MR contrast agent. Magn.
Reson. Med. 37, 769–777.
Spetzler, R., Wilson, C., Weinstein, P., Mehdorn, M., Townsend, J., Telles, D., 1977. Normal
perfusion pressure breakthrough theory. Clin. Neurosurg. 25, 651–672.
Toyama, H., Takeshita, G., Takeuchi, A., Anno, H., Ejiri, K., Maeda, H., Katada, K., Koga, S.,
Ishiyama, N., Kanno, T., 1990. Cerebral hemodynamics in patients with chronic obstructive carotid disease by rCBF, rCBV, and rCBV/rCBF ratio using SPECT. J. Nucl.
Med. (U. S. A.) 31.
Tropres, I., Grimault, S., Vaeth, A., Grillon, E., Julien, C., Payen, J.F., Lamalle, L., Decorps, M.,
2001. Vessel size imaging. Magn. Reson. Med. 45, 397–408.
Warach, S., Li, W., Ronthal, M., Edelman, R., 1992. Acute cerebral ischemia: evaluation
with dynamic contrast-enhanced MR imaging and MR angiography. Radiology 182,
41–47.
Weissleder, R., Elizondo, G., Wittenberg, J., Rabito, C., Bengele, H., Josephson, L., 1990. Ultrasmall superparamagnetic iron oxide: characterization of a new class of contrast
agents for MR imaging. Radiology 175, 489–493.
Wu, E.X., Tang, H., Jensen, J.H., 2004a. Applications of ultrasmall superparamagnetic iron
oxide contrast agents in the MR study of animal models. NMR Biomed. 17, 478–483.
Wu, E.X., Tang, H., Asai, T., Yan, S.D., 2004b. Regional cerebral blood volume reduction in
transgenic mutant APP (V717F, K670N/M671L) mice. Neurosci. Lett. 365, 223–227.
Yablonskiy, D.A., Haacke, E.M., 1994. Theory of NMR signal behavior in magnetically inhomogeneous tissues: the static dephasing regime. Magn. Reson. Med. 32, 749–763.