|
|
||||||||
Neuroradiology |
1 From the Imaging Research Laboratories, John P Robarts Research Institute, PO Box 5015, 100 Perth Drive, London, Ontario, N6A 5K8 Canada (D.G.N., A.C., T.Y.L.); Dept of Radiology and Lawson Research Institute, St Joseph's Health Centre, London, Ontario (A.C., J.D.B., R.K., T.Y.L.); Dept of Anaesthesia (R.A.C., A.W.G.), London Health Sciences Centre, University Campus, Ontario; and Dept of Neurology, Westfälische Wilhelms-Universität, Munster, Germany (D.G.N.). Received Aug 26, 1998; revision requested Oct 22; final revision received Jan 14, 1999; accepted Feb 15. Supported in part by Medical Research Council of Canada, Heart and Stroke Foundation of Canada, and GE Medical Systems. D.G.N. supported by a research grant of the Deutsche Forschungsgemeinschaft. Address reprint requests to T.Y.L. (e-mail: tlee@irus.rri.on.ca).
| Abstract |
|---|
|
|
|---|
MATERIALS AND METHODS: Twelve dynamic contrast materialenhanced CT studies were performed in beagles. CBV, CBF, and mean transit time (MTT) values were calculated by using an internal carotid artery (ICA) and a noncarotid artery as the input artery to the brain. Patient studies with use of the radial artery as the input were performed (a) repetitively in two patients after subarachnoid hemorrhage, (b) in a patient with a symptomatic ICA occlusion before and after the intravenous injection of 1 g of acetazolamide, and (c) in a patient with a malignant brain tumor.
RESULTS: Linear regression analyses revealed highly significant correlations (P < .001) between CBV (r, 0.98; slope, 0.96), CBF (r, 0.89; slope, 0.87), and MTT (r, 0.80; slope, 0.76) values calculated with the ICA and the noncarotid inputs. The CT-derived patient data correlated well with ancillary clinical and neuroradiologic findings.
CONCLUSION: Dynamic single-section CT scanning to measure CBV and CBF on the basis of a noncarotid input is a highly accessible and cost-effective blood flow measurement technique.
Index terms: Blood vessels, stenosis or obstruction, 17.72 Brain, blood flow, 10.919 Brain, CT, 10.12113, 10.12115 Brain, hemorrhage, 10.2521 Brain, perfusion, 10.919 Cerebral blood vessels, flow dynamics, 17.12112, 17.919 Computed tomography (CT), perfusion study, 10.919
| Introduction |
|---|
|
|
|---|
Several in vivo methods, such as positron emission tomography (PET) (3), single photon emission computed tomography (SPECT) (4), and magnetic resonance (MR) imaging (5,6) are presently available to measure CBF and cerebral blood volume (CBV). However, high costs together with the limited accessibility of and the difficulty in obtaining quantitative values by using SPECT and MR imaging (7,8) have restricted the widespread clinical use of these imaging modalities. Because computed tomographic (CT) scanners are available in most clinical facilities, a CT-based measurement technique has the advantage of easy accessibility for clinical and experimental studies.
As described nearly 20 years ago (9), by using injection of a bolus of contrast material and subsequent dynamic CT imaging to follow the contrast enhancement in an input artery and the brain itself, CBF and CBV values can be obtained by application of the "central volume principle" and the technique of deconvolution. However, due to the limited scanning frequency of CT scanners at the time, this method was not used (9,10).
Since the publication of Axel's work (9,10), several alternative approaches have been proposed to evaluate the mean transit time (MTT) with CT, with promising experimental and clinical results (1117). Compared with the deconvolution method, however, these techniques require additional assumptions and curve fitting procedures, each representing a potential source of variability. In a previous study (18), we validated our CT technique to measure CBF in a rabbit model against the reference standard of microspheres. When cerebral perfusion in humans at the level of the basal ganglia or higher is measured, however, the true supplying arteries (ie, the internal carotid artery [ICA]) usually are not depicted on the CT image. Therefore, other arteries must be selected to obtain the arterial contrast enhancement curve required for the deconvolution process. One suggestion was to use intraparenchymal cerebral arteries for this purpose (9,10), but they are very small and of variable appearance on the image, which reduces the reliability and reproducibility of this approach.
Alternatively, larger extracranial noncarotid arteries, such as the radial artery, can be scanned simultaneously together with the brain tissue, which provides a more consistent arterial input curve without partial volume averaging. It is known, however, that the vascular resistance of peripheral tissues is much higher than that of the brain parenchyma (19,20). Thus, differences in the shape of the contrast enhancement curves for normal ICA and normal noncarotid artery can be expected. To our knowledge, the effect of these differences on the calculation of CBF and CBV values is yet unknown.
We undertook this study to validate experimentally our CT method to measure CBF and CBV on the basis of a noncarotid input in beagles under normal physiologic conditions. In the second part of the study, the clinical feasibility of this technique and the spectrum of potential diagnostic applications is illustrated.
The theoretic basis of our measurement technique is the central volume principle (21). This principle relates CBF, CBV, and MTT values in the simple relationship
denotes the convolution operator and R(t) is the impulse residue function (Appendix, Fig 1b). R(t) is calculated by means of deconvolution (Appendix), and MTT is calculated by means of the area-over-height formula (10):
|
|
|
| MATERIALS AND METHODS |
|---|
|
|
|---|
After the surgical procedures were completed, each dog was placed in the prone position on the couch of the CT scanner, and the head was fixed in a conventional CT head holder to prevent movement during the subsequent experiments. Vecuronium bromide ([0.2 mg per kilogram of body weight] Norcuron; Organon Canada, Scarborough, Ontario, Canada) was administered intravenously every hour to stop spontaneous breathing. At completion of the experiments, the animals were sacrificed by means of barbiturate overdose, which was intravenously administered sodium pentobarbital ([162 mg/kg] Euthanyl Forte; MTC Pharmaceuticals, Cambridge, Ontario, Canada).
Patient Studies
The baseline characteristics of the four patients examined with the dynamic CT method are shown in Table 1. All patient studies were performed, after written informed consent was obtained, for diagnostic reasons to assess brain perfusion. In the two patients with subarachnoid hemorrhage, serial CT studies were performed twice per week within the first 3 weeks after subarachnoid hemorrhage and surgical clipping of the aneurysm. In the patient who experienced repetitive transient ischemic attacks owing to a known right ICA occlusion, digital subtraction angiograms revealed evidence of insufficient intracranial collateral flow with hypoplasia of the anterior communicating artery and very limited leptomeningeal collateral flow. One gram of acetazolamide (Diamox; Wyeth-Ayerst Canada, Saint-Laurent, Quebec, Canada) was injected intravenously to evaluate the reserve capacity of cerebral autoregulation. In patient 4, a primary malignant brain tumor within the left frontal cortical area was diagnosed by means of CT and MR imaging and was later confirmed by means of histologic findings.
|
|
CT scanning was initiated 5 seconds prior to the intravenous injection of 1.5 mL/kg iopromide by means of an automated injector (Mark IV; Medrad, Pittsburgh, Pa). This delay in injection of contrast material allowed for the acquisition of nonenhanced baseline images. To obtain a large range of CBF values, CT studies were performed with PaCO2 levels between 30 and 55 mm Hg. An interstudy delay of at least 20 minutes was chosen to ensure that most of the contrast material from previous injections had left the circulatory system.
CT Imaging Protocol for the Human Studies
Owing to the longer circulation time in humans, cine scanning of 90-second duration was performed with the same radiologic parameters as in the experimental studies. In the patient with the brain tumor, the section with the largest tumor cross section was chosen. In the remaining patients, the section containing the upper parts of the lateral ventricles was selected (Fig 3). The right forearm was placed on a specially designed support above the head fixed with Velcro to the head holder (Fig 3). This enabled axial scanning of the arteries of the forearm simultaneously with the head. Image reconstruction was performed as described earlier. A dose of 1 mL/kg iopromide was injected at a rate of 1.5 mL/sec.
|
In the human studies, a standardized set of ROIs was used as illustrated in Figure 3. A circular ROI was drawn (D.G.N.) to encompass the entire brain. A smaller copy of this ROI was then automatically generated and placed at the border between white and gray matters. The brain area was then automatically partitioned by radial lines at equal angular increments from a starting radial line chosen by the operator. In the patient with the brain tumor, additional ROIs were drawn that covered the tumor and the surrounding tissue.
The contrast enhancement curves were obtained by subtracting the regional mean baseline CT number (Hounsfield units) in the images obtained before the administration of contrast material from the mean CT number in serial contrast-enhanced images. Correction for partial volume averaging for the ICA and noncarotid contrast enhancement curves was performed as described previously (18). In brief, Gaussian curves were fitted to either the horizontal or vertical profile through the center of each arterial ROI. As has been shown experimentally, the SD of the fitted Gaussian curve is a good estimate of the internal diameter of the artery (18). The partial volume averaging correction factor can then be obtained from an experimentally determined calibration curve relating the correction factor to the internal diameter.
Statistics
Statistical analysis was performed by using the SPSS Package (SPSS, Chicago, Ill). The Student t test and Mann-Whitney U test were used to compare normally and nonnormally distributed data, respectively. Linear regression analysis and the Pearson product moment test were used to analyze the correlation of data sets. Significance was declared at P values less than .05.
| RESULTS |
|---|
|
|
|---|
For both calculation modes (ICA or noncarotid input), significantly higher CBV and CBF and lower MTT values were found for gray matter ROIs as compared with those for white matter ROIs (all P < .001, Mann-Whitney U test). A trend toward higher values for CBV and CBF was found when using the noncarotid artery as compared with the ICA to obtain the input curve (all P > .1, Mann-Whitney U test). This was found consistently for all subgroup analyses according to tissue type and size of ROI, although no significance was found (all P > .1, Mann-Whitney U test) (Table 2). Notably, the MTT showed an opposite trend toward lower values for the noncarotid- as compared with the ICA-based calculations (all P > .1), which reflects subtle differences in the shapes of the curves. The difference in the ICA- and noncarotid-derived MTT among the ROI subgroups approximated 0.2 second (Table 2).
|
|
|
|
|
Malignant brain tumor.Within the tumor, a pathologic increase in both CBF (155 mL/100 g/min) and CBV (5.5 mL/100 g) values was found. These were markedly higher than both the mean CBF and mean CBV values (44.3 mL/100 g/min and 3.3 mL/100 g) of all the other ROIs, as well as the mean values (52.1 mL/100 g/min and 3.9 mL/100 g) for the gray matter ROIs. Owing to the mass effect of the tumor, the two cortical ROIs adjacent to the tumor had much lower CBF (17.5 mL/100 g/min) and CBV (2.0 mL/100 g) values than those of the average of the remaining cortical ROI (52.1 mL/100 g/min and 3.9 mL/100 g).
| DISCUSSION |
|---|
|
|
|---|
The MTT and the CBF measurements likewise showed a strong correlation between the ICA- and the noncarotid-derived values (r, 0.80 and r, 0.89, respectively). However, the slopes of their regression lines, especially for the MTT, were not as good as the CBV results. In contrast to the CBV calculation, the initial slope of the arterial enhancement curve determines the results of the deconvolution-based MTT and CBF values. Our results suggest that minor differences between the shape of the ICA contrast enhancement curve and the noncarotid contrast enhancement curve existed, which affected the deconvolution process (Fig 1a). With lower MTT values (eg, in gray matter ROIs), an increasing effect on the flow measurements is expected, even with subtle differences in the ascending slopes of two different arterial curves. The latter is supported by the observation that the largest differences between the ICA and noncarotid CBF measurements (>20 mL/100 g/min) mostly occurred for tissue flow values higher than 80 mL/100 g/min (Figs 4b, 5). In human subjects, however, the range in which accurate CBF measurements are most critical is 080 mL/100 g/min. Within this range, the absolute measurement error with use of the surrogate noncarotid input was acceptable for both CBV and CBF values (Figs 4, 5).
In the second part of the study, the clinical feasibility of our dynamic single-section CT technique was demonstrated. By using a simple support for one of the patient's forearms, which can be easily attached to the scanner's head holder, the radial artery can be scanned simultaneously with the patient's head, which provides a consistent and reliable arterial input for the CT measurements. Similar to the extracranial neck arteries used in the beagle experiments, the radial artery also supplies muscle and connective tissue. Thus, similar vascular resistances in the vascular beds downstream from both arteries can be expected. However, in the beagle studies, both ICA and noncarotid enhancement occurred simultaneously in the CT images. In contrast, in the patient studies, radial arterial enhancement appeared with a short delay after the onset of brain enhancement owing to the longer travel time to the forearm.
To demonstrate the validity of our radial artery approach in humans, it has to be shown that this longer travel time has no substantial effect on the deconvolution-based calculations. The delay itself can be corrected for by shifting the arterial enhancement curve to the left of the brain enhancement curve and letting the deconvolution algorithm determine the optimal separation (in the least squares sense) between the two curves. However, it generally is believed that a longer arterial travel time increases the dispersion of a given tracer (25). A more pronounced dispersion, as compared with that of the ICA curve, would make the radial artery invalid as a surrogate for the arterial input for the assessment of cerebral hemodynamics.
We have shown recently by means of a theoretic approach, as well as in vitro experiments in which vessel phantoms were used, that arterial dispersion is not a substantial point of concern for normal flow conditions (23). Thus, with the normal hemodynamic conditions present in the circulation system, an intravenously introduced tracer will appear with nearly identical dispersion in all the major vessels of the arterial tree. In addition, pathologically increased dispersion due to a high-grade arterial stenosis along the path to the forearm is a potential source of errors. However, peripheral arterial disease mostly occurs in the legs, and high-grade arterial obstruction in the upper extremities represents less than 5% of all clinical cases (26). Furthermore, substantial arteriopathy of the arms can be excluded by routine clinical history and the presence of regular radial arterial pulses (27).
Thus, alterations of the shape of the radial artery curve, owing to either the longer path to the forearm, as compared with the path to the brain, or arterial stenoses are not substantial points of concern. Hence, the radial artery's contrast enhancement curve can serve as a surrogate of the arterial input to the brain to measure CBV and CBF values. We are now investigating the CBF in healthy human subjects to obtain normal values for CBV and CBF by using this method.
However, the case in which a normal radial artery is used as the input function for brain tissue supplied by a severely obstructed ICA has to be discussed. For the core of an ischemic region, it is reasonable to expect that the arterial input to this type of tissue will be subjected to more dispersion than regions with normal or slight decrease in blood flow owing to the increase in the length of the collateral pathways supplying that region. Østergaard et al (28) have shown preliminary evidence of this phenomenon in a patient with cerebrovascular disorder.
If we assume that the additional dispersion can be characterized by the function h(t), then instead of the true impulse residue function, R(t), our radial arterybased deconvolution algorithm will calculate the convolution of h(t) and R(t) (ie, h(t)
R(t)). In the case of insubstantial dispersion, the width of h(t) is equal to or less than the plateau width of R(t). Then the plateau height of R(t), and hence, CBF calculation, is not affected. In addition, because the area of both R(t) and h(t)
R(t) will be identical, MTT will not be affected. However, in case of substantial dispersion (ie, the width of h[t] is greater than the plateau width of R[t]), then the maximum height of h(t)
R(t) will be less than that of R(t). As a result, CBF will be underestimated. Even in this situation, the area of h(t)
R(t) will be the same as R(t), so that the CBV value calculated by using the radial arterial input will not be affected.
The magnitude of the underestimation in CBF and overestimation in MTT cannot be quantified without detailed knowledge of h(t). Greitz (29) has shown that intracranial dispersion in the brain is generally insubstantial both in healthy subjects and in patients with obstructed input pathways. However, this was an old study based on nondigital angiography with limited time resolution. Owing to the importance of this assumption for our method, confirmation of Greitz's results by using modern digital subtraction angiographic equipment capable of higher framing rates is desirable and being planned.
Results in our pilot series were in good agreement with the current understanding of pathophysiology. In all patients, gray matter CBF (~4550 mL/100 g/min) and CBV (~4 mL/100 g) values were approximately two times higher than the respective white matter values (3,4). During the follow-up in two patients with subarachnoid hemorrhage, a marked reduction of the CBF was observed in the patient with vasospasm in contrast to normal flow rates in the other patient. In the patient with vasospasm, the time course and amount of flow changes for white and gray matters correlated well with observations from a large-scale study in which the xenon 133 inhalation method was used (30).
In the patient with symptomatic ICA occlusion, evidence of activated autoregulation at rest and exhausted vasomotor reactivity after injection of acetazolamide were found for the affected hemisphere, in accordance with clinical and angiographic findings. The contralateral increases in CBF and CBV values in this patient (~20%) were notably lower than the fractional increases of 30%40% reported for healthy subjects after injection of acetazolamide (31). This could be explained by underlying arteriopathy on the clinically unaffected side or by a steal effect due to collateral flow toward the ischemic hemisphere. The striking increase in CBF (150 mL/100 g/min) and CBV (5.5 mL/100 g) values measured within the tumor, and substantial ischemia present in the adjacent cortical tissue, is in accordance with knowledge of the pathophysiology and hemodynamics of brain tumors (32). Thus, this pilot series demonstrates the feasibility of this technical approach and illustrates the broad clinical spectrum of potential applications of this CT method. Nevertheless, controlled clinical studies are necessary to delineate more accurately the clinical capability and limitations of this measurement technique.
Further methodological issues have to be discussed. First, the radiation dose of this CT technique has to be defined. Although the number of sections required for a dynamic CT study in humans is about six to seven times higher compared with that required for a regular diagnostic scan of the brain (90 vs 1014 scans), much lower radiologic parameters are used (80 kVp and 80 mAs vs 120 kVp and 340 mAs). Thus, the overall effective dose equivalent required for a dynamic CT study (~2.0 mSv) is not much higher than that required for a routine CT scan (~1.5 mSv) (33). Moreover, this dose is below the dose equivalent delivered with other blood flow measurement techniques such as PET (34) and SPECT (35).
Second, our CT measurements are restricted to a single anatomic level of the brain. Although measurements at different brain levels can be obtained sequentially, this would lead to a very low time resolution of the measured contrast enhancement curves. The latter would result in less than optimal data interpolation and curve fitting procedures and lead to an unacceptably low accuracy of CBF measurements. The introduction of multisection CT scanners with the ability to scan simultaneously at different levels may overcome this limitation in the future.
In summary, we have validated a dynamic single-section CT technique to measure CBV and CBF values by using a noncarotid artery as a surrogate arterial input for deconvolution-based calculations and have demonstrated its feasibility in patient studies. Because the use of the radial artery in humans is easy to apply and provides consistent enhancement curves without substantial partial volume averaging, this CT approach enables accurate, readily available, and cost-effective measurements of CBV and CBF values in a variety of clinical and experimental settings.
| Appendix 1 |
|---|
|
|
|---|
Correction for Extravasation of Contrast Material
When the capillary endothelium leaks contrast material, we have shown both by heuristic arguments (22) and experimental validation (23) that R(t), in this case, can be expressed as:
Correction for Difference in Tissue and Large-Vessel Hematocrit Values
Contrast material is confined to the plasma phase of blood. As long as the hematocrit values of blood remain the same both in peripheral (large) blood vessels and in tissue capillaries, then the enhancement measured in tissue and blood are equivalent. To correct the tissue enhancement for the difference in large-vessel and small-vessel (tissue) hematocrit values, Q(t) has to be multiplied by the factor
= (1 - H)/(1 - rH) = 0.88, where H (0.4) is the hematocrit value of blood in large vessels and r (0.80.85) is the ratio of small- to large-vessel hematocrit values (4).
| Acknowledgments |
|---|
| Footnotes |
|---|
Author contributions: Guarantors of integrity of entire study, D.G.N., T.Y.L.; study concepts, T.Y.L., A.W.G.; study design, A.C., R.A.C.; definition of intellectual content, D.G.N., T.Y.L., A.W.G.; literature research, D.G.N., A.C.; clinical studies, J.D.B., R.K., T.Y.L.; experimental studies, J.D.B., R.K., D.G.N., A.C., R.A.C., T.Y.L.; data acquisition, A.C., R.A.C.; data analysis, D.G.N., A.C.; statistical analysis, D.G.N., T.Y.L., R.A.C.; manuscript preparation, D.G.N., T.Y.L.; manuscript editing and review, all authors
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
K. Murayama, K. Katada, M. Nakane, H. Toyama, H. Anno, M. Hayakawa, D. San Millan Ruiz, and K. J. Murphy Whole-Brain Perfusion CT Performed with a Prototype 256-Detector Row CT System: Initial Experience Radiology, January 1, 2009; 250(1): 202 - 211. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. G. Sheiman and A. Sitek Feasibility of Measurement of Pancreatic Perfusion Parameters with Single-Compartment Kinetic Model Applied to Dynamic Contrast-enhanced CT Images Radiology, December 1, 2008; 249(3): 878 - 882. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. Makkat, R. Luypaert, T. Stadnik, C. Bourgain, S. Sourbron, M. Dujardin, J. De Greve, and J. De Mey Deconvolution-based Dynamic Contrast-enhanced MR Imaging of Breast Tumors: Correlation of Tumor Blood Flow with Human Epidermal Growth Factor Receptor 2 Status and Clinicopathologic Findings--Preliminary Results Radiology, November 1, 2008; 249(2): 471 - 482. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. Sabir, R. Schor-Bardach, C. J. Wilcox, S. Rahmanuddin, M. B. Atkins, J. B. Kruskal, S. Signoretti, V. D. Raptopoulos, and S. N. Goldberg Perfusion MDCT Enables Early Detection of Therapeutic Response to Antiangiogenic Therapy Am. J. Roentgenol., July 1, 2008; 191(1): 133 - 139. [Abstract] [Full Text] [PDF] |
||||
![]() |
H.M. Silvennoinen, L.M. Hamberg, P.J. Lindsberg, L. Valanne, and G.J. Hunter CT Perfusion Identifies Increased Salvage of Tissue in Patients Receiving Intravenous Recombinant Tissue Plasminogen Activator within 3 Hours of Stroke Onset AJNR Am. J. Neuroradiol., June 1, 2008; 29(6): 1118 - 1123. [Abstract] [Full Text] [PDF] |
||||
![]() |
F. Gaudiello, V. Colangelo, F. Bolacchi, M. Melis, R. Gandini, F.G. Garaci, V. Cozzolino, R. Floris, and G. Simonetti Sixty-Four-Section CT Cerebral Perfusion Evaluation in Patients with Carotid Artery Stenosis before and after Stenting with a Cerebral Protection Device AJNR Am. J. Neuroradiol., May 1, 2008; 29(5): 919 - 923. [Abstract] [Full Text] [PDF] |
||||
![]() |
S.K. Ellika, R. Jain, S.C. Patel, L. Scarpace, L.R. Schultz, J.P. Rock, and T. Mikkelsen Role of Perfusion CT in Glioma Grading and Comparison with Conventional MR Imaging Features AJNR Am. J. Neuroradiol., November 1, 2007; 28(10): 1981 - 1987. [Abstract] [Full Text] [PDF] |
||||
![]() |
H.M. Silvennoinen, L.M. Hamberg, L. Valanne, and G.J. Hunter Increasing Contrast Agent Concentration Improves Enhancement in First-Pass CT Perfusion AJNR Am. J. Neuroradiol., August 1, 2007; 28(7): 1299 - 1303. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. C. Faria, C. S. Ng, K. R. Hess, S. Phongkitkarun, J. Szejnfeld, D. Daliani, and C. Charnsangavej CT Quantification of Effects of Thalidomide in Patients with Metastatic Renal Cell Carcinoma Am. J. Roentgenol., August 1, 2007; 189(2): 378 - 385. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. P. Coles Imaging after brain injury Br. J. Anaesth., July 1, 2007; 99(1): 49 - 60. [Abstract] [Full Text] [PDF] |