Indirect Proton MR Imaging and Kinetic Analysis of 17O-Labeled Water Tracer in the Brain.

PURPOSE: The feasibility of steady-state sequences for 17O imaging was evaluated based on a kinetic analysis of the brain parenchyma and cerebrospinal fluid (CSF).
MATERIALS AND METHODS: The institutional review board approved this prospective study with written informed consent. Dynamic 2D or 3D steady-state sequences were performed in five and nine participants, respectively, with different parameters using a 3T scanner. During two consecutive dynamic scans, saline was intravenously administered for control purposes in the first scan, and 20% 17O-labeled water (1 mL/Kg) was administered in the second scan. Signal changes relative to the baseline were calculated, and kinetic analyses of the curves were conducted for all voxels. Region of interest analysis was performed in the brain parenchyma, choroid plexus, and CSF spaces.
RESULTS: Average signal drops were significantly larger in the 17O group than in the controls for most of the imaging parameters. Different kinetic parameters were observed between the brain parenchyma and CSF spaces. Average and maximum signal drops were significantly larger in the CSF spaces and choroid plexus than in the brain parenchyma. Bolus arrival, time to peak, and the first moment of dynamic curves of 17O in the CSF space were delayed compared to that in the brain parenchyma. Significant differences between the ventricle and subarachnoid space were also noted.
CONCLUSION: Steady-state sequences are feasible for indirect 17O imaging with reasonable temporal resolution; this result is potentially important for the analysis of water kinetics and aquaporin function for several disorders.

Introduction

The 17O isotope is the only stable isotope of oxygen that is capable of producing an MRI signal, with a nuclear magnetic moment γ of 5.7 MHz/Tesla and short T1/T2 values[1,2]. As the natural abundance of 17O is only 0.037%, enrichment of 17O is required for MR imaging; however, the enrichment process is expensive. Previous studies have utilized both direct and indirect MRI methods for 17O imaging. The direct method relies on the radiofrequency (RF) of 17O, and both human[3,4] and animal[5,6] studies have been reported; however, it requires specific hardware, such as an amplifier, exciter, and RF coils for 17O.

In contrast, the indirect method relies on the 17O-induced changes in the 1H signal. Therefore, standard clinical MRI machines can be used without any additional hardware. This method utilizes T2 shortening of 17O-induced proton signals based on scalar coupling between 17O and 1H2. For example, to obtain image contrasts using T2-weighted imaging, a long TE is required in fast spin echo sequences. However, the brain parenchyma becomes too dark with a long TE, making it difficult to determine the signal change. Therefore, indirect 17O imaging has been mostly conducted using T1-ρ or decoupling sequences[7-10], although such methods remain relatively uncommon, and the decoupling sequence requires a specific RF coil for 17O, similar to the direct method. Further, no human research has been conducted using the indirect method.

Steady-state sequences have been used in conjunction with long scan times for high-resolution cisternography based on their strong T2 contrast. We hypothesized that modifying imaging parameters would allow for dynamic steady-state indirect imaging of 17O with reasonable temporal resolution.

The purpose of the present study was to evaluate the feasibility of steady-state sequences for 17O imaging via kinetic analysis of the brain parenchyma and cerebrospinal fluid (CSF) following administration of a 17O water tracer.

Materials and Methods

Participants

After the Institutional Review Board of our institution approved our study protocol, we recruited male volunteers with no past history of neurological disease from December 2012 to February 2015. The inclusion criteria included an age between 20 and 65 years old, and a body weight below 60 kg. Fourteen male volunteers (mean age: 28.6 years; age range: 21–39 years) were eligible. Written informed consent was obtained from all participants prior to enrollment in our study.

Magnetic resonance imaging

MRI scans were conducted using a 3.0 Tesla scanner (MR750; GE Healthcare, Milwaukee, WI, USA) with an 8-channel head coil. After the scout scan, dynamic 2D or 3D scans were performed with fast imaging employing steady state acquisition (FIESTA) sequences with variable parameters (Table 1). Five and nine participants underwent 2D and 3D imaging, respectively.

