A new non-invasive graphical method for quantification of cerebral blood flow with[[Formula: see text]] IMP.

OBJECTIVE: [[Formula: see text]] N-isopropyl-p-iodoamphetamine ([Formula: see text]-IMP) is an ideal perfusion tracer for single photon emission computed tomography, which shows good linearity between cerebral blood flow (CBF) and accumulation. However, quantification of CBF using [Formula: see text]-IMP without arterial blood sampling has been challenging, with previous methods requiring empirically obtained regression formulae to estimate CBF. Furthermore, the CBF value obtained via some of the previous methods would be affected by the clearance rate of [Formula: see text]-IMP from the lungs. This paper introduces a new non-invasive quantification method for CBF using [Formula: see text]-IMP and dynamic planar images.
METHODS: We have developed a theory based on Microsphere model. This method does not involve regression formulae for estimation and allows for direct measurement of CBF, considering the clearance rate of [Formula: see text]-IMP from the lungs. The study method is executed as easily as conventional Graph-Plot method. We compared the CBF values obtained by our study method and the established autoradiograph (ARG) method.
RESULTS: CBF values obtained using the new method demonstrated significant correlation with values determined using ARG method.
CONCLUSIONS: The novel method described presents a reliable and more simple way of determining CBF when compared to current methods.

Introduction

For quantification of cerebral blood flow (CBF) with single photon emission computed tomography (SPECT), [] N-isopropyl-p-iodoamphetamine (-IMP) is an ideal tracer as it has an advantage over other brain perfusion SPECT tracers, such as Tc-hexamethylpropyleneamine oxime (Tc-HMPAO) and Tc-ethyl cysteinate dimer (Tc-ECD), in that -IMP does not suffer underestimation of regional cerebral blood flow (rCBF) at high flow values [1]. -IMP shows retention in brain tissue by nonspecific, high-capacity, low-affinity amine binding sites within the brain after the initial uptake in the lung endothelium [2, 3].

The Microsphere model [4-6] is the simplest model for CBF estimation. However, the model assumes that the first-pass extraction (E) of tracer is high and efflux from the brain tissue is negligible. Fortunately, the E of -IMP is extremely high [2] and its volume of distribution (=) is around 30 [7]. As -IMP concentration in brain tissue () is limited during the early period, efflux from the brain tissue () is negligible. Thus, -IMP fulfills the assumptions of the Microsphere method if limited to the early phase of examination. However, the Microsphere method requires continuous or frequent arterial blood sampling to obtain input function. The Table Look Up method [7] and its simplification, the Autoradiograph (ARG) method [8] based on the one-tissue two-compartment model were developed to reduce the number of arterial blood samples to only one using standard input function, which was generated by combining the input functions from 12 subjects. Currently, the ARG method is the standard CBF quantification method using -IMP.

The Gjedde–Patlak–Matsuda method, a non-invasive application of the Gjedde–Patlak plot [9-11] with Tc-HMPAO [12, 13] or Tc-ECD [14], which considers tracer activity in the aortic arch with planar acquisition as an input function, inspired investigation of a non-invasive CBF quantification method with -IMP. The simple application of Gjedde–Patlak–Matsuda method to -IMP was not successful owing to scatter radiation from the lungs into the region in aortic arch. The non-invasive microsphere (NIMS) method [15, 16] utilized pulmonary clearance of the tracer for the estimation of input function. However, this method was not pursued due to its requirement for complicated procedures such as measurement of right ventricular cardiac output.

The Graph-Plot (GP) method [17] is relatively simple to carry out. However, it contains several critical methodological flaws. The GP method is also an application of Gjedde–Patlak plot, which regards activity in the pulmonary artery trunk as an input function, instead of activity in the aortic arch as with the Gjedde–Patlak–Matsuda method. Input function for brain response function must be the activity in the artery of systemic circulation. Thus, use of pulmonary artery measurement as an input function is not appropriate. Though the cause has yet to be elucidated, the lung clearance of -IMP showed great variation among patients [18]. Thus, the CBF value obtained via the GP method would be affected by the clearance rate of -IMP from the lungs. For example, slow clearance from the lungs would reduce tracer activity in the artery (), as it would in the brain (), but would not affect tracer activity in the pulmonary artery () (Fig. 1).

Fig. 1

Schema of time–activity curve of pulmonary artery , lung L(t) and brain . Black line is , blue line is L(t) and red line is . light blue line is L(t) and pink line is when -IMP washout from lungs is slow. Note that even if the washout from lungs is slow, does not change

Furthermore, the GP method requires the application of a regression formula obtained via the ARG method to calculate an absolute value. Similarly, the Gjedde–Patlak–Matsuda method also involves a regression formula determined by the Kanno–Lassen method [19]. The simple non-invasive -IMP quantification (SIMS) method that reduces the burden of NIMS also requires a regression formula to translate their SIMS index to input function from arterial blood sampling [20].

