Quantitative Assessment of Water Pools by T 1 p and T 2 p MRI in Acute Cerebral Ischemia of the Rat

Abstract

The rotating frame longitudinal relaxation magnetic resonance imaging (MRI) contrast, T_1ρ, obtained with on-resonance continuous wave (CW) spin-lock field is a sensitive indicator of tissue changes associated with hyperacute stroke. Here, the rotating frame relaxation concept was extended by acquiring both T_1ρ and transverse rotating frame (T_2ρ) MRI data using both CW and adiabatic hyperbolic secant (HSn; n = 1, 4, or 8) pulses in a rat stroke model of middle cerebral artery occlusion. The results show differences in the sensitivity of spinlock T_1ρ and T_2ρ MRI to detect hyperacute ischemia. The most sensitive techniques were CW-T_1ρ and T_1ρ using HS4 or HS8 pulses. Fitting a two-pool exchange model to the T_1ρ and T_2ρ MRI data acquired from the infarcting brain indicated time-dependent increase in free water fraction, decrease in the correlation time of water fraction associated with macromolecules, and increase in the exchange correlation time. These findings are consistent with known pathology in acute stroke, including vasogenic edema, destructive processes, and tissue acidification. Our results show that the sensitivity of the spinlock MRI contrast in vivo can be modified using different spinlock preparation blocks, and that physicochemical models of the rotating frame relaxation may provide insight into progression of ischemia in vivo.

Keywords

adiabatic pulses brain stroke

Introduction

Magnetic resonance imaging (MRI) is the modality of choice to diagnose acute cerebral ischemia owing to the fact that several endogenous contrasts develop on ischemia. The apparent diffusion coefficient of brain water rapidly decreases (Moseley et al, 1990) after compromised cerebral blood flow (CBF). It has been shown that apparent diffusion coefficient decrease has two CBF thresholds: a small decrease occurs close to CBF of 40 to 50 mL per 100 g/mins and an extensive one immediately below CBF of 20 mL per 100 g/mins (Gröhn et al, 2000b; Hossmann, 1994). The latter CBF value is close to the threshold for ischemic energy failure (Crockard et al, 1987; Naritomi et al, 1988). Diffusion MRI is part of stroke MRI protocols in clinical settings and provides excellent sensitivity and specificity to acute ischemia (Lythgoe et al, 1997; Mullins et al, 2002). It is a common observation both in stroke patients (Baird et al, 1997) and experimental models for stroke (Meng et al, 2004; Shen et al, 2003) that the volume of hypoperfusion, as determined by perfusion MRI, is larger than the volume displaying lowered diffusion. The so-called diffusion-perfusion mismatch is suggested to predict enlargement of an ischemic lesion.

Detection of acute ischemia by MRI can be accomplished by several other endogenous contrasts than diffusion and perfusion, including T₁, T₂, T^2*₂, and longitudinal relaxation time in the rotating frame (T1ρ). Factors leading to hyperacute T₁ prolongation in ischemia have been shown to include both CBF-dependent and -independent components (Calamante et al, 1999; Kettunen et al, 2000). Build up of deoxyhemoglobin shortens T₂ and TJ in the hyperacute phase of cerebral ischemia (Calamante et al, 1999; Grohn et al, 2008). T_1ρ obtained by on-resonance continuous wave (CW) spin-lock (SL) field has been shown to be a specifically interesting MRI contrast to assess hyperacute cerebral ischemia (Grohn et al, 2008). This is because of the ability of T1ρ MRI to indicate ischemia as follows: first, it reveals tissue response to compromised blood flow almost immediately after a drop of CBF below 20 mL per 100 g/mins (Grohn et al, 2008) and second, the increase in absolute parenchymal T1ρ correlates with long-term histologic outcome (Grohn et al, 2008). These factors make T^1ρ a very attractive MRI contrast to reveal acute stroke, as it predicts the tissue outcome in the acute phase when therapies based on clot removal and/or neuroprotection are most likely to be beneficial. Despite the limitations imposed by specific absorption rate (SAR) of radiofrequency (RF) energy, T_1ρ techniques have been applied to human brain imaging (Grohn et al, 2005; Michaeli et al, 2006; Wheaton et al, 2004).

Recently, the concept of rotating frame relaxation contrasts was extended by exploiting relaxation during adiabatic RF pulses. It was shown that in biologic systems rotating frame relaxation time constants obtained with adiabatic pulses enable assessment of slow molecular motion with an extended range of correlation times (Michaeli et al, 2005, 2006). Sensitivity to slow motion occurs because the rotating frame relaxation time constants are sensitive to magnetic field fluctuations caused by tumbling dipoles at or close to the effective frequency (Ω_eff) in the rotating frame (the Rabi frequency). The sensitivity of T1_ρ to molecular dynamics in the kHz range makes it a practical tool for studying dipolar fluctuations in vivo, to gain information about water dynamics and interactions with endogenous macromolecules (Liepinsh and Otting, 1996). Similarly to T_1ρ, the transverse relaxation time constant in the rotating frame, T_2ρ, is a function of dipolar interactions and is additionally affected by the exchange of protons between different magnetic sites. Adiabatic rotating frame time constants are dependent on the choice of amplitude- and frequency-modulation functions of the adiabatic full passage (AFP) pulse (Michaeli et al, 2005). Therefore, adiabatic rotating frame relaxation rates can be altered by using different modulation functions, which provides a basis to assess changes in different intrinsic tissue parameters probing water dynamics.

The aim of this work was twofold as follows: first, to investigate abilities of various T_ρ1 and Tρ magnetic resonance (MR) approaches with inherently differing SARs to reveal hyperacute ischemia in a rat model of middle cerebral artery occlusion (MCAo), and second, to evaluate alterations in water spin dynamics during evolving cerebral ischemia using the MR relaxation data in the context of an equilibrium two-site exchange (2SX) model (Michaeli et al, 2006).