Table 1

Imaging parameters

Imaging sequence			No.	FA	TR (ms)	TE (ms)	Matrix	NEX	No. of section	Section thickness (mm)	No. of dynamic phase	Temporal resolution (sec)	FOV (mm)	Phase FOV	Total scan time (sec)
2D			1	90	32.6	30.0	256 × 128	1	1	5	77	4.7	240 × 180	0.75	363
			2	90	42.6	40.0	256 × 128	1	1	5	54	6.8	240 × 180	0.75	369
			3	90	42.6	40.0	256 × 128	1	1	5	62	5.8	240 × 180	0.75	362
			4	90	42.3	40.0	256 × 128	2	1	5	28	13.4	240 × 180	0.75	376
			5	90	42.3	40.0	256 × 128	4	1	5	15	24.3	240 × 180	0.75	365
3D	Low	Short	6	5	2.8	1.4	128 × 128	1	12	5	110	3.3	240 × 180	0.75	362
	FA	TR/TE	7	20	2.4	0.9	128 × 128	1	28	5	74	4.8	240 × 180	0.75	354
			8	25	2.8	1.4	128 × 128	1	12	5	110	3.3	240 × 180	0.75	362
	High	Short	9	70	3.2	1.5	128 × 90	0.75	40	4	120	3.0	240 × 180	0.75	356
	FA	TR/TE	10	70	3.2	1.5	128 × 128	1	20	4	24	15.1	240 × 180	0.75	362
			11	90	11.1	5.4	512 × 128	1	4	14	20	19.0	240 × 192	0.8	379
		Long	12	90	21.3	20.0	128 × 128	1	4	10	23	16.0	320 × 160	0.5	367
		TR/TE	13	90	31.4	30.0	128 × 64	0.75	4	14	22	16.5	240 × 180	0.75	362
			14	90	51.3	50.0	128 × 64	1	4	14	11	35.0	300 × 180	0.6	385

FA, flip angle; NEX, number of excitation.

Single section acquisition with high flip angle (flip angle [FA], 90 degrees) and long repetition time (TR, 32.6 to 42.6 ms)/echo time (TE, 30 to 40 ms) were used for 2D imaging, with the section location at the centrum semiovale. For 3D imaging, low FA and short TR/TE, high FA and short TR/TE, and high FA and long TR/TE were respectively applied to three participants. The middle of the scan slab was set at the level of the basal ganglia. The total scan time was approximately 6 minutes.

Two consecutive dynamic scans were performed for each participant. In the first dynamic scan, saline (1 mL/kg weight) was intravenously administered 60 or 120 s after the start of scanning for control purposes via the antecubital vein at a rate of 3 mL/s. In the second dynamic scan, 17O-labeled water (enriched to 20%; Taiyo Nippon Sanso Corp, Tokyo, Japan) was administered at the same dose, timing, and rate. This water tracer was made with good manufacturing practice (GMP) standards for intravenous injection in human participants.

Data analysis

The acquired image data were post-processed with in-house software (Perfusion Mismatch Analyzer [PMA] Research Version 5.2.0.1; Acute Stroke Imaging Standardization Group, Tokyo, Japan). Relative signals to the average baseline were calculated for all voxels, and the dynamic curves were analyzed. To exclude the unstable signals, the first 15 s data was discarded and the baseline was defined between 15 s after the scan start and the injection point of saline or 17O-labeled water. Gaussian smoothing was applied for the temporal axis (the sigma of which was 10 s), and data interpolation was not used. No spatial smoothing was applied to avoid the contamination of signal changes from adjacent voxels.

The average signal drop (ASD) was calculated within the period from the injection start of saline or 17O-labeled water to the end of the scan. Other kinetic parameters such as bolus arrival time (BAT), time to peak (TTP), first moment (FM), and maximum signal drop (MSD) were also calculated voxel by voxel. BAT was defined as the time at which a signal drop larger than 5% of maximum occurred. TTP was defined as the time of bolus peak, while FM was defined as the first moment of bolus area (gravity center).

ROI analysis was performed using the same software (Fig. 1). Ten square ROIs were manually defined in the cerebral cortex (CC), basal ganglia/thalamus (BG/TM), white matter (WM), choroid plexus (CP), ventricles, and subarachnoid space (SAS), respectively. To avoid the contamination of signal changes between CC and SAS, and between CP and ventricle, small size of ROI (4×4 mm) was used. In addition, 10 ROIs were evenly placed and banding artifact of FIESTA sequence was carefully avoided. ASD, BAT, TTP, FM, and MSD were measured for each of these ROIs. However, these values were not measured in the BG/TM, CP, and ventricles of five, four, and two participants who underwent 2D imaging, respectively, because those structures were not included within the section.

Fig. 1

Regions of interest. Ten square ROIs (4 mm × 4 mm) are defined in the cerebral cortex (red), basal ganglia/thalamus (pink), white matter (yellow), choroid plexus (purple), ventricles (blue), and subarachnoid space (green).

Statistical analysis

The ASDs of the 17O-enriched water and saline control conditions were compared using paired t-tests for each location in each participant to assess whether a significant signal drop caused by T2 shortening of 17O-enriched water had occurred. For participants in whom significant differences in the ASD were noted, we then compared the kinetic parameters of 17O-enriched water in each location using an analysis of variance followed by multiple comparisons using the Holm-Sidak method. A P-value less than 0.05 was considered significant in all statistical tests.

Results

All participants completed the study protocol, and no adverse effects were noted following the intravenous administration of 17O-enriched water.

In 2D imaging, the ASD of the CC, CP, ventricles, and SAS were significantly larger in the 17O than control condition for all participants (Table 2). Four of five participants also exhibited significant differences in ASD of WM. In 3D imaging, there were no significant differences in ASD at any location for any participant of the low FA and short TR/TE value group. In the high FA and short TR/TE group, significant differences in ASD were noted in two, one, two, three, and three participants for the CC, WM, CP, ventricle, and SAS, respectively. No significant differences were observed in the BG/TM for any of the three participants. In the high FA and long TR/TE group, significant differences were observed in all locations in all participants.

Table 2

Difference of ASD between 17O and control

			No.	Cerebral cortex			Basal ganglia/thalamus			White matter			Choroid plexus			Ventricle			Subarachnoid space
				17O	Control	P value	17O	Control	P value	17O	Control	P value	17O	Control	P value	17O	Control	P value	17O	Control	P value
2D			1	1.9 ± 1.1	−0.3 ± 0.9	< 0.001	N.A.	N.A.	N.A.	0.4 ± 0.4	−0.3 ± 0.3	< 0.01	N.A.	N.A.	N.A.	N.A.	N.A.	N.A.	4.3 ± 1.6	−0.4 ± 0.6	< 0.001
			2	2.5 ± 1.5	−0.1 ± 1.1	< 0.001	N.A.	N.A.	N.A.	2.3 ± 2.0	−0.2 ± 0.7	< 0.05	N.A.	N.A.	N.A.	N.A.	N.A.	N.A.	8.1 ± 2.3	0.3 ± 1.0	< 0.001
			3	4.3 ± 3.6	0.3 ± 1.8	< 0.01	N.A.	N.A.	N.A.	0.6 ± 0.8	−0.8 ± 0.6	< 0.05	7.3 ± 2.8	0.6 ± 0.5	< 0.001	5.7 ± 2.0	0.1 ± 1.5	< 0.001	11.0 ± 3.0	−0.2 ± 1.5	< 0.001
			4	9.4 ± 3.0	1.1 ± 5.3	< 0.01	N.A.	N.A.	N.A.	0.6 ± 1.4	−1.6 ± 1.3	< 0.01	N.A.	N.A.	N.A.	2.6 ± 1.6	−0.6 ± 1.1	< 0.01	8.5 ± 2.8	−2.3 ± 4.1	< 0.001
			5	3.1 ± 1.7	−0.7 ± 1.1	< 0.001	N.A.	N.A.	N.A.	1.4 ± 0.7	0.3 ± 1.5	N.S.	N.A.	N.A.	N.A.	2.0 ± 1.8	−0.5 ± 1.2	< 0.01	5.5 ± 2.0	−0.3 ± 0.6	< 0.001
3D	Low	Short	6	0.5 ± 0.5	0.4 ± 0.4	N.S.	0.4 ± 0.2	0.4 ± 0.2	N.S.	0.4 ± 0.3	0.3 ± 0.3	N.S.	0.5 ± 0.2	0.3 ± 0.3	N.S.	0.4 ± 0.3	0.1 ± 0.4	N.S.	0.3 ± 0.4	0.2 ± 0.3	N.S.
	FA	TR/TE	7	0.8 ± 1.3	0.7 ± 0.5	N.S.	0.2 ± 0.6	0.5 ± 0.6	N.S.	0.2 ± 0.6	0.1 ± 0.7	N.S.	0.4 ± 0.9	1.0 ± 1.0	N.S.	1.3 ± 1.3	0.8 ± 1.0	N.S.	0.7 ± 1.0	0.7 ± 1.0	N.S.
			8	0.1 ± 0.5	−0.2 ± 0.4	N.S.	0.7 ± 0.5	0.7 ± 0.3	N.S.	0.3 ± 0.4	−0.2 ± 0.8	N.S.	0.1 ± 2.1	−0.4 ± 2.0	N.S.	−1.1 ± 3.5	0.8 ± 3.0	N.S.	0.6 ± 1.0	1.2 ± 0.9	N.S.
	High	Short	9	−0.5 ± 0.5	−0.3 ± 0.6	N.S.	−0.1 ± 2.5	−2.4 ± 3.1	N.S.	−0.6 ± 1.2	−0.8 ± 0.9	N.S.	0.1 ± 6.2	3.5 ± 6.6	N.S.	6.1 ± 5.7	−4.5 ± 6.1	< 0.05	4.1 ± 1.5	−0.6 ± 3.1	< 0.001
	FA	TR/TE	10	1.6 ± 1.6	−0.9 ± 0.6	< 0.001	1.8 ± 4.1	1.9 ± 3.0	N.S.	2.8 ± 2.4	0.1 ± 1.7	N.S.	3.7 ± 2.4	−0.1 ± 0.6	< 0.001	5.0 ± 1.5	−0.2 ± 0.7	< 0.001	4.6 ± 1.9	0.5 ± 0.5	< 0.001
			11	1.9 ± 1.2	−0.3 ± 0.5	< 0.001	1.9 ± 1.5	1.0 ± 2.1	N.S.	1.8 ± 1.0	−0.1 ± 1.2	< 0.01	8.1 ± 1.6	0.3 ± 1.3	< 0.001	6.6 ± 1.7	0.6 ± 1.4	< 0.001	8.2 ± 1.5	−0.8 ± 0.7	< 0.001
		Long	12	5.9 ± 2.0	0.0 ± 0.7	< 0.001	1.2 ± 0.9	−0.3 ± .9	< 0.01	1.2 ± 0.8	−0.3 ± 1.4	< 0.05	8.7 ± 2.3	0.1 ± 0.6	< 0.001	5.0 ± 1.1	0.0 ± 0.8	< 0.001	12.6 ± 2.6	−0.1 ± 1.0	< 0.001
		TR/TE	13	4.2 ± 1.2	−0.6 ± 0.7	< 0.001	1.9 ± 1.5	0.1 ± 2.5	< 0.05	1.3 ± 2.0	−1.2 ± 2.0	< 0.05	8.5 ± 3.7	0.5 ± 0.7	< 0.001	3.6 ± 1.3	0.7 ± 1.2	< 0.001	9.4 ± 1.2	−0.4 ± 0.9	< 0.001
			14	5.9 ± 4.6	−2.5 ± 3.7	< 0.001	7.4 ± 3.0	0.8 ± 3.6	< 0.001	5.0 ± 4.0	0.6 ± 3.8	< 0.05	9.8 ± 5.5	0.7 ± 1.8	< 0.001	6.4 ± 2.8	−1.7 ± 2.6	< 0.001	12.0 ± 1.8	−1.0 ± 2.2	< 0.001

N.A., not applicable; N.S., not significant.

As no significant differences were noted in any participant of the low FA and short TR/TE group, these three participants were excluded from the following analysis of kinetic parameters. Examples of parameter maps and dynamic curves in a single participant are shown in Figs. 2 and 3, respectively.

Fig. 2

Parameter maps. Examples of parameter maps (3D imaging with a high FA and a long TR/TE) are shown. In 17O maps, higher ASD and MSD values are observed in the subarachnoid space, choroid plexus, and cerebral cortex when compared to the saline control condition. Time-related parameters such as BAT, TTP, and FM are noisy; however, some anatomical structures are more easily visualized following 17O injection when compared to saline control conditions. BAT, bolus arrival time; TTP, time to peak; FM, first moment; MSD, maximum signal drop; ASD, average signal drop; MS, maximum slope; FA, flip angle.

Fig. 3

Dynamic curves of relative signals. Examples of dynamic curves in a single participant (same as Fig. 2) are shown. The bolus injection of saline or 17O is initiated 120 s after the scan start. Significant signal drops are observed in the cerebral cortex (a), choroid plexus (d), ventricle (e), and subarachnoid space (f) in the 17O condition, while slight decreases in the signals are noted in the basal ganglia/thalamus (b) and white matter (c). In the cerebral cortex, choroid plexus, and subarachnoid space, the curves are characterized by two phases, i.e., an initial signal drop phase and a plateau phase. In the ventricle, the signals gradually and continuously decreased during the acquisition time window.

The kinetic parameters of eleven participants are summarized in Fig. 4. Significantly larger BAT values were observed for the ventricles and SAS relative to other locations, indicating that bolus arrivals of 17O in the CSF space were delayed compared to the brain parenchyma. Similarly, TTP and FM of the ventricles and SAS were significantly larger than values in other locations. In addition, TTP values of the CP were significantly larger than those of the CC and BG/TM. WM exhibited slightly but significantly larger TTP values than the BG/TM. FM values of the BG/TM were significantly smaller than those in other locations, with the exception of the CC.

Fig. 4

Comparisons of kinetic parameters among locations. The kinetic parameters of 11 participants are summarized. Three participants are excluded due to no significant signal drops following injection of 17O. Significant differences (P < 0.05) are indicated with an asterisk (*) and a bold line as a reference. For example, significant differences in BAT are observed for the SAS when compared with the cerebral cortex, basal ganglia/thalamus, choroid plexus, and ventricles, although no such differences are noted when this region is compared with the white matter. BAT (a), TTP (b), and FM (c) of the ventricle and SAS are significantly larger than the brain parenchyma. Similarly, MSD (d) and ASD (e) of the choroid plexus, ventricle, and SAS are significantly larger than the brain parenchyma. MS (f) of the SAS is significantly larger than the ventricle, indicating rapid signal change. BG, basal ganglia; CC, cerebral cortex; CP, choroid plexus; SAS, subarachnoid space; TM, thalamus; WM, white matter.

MSD and ASD values in the SAS, ventricles, and CP were significantly larger than those in other locations. Among these, values for the SAS were significantly higher than those for the CP and ventricles, while values for the CP were significantly higher than those in the ventricles. ASD and MSD values were significantly smaller in the WM than in the CC. Significantly larger MSD values were noted in the SAS when compared to those of the ventricles and BG/TM. Furthermore, the ventricle exhibited significantly smaller MSD values than the CC.

Discussion

In the present study, we applied steady-state sequences for the indirect imaging of 17O. Steady-state imaging is usually applied for high-resolution hydrography; however, we reduced the spatial resolution to achieve a reasonable temporal resolution for the evaluation of dynamic signal changes following the administration of 17O water tracer. Although dependencies in scan parameters existed, significant 17O-induced signal drops were observed using both 2D and 3D steady-state imaging.

Steady-state sequences rely on T2/T1 contrast and are usually used with a short TR/TE to reduce the influence of banding artifacts[11]. Better temporal resolution (3.3–4.8 s) was obtained with short TR/TE in 3D imaging; however, no significant signal drops of 17O were observed. In contrast, significant signal drops in 17O were observed using a long TR/TE in conjunction with a high FA. Although temporal resolution became longer with a longer TR, this combination of parameters could be effectively used for indirect 17O imaging. Other sequences such as fast spin echo, echo planar imaging (EPI), and gradient echo are potential candidates for T2-weighted imaging of 17O. However, a longer TR/TE may be required to obtain T2 contrast of 17O, which may result in a significant signal drop before the administration of 17O in the brain parenchyma, a lower signal-to-noise ratio, and a decreased temporal resolution.

In the present study, signal changes (ASD and MSD) in CC were significantly larger than in WM. Although quantification has not been conducted, these differences in signal change may reflect the difference of cerebral blood flow and/or volume. In CT and MR perfusion study, Cerebral blood flow can be measured by maximum slope model; however, there was no significant differences in MS between CC and WM. This was probably because we used diffusible tracer which is different from intravascular tracer of CT perfusion (iodine contrast) and MR perfusion (Gd contrast). In addition, MS calculation might be susceptible to image noise as the signal change is too small to measure, while ASD and MSD were relatively robust to noise as these parameters were simply calculated as the average or maximum of signals.

Signal changes in the CSF were larger than brain parenchyma, likely because 17O is a water tracer, and difference in relaxation time between brain parenchyma and CSF. Such a result has not been observed using 15O positron emission tomography (PET), probably because of the larger signal change with 17O MRI in CSF, or due to the short half-life of 15O which leaks out to the CSF. In contrast, 17O is a stable isotope, allowing for the visualization of a relatively slower water leakage from vessels or brain parenchyma into the CSF space, such as the ventricles and SAS, compared to a faster blood flow. In the classic theory, the choroid plexus produces CSF; however, in our study, earlier and more rapid signal changes were observed in the SAS compared to the ventricle, indicating fluid leakage from vessels or the surface of the brain. This observation is compatible with the findings of previous studies that used 17O as a water tracer[12]; therefore, we conclude that 17O can be used as a tracer in the analysis of CSF kinetics.

Previous studies have indicated that aquaporin (AQP) may be associated with the pathogenesis of several diseases. For example, neuromyelitis optica (NMO) is associated with antibodies to AQP[13]. Additional studies have reported that Alzheimer’s disease and Parkinson’s disease are associated with changes in the expression of AQP[14-16]. Other brain disorders such as stroke, intracranial hematoma, and edema may also be related to AQP function[17-20]. In addition to the central nervous system, altered expression of AQP has been noted in several diseases in the lungs, kidneys, liver, and other organs[21-24]. Although the distribution of AQP can be histologically analyzed, the dynamic function of AQP in vivo cannot be evaluated without the use of a water tracer. Currently, no tracers have been useful for AQP kinetics analysis; however, our findings indicate that 17O may be one such candidate for use in MR imaging. The dysfunction of AQP channel might be imaged with 17O-water as decreased signal change of impaired lesion.

The present study has several limitations. First, the temporal and spatial resolution was not high enough, especially for a longer TR/TE, and a much higher resolution may be required for a detailed examination of water kinetics in small structures of the brain. Second, a longer TR/TE may produce strong banding artifacts. Therefore, improved homogeneity of the static magnetic field may be required. Third, contamination of signal change between CC and SAS, and between CP and ventricle might be occurred even we used small ROI. Although we carefully avoided the location of banding artifact of FIESTA sequence, several ROIs may suffer signal decrease derived from banding artifact Furthermore, although we observed significant changes in relative signals, we did not quantify the absolute amount of 17O. Further studies are required for such quantification, although our results have revealed significant differences in the kinetics of dynamic curves using relative signal changes.

Conclusion

Our findings indicate that steady-state sequences are feasible for indirect 17O imaging with reasonable temporal resolution; this result is potentially important for the analysis of water kinetics and AQP function for several disorders.