Distinction Between Cerebral Blood Flow and Vascular Activity by PET Tracer Kinetic Analysis: Practical and Physiologic Considerations ====================================================================================================================================== * Satoshi Minoshima * Robert A. Koeppe Quantitative measurement of cerebral blood flow using PET dates back more than 20 y. The emergence of PET as a clinical experimental tool permitted 3-dimensional measurement of radioactivity within the brain and applications of tracer kinetic analyses to estimate quantitatively various physiologic parameters of the brain. The estimation of cerebral blood flow was one of the earliest and most extensively researched applications by pioneers of PET imaging. Early investigations of cerebral blood flow measurements commonly used Kety’s (1) 1-compartment model, in which radioactive water was treated as a freely diffusible tracer between blood and brain. However, the diffusibility of water across the blood-brain barrier was questioned (2–5), and 2-compartment models including a vascular component subsequently were developed (6, 7). Ohta et al.8) from the Montreal Neurological Institute examined further influences of residual vascular activity in the estimation of cerebral blood flow and introduced a 2-compartment model that included an “apparent vascular distribution volume,” the parameters of which could be estimated efficiently using a weighted integration method (9). The study presented by Okazawa and Vafaee (10) used the 2-compartment model and reexamined influences of the vascular activity in the estimation of cerebral blood flow under basal conditions and visual stimulation using statistical parametric mapping. As has been reported (8), a conventional 1-compartment analysis resulted in overestimation (though not significant in this study) of true water clearance in brain regions with greater vascularity. The distribution of such overestimation is better appreciated in 3 dimensions by means of statistical parametric mapping (Fig. 4). The authors also confirmed the work of Koeppe et al. (5), which showed that eliminating initial dynamic frames in a 1-compartment autoradiographic approach could reduce overestimation inherent from vascular activity. However, even when optimal elimination of initial dynamic frames was conducted, there still remained overestimation of cerebral blood flow in comparison with estimation by the 2-compartment model (Table 1). These findings indicate that the use of the 2-compartment model in conjunction with dynamic PET imaging has an advantage in estimating absolute levels of cerebral blood flow with minimal influences from residual vascular activity because of the limited diffusion of water. Cerebral blood flow and vascular distribution volume were also examined within brain regions where transient hemodynamic changes were evoked by sensory stimulation. The authors used visual stimulation delivered by a reversing checkerboard pattern. When comparing 1-compartment and 2-compartment models, 2 important observations emerged, as was shown previously with vibrotactile stimulation (11). First, increased neuronal activity in the occipital cortex was associated with increased cerebral blood flow and increased vascular distribution volume, which could be demonstrated by the 2-compartment model. Because the estimated vascular distribution volume by the 2-compartment model is reported to reflect radioactivity in the arteries, arterioles, and a small fraction of the capillary (12), this finding indicates the presence of arterial hyperemia at activation foci. This is an interesting contrast to functional MRI (fMRI), which detects blood oxygenation level dependent signals primarily arising from draining venules (and parenchyma, particularly with a high magnetic filed scanner) near activation foci (13). Second, peak locations of changes in cerebral blood flow and vascular distribution volume revealed by the 2-compartment model were not identical. The previous (11) and current study (10) showed differences of several millimeters between peak changes of cerebral blood flow and vascular distribution volume. The peak changes of vascular distribution volume were located closer to arterial vascular trees (11). The peak change detected by a conventional 1-compartment model resides between these 2 peaks. Similarly, the nonquantitative 1-compartment autoradiographic approach commonly used in brain activation studies suffers from a mixture of changes in blood flow and arterial vascular volume, which makes comparisons between PET and fMRI signal localization in relation to local neuronal activity even more difficult (14). The use of the 2-compartment model improves spatial accuracy in identifying regional cerebral blood flow changes that are associated with evoked local neuronal activity in brain activation studies. Better quantification of cerebral blood flow and better localization of brain activation can be achieved by the proposed 2-compartment model. Should this technique routinely be used in cerebral blood flow measurements? For the purpose of better quantifying cerebral blood flow in certain brain disorders, such as strokes, the 2-compartment model minimizes potential overestimation of cerebral blood flow caused by vascular activities. This is advantageous if the primary goal of the investigation or clinical application is to measure quantitatively local cerebral blood flow. For brain activation studies, some practical and physiologic issues need to be addressed. One of the practical issues is the requirement of arterial blood sampling. Brain activation studies typically require repeat imaging with different types of stimulation or task paradigms. In addition to the somewhat invasive procedure of arterial catheter placement, the total amount of blood withdrawn during the study cannot be ignored. Although peripheral arterial sampling is nothing extraordinary to brain investigators who perform much more invasive and complex interventions to research subjects, frequent and repeat blood sampling introduces additional uncertainty and increases the potential for technical error. One assumption of brain activation studies is the physiologic coupling between local neuronal activity and blood flow. However, absolute levels of cerebral blood flow can be modulated by factors independent of local neuronal activities. For example, the plasma carbon dioxide level is a significant factor that influences global and local cerebral blood flow (15). In this regard, an absolute level of cerebral blood flow is not a sensitive or specific marker of neuronal activity (16), and proper data normalization procedures can identify regional cerebral blood flow changes associated with neuronal activities sensitively and semiquantitatively (17,18). Attempts to measure absolute levels of cerebral blood flow changes quantitatively have been largely abandoned in brain activation studies for this reason. Another factor that can potentially bias cerebral blood flow measurement that is often overlooked in brain activation studies is the assumption of a physiologic “steady state” of neural activities and cerebral blood flow changes evoked by task or stimulation. For example, when constant repetitive nociceptive stimuli are given to subjects, neuronal activities within the brain gradually can change within a minute, partly because of spatial and temporal summation (19). Neuronal activities relating to attention, anxiety, tolerance, and fatigue during a single imaging session can potentially fluctuate and influence local neuronal activities and patterns of activities. These changes can be difficult to appreciate unless they are specifically investigated. Changes in neuronal activities during image acquisition were shown to bias significantly the final estimation of cerebral blood flow in a conventional autoradiographic approach (20). In the kinetic models presented in this article, regional cerebral blood flow is assumed to be constant (Fig. 1). In addition, the dynamic imaging method requires longer image acquisition (3 min), compared with a generally shorter acquisition (1 min) for the nonquantitative autoradiographic method that is often used in conventional brain activation studies. The longer scan duration increases likelihood for violation of the physiologic steady state assumption. The magnitude of potential bias introduced by such violations can be estimated for the proposed 2-compartment model. This is an important tradeoff between physics (more imaging data, better accuracy) and physiology (long scan duration, more difficult to maintain physiologic steady state) in the optimization of imaging protocols. Finally, it is interesting to note that statistical significance in the detection of brain activation was consistently greater with the 1-compartment model compared with the cerebral blood flow maps generated by the 2-compartment model (10,11). There are 2 probable reasons for this observation. First, the additive effects of increases in cerebral blood flow and vascular volume associated with neuronal excitation produce greater signals at activation foci, despite small blurring effects caused by mismatching of peaks. Second, the precision of parameter estimation generally improves when fewer parameters are estimated. Thus, increased signals and smaller variances at activation foci probably contribute to the greater statistical significance. This result indicates that the widely used nonquantitative autoradiographic method with proper data normalization is a sensitive way to detect the presence of local neuronal activation evoked by task or stimulation, but with a small compromise in spatial accuracy. Advanced and complex image analysis improves our understanding of local brain physiology and the nature of signals measured by in vivo imaging, but it often comes back to the simplest form when the methodology matures for more general applications. ## Acknowledgments The authors thank Donna J. Cross for her assistance in preparing this manuscript. ## Footnotes * Received Mar. 9, 2001; revision accepted Mar. 9, 2001. For correspondence and reprints contact: Satoshi Minoshima, MD, PhD, Department of Radiology, University of Washington School of Medicine, Health Sciences Building, NW040J, 1959 N.E. Pacific St., Box 356004, Seattle, WA 98195-6004. ## References 1. Kety SS. Measurement of local blood flow by the exchange of an inert, diffusible substance. Methods Med Res 1960;8:228–236. 2. Eichling JO, Raichle ME, Grubb RL, Ter-Pogossian MM. Evidence of the limitations of water as a freely diffusible tracer in brain of the rhesus monkey. Circ Res 1974;35:358–364. [Abstract/FREE Full Text](http://jnm.snmjournals.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTA6ImNpcmNyZXNhaGEiO3M6NToicmVzaWQiO3M6ODoiMzUvMy8zNTgiO3M6NDoiYXRvbSI7czoyMjoiL2pudW1lZC80Mi83LzEwNDAuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 3. Gjedde A, Andersson J, Eklof B. Brain uptake of lactate, antipyrine, water and ethanol. Acta Physiol Scand 1975;93:145–149. [PubMed](http://jnm.snmjournals.org/lookup/external-ref?access_num=1146566&link_type=MED&atom=%2Fjnumed%2F42%2F7%2F1040.atom) 4. Herscovitch P, Markham J, Raichle ME. Brain blood flow measured with intravenous H215O. I. Theory and error analysis. J Nucl Med 1983;24:782–789. [Abstract/FREE Full Text](http://jnm.snmjournals.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Njoiam51bWVkIjtzOjU6InJlc2lkIjtzOjg6IjI0LzkvNzgyIjtzOjQ6ImF0b20iO3M6MjI6Ii9qbnVtZWQvNDIvNy8xMDQwLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 5. Koeppe RA, Hutchins GD, Rothley JM, Hichwa RD. Examination of assumptions for local cerebral blood flow studies in PET. J Nucl Med 1987;28:1695–1703. [Abstract/FREE Full Text](http://jnm.snmjournals.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Njoiam51bWVkIjtzOjU6InJlc2lkIjtzOjEwOiIyOC8xMS8xNjk1IjtzOjQ6ImF0b20iO3M6MjI6Ii9qbnVtZWQvNDIvNy8xMDQwLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 6. Gambhir SS, Huang SC, Hawkins RA, Phelps ME. A study of the single compartment tracer kinetic model for the measurement of local cerebral blood flow using 15O-water and positron emission tomography. J Cereb Blood Flow Metab 1987;7:13–20. [PubMed](http://jnm.snmjournals.org/lookup/external-ref?access_num=3492506&link_type=MED&atom=%2Fjnumed%2F42%2F7%2F1040.atom) [Web of Science](http://jnm.snmjournals.org/lookup/external-ref?access_num=A1987G010000003&link_type=ISI) 7. Larson KB, Markham J, Raichle ME. Tracer-kinetic models for measuring cerebral blood flow using externally detected radiotracers. J Cereb Blood Flow Metab 1987;7:443–463. [PubMed](http://jnm.snmjournals.org/lookup/external-ref?access_num=3611204&link_type=MED&atom=%2Fjnumed%2F42%2F7%2F1040.atom) 8. Ohta S, Meyer E, Fujita H, Reutens DC, Evans A, Gjedde A. Cerebral [15O]water clearance in humans determined by PET: I. Theory and normal values. J Cereb Blood Flow Metab 1996;16:765–780. [CrossRef](http://jnm.snmjournals.org/lookup/external-ref?access_num=10.1097/00004647-199609000-00002&link_type=DOI) [PubMed](http://jnm.snmjournals.org/lookup/external-ref?access_num=8784222&link_type=MED&atom=%2Fjnumed%2F42%2F7%2F1040.atom) [Web of Science](http://jnm.snmjournals.org/lookup/external-ref?access_num=A1996VE55500002&link_type=ISI) 9. Alpert NM, Eriksson L, Chang JY, et al. Strategy for the measurement of regional cerebral blood flow using short-lived tracers and emission tomography. J Cereb Blood Flow Metab 1984;4:28–34. [PubMed](http://jnm.snmjournals.org/lookup/external-ref?access_num=6607259&link_type=MED&atom=%2Fjnumed%2F42%2F7%2F1040.atom) [Web of Science](http://jnm.snmjournals.org/lookup/external-ref?access_num=A1984SF18600004&link_type=ISI) 10. Okazawa H, Vafaee M. Effect of vascular radioactivity on regional values of cerebral blood flow: evaluation of methods for H215O PET to distinguish cerebral perfusion from blood volume. J Nucl Med 2001;42:1032–1039. [Abstract/FREE Full Text](http://jnm.snmjournals.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Njoiam51bWVkIjtzOjU6InJlc2lkIjtzOjk6IjQyLzcvMTAzMiI7czo0OiJhdG9tIjtzOjIyOiIvam51bWVkLzQyLzcvMTA0MC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 11. Fujita H, Meyer E, Reutens DC, Kuwabara H, Evans AC, Gjedde A. Cerebral [15O] water clearance in humans determined by positron emission tomography: II. Vascular responses to vibrotactile stimulation. J Cereb Blood Flow Metab 1997;17:73–79. [PubMed](http://jnm.snmjournals.org/lookup/external-ref?access_num=8978389&link_type=MED&atom=%2Fjnumed%2F42%2F7%2F1040.atom) [Web of Science](http://jnm.snmjournals.org/lookup/external-ref?access_num=A1997VZ02800010&link_type=ISI) 12. Fujita H, Meyer E, Kuwabara H, Marrett S, Gjedde A. Comparison of apparent V images for water, oxygen and CO in PET. J Nucl Med 1993;34.(suppl):S199. 13. Chen W, Zhu XH, Kato T, Andersen P, Ugurbil K. Spatial and temporal differentiation of fMRI BOLD response in primary visual cortex of human brain during sustained visual simulation. Magn Reson Med 1998;39:520–527. [PubMed](http://jnm.snmjournals.org/lookup/external-ref?access_num=9543413&link_type=MED&atom=%2Fjnumed%2F42%2F7%2F1040.atom) [Web of Science](http://jnm.snmjournals.org/lookup/external-ref?access_num=000072700900003&link_type=ISI) 14. Kinahan PE, Noll DC. A direct comparison between whole-brain PET and BOLD fMRI measurements of single-subject activation response. Neuroimage 1999;9:430–438. [CrossRef](http://jnm.snmjournals.org/lookup/external-ref?access_num=10.1006/nimg.1998.0420&link_type=DOI) [PubMed](http://jnm.snmjournals.org/lookup/external-ref?access_num=10191171&link_type=MED&atom=%2Fjnumed%2F42%2F7%2F1040.atom) [Web of Science](http://jnm.snmjournals.org/lookup/external-ref?access_num=000079863600007&link_type=ISI) 15. Minoshima S, Cross DJ, Koeppe RA, Casey KL. Brain activation studies using PET and SPECT: execution and analysis. In: Casey KL, Bushnell MC, eds. Pain Imaging Seattle, WA: IASP Press; 2000:95–121. 16. Kanno I, Shimosegawa E, Fujita H, Hatazawa J. Uncoupling of absolute CBF to neural activity. Adv Exp Med Biol 1997;413:209–214. [PubMed](http://jnm.snmjournals.org/lookup/external-ref?access_num=9238502&link_type=MED&atom=%2Fjnumed%2F42%2F7%2F1040.atom) 17. Fox PT, Mintun MA, Raichle ME, Herscovitch P. A noninvasive approach to quantitative functional brain mapping with H215O and positron emission tomography. J Cereb Blood Flow Metab 1984;4:329–333. [PubMed](http://jnm.snmjournals.org/lookup/external-ref?access_num=6470051&link_type=MED&atom=%2Fjnumed%2F42%2F7%2F1040.atom) [Web of Science](http://jnm.snmjournals.org/lookup/external-ref?access_num=A1984TG21300004&link_type=ISI) 18. Friston KJ, Frith CD, Liddle PF, Dolan RJ, Lammertsma AA, Frackowiak RS. The relationship between global and local changes in PET scans. J Cereb Blood Flow Metab 1990;10:458–466. [PubMed](http://jnm.snmjournals.org/lookup/external-ref?access_num=2347879&link_type=MED&atom=%2Fjnumed%2F42%2F7%2F1040.atom) [Web of Science](http://jnm.snmjournals.org/lookup/external-ref?access_num=A1990DK16400002&link_type=ISI) 19. Casey KL, Morrow TJ, Lorenz J, Minoshima S. Temporal and spatial dynamics of human forebrain activity during heat pain: analysis by positron emission tomography. J Neurophysiol 2001;85:951–959. [Abstract/FREE Full Text](http://jnm.snmjournals.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mjoiam4iO3M6NToicmVzaWQiO3M6ODoiODUvMi85NTEiO3M6NDoiYXRvbSI7czoyMjoiL2pudW1lZC80Mi83LzEwNDAuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 20. Cherry SR, Woods RP, Doshi NK, Banerjee PK, Mazziotta JC. Improved signal-to-noise in PET activation studies using switched paradigms. J Nucl Med 1995;36:307–314. [Abstract/FREE Full Text](http://jnm.snmjournals.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Njoiam51bWVkIjtzOjU6InJlc2lkIjtzOjg6IjM2LzIvMzA3IjtzOjQ6ImF0b20iO3M6MjI6Ii9qbnVtZWQvNDIvNy8xMDQwLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==)