Methods

Animals

Male Wistar rats (weight ranging from 280 to 320 g) were anesthetized with isoflurane (2.5% during surgery and 1.5% to 2% during MRI) with a constant flow of 70:30 N₂O/O₂ through a nose cone. Temporary MCAo (n = 8) for 90 mins was induced by the intraluminal thread method (Longa et al, 1989). A satellite group of male Wistar rats (n = 3) were operated similarly as those for temporary MCAo, but leaving the occluder thread in place for permanent MCAo. A shamoperated rat was used as a control. Arterial blood gases and pH were analyzed immediately before and during MR data acquisition (i-Stat Co, East Windsor, NJ, USA). The core temperature was monitored online and maintained close to 37° Cby warm water circulation in a heating pad placed under the torso. Blood pressure was monitored online during the scanning (CardioCap II, Datex, Helsinki, Finland). All animal experiments were approved by the ethics committee of the National Laboratory Animal Center of the University of Kuopio.

Magnetic Resonance Imaging

Magnetic resonance experiments were performed in a horizontal 4.7 T magnet (Magnex Scientific Ltd, Abingdon, UK) interfaced to a Varian Inova console (Varian, Palo Alto, CA, USA). A quadrature half-volume coil consisting of two loops with 20 mm diameters (HF Imaging LLC, Minneapolis, MN, USA) was used to transmit RF pulses and to receive MRI signals.

The transversal imaging plane (field of view 2.56 × 2.56cm², 128 × 64 points) was positioned so that the center of the slice (thickness 1.5 mm) was 4 mm from the surface of the brain. In the transient MCAo group the MR data were measured once during ischemia (between 60 and 90 mins after onset of MCAo), followed by retraction of the occluder thread at 90 mins without moving the animal in the magnet, and three more MR data sets were acquired between 0 and 30, 30 and 60, and finally, 60 and 90 mins after the onset of reperfusion. In the permanent MCAo group, MRI scans were acquired at two time points; the first set between 120 to 150 mins and the second between 230 and 260 mins from induction of ischemia.

Data for the trace of the diffusion tensor (D_av = 1/3* Trace(D)) map were first collected and used to localize the ischemic tissue (Kettunen et al, 2000). The D_av was quantified using a spin echo MRI sequence incorporating four bipolar gradients along each axis (Mori and van Zijl, 1995) with four b-values ranging from 0 to 1,370 secs/mm² (time to repetition (TR) 1.5 sees, time to echo (TE) 55 ms).

The MRI data for on-resonance T_1ρ were collected with a conventional CW-type technique consisting of the pulse sequence module: adiabatic half passage (AHP)-SL-reverse AHP (Figure 1A). This module was immediately followed by a fast spin echo readout sequence (echo spacing 10 ms, 16 echoes, center-out k-space filling, TR = 2.5 secs). The SL period of the preparation module (Grohn et al, 2005) ranged from 8 to 64 ms with a target SL amplitude B_1,SL = 0.8G. CW-type T_2ρ data were acquired with the following preparation module: AHP-rotary echo—reverse AHP (Solomon, 1959; Figure 1B). In this CW-T_2ρ preparation module, the first AHP pulse created transverse magnetization along the x' axis of the rotating frame (phase = 0°), and the phases of the subsequent rotary echo pulses alternated between +90° and + 270° every 2 ms. A final reverse AHP returned the magnetization to the longitudinal (z') axis and its magnitude was measured using the fast spin echo readout sequence. Spin-lock lengths and B₁ were identical to those used in the CW-T_1ρ experiments.

Rotating frame relaxation using adiabatic preparation pulses was measured in a similar manner, except the amplitude and frequency of the SL and rotary echo pulses were modulated under the adiabatic condition. The modulations were based on the hyperbolic secant family of pulses (HSn), using conditions and parameters that achieved magnetization inversion (i.e., AFP; Garwood and DelaBarre, 2001; Silver et al, 1984). Under this so-called ‘adiabatic condition’, the magnetization follows the trajectory of the effective magnetic field vector, B_eff. For convenience, the effective field is often expressed in units of angular frequency (rad per sec), and as such is given by Ω_eff = γB_eff, where g is the gyromagnetic ratio. The amplitude and orientation of oeff are time dependent and these in turn are functions of the amplitude modulations and frequency modulations used for the pulse. One way to alter the modulations is to change the stretching factor, n Accordingly, the amplitude modulation function of the HSn pulse becomes flatter as the factor n increases (Figure 1E), as given by

ω_{1} (τ) = ω_{1}^{max} sec h (β τ^{n})

(1)

where Ω^max₁ is the maximum pulse amplitude in rad per sec and β is a constant. When n = 1, the pulse is the original hyperbolic secant pulse (Silver et al, 1984), whereas when n > 1, a stretched version of the pulse is produced (Tannús and Garwood, 1996). In the experiments performed here, all HSn pulse parameters were held constant except for the stretching factor (n = 1, 4, or 8). With these adiabatic pulses, the rotating frame relaxation times (i.e., adiabatic T^1ρ and adiabatic T_2ρ) were measured using a preparation module consisting of consecutive identical HSn pulses, with no delays between pulses (i.e., negligible ‘RF-free’ intervals; Figure 1C). To obtain data for fitting, the number of HSn pulses was increased in separate experiments. Specifically, the train length increased from 4 to 32 pulses, each having a duration of 2 ms. The phases of the HSn pulses were prescribed according to MLEV4 (Levitt et al, 1982). For adiabatic T_2ρ measurements, AHP and reverse AHP pulses were added before and after the HSn pulse train (Figure 1D). These preparation modules preceded a fast spin echo readout sequence.

Figure 1