Many non-invasive methods for the calculation of CBF make use of regression formulae for estimation (Table 1). However, regression formulae obtained empirically generally amplify errors and reduce reliability of the system, because they contain dispersion and sometimes the formulae themselves are not appropriate. For example, the linear regression equation between HbA1c and average glucose by A1c-Derived Average Glucose (ADAG) study, [21], is empirical and contains dispersion which is mainly due to between-patient variation in derived mean RBC age. We have just proposed a theoretical approach including mean RBC age, [22]. Thus, equations obtained empirically using regression analysis are less reliable, and there is no evidence that the true relationship is linear.

Table 1

Summary of quantification methods using -IMP

Method	Invasion	Empirical regression formulae	Model	Remarks
Kameyama (2018)	None	None	Microsphere	This study
Microsphere (1982) [4]	Continuous ABS	None	Microsphere
Yokoi plot (1993) [23]	Continuous ABS	None	1 tissue	Dynamic SPECT
Table lookup (1994) [7]	1 ABS	None	1 tissue	2 × SPECT, standard IF
ARG (1994) [8]	1 ABS	None	1 tissue	Fixed DV, standard IF
NIMS (1997) [15]	None	Cardiac output	Microsphere	Too complicated
GP (2002) [17]	None	CBF from ARG	Gjedde–Patlak	IF: pulmonary artery
SIMS (2016) [20]	None	Continuous ABS	Microsphere	Simplified NIMS
SIARG (2018) [24]	None	Continuous ABS	1 tissue	SIMS derivative

ABS arterial blood sampling, IF input function, DV distribution volume ()

Therefore, a new, non-invasive CBF measurement method using -IMP but not requiring regression analysis has been long sought for, a problem solved by the appropriate application of mathematics in the method proposed in this article.

Materials and methods

Mathematical theory

The Microsphere method assumes that the brain absorbs all tracer from the blood, so that first-pass extraction (E) is 1 and , the rate constant which represents tracer movement from brain tissue to arterial plasma, can be ignored.where (Bq/g) is -IMP radioactivity in the brain, (Bq/mL) is concentration of -IMP in arterial blood, F (mL/g/min) is CBF, and is the lipophilic fraction of -IMP, considering the metabolism of -IMP. is assumed to be 0.8 according to the measurement of octanol extraction [7].

Fick’s principle [25, 26] yields the following equation:where L(t) (Bq) is the total amount of -IMP in the lungs, P (mL/min) is cardiac output from right ventricle, and (Bq/mL) is concentration of -IMP in the pulmonary artery.

Substituting Eq. (2) into Eq. (1) yieldsEquation (3) divided by yieldsWhen is plotted on the y-axis and is plotted on the x-axis, the slope of the line will represent F.

Equations (5) and (6) illustrate that imaging of a portion of the lungs is sufficient to calculate F.where l(t) (Bq) denotes the amount of the tracer in the visible portion of the lungs.

Data acquisition

Fifty-nine patients at Sapporo Shuyukai Hospital, Japan, underwent cerebral blood flow SPECT for clinical reasons in September 2013 and from January 2015 to October 2015 and CBF quantification was undertaken using both ARG and GP methods. GP method was executed by a software package (AZE VirtualPlace Hayabusa®, AZE, Tokyo, Japan). ARG method was executed by software on a SPECT console.

The data were retrospectively analyzed. The committee in the hospital approved the use of de-identified data from these patients including dynamic planar images and the results of ARG and GP methods. The hospital provided the data under the contract between the hospital and Nihon Medi-Physics. For this type of retrospective study, written consents were waived. Dynamic front planar scans (2 s/frame, 60 frames, 4.4181 mm/pixel) were acquired by a camera (GE Infinia 3, GE Healthcare, Waukesha, WI, USA) with an ELEGP collimator from the point of intra-venous bolus administration of 222 MBq -IMP. Arterial blood sampling occurred 10 min after the -IMP injection. SPECT images were acquired for 20 min commencing 20 min after the injection.

Data analysis

Data were analyzed on a standard spreadsheet software (Excel® 2016, Microsoft Corporation, Redmond, WA, USA). Four regions of interest (ROIs) were selected on planar images. One on the pulmonary artery trunk, one on the brain and one on both lungs (Fig. 2). The GP method also utilizes ROIs in the brain and pulmonary artery trunk [17]. Time–activity curves (TACs) were derived from these four ROIs (Fig. 3).

Fig. 2

An example of the regions of interest; pulmonary artery trunk (a red), brain (b yellow) and lungs (b blue)

Fig. 3

An example of time–activity curves of the three regions of interests; pulmonary artery trunk (), brain () and lungs (l(t))

We have obtained the mCBF value of ARG method from the same slice of CBF map as GP method and the study method.

The TACs calculated above were smoothed by low-pass filter to reduce statistical noise. Time zero-adjustment was then performed for time points where activity rose sharply [12, 17]. was plotted and the slope of the linear section was determined (Fig. 4).

Fig. 4

An example of plotting

Results

A significant correlation was observed between the CBF values obtained by the ARG method and the study method (, ). Similarly, a significant correlation was found between figures derived via the ARG and GP methods (, ) (Fig. 5).

Fig. 5

Relationship between ARG method and the study method/GP method. a Scatter plot of ARG and the study methods. b Scatter plot of ARG and GP methods. Green lines are identity lines

However, CBF values determined by the GP method were generally higher than those calculated by the ARG method.

Discussion

The aim of this study was to develop and determine the validity of a new non-invasive method for CBF quantification using -IMP that would remove the need for arterial blood sampling and regression formulae for estimation. CBF values obtained using the novel method showed significant correlation with the ARG method, suggesting that this method would be of practical use.

We can estimate regional CBF values by distributing the acquired mCBF value according to mean brain counts and constructed perfusion images, like Gjedde–Patlak–Matsuda method [13] and the GP method [27].

Study methods and NIMS

The novel method developed has many similarities with the NIMS method. However, modification of the original NIMS equation has allowed measurement of CBF without the need for measuring cardiac output from right ventricle (P) nor activity within entire lungs (L(t)).

Equations (7), (8) and (9) show the mathematical similarity between the study method and NIMS. Accumulation of the entire tracer dose in the lungs will occur within several seconds of injection ( (min)).where D (Bq) is the injected dose and (Bq) is the peak L(t). Due to clearance from the lungs, should be smaller than D. Note that . Equation (7) is substituted into Eqs. (2) and (4)The NIMS equation is as follows:The study Eq. (8) and the NIMS Eq. (9) differ only in that D is substituted for . It has been reported that [right side of Eq. (9)] consistently underestimates arterial blood sampling data [left side of the Eq. (9)] [16]. This is likely due to the use of instead of D. SIMS method [20] and SIARG method [24] suffer the same issue.

Study method and the GP method

While the study equation and the GP equation appear similar, the underlying principles differ greatly.where is the index value used in the GP equation to calculate CBF with the aid of a regression formula. The GP method is, in fact, an application of the Gjedde–Patlak plot. This method is flawed in that it does not consider tracer pharmacodynamics within the lungs and K is not a constant. Substituting for is not appropriate. Certainly we used , but it is not an input function in the study method.

Furthermore, the use of a regression formula in itself can be a source of error. CBF as determined by the GP method was higher than the values obtained via ARG, and this was especially so in instances of low flow, a result of regression analysis (Fig. 5). Also, the regression formula for each camera needs to be determined separately as one regression formula cannot be applied to all. Despite these flaws, figures calculated using the GP equation were significantly correlated with values calculated using the more reliable ARG method. It is possible that, despite previous reports [18], the variability in clearance of -IMP from the lungs was not great enough to influence the outcome of the equation. It should be noted that the ARG method also suffers over-simplification as it uses standard input function, which neglects variability in input functions.

The equations below aim to elucidate the relationship between and F. The overall tracer uptake rate constant for the system has been shown to be [11]: Here, Eqs. (7) and (11) yieldEquations (7) and (2) and yieldEquation (14) makes Eq. (12) as follows:These Eqs. (13) and (15) yieldwhere is the fraction of tracer in the lungs. This shows mathematically that the GP equation overestimates input function due to its lack of consideration of tracer retained in the lungs, leading to underestimation of . However, as mentioned earlier, if the fraction of tracer retained in the lungs remains largely constant, the GP method will provide valid results.

Limitations

The study method is based on the Microsphere model which assumes first-pass extraction (E) is 1 and neglects to consider back-diffusion. Certainly, extraction of -IMP approaches 1 and of -IMP is small; however, E is not 1 and is not 0, thus underestimation of F is possible.

The authors are aware that measurements taken from the pulmonary artery trunk include scatter radiation from the lungs. While subtraction of the scatter radiation would yield more accurate measurements, a reliable technique to correct TAC in this way could not be devised. Therefore, overestimation of tracer activity in the pulmonary artery trunk due to scatter radiation artifact is likely to cause underestimation of F.

We treated lipophilic fraction, as a constant to make calculation simple, which is in fact a variable of time, reflecting the process of the metabolism of 123I-IMP. However, after quick conversion, lipophilic fraction seems relatively stable [7], and NIMS method also adopted a constant lipophilic fraction [16].

We assumed that lung tissue is everywhere the same (Eq. 5). We believe that the assumption would be fulfilled if the patient does not have an evident heterogeneous lung disease. Moreover, we have in fact confirmed the robustness of the study method by various types of lung ROIs (data not shown). However, we admit that patients with lung disease (e.g., a patient with her right lung damaged by tuberculosis) would not fulfill the assumption, although other methods that use standard input function, including ARG method cannot be applied to the patients with lung disease, either.

Finally, further study into the attenuation of planar images in the head versus the pulmonary artery trunk is required to correct for differences in attenuation.