Schematic illustration of the preparation pulse sequences: conventional spin-lock (SL) T_1ρ (A), rotary echo T_2ρ (B), T_1&rH; using HSn pulses (C), and T_2ρ using HSn pulses (D). Spin-lock time in each is 8 ms. (E) Amplitude and frequency modulation functions for AFP pulses HS1, HS4, and HS8. AHP, adiabatic half passage.

Data for T₂ maps were collected with the same fast spin echo sequence with a T₂ preparation block consisting of AHP-TE/4-AFP-TE/2-AFP-TE/4-reverse AHP. To map the B₁ field, data were acquired using a variable length square preparation pulse and a crusher gradient in front of a FLASH pulse sequence (TR 4.5 ms, TE 2.2 ms), and a cosine function was fitted to the measured signal intensity (SI) oscillation.

For fitting the 2SX model, more complete T_1ρ and T_2ρ relaxation data sets using eight different HSn pulses (n =1, 2, …, 8) were acquired from a satellite group of permanent MCAo rats at two time points between. A sham-operated animal was used as control to the permanent ischemia group.

Data Analysis

All data analyses were performed with Matlab (The Math Works, Natick, MA, USA). Regions of interest (ROIs) were positioned to cover the ipsilateral striatum, which represents the core of the ischemic region, and the corresponding contralateral area. To evaluate whether the MRI contrasts used in this study could detect spatiotemporal progression during MCAo and reperfusion, a threshold-based analysis was performed. The threshold was set to be 15% above (or below for diffusion) the mean value in contralateral ROI consisting of entire striatum, and cortex and number of pixels above (or below for diffusion) this were calculated both in contralateral and corresponding ipsilateral ROIs. The difference (ipsi—contra) in the number of pixels passing threshold in ROIs was calculated and assumed to represent the amount of abnormal pixels. All results are presented as mean±standard deviation (s.d.). Paired Student's t-test was used for statistical analysis of data obtained from ipsilateral and contralateral hemispheres.

Under the intermediate and slow exchange conditions nonmonoexponential MR signal intensity decay is detected. In our experiments mono exponential signal intensity decay was observed, within the range of SL times used, and therefore fast exchange condition was assumed. For the interpretation of nuclear magnetic resonance data the dynamic processes were modeled with an equilibrium 2SX system approximation, comprising of two water populations coupled by an equilibrium exchange as follows: ‘A’ is the water population associated with macromolecules, and ‘B’ is the free water population. These two sites have relaxation times T_1ρ,A and T_1ρ,B that maintain equilibrium during RF pulse. At the two exchanging sites, the relaxation rate constants are modulated as oeff is time dependent. Thus, the sensitivity to different correlation times changes during the RF pulse. The formalism is detailed in the Appendix.

The fundamental MR parameters, that is, correlation times at specific sites A and B (τc,A and τc,B), the relative spin populations (P_A and P_B), and the exchange correlation time (τe,x), were estimated by fitting the 2SX model into the R_1ρ,obs and R_2ρ,obs measurements from each animal using the Levenberg—Marquardt nonlinear least squares regression algorithm. The fittings were performed using combined and data to find a set of parameters that optimally fit the measured data. To reduce the number of fitted parameters and to keep fitting unambiguous, the correlation time of free water τc,A and the chemical shift difference δΩ between exchanging sites were assumed to be 10⁻¹² secs and 125 Hz, respectively (Sierra et al, 2008). The resulting parameter estimates were averaged and tested for statistical significance using two-tailed two-sample Student's t-test.

Results

During the MCAo experiments physiologic parameters were within the physiologic range except pCO₂, which was slightly elevated (Table 1). Mild hypercapnia is commonly observed in spontaneously breathing animals under prolonged isoflurane anesthesia.

Table 1

Physiological parameters for sham-operated and stoke rats during MCAo and 30 mins after occlusion

	pCO₂ (mm Hg)	pO₂ (mm Hg)	pH	BP_mean (mmHg)	Temp (°)
MCAo	49±9	99±21	7.35±0.05	117±13	36.2±1.1
After occlusion	52±6	121±27	7.33±0.04	113±13	36.7±0.6
occlusion

BP blood pressure; MCAo, middle cerebral artery occlusion. Values are mean±s.d.

D_av decreased by approximately 20% by 60 mins of MCAo ipsilaterally consistent with severe ischemia and it did not recover on retraction of the MCA occluder (Figure 2B). T₂ was not significantly different between ipsilateral and contralateral brain by 60 mins (Figure 2A), but showed prolonged values at 90 mins of MCAo and T₂ continued to rise despite the retraction of the occluder.

Figure 2

T₂ (A) and diffusion (D_av) (B) values for ischemic (‘ipsi’) and contralateral (‘contra’) sides of the brain. Data are shown as mean±s.d. Vertical dotted line denotes the time of occluder removal.

T_1p measured with HS_n pulses from normal brain decreased with increasing n, whereas T_2ρ behaved in the opposite manner (Figure 3A). All T_1p relaxation time constants obtained with HSn pulses (values ranging from 100 to 180 ms) were longer than those with CW measurements (ranging from 80 to 100 ms). In contrast, the T_2ρ values acquired with HSn pulses (T_2ρ,HSn) were shorter (ranging from 60 to 85 ms; Figure 3A) than those obtained with the CW technique (ranging from 80 to 115 ms; Figure 3B).

Figure 3

T_1ρ and T_2ρ values measured using HSn pulses (A) and CW pulses (B) for ischemic (‘ipsi’) and contralateral (‘contra’) brain. Data are shown as mean±s.d. and obtained from animals with transient MCA occlusion. R_1ρ and R_2p values (C) measured with HSn (n = 1 to 8) for ischemic (‘ipsi’) and contralateral (‘contra’) brain from all three animals with permanent MCA occlusion and a control rat. The fitted lines for two time points (Tp 1 = 120 to 150 mins after MCAo, Tp 2 = 230 to 260 mins after MCAo) were calculated using the measured T_1ρ and T_2ρ values as described in Methods.

By 60 mins of MCAo, all ipsilateral T^1ρ and T^2ρ values increased relative to the contralateral ones and T^2p remained constant throughout the MRI observation time (Figures 3A and 3B). All rotating frame relaxation times significantly differed between ipsilateral and contralateral sides (P < 0.005) at all time points during stroke. The rotating frame relaxation times most sensitive to ischemic changes in the core of the lesion were T^1ρ,cw, T^1ρ,HS4, and T_1ρHS8, all showing increased values by 9% to 10% by 60 mins and by 13% to 15% by 150 mins after induction of MCAo, respectively (Table 2). Figure 4 shows an example of the calculated parametric maps acquired at 60 to 90 mins after MCAo and at 90 to 120 mins after reperfusion began.

Figure 4

Representative relaxation and diffusion maps acquired at 60 to 90 mins during MCA occlusion and 60 to 120 mins after beginning of reperfusion (after 90 mins of occlusion) showing the difference in contrast between the contralateral (left) and ischemic (right) rat striatum. Color scaling is adjusted to ± 20% of the contralateral striatum in all relaxation maps for easier comparison. Schematic brain outline adapted from Paxinos and Watson (1998).

Table 2

Average difference between the ischemic and contralateral side of ΔT₁ and ΔT₂ relaxation times, diffusion (Δ_av) and ΔT₂

	ΔT₁(%)	ΔT₂(%)	ΔT₁(%)
	CW	CW	HS1	HS4	HS8
60-90 mins	9±4**	6±3**	4±3**	10±3**	9±4**
90-120 mins	11±5**	7±5**	7±3**	10±5**	11±6**
120-150 mins	13±7**	9±6**	8±2**	12±7**	12±7**
150-180 mins	15±8**	11±5**	9±4**	13±7**	15±7**
	ΔT₂(%)			ΔD_av(%)	ΔT₂(%)
	HS1	HS4	HS8

60-90 mins	4±3**	4±2**	5±2**	−9±12**	0±3
90-120 mins	7±3**	7±4**	8±3**	−20±8**	6±3**
120-150 mins	10±4*	10±4**	11±3**	−22±11**	8±3**
150-180 mins	11±5**	12±5**	11±5**	−19±16*	9±5**

Values are in percent (mean±s.d.).

Significance was calculated by Student's t-test

P<0.05

P<0.01

When ROI analysis focused on the core of the lesion, threshold-based analysis was also performed to assess spatial progression or recovery in the periphery of the lesion (Figure 5). Diffusion and all measured relaxation times, except T₂ in the first measurement point, showed increased number of abnormal pixels in the ipsilateral side compared with the contralateral side. Importantly, spatiotemporal progression was evident from the increasing number of pixels in T_1ρ,cw, T_2ρ,cw, T_2ρ,HSn, and in T₂ contrasts whereas diffusion showed significantly decreased lesion area during reperfusion. These data show the complementary nature of these MR parameters to evaluate postischemic brain tissue.

Figure 5

Threshold analysis showing the progression of the lesion area at four time points after MCAo. Data are presented as difference (ipsi—contra) in the number of pixels in striatum and cortex with values over 115% (or under 85% for diffusion) of the mean value in contralateral side (for details see Methods). Data are shown as mean±s.d. Significance is shown as follows: *P < 0.05, **P < 0.01, different from zero; + P < 0.05, + + P < 0.01, different from the first time point.

For the quality control, possible spatial variation of the B₁ field within the field of view was determined by means of B₁ maps. The average γB₁ measured with the same transmitter power settings as used for rotating frame relaxation measurements was 3130±40 Hz (i.e., B₁ = 0.74±0.01 G) in ROIs covering the ischemic lesion. The mean difference between ipsilateral and contralateral sides was −1.1±0.4%, demonstrating insignificant B₁ variation between the hemispheres.

Magnetic resonance imaging data from the permanent MCAo rats were fitted to the equilibrium 2SX model (Michaeli et al, 2006) to assess changes in water dynamics during evolution of ischemia (Figure 3C). Using the combined HSn (n = 1 to 8) T1ρ and T^2ρ MR data as input, the 2SX model revealed increasing difference between ipsilateral and contralateral sides in the fitted parameters during progression of ischemia (Table 3). The relative size of the free water pool (site B) increased significantly during MCAo. At the same time, the rotational correlation time τ_c,A showed drastic decrease indicating increased mobility of the bound water fraction. Conversely, τ_ex increased in the brain parenchyma after induction of ischemia.

Table 3

Averaged fitted model parameters of the adiabatic T₁ and T₂ of water pools (averaged parameter values±s.d.) for contralateral and ipsilateral sides at two time points (Tp 1 = 120 to 150 and Tp 2 = 230 to 260 mins after MCAo)

	P _A	P _B	τ_c,A (sec)	τ_ex(sec)
Contra	0.38±0.02	0.62±0.02	(6.4±1.3)×10⁻¹⁰	(6.6±0.2)×10⁻⁵
Ipsi Tp 1	0.37±0.02	0.63±0.02	(4.0±0.9)×10⁻¹⁰*	(7.3±0.1)×10⁻⁵*
Ipsi Tp 2	0.34±0.02*	0.66±0.02*	(4.3±1.1)×10⁻¹⁰*	(7.01±0.03)×10⁻⁵*

Statistically significant changes were observed between contralateral and ipsilateral sides

P × 0.05, two-sample two-tailed Student's t-test).

Discussion

We examined the abilities of the adiabatic modulation functions and conventional CW-type RF pulses to generate rotating frame relaxation MR contrast for detection of hyperacute brain ischemia. It was observed that the adiabatic HS pulses used provide T_1p or T_2p relaxation time constants with inherent differences to detect time-dependent changes in acute stroke. Fitting the equilibrium 2SX model into the rotating frame relaxation data indicated substantial changes in the pool sizes, mobility, and exchange times of cerebral water pools during the evolution of ischemia. These data are taken to provide an insight into the physicochemical mechanisms underpinning the observed MR relaxation changes occurring in stroke.

In normal brain parenchyma, CW-T_1ρ was found to be shorter than T_1ρ acquired with HSn pulses. T_1ρ relaxation times measured with HSn pulses became shorter with increasing number of HS pulses in contrast to T_2ρ relaxation times, which behaved in the opposite fashion. These observations are consistent with the previous results from human brain and water/ethanol studies by Michaeli et al (2004, 2006). Our data show that rotating frame relaxation in vivo is indeed influenced by the choice of pulse modulation functions. Specifically, the effective SL frequency with the HSn pulses (ω_eff(t)) is off-resonance during a significant portion of the pulse duration, so the method therefore probes more of the shorter correlation time regime of molecular motion than the on-resonance CW-type technique (Gröhn et al, 2003).

Our results show that CW-T^1ρ and T_1ρ measured with HS4 and HS8 pulses provide a sensitive means to reveal brain ischemia. Our earlier work shows greater sensitivity of the on-resonance SL than the off-resonance T_1ρ to the tissue changes caused by ischemia during MCAo (Gröhn et al, 2003). In this regard the current data are in agreement with the previous observations (Gröhn et al, 2003). Especially, during adiabatic pulse modulation using HSn pulses, as n decreases in the HSn-type technique, the spins experience an effective field that is off-resonance during a greater portion of each pulse, which creates a less sensitive situation to detect acute stroke than what is obtained with CW-T^1ρ.

The sensitivity difference between T_1ρ and T_2ρ MRI to acute stroke could be considered being analogous to the established one between T_1ρ and T₂ T_1ρ relaxation time is elevated in the early minutes of ischemia, whereas increase inT₂ lags much behind. T_1ρ is predominantly sensitive to dipolar interactions at the specific sites and isochronous exchange, that is, the exchange of spins between sites with the same chemical shift, but with relaxation time constants that differ between sites: T_1ρ,A ≠ T_1ρ.B. As compared with T_2ρ, the value of T_1ρ is less influenced by anisochronous mechanisms that are thought to arise from dynamic processes on a microscopic scale, such as water diffusing between sites having different magnetic susceptibilities and/or exchange of spins between chemical environments with different chemical shifts. With both T_1ρ techniques (CW and HSn type), the influence of the latter types of frequency shifts on the measured T_1ρ is at least partially ‘suppressed’ by the effective locking field. This may be the main reason for the sensitivity of T_1ρ to water—protein interactions and to isochronous exchange. Because of suppression of magnetic susceptibility effects from the T_1ρ relaxation, the sensitivity of T_1ρ to blood oxygenation level dependent (Gröhn et al, 1998) is small. However, the anisochronous mechanisms significantly contribute to the T_2ρ relaxation. It has been recently shown that T_2ρ can be used to probe iron accumulation in the substantia nigra in Parkinson's disease (Michaeli et al, 2007) and potentially is able to differentiate patients with Parkinson's disease from age-matched controls on an individual basis. The work involving iron content in the brain implies that T_2ρ is sensitive to susceptibility effects in vivo.

Inherent properties of T_1ρ and T_2ρ together with the physicochemistry of ischemic tissue may reconcile their differences in revealing acute stroke. An almost immediate drop of temperature by 2°C to 3°C in the ischemic core (Gröhn et al, 1999) and acidification of tissue is expected to decrease contributions of the anisochronous mechanisms to the relaxation. Anisochronous effects may not be significant for the T_1ρ relaxation (Abragam, 1961; Michaeli et al, 2006), but they may contribute to the T_2ρ relaxation.

One of the motivations behind this study was to search for alternative approaches to the CW-T1ρ technique to lower SAR, a factor that limits the clinical applicability of the rotating frame MR techniques (Wheaton et al, 2004). We used a relatively high SL amplitude of 0.8 G, as contrast between ischemic and normal brain tissue with CW-T_1ρ has been shown to be a function of the SL amplitude (Gröhn et al, 2000a; Kettunen et al, 2005). Recently, SL amplitudes of ~ 0.2 G in 4 T (Gröhn et al, 2005) and 0.12 G in 1.5 T (Borthakur et al, 2004) have been used in human head MRI. With the higher SL power it is difficult to reach full-brain coverage required for clinical MRI. Evidently, the HSn-type method produces lower SAR values than the equivalent CW-type method. We calculate that compared with the CW-T1ρ method the SAR was lower by 81%, 37%, and 21% when scans were acquired with HS1, HS4, or HS8 pulses, respectively, with the same Ω₁^max. To our dissatisfaction, HSn pulses with high numbers are needed to produce MR contrast able to detect acute ischemia, partially counteracting the SAR advantage. In practice, the difference between SARs in the CW and HSn techniques is slightly smaller, when longer relaxation times with the HSn technique compared with the CW approach are accounted for, as the HSn technique requires longer SL times to obtain comparable contrast.

The MRI contrast detected between ischemic and normal tissue with the most sensitive rotating frame relaxation approach at 1 to 3 h after onset of stroke was comparable in size to the decrease in diffusion coefficient. It should be noted that this is typically the earliest time point when stroke patients are examined with MRI and often shows significant diffusion-perfusion mismatch (Baird et al, 1997; Meng et al, 2004; Shen et al, 2003) as a sign of tissue at risk of infarction. It has been shown both in animal models (Gröhn et al, 1999; Minematsu et al, 1992) and in human stroke patients (Doege et al, 2000) that the original diffusion abnormality can transiently recover, notwithstanding the evolving lesion. Also, T_1ρ measured during the early hours of ischemia has been shown to better correlate with cell death (Gröhn et al, 1999). Furthermore, T_1ρ progressively increases during the first hours of stroke, whereas diffusion remains at a constant low value, showing that these two MR contrasts probe different aspects of ischemia. These known features of diffusion and rotating frame relaxation data were supplemented by results from threshold analysis in the present study: decreased area of diffusion abnormality was seen during reperfusion whereas many of the rotating frame relaxation contrasts used as well asT₂ were able to highlight progression of the lesion. Data from this and other experimental stroke studies (Gröhn et al, 1999, 2000a; Makela et al, 2002) suggest that rotating frame relaxation contrasts could be exploited to estimate the age and severity of the ischemia in the brain parenchyma, and thus, may provide complementary information to the current clinically used MRI techniques. In this way, rotating frame relaxation contrasts can potentially facilitate attempts to predict outcome after brain ischemia (Shen et al, 2005; Wu et al, 2007).

The 2SX model was used to extract information from different water pools during evolution of acute ischemic stroke. It should be realized that such models created to describe water dynamics in tissue are often oversimplified and thus, caution should be exercised when interpreting the results in terms of true cerebral water pools. Two-pool model has often been used in the context of magnetization transfer (Henkelman et al, 1993; McConnell, 1958). In this model, the immobile fraction is named ‘semisolid pool’ and it includes protons in macromolecules within tissue that have aT₂ relaxation time in the order of 10 ms (Henkelman et al, 1993) and a typical t_c in the range of 10⁻³ to 10⁻⁶ secs. In the present study the more immobile fraction is named ‘water associated with macromolecules' and has an effective t_c in the order of 10⁻⁹ secs indicating much higher mobility. Accordingly, the pool size for semisolid pool in MT-based models is typically between 3% and 15% (Makela et al, 2002; Portnoy and Stanisz, 2007; Sled and Pike, 2001) whereas we obtained 34% to 38% for this pool. We initially tried to constrain the pool sizes in the 2SX model fitting to be more in line with MT literature, but this resulted in a physically impossible (negative) parameter estimate for the exchange correlation time. Therefore, it is evident that these two models describe behavior of different water pools in vivo and results obtained by modeling of Z-spectra data from MT experiments cannot be directly compared with those of the present study.

The results obtained from the 2SX model agree with pathophysiology of ischemia, justifying discussion concerning the mechanisms contributing to these early relaxation changes. Recently, acute T₁ relaxation changes were linked to increased water content (increase by 2.3% by 2 h of ischemia) in tissue during acute cerebral ischemia (Barbier et al, 2005). Similarly, an increase in the amount of MT has been associated with elevated water content in ischemic tissue (Makela et al, 2002). Interestingly, we observed an increase by some 2% in relative size of free water pool about 2 h after MCAo. It should be stressed, however, that our results indicate that not only water content, but also other factors may contribute to the prolongation of rotating frame relaxation times in ischemic tissue. Our results show a very pronounced decrease (38%) of the correlation time of the water pool associated with macromolecules and an increase (11%) of the exchange correlation time. It is well known that the collapse of cerebral energy state leads to complex biochemical cascades involving the activation of proteases and other destructive processes. It seems reasonable to speculate that decreased correlation time is associated with the accumulation of protein fragments that arise during proteolysis (Ewing et al, 1999). The increase in exchange time may be caused by the concerted actions of temperature decline and pH decline. Proton exchange rate decreases in acute ischemia for the majority of exchangeable groups in proteins when pH drops from 7.1 to close to 6 (Jokivarsi et al, 2007; Liepinsh and Otting, 1996). However, it should be noted that the exchange term in the present model includes not only chemical exchange but also a broad variety of different exchange phenomena. Furthermore, a fixed value for δΩ was used in the fitting. However, it cannot be ruled out that a possible change in δΩ during acute ischemia could have contributed to the fitting results.

To summarize, T_1ρ measured using either a CW on-resonance pulse or HSn pulses with high n proved to be the most sensitive for acute ischemic stroke. The relaxation changes are likely to be associated with increased tissue water content, an increase in the water pool associated with macromolecules, and a change in exchange behavior of water molecules. Sensitivity of T_1ρ MRI techniques is comparable to diffusion MRI at clinically relevant time points after onset of ischemia and the combined use of these MRI techniques may facilitate evaluation of the time of the onset of stroke and reversibility of the tissue damage (Shen et al, 2005). Recently, T_1p and T_2p MRI techniques have shown to bear great potential for diagnostic assessment of chronic neurodegenerative disorders, such as Parkinson's disease (Michaeli et al, 2007) and Alzheimer's disease (Borthakur et al, 2006).

Footnotes

Acknowledgements

We thank Ms Maarit Pulkkinen for expert technical assistance and Mr Nick Hayward for proofreading and editing the paper.

Appendix

During an adiabatic pulse the orientation of the effective frequency vector Ω_eff(t) changes, as described by the angle α(t) between Ω_eff(t) and the Z‘ axis,

where Ω₁(t) is the amplitude modulation function of the RF pulse, Ω_RF(t) the time-dependent frequency of the RF pulse, and Ω₀ the Larmor precession frequency. The time-dependent magnitude of the effective frequency vector is given by

In accordance with the 2SX model, with two water populations (P_A and P_B, where P_A + P_B = 1) coupled by equilibrium exchange, the rate constant describing longitudinal rotating frame relaxation (R_1p,obs = (T_1ρ,obs)⁻¹) during an adiabatic pulse is given by (Michaeli et al, 2004; Sorce et al, 2006)

where T_P is the length of individual pulses. Contributions from dipolar fluctuations (first and second terms of Equation (4)) can be calculated from

where

and exchange-induced relaxation (third term in Equation (4)) is given by

For transverse rotating frame relaxation, the rate constant (R_1ρ,obs = 1/T_2ρ,obs) can be written as (Michaeli et al, 2004)

where

and the exchange relaxation rate constant is given by:

Here δΩ is the chemical shift difference between sites A and B, τ_c the correlation time for sites A and B(t_c,A and τ_c,B, respectively), and τ_ex the exchange correlation time (the mean lifetime of the exchanging species at the two sites), r is the internuclear distance (1.58 Å for water), I and γ are the equivalent spin and gyromagnetic ratio of the nucleus, respectively, and h is Planck's constant (Blicharski, 1972; Michaeli et al, 2005).

For the spin-lock steady-state relaxation (or instantaneous approximation for adiabatic rotation), the time dependence from the above equations can be ignored. Furthermore, for the on-resonance case (i.e., α = π/2), the R_1ρ,obs and R_2ρ,obs can be written as

Where

and the exchange-induced relaxation rate constants are given by

References

Abragam

(1961) The Principles of Nuclear Magnetism. Oxford: Clarendon Press.

Baird

Benfield

Schlaug

Siewert

Lovblad

K-O

Edelman

Warach

(1997) Enlargement of human cerebral ischemic lesion volumes measured by diffusion-weighted magnetic resonance imaging. Ann Neurol 41:581–9.

Barbier

Liu

Grillon

Payen

Lebas

Segebarth

Remy

(2005) Focal brain ischemia in rat: acute changes in brain tissue T1 reflect acute increase in brain tissue water content. NMR Biomed 18:499–506.

Blicharski

(1972) Nuclear magnetic relaxation in rotating frame. Acta Phys Pol A 41:223–36.

Borthakur

Gur

Wheaton

Corbo

Trojanowski

Lee

Reddy

(2006) In vivo measurement of plaque burden in a mouse model of Alzheimer's disease. J Magn Reson Imaging 24:1011–7.

Borthakur

Wheaton

Gougoutas

Akella

Regatte

Charagundla

Reddy

(2004) In vivo measurement of T1rho dispersion in the human brain at 1.5 Tesla. J Magn Reson Imaging 19:403–9.

Calamante

Lythgoe

Pell

Thomas

King

Busza

Sotak

Williams

Ordidge

Gadian

(1999) Early changes in water diffusion, perfusion, T1, and T2 during focal cerebral ischemia in the rat studied at 8.5 T. Magn Reson Med 41:479–85.

Crockard

Gadian

Frackowiak

RSJ

Proctor

Allen

Williams

Russell

RWR

(1987) Acute cerebral ischaemia: concurrent changes in cerebral blood flow, energy metabolites, pH, and lactate measured with hydrogen clearance and 31P and 1H nuclear magnetic resonance spectroscopy. II. Changes during ischaemia. J Cereb Blood Flow Metab 7:394–402.

Doege

Kerskens

Romero

Brunecker

Junge-Hulsing

Muller

Villringer

(2000) MRI of small human stroke shows reversible diffusion changes in subcortical gray matter. Neuroreport 11:2021–4.

10.

Ewing

Jiang

Boska

Zhang

Brown

Divine

Chopp

(1999) T1 and magnetization transfer at 7 Tesla in acute ischemic infarct in the rat. Magn Reson Med 41:696–705.

11.

Garwood

DelaBarre

(2001) The return of the frequency sweep: designing adiabatic pulses for contemporary NMR. J Magn Reson 153:155–77.

12.

Gröhn

Michaeli

Garwood

Kauppinen

Grohn

(2005) Quantitative T(1rho) and adiabatic Carr-Purcell T2 magnetic resonance imaging of human occipital lobe at 4 T. Magn Reson Med 54:14–9.

13.

Gröhn

Mäkelä

Lukkarinen

DelaBarre

Lin

Garwood

Kauppinen

(2003) On- and off-resonance T(1rho) MRI in acute cerebral ischemia of the rat. Magn Reson Med 49:172–6.

14.

Gröhn

Lukkarinen

Oja

van Zijl

Ulatowski

Traystman

Kauppinen

(1998) Noninvasive detection of cerebral hypoperfusion and reversible ischemia from reductions in the magnetic resonance imaging relaxation time, T2. J Cereb Blood Flow Metab 18:911–20.

15.

Gröhn

Lukkarinen

Silvennoinen

Pitkanen

van Zijl

Kauppinen

(1999) Quantitative magnetic resonance imaging assessment of cerebral ischemia in rat using on-resonance T1 in the rotating frame. Magn Reson Med 42:268–76.

16.

Gröhn

OHJ

Kettunen

Makela

Penttonen

Pitkanen

Lukkarinen

Kauppinen

RAM

(2000a) Early detection of irreversible cerebral ischemia in the rat using dispersion of the MRI relaxation time, T1ρ. J Cereb Blood Flow Metab 20:1457–66.

17.

Gröhn

OHJ

Kettunen

Penttonen

Oja

JME

van Zijl

PCM

Kauppinen

(2000b) Graded reduction of cerebral blood flow in rat as detected by the nuclear magnetic resonance relaxation time T2: a theoretical and experimental approach. J Cereb Blood Flow Metab 20:316–26.

18.

Henkelman

Huang

Xiang

Q-S

Stanisz

Swanson

Bronskill

(1993) Quantitative interpretation of magnetization transfer. Magn Reson Med 29:759–66.

19.

Hossmann

(1994) Viability thresholds and the penumbra of focal ischemia. Ann Neurol 36:557–65.

20.

Jokivarsi

Grohn

Kauppinen

(2007) Proton transfer ratio, lactate, and intracellular pH in acute cerebral ischemia. Magn Reson Med 57:647–53.

21.

Kettunen

Gröhn

OHJ

Lukkarinen

Vainio

Silvennoinen

Kauppinen

(2000) Interrelations of T1 and diffusion of water in acute cerebral ischemia of the rat. Magn Reson Med 44:833–9.

22.

Kettunen

Kauppinen

Gröhn

OHJ

(2005) Dispersion of cerebral on-resonance T1 in the rotating frame (T1 rho) in global ischaemia. Appl Magn Reson 29:89–106.

23.

Levitt

Freeman

Frenkiel

(1982) Broadband heteronuclear decoupling. J Magn Reson 47:328–30.

24.

Liepinsh

Otting

(1996) Proton exchange rates from amino acid side chains-implications for image contrast. Magn Reson Med 35:30–42.

25.

Longa

Weinstein

Carlson

Cummins

(1989) Reversible middle cerebral artery occlusion without craniectomy in rats. Stroke 20:84–91.

26.

Lythgoe

Busza

Calamante

Sotak

King

Bingham

Williams

Gadian

(1997) Effects of diffusion anisotropy on lesion delineation in a rat model of cerebral ischemia. Magn Reson Med 38:662–8.

27.

Makela

Kettunen

Gröhn

Kauppinen

(2002) Quantitative T(1rho) and magnetization transfer magnetic resonance imaging of acute cerebral ischemia in the rat. J Cereb Blood Flow Metab 22:547–58.

28.

McConnell

(1958) Reaction rates by nuclear magnetic resonance. J Chem Phys 28:430–1.

29.

Meng

Fisher

Shen

Sotak

Duong

(2004) Characterizing the diffusion/perfusion mismatch in experimental focal cerebral ischemia. Ann Neurol 55:207–12.

30.

Michaeli

Grohn

Sorce

Kauppinen

Springer

Jr Ugurbil

Garwood

(2005) Exchange-influenced T2rho contrast in human brain images measured with adiabatic radio frequency pulses. Magn Reson Med 53:823–9.

31.

Michaeli

Sorce

Garwood

Ugurbil

Majestic

Tuite

(2007) Assessment of brain iron and neuronal integrity in patients with Parkinson's disease using novel MRI contrasts. Mov Disord 22: 334–40.

32.

Michaeli

Sorce

Idiyatullin

Ugurbil

Garwood

(2004) Transverse relaxation in the rotating frame induced by chemical exchange. J Magn Reson 169:293–9.

33.

Michaeli

Sorce

Springer

Jr Ugurbil

Garwood

(2006) T1rho MRI contrast in the human brain: modulation of the longitudinal rotating frame relaxation shutter-speed during an adiabatic RF pulse. J Magn Reson 181:135–47.

34.

Minematsu

Sotak

Davis

Fisher

(1992) Reversible focal ischemic injury demonstrated by diffusion-weighted magnetic resonance imaging in rats. Stroke 23:1304–11.

35.

Mori

van Zijl

PCM

(1995) Diffusion weighting by the trace of the diffusion tensor within a single scan. Magn Reson Med 33:41–52.

36.

Moseley

Cohen

Mintorovitch

Chileuitt

Shimizu

Kucharczyk

Wendland

Weinstein

(1990) Early detection of regional cerebral ischemia in cats: comparison of diffusion- and T2-weighted MRI and spectroscopy. Magn Reson Med 14:330–46.

37.

Mullins

Schaefer

Sorensen

Halpern

Koroshetz

Gonzalez

(2002) CT and conventional and diffusion-weighted MR imaging in acute stroke: study in 691 patients at presentation to the emergency department. Radiology 224:353–60.

38.

Naritomi

Sasaki

Kanashiro

Kitani

Sawada

(1988) Flow thresholds for cerebral energy disturbance and Na+ pump failure as studied by 31P and 23Na NMR spectroscopy. J Cereb Blood Flow Metab 8:16–23.

39.

Paxinos

Watson

(1998) The rat brain in stereotaxic coordinates. London: Academic Press.

40.

Portnoy

Stanisz

(2007) Modeling pulsed magnetization transfer. Magn Reson Med 58:144–55.

41.

Shen

Meng

Fisher

Sotak

Duong

(2003) Pixel-by-pixel spatiotemporal progression of focal ischemia derived using quantitative perfusion and diffusion imaging. J Cereb Blood Flow Metab 23: 1479–88.

42.

Shen

Ren

Fisher

Duong

(2005) Statistical prediction of tissue fate in acute ischemic brain injury. J Cereb Blood Flow Metab 25:1336–45.

43.

Sierra

Michaeli

Niskanen

J-P

Valonen

Grohn

Ylä-Herttuala

Garwood

Grohn

(2008) Water spin dynamics during apoptotic cell death in glioma gene therapy probed by T1rho and T2rho. Magn Reson Med 59:1311–9.

44.

Silver

Joseph

Hoult

(1984) Highly selective π/2 and π pulse generation. J Magn Reson 59:347–51.

45.

Sled

Pike

(2001) Quantitative imaging of magnetization transfer exchange and relaxation properties in vivo using MRI. Magn Reson Med 46:923–31.

46.

Solomon

(1959) Rotary spin echoes. Phys Rev Lett 2: 301–2.

47.

Sorce

Michaeli

Garwood

(2006) The time-dependence of exchange-induced relaxation during modulated radio frequency pulses. J Magn Reson 179:136–9.

48.

Tannuüs

Garwood

(1996) Improved performance of frequency-swept pulses using offset-independent adiabaticity. J Magn Reson A 120:133–7.

49.

Wheaton

Borthakur

Corbo

Charagundla

Reddy

(2004) Method for reduced SAR T1rho- weighted MRI. Magn Reson Med 51:1096–102.

50.

Sumii

Asahi

Sasamata

Ostergaard

Rosen

Dijkhuizen

(2007) Infarct prediction and treatment assessment with MRI-based algorithms in experimental stroke models. J Cereb Blood Flow Metab 27:196–204.