Radiology
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


Published online before print June 11, 2003, 10.1148/radiol.2282020699
This Article
Right arrow Abstract Freely available
Right arrow Figures Only
Right arrow Full Text (PDF)
Right arrow All Versions of this Article:
2282020699v1
228/2/473    most recent
Right arrow Submit a response
Right arrow Alert me when this article is cited
Right arrow Alert me when eLetters are posted
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Lucidarme, O.
Right arrow Articles by Berger, G.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Lucidarme, O.
Right arrow Articles by Berger, G.
(Radiology 2003;228:473-479.)
© RSNA, 2003


Experimental Studies

Blood Flow Quantification with Contrast-enhanced US: "Entrance in the Section" Phenomenon—Phantom and Rabbit Study1

Olivier Lucidarme, MD, Stephanie Franchi-Abella, MD, Jean-Michel Correas, MD, PhD, S. Lori Bridal, PhD, Erol Kurtisovski, PhD and Geneviève Berger, MD, PhD

1 From the Parametric Imaging Laboratory, UMR 7623 CNRS and Paris University VI, 15 rue de l’École de Médecine, 75006 Paris, France (O.L., S.F.A., J.M.C., S.L.B., E.K., G.B.); Department of Radiology, Necker Hospital, AP-HP, Paris, France (O.L., S.F.A., J.M.C.); and Department of Radiology, Pitié-Salpêtrière Hospital, AP-HP, Paris, France (O.L.). Received June 13, 2002; revision requested August 8; revision received October 17; accepted December 10. Address correspondence to S.L.B. (e-mail: bridal@lip.bhdc.jussieu.fr).


    ABSTRACT
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 APPENDIX
 REFERENCES
 
PURPOSE: To investigate changes in destruction-replenishment curves (in vitro and in vivo) that result from microbubble destruction in feeding vessels that pass through the imaging plane before microbubbles enter the region of interest (ROI).

MATERIALS AND METHODS: During continuous injections of an ultrasonographic contrast agent, nonlinear gray-scale images were obtained in vitro in the longitudinal plane of a renal dialysis cartridge flow phantom (flow rates of 100, 200, and 400 mL/min) and in vivo in the coronal plane of the left kidneys of two rabbits (two kidneys). Destruction-replenishment curves were obtained for the dialysis cartridge in ROIs located immediately after the entrance of the microbubbles into the image plane and further from the entrance, after microbubbles had traveled across the complete length of the imaging plane. Replenishment curves were also obtained from ROIs in the rabbit kidneys at the level of segmental arteries, distal interlobar arteries, and the cortex.

RESULTS: The ROIs immediately after the entrance of the microbubbles in the image plane of the dialysis cartridge and in the segmental artery of the kidney followed a typical exponential function, A(1 - e-{alpha}t). Early portions of curves obtained in ROIs filled with microbubbles that had already passed through the image plane of the dialysis cartridge or in the renal cortex were not well described by such a function. The shape of the curve and the variations as a function of flow rate can be explained by means of a mathematical model based on indicator-dilution theory.

CONCLUSION: When the feeding vessels of an ROI travel across the ultrasound field before they reach the measurement region, the typical shape of the replenishment curve is modified (reduced velocity parameter and plateau).

© RSNA, 2003

Index terms: Blood, flow dynamics, 81.12988 • Kidney, US, 81.12984, 81.12988 • Microbubbles, 81.12988 • Phantoms, 81.12984, 81.12988 • Test objects, 81.12984, 81.12988 • Ultrasound (US), contrast media 81.12984, 81.12988


    INTRODUCTION
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 APPENDIX
 REFERENCES
 
Findings in previous studies suggest that intravenous injection of microbubble-based ultrasonographic (US) contrast agents allows estimation of relative blood flow and fractional blood volume of the microvasculature in a region of interest (ROI). Because US can destroy microbubbles (1,2) the change in echogenicity after tissue has been cleared of microbubbles can be followed (destruction-replenishment technique). Assuming a steady state for the microbubble concentration in the circulation and a linear relationship between microbubble concentration in the studied organ and the backscattered signals, such destruction-replenishment curves are described by an exponential function that asymptotically approaches a maximum (3). In this case, the initial slope of the curve is related directly to the velocity at which microbubbles replenish the ROI that represents the blood velocity, while the plateau can be related to the fractional blood volume (3,4). This theory has been used with intermittent high transmit power or continuous low transmit power imaging (interspersed with destructive pulses) to assess the perfusion of myocardium (5,6), brain (7), and kidney (8). This model, however, is based on a continuous infusion of contrast material into a single open compartment and is valid only if we assume that a constant number of microbubbles enters the ROI per unit time. Specifically, microbubble destruction that occurs in vessels before their entrance in the ROI must be negligible.

The purpose of this study was to investigate changes in destruction-replenishment curves (in vitro and in vivo) that result from microbubble destruction in feeding vessels that pass through the imaging plane before microbubbles enter the ROI.


    MATERIALS AND METHODS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 APPENDIX
 REFERENCES
 
An indicator-dilution model was developed that takes into account the effect of microbubble destruction located in the imaging plane during the destruction pulse on the form of the destruction-reperfusion curve for a selected ROI. Details of this theoretic model are given in the Appendix.

In Vitro Experiment
All in vitro measurements were performed with a nonrecirculating flow phantom consisting of a renal dialysis cartridge (Hémoflow model F60S; Fresenius Medical Care, Bad Homburg, Germany) with an internal diameter of 47 mm, which contained approximately 9,000 parallel polysulfone hollow fibers with a diameter of 200 µm, length of 25 cm, and total fiber volume of 70 mL (Fig 1). A dialysis cartridge was chosen because of the unidirectional alignment of the hollow fibers and because flow rate and the fractional circulating volume are constant throughout. The dialysis cartridge was connected to a 10-L reservoir positioned 2 m above the dialysis cartridge. The distal tubing was connected to a peristaltic pump (Masterflex, model 7523-27; Cole-Parmer Instruments, Vernon Hills, Ill). An electric syringe (Pulsar Ultrasound Injection System; Medrad, Indianola, Pa) was connected to the proximal side of the tubing system to inject the US contrast agent in the circuit. The direction of the flow at measurement sites within the circuit was oriented upward through vertical segments to minimize microbubble distribution gradients across the diameter of the filter and tubing that could be caused by upward movement of microbubbles that results from buoyant forces.



View larger version (19K):
[in this window]
[in a new window]
[Download PPT slide]
 
Figure 1. Diagram of flow phantom. Dialysis filter and tubing were positioned so that the solution was pumped up a vertical segment of tubing at the measurement sites. Parts of phantom are reservoir located 2 m above the dialysis cartridge (1), automatic injector filled with US contrast agent (2), continuous wave (CW) Doppler measurement system (3), dialysis cartridge with 9,000 200-µm-diameter hollow fibers (4), US machine with a sector transducer (5), and peristaltic pump (6).

 
The fluid used in the experiment was nondegassed water at 24°C that was circulated at three flow rates: 100, 200, and 400 mL/min. Intensity of the Doppler audio signal is proportional to the number of moving targets; therefore, the microbubble concentration in the main tubing was monitored with Doppler intensitometry (continuous wave Doppler) (model Ethydop 100E; Ethypharm, Saint-Cloud, France) (9). The Doppler signals were recorded in quadrature at low acoustic power with an 8-MHz pencil probe stabilized with a mechanical device (O.L.). Acoustic pressure at the measurement site was 95 kPa as evaluated (S.L.B., O.L.) with a calibrated 1-mm-diameter hydrophone (PVDF; Precision Acoustic, Dorchester, England). The Doppler signal was digitized with a sound card and analyzed with software (MatLab; Mathworks, Natick, Mass) to calculate the mean Doppler powers (in decibels) for each 0.14-second period from the Doppler spectrum (10).

Contrast agent injection.—The US contrast agent used in this study (SonoVue; Bracco, Milan, Italy) consists of sulphur hexafluoride microbubbles obtained by means of reconstitution of a lyophilized powder with 5 mL of 0.9% saline solution. One milliliter of solution contains approximately 2 x 108 microbubbles stabilized by a layer of phospholipids (11). The mean microbubble diameter is 2.5 µm with a narrowed distribution (90% with a diameter of less than 8 µm) (11). Concentration of the microbubbles was kept constant at 5 mL of the US contrast agent per liter of fluid in the phantom; thus, the infusion rates of the US contrast agent were 0.5, 1.0, and 2.0 mL/min, respectively, for peristaltic pump flow rates of 100, 200, and 400 mL/min. No attempt was made to agitate the contrast agent in the syringe during the infusion. Duration of each infusion was less than 120 seconds.

Image acquisition.—Once the microbubble concentration in the circulation had stabilized, the dialysis cartridge was imaged (O.L.) with a clinical US unit (HDI 5000; ATL-Philips, Bothell, Wash) with a curvilinear transducer (C5-2 probe; ATL-Philips). The image plane was chosen parallel to the hollow fibers (sagittal plane). Nonlinear gray-scale mode was obtained with a pulse-inversion technique. To eliminate all microbubbles from the region, we acquired a set of 20 frames at maximum power (mechanical index of 1.1) for 0.8 second (24 frames per second). This destruction set was followed immediately by acquisition of 408 observation frames at low output power (mechanical index of 0.08) for 17 seconds (24 frames per second). Persistence of the US machine image display was set to zero. Focal zone was positioned in the first third (1.5 cm deep) of the dialysis cartridge. The destruction-replenishment experiments were repeated three times—30, 60, and 90 seconds after the beginning of the infusion—to allow averaging of the data.

In Vivo Experiment
Animal model.—After the institutional animal care committee at Necker Hospital approved the experimental procedure, two 1.5-month-old New Zealand male rabbits, which weighed 1.7 and 1.8 kg, were studied with general anesthesia (atropine; Aguettant, Lyon, France), 20 µg; ketamine (Ketalar; Parke Davis, Courbevoie, France), 35 mg per kilogram of body weight; and xylazine (Paxman; Vibrac Production Animale, Carros, France), 5 mg/kg. No paralytic agents were used, and oxygen (4 L/min) was delivered during the procedure. The animals were placed in the supine position, and the left kidney was imaged. The study was focused on the kidney because it presents a well-organized arterial tree with the feeding vessel from the hilum to the cortex included in a coronal section. The US contrast agent was administered in the right femoral vein at a continuous infusion rate of 0.5 mL/min (total dose of 3 mL). As in the flow phantom experiment, Doppler intensitometry was performed to detect the steady state of microbubble blood concentration. The 8-MHz pencil probe was stabilized with a mechanical device (S.F.A.) to measure the Doppler intensity in the left femoral vein.

The imaging sequence was performed as soon as the blood concentration of microbubbles reached equilibrium. Nonlinear gray-scale imaging with the pulse-inversion technique was performed with a low mechanical index (0.08) by two observers (S.F.A., J.M.C.) by using a linear transducer (model L7-4; ATL-Philips). The left kidney of each rabbit was imaged in the coronal plane (two kidneys). After four destructive frames were acquired at maximum power, 225 frames were acquired at a rate of 15 frames per second for 15 seconds. Persistence of the pulse sequence was set to zero. The focal zone was positioned in the middle of the kidney.

Imaging Analysis
Digital cine loops were transferred to a personal computer. Gray-scale enhancement was measured with the manufacturer software (HDI Lab, version 1.90a; ATL-Philips), which allows quantification of pixel intensity before logarithmic compression for video display.

For the in vitro experiment, two ROIs of 5 x 5 mm were positioned (O.L.) at the proximal and distal ends of the field of view. The proximal ROI included microbubbles that were just entering the ultrasound field. The distal ROI was filled with microbubbles that had traveled across the complete length of the ultrasound field. The proximal and distal ROIs were both located 0.5 cm deep within the phantom to minimize effects of shadowing. They were positioned symmetrically with respect to the central axis of the probe to ensure comparable beam characteristics (which influence backscatter amplitude and regularity) in the two zones.

For the in vivo experiment, four contiguous or slightly overlapping ROIs (1 x 1 mm) were positioned (O.L.) on the segmental arteries near the hilum (segmental artery ROI), on the distal part of the corresponding interlobar artery (interlobar artery ROI), and on the renal microcirculation of the cortex fed by this artery (cortex ROI) (Fig 2). For each location, the signal intensities in the four ROIs were averaged to reduce the signal variations and noise within each ROI. The segmental artery ROI was the first ROI to intersect with flowing contrast agent (the renal artery was not included in the imaging plane). In the interlobar and cortex ROIs, the signal intensity was measured in ROIs filled with microbubbles that had traversed the ultrasound field via the segmental artery or other portions of the microvasculature. In addition, we measured the enhancement of the cortex of each left kidney in a large cortex ROI, which surrounded the cortical area (Fig 2), to obtain maximum smoothing of the replenishment curve of the cortex.



View larger version (159K):
[in this window]
[in a new window]
[Download PPT slide]
 
Figure 2. Coronal US scan in left kidney of rabbit 2 shows positions of the large cortex ROI (ROI Large), the three groups of four ROIs positioned on a segmental artery (ROIs Seg), the corresponding interlobar artery ROI (ROIs Ila), and the cortex ROI (ROIs Cort) fed by this artery.

 
The intensity of contrast enhancement was defined as the difference between the signal intensity (in linear units) after injection and the mean baseline value. Experimental time-intensity curves were plotted and compared with those calculated with the theoretic model described in the Appendix.


    RESULTS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 APPENDIX
 REFERENCES
 
In Vitro Experiment
The circulating microbubble concentration, monitored with Doppler intensitometry, always reached equilibrium less than 30 seconds after the beginning of infusion. Between 30 and 120 seconds, the Doppler intensity in the tubing immediately before entrance into the dialysis cartridge was 19.5 dB ± 0.5 (mean ± SD), 19.6 dB ± 0.4, and 19.0 dB ± 0.2 for infusion rates and flow rates of 0.5 and 100 mL/min, 1.0 and 200 mL/min, and 2.0 and 400 mL/min, respectively. The signal intensities observed in the dialysis cartridge as a function of time and flow rate for proximal and distal ROIs are plotted in Figure 3a and 3b. The corresponding values calculated for the concentration of the microbubbles as a function of time, given by the model developed in the present study with a fixed arbitrary value of the destruction coefficient and an arbitrary base value of the transit time, adjusted for the three curves by means of ratios that correspond to those between the three flow rates, are plotted in Figure 3c and 3d. The destruction-reperfusion curves predicted by the model present variations as a function of ROI position and flow rate that are similar to those for experimentally measured curves.



View larger version (40K):
[in this window]
[in a new window]
[Download PPT slide]
 
Figure 3a. Scatterplots for the (a) proximal and (b) distal ROIs depict signal intensities in the dialysis cartridge as a function of time and flow rate. (c, d) Scatterplots depict the corresponding value calculated for the concentration of microbubbles as a function of time as given by the model in the first (c) and third (d) subvolumes. Destruction coefficient {lambda} was fixed arbitrarily at 0.5, with three values of {tau} (related with the flow rate by the relation 1/{tau} = F/Vb) (Appendix). Values of {tau} were decreased each time by a factor of 2 to simulate a doubling of the flow rate. Calculated destruction-reperfusion curves in c and d depict results that correspond well to the experimental flows in a and b, respectively.

 


View larger version (46K):
[in this window]
[in a new window]
[Download PPT slide]
 
Figure 3b. Scatterplots for the (a) proximal and (b) distal ROIs depict signal intensities in the dialysis cartridge as a function of time and flow rate. (c, d) Scatterplots depict the corresponding value calculated for the concentration of microbubbles as a function of time as given by the model in the first (c) and third (d) subvolumes. Destruction coefficient {lambda} was fixed arbitrarily at 0.5, with three values of {tau} (related with the flow rate by the relation 1/{tau} = F/Vb) (Appendix). Values of {tau} were decreased each time by a factor of 2 to simulate a doubling of the flow rate. Calculated destruction-reperfusion curves in c and d depict results that correspond well to the experimental flows in a and b, respectively.

 


View larger version (22K):
[in this window]
[in a new window]
[Download PPT slide]
 
Figure 3c. Scatterplots for the (a) proximal and (b) distal ROIs depict signal intensities in the dialysis cartridge as a function of time and flow rate. (c, d) Scatterplots depict the corresponding value calculated for the concentration of microbubbles as a function of time as given by the model in the first (c) and third (d) subvolumes. Destruction coefficient {lambda} was fixed arbitrarily at 0.5, with three values of {tau} (related with the flow rate by the relation 1/{tau} = F/Vb) (Appendix). Values of {tau} were decreased each time by a factor of 2 to simulate a doubling of the flow rate. Calculated destruction-reperfusion curves in c and d depict results that correspond well to the experimental flows in a and b, respectively.

 


View larger version (22K):
[in this window]
[in a new window]
[Download PPT slide]
 
Figure 3d. Scatterplots for the (a) proximal and (b) distal ROIs depict signal intensities in the dialysis cartridge as a function of time and flow rate. (c, d) Scatterplots depict the corresponding value calculated for the concentration of microbubbles as a function of time as given by the model in the first (c) and third (d) subvolumes. Destruction coefficient {lambda} was fixed arbitrarily at 0.5, with three values of {tau} (related with the flow rate by the relation 1/{tau} = F/Vb) (Appendix). Values of {tau} were decreased each time by a factor of 2 to simulate a doubling of the flow rate. Calculated destruction-reperfusion curves in c and d depict results that correspond well to the experimental flows in a and b, respectively.

 
In Vivo Experiment
Microbubble concentration in the left femoral vein of each rabbit, which was monitored with Doppler intensitometry, reached an equilibrium 80 seconds after the beginning of infusion. Between 80 and 100 seconds, mean Doppler intensity was 21.8 dB ± 0.8 and 19.4 dB ± 0.2 for rabbits 1 and 2. Evolution of the signal intensity measurements obtained in segmental artery, interlobar artery, and cortex ROIs in the left kidney of rabbits 1 and 2 are shown in Figure 4a and 4b. Figure 5 shows the destruction-replenishment curve for the large cortex ROIs in the left kidney of the two rabbits. Replenishment curves obtained for segmental artery ROIs were comparable for the two rabbits, with the same mean signal intensity values at the plateau (mean gray level of approximately 90). In the interlobar artery ROIs and mainly in the cortex ROIs, the replenishment curves exhibited a longer initial low slope portion and a lower plateau, particularly in rabbit 2.



View larger version (34K):
[in this window]
[in a new window]
[Download PPT slide]
 
Figure 4a. Mean destruction-reperfusion curves in the left kidney of (a) rabbit 1 and (b) rabbit 2 were obtained for segmental artery ROIs (top), the corresponding interlobar artery ROIs (middle), and the cortex ROIs (bottom). Replenishment curves begin with a lower slope than would be described by the exponential form A(1 - e-{alpha}t). This sigmoid aspect is weakest for the segmental artery ROIs and strongest for the cortex ROIs.

 


View larger version (37K):
[in this window]
[in a new window]
[Download PPT slide]
 
Figure 4b. Mean destruction-reperfusion curves in the left kidney of (a) rabbit 1 and (b) rabbit 2 were obtained for segmental artery ROIs (top), the corresponding interlobar artery ROIs (middle), and the cortex ROIs (bottom). Replenishment curves begin with a lower slope than would be described by the exponential form A(1 - e-{alpha}t). This sigmoid aspect is weakest for the segmental artery ROIs and strongest for the cortex ROIs.

 


View larger version (22K):
[in this window]
[in a new window]
[Download PPT slide]
 
Figure 5. Destruction-reperfusion curves obtained in the large cortex ROIs in rabbit 1 ({bullet}) and rabbit 2 ({circ}). Averaging of the data from all pixels in these ROIs leads to smoothing of the curves, which emphasizes the sigmoid aspect of the initial portion of the replenishment curves obtained in the parenchyma when the imaging plane intersects the feeding arteries.

 

    DISCUSSION
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 APPENDIX
 REFERENCES
 
An exponential function approaching a maximum A(1 - e-{alpha}t) based on a single-tank model has generally been proposed to fit destruction-replenishment curves obtained in vivo in myocardium, brain, and kidney (5,6,8,12,13) during continuous perfusion of US contrast agents. According to this equation, microbubble velocity in the microcirculation is related directly to {alpha}, and the plateau echo enhancement A correlates with the fractional blood volume when the microbubble destruction induced by the ultrasound field is null. The most critical condition for this type of evaluation is to know the function that relates image intensity to microbubble concentration. Since such a function remains difficult to characterize, often only a relative flow assessment can be achieved (same subject and same machine settings). This approach can still be very useful, for example, to study the effectiveness of angiogenic drugs or reperfusion therapy (5).

This exponential model, however, is valid only if the concentration of microbubbles that enter the ROI immediately after the destruction pulse is constant. This supposes that two hypotheses must be fulfilled: First, the injection rate of the contrast media remains constant; and second, the blood vessels that enter the ROI have not already been submitted to the ultrasound beam. Findings in our in vitro study show that when vessels that feed the ROI have previously been emptied of microbubbles, the replenishment curve is no longer an exponential that approaches a maximum but is a sigmoid; we describe this theoretically with an indicator-dilution model (Appendix). In this case, the initial slope of the replenishment curve is theoretically equal to zero. For example, the slope measured at the inflexion point, even with a constant flow rate and a constant microbubble concentration in the circulation, decreases when the length of the path followed by the microbubbles in the ultrasound field increases before they reach the ROI (Fig 3). Consequently, microbubble velocity in the ROI estimated from the initial slope ({alpha}) of the refill curve, in the classic model described by an exponential function that approaches a maximum (3), is not correct if microbubble destruction occurs in the feeding vessels.

In the same way, the plateau of enhancement, which classically gives information about the fractional blood volume, is related to the value C0/(1 + {tau}{lambda})n given by our model (Appendix). When the velocity of the microbubbles is high (ie, when {tau} is very small) compared with the destruction rate or when the destruction rate is close to zero, {tau}{lambda} is negligible, and the measured signal intensity reflects the fractional blood volume. However, findings in our in vitro experiment (Fig 3) showed a decrease in the signal intensity measured at the plateau when the flow rates decreased (when {tau} increased). This effect becomes more visible at the distal extremity of the field of view (distal ROI). This can be explained by the effect of the number of subvolumes n, which reflects the length of the path followed by the microbubbles in the ultrasound field before they reach the ROI. The higher the value of n, the more the microbubble destruction coefficient becomes perceptible. In our study, even with the mechanical index set at 0.08, the effect of destruction induced by the continuous observation on the level of the plateau was clearly perceptible at a flow rate of 100 mL/min, particularly in the distal part of the ultrasound field (Fig 3).

Consequently, the plateau of enhancement depends on the fractional blood volume in the ROI, the length of the feeding vessels subjected to the ultrasound beam before they reach this region, and the flow rate. Thus, when possible, to relate the plateau level to the fractional blood volume, the signal intensity in tissue should be normalized by that in a blood vessel within the ROI (mean intensity in the lumen of a vessel filled by the same feeding vessel and at the same location and depth in the imaging plane). In addition to these factors, the total time needed to reach the plateau increases as the length of feeding vessels exposed to the initial ultrasound destruction pulse increases. This in turn increases the necessary observation duration. The parameter {lambda} in our model describes the fraction of microbubbles destroyed by ultrasound per unit time. The more stable the microbubbles are, the smaller this parameter becomes. Thus, the use of microbubbles that are less sensitive to ultrasound should reduce the sigmoid distortion of the curve and increase the rapidity of arrival time at the plateau.

In our in vivo study, the cortical territories fed by the segmental arteries already included in the imaging and destruction plane exhibited replenishment curves that were comparable to those obtained in the distal ROI of the dialysis cartridge (Figs 4, 5). Since the two rabbits were healthy and were of the same age and weight, the fractional blood volume in the renal cortex of both animals should be comparable. The distances between each group of ROIs were similar in the two rabbits. The initial concentrations of the arriving contrast material were similar as assessed on the basis of the signal intensity at the plateau in the lumen of the segmental artery ROI. Consequently, the higher plateau in the interlobar arterial ROI and the cortex ROI curves for rabbit 1 may indicate a higher flow rate than that for rabbit 2. This would reduce the destruction that occurred in feeding vessels in rabbit 1 and, according to the model and results observed in vitro, should also lead to shortening of the low slope portion of the sigmoid curve in the interlobar artery ROI. Such shortening was observed (Fig 4), but no attempt was made in this study to correlate results with a physical measurement of renal blood flow or with pathologic findings in the kidneys.

Recently, Schlosser et al showed in an in vitro kidney perfusion model that different forms (comparable to our results) are observed for the replenishment curves from hilum vessels and the renal cortex (8). The hilum replenishment exhibited an exponential that asymptotically approached a maximum, but the cortical replenishment showed a sigmoid curve. We consider, as do these authors, that the important difference between the slope they calculated in the macrocirculation (hilum) and the microcirculation (cortex) is mainly a result of physiologic differences of local flow velocities. However, we also believe that the phenomenon of "entrance in the section" described in our study plays an important role and explains the delay and the shape of the curves observed by the authors and by ourselves for cortical replenishment. Furthermore, a voxel in the cortex may be supplied by several in-flow circuits that may bring microbubbles to the voxel at different velocities. Variations in microbubble arrival time due to effects such as these (unrelated to microbubble destruction) could also contribute to the sigmoid-shaped destruction-reperfusion curve. Similar results were obtained in an in vivo brain perfusion study in dogs (5,6). The replenishment curves the authors obtained in the brain cortex showed a sigmoid shape rather than an exponential shape that asymptotically approached a maximum. We believe that feeding vessels may have been included in the imaging plane.

Some limitations of the study should be mentioned. One concern is that if the microbubbles are layered in the electric syringe during injection, distribution of the contrast agent microbubbles may vary from one injection to another. To limit potential consequences of layering, in vitro injection times did not exceed 120 seconds. Furthermore, during each injection, three sets of data were obtained at fixed times (30, 60, and 90 seconds) after the beginning of injection and averaged. In vivo, the replenishment curves were acquired for both rabbits immediately after Doppler intensity equilibrium was reached, which occurred approximately 80 seconds after the beginning of the injection. Thus, because replenishment curves were measured at comparable times after the beginning of injection, the distribution of contrast agent in the syringe should have been comparable. Another consideration is that the measurements in selected ROIs in the dialysis cartridge were not compensated for beam characteristics or attenuation. For estimation of sets of compared curves, ROIs were positioned as similarly as possible to minimize potential effects of attenuation or field differences.

Additional in vivo studies, particularly in organs with well-organized terminal vascularization, such as the myocardium and the kidney, or on the contrary in a disorganized vascular lesion, such as a tumor, must be performed to better appreciate the consequences of the entrance in the section effect on blood flow and fractional blood volume quantification.

In conclusion, we highlighted and modeled an important and previously unexplained aspect of destruction-reperfusion curves. Results of our in vitro experiments show the effect on the initial slope and the plateau of the replenishment curve when the vessels that feed the ROI travel across the ultrasound field. This effect is predicted by our model, which was developed on the basis of indicator-dilution behavior. In vivo, this effect can explain the aspect of replenishment curves observed in the renal cortex. It also adds another complicating factor that must be taken into account for contrast material–enhanced functional US of an organ, especially one with a well-organized vasculature.

Practical application: Findings in this study suggest that it is best to minimize insonification of feeding arteries during destruction-reperfusion measurements in parenchyma and to limit the zone of microbubble destruction to only the studied zone. Such precautions should help determine reliable and reproducible estimates of both the flow rate and the fractional blood volume of a selected organ.


    APPENDIX
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 APPENDIX
 REFERENCES
 
During constant infusion of a US contrast agent, if observation time (t) is short compared with the spontaneous half-life of microbubbles, the refilling evolution of the microbubble concentration C(t) in an ROI can be described by the following equations according to classic indicator-dilution theory with continuous infusion in a single mixed tank:

where Vb represents the volume of blood in the ROI (in liters); F represents the rate of the inflow, which equals that of the outflow (in liters of blood per second); and {lambda} represents the fraction of microbubbles destroyed by the ultrasound beam per unit time (per second), which is assumed to be constant.

C0 is the concentration of microbubbles in blood vessels that enter the ROI (number of microbubbles per liter), which is assumed to be constant with time. This supposes that the blood vessels that enter the ROI were not themselves submitted to the destruction pulse. Consider now several subvolumes in the insonified zone that are oriented serially with respect to the direction of flow. When the microbubbles that come into the subvolume n have already traveled through the ultrasound field in the subvolume n - 1, the number of microbubbles that enter the subvolume n is itself a function of the output function of the subvolume n - 1 (Fig A1a). Evolution of the microbubble concentration in the subvolume n, Cn(t), is given by the following equation:

where 1/{tau} = F/Vb. Cn-1 is the microbubble concentration in subvolume n - 1.



View larger version (6K):
[in this window]
[in a new window]
[Download PPT slide]
 
Figure A1a. (a) Diagram of three subvolumes in the ultrasound field placed serially with respect to the flow direction. Cbo represents the constant blood concentration of the microbubbles in the equilibrium phase that enter the first subvolume. Cn(t) is the time-varying concentration of the microbubbles in the subvolume n. Vb is the volume of blood in each subvolume. (b) Theoretic curves derived from our mathematical model describe the evolution of the microbubble concentration Cbn(t) in the first subvolume (n = 1, top solid and dashed-dotted lines) and third subvolume (n = 3, bottom dashed and dashed-dotted lines). Solid and dashed lines represent the cases where destruction of the microbubbles induced by the observation ultrasound field is null. The dashed-dotted lines show the results when the destruction coefficient {lambda} is 0.25. Other parameters used for calculation were C0 = 500 microbubbles per milliliter, inflow = outflow = 1 mL/sec, and Vb = 1 mL.

 
Consider that at t = 0, Cn = 0 for all values of n except n = 0, where C0 represents the concentration of the microbubbles in blood entering the first subvolume in the ultrasound field. Then the solution of Equation (A2) is the following:

with ß = (1 + {tau}{lambda})/{tau}.

Finally, for the first subvolume, Equation (A3) gives the following equation:

When the destruction induced by the ultrasound between destruction pulses is null ({lambda} = 0, case of intermittent imaging) or considered negligible (case of continuous imaging with low transmit power), Equation (A4) corresponds to the classic exponential function, which asymptotically approaches a maximum, that results from a single-tank model; that model was initially proposed by Wei et al (3).

For subvolumes other than the first (when n > 1),

Equation (A5) describes a curve with a low initial slope that increases secondarily to an inflexion point (sigmoid). Figure A1b shows the microbubble concentration as a function of time after destruction (calculated with Eqq [A4, A5]) in the first subvolume crossed by the microbubbles that entered into the ultrasound field (n = 1) and into the third subvolume (n = 3), which was filled at the same flow rate but was fed by portions of vessels that were emptied by means of the destruction pulse and refilled slowly. This model does not describe the effect of microbubble destruction and flow perpendicular to the image plane.



View larger version (16K):
[in this window]
[in a new window]
[Download PPT slide]
 
Figure A1b. (a) Diagram of three subvolumes in the ultrasound field placed serially with respect to the flow direction. Cbo represents the constant blood concentration of the microbubbles in the equilibrium phase that enter the first subvolume. Cn(t) is the time-varying concentration of the microbubbles in the subvolume n. Vb is the volume of blood in each subvolume. (b) Theoretic curves derived from our mathematical model describe the evolution of the microbubble concentration Cbn(t) in the first subvolume (n = 1, top solid and dashed-dotted lines) and third subvolume (n = 3, bottom dashed and dashed-dotted lines). Solid and dashed lines represent the cases where destruction of the microbubbles induced by the observation ultrasound field is null. The dashed-dotted lines show the results when the destruction coefficient {lambda} is 0.25. Other parameters used for calculation were C0 = 500 microbubbles per milliliter, inflow = outflow = 1 mL/sec, and Vb = 1 mL.

 


    ACKNOWLEDGMENTS
 
We thank Bracco Pharmaceuticals for providing the US contrast material for this study.


    FOOTNOTES
 
See also Science to Practice in this issue.

Abbreviation: ROI = region of interest

Author contributions: Guarantors of integrity of entire study, O.L., S.F.A.; study concepts, O.L., S.F.A., J.M.C., G.B.; study design, O.L., S.F.A., J.M.C.; literature research, O.L., S.F.A.; experimental studies, O.L., S.F.A., J.M.C.; data acquisition, O.L., S.F.A.; data analysis/interpretation, O.L., S.F.A., E.K., S.L.B.; manuscript preparation, O.L., S.L.B.; manuscript definition of intellectual content, O.L., S.L.B., E.K.; manuscript editing, O.L.; manuscript revision/review, O.L., J.M.C., E.K., G.B.; manuscript final version approval, S.F.A., J.M.C., S.L.B., E.K., G.B.


    REFERENCES
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 APPENDIX
 REFERENCES
 

  1. Sirlin CB, Girard MS, Baker KG, Steinbach GC, Deiranieh LH, Mattrey RF. Effect of acquisition rate on liver and portal vein enhancement with microbubble contrast. Ultrasound Med Biol 1999; 25:331-338.[CrossRef][Medline]
  2. Wei K, Skyba DM, Firschke C, Jayaweera AR, Lindner JR, Kaul S. Interactions between microbubbles and ultrasound: in vitro and in vivo observations. J Am Coll Cardiol 1997; 29:1081-1088.[Abstract]
  3. Wei K, Jayaweera AR, Firoozan S, Linka A, Skyba DM, Kaul S. Quantification of myocardial blood flow with ultrasound-induced destruction of microbubbles administered as a constant venous infusion. Circulation 1998; 97:473-483.[Abstract/Free Full Text]
  4. Cosgrove D, Eckersley R, Blomley M, Harvey C. Quantification of blood flow. Eur Radiol 2001; 11:1338-1344.[CrossRef][Medline]
  5. Villanueva FS, Abraham JA, Schreiner GF, et al. Myocardial contrast echocardiography can be used to assess the microvascular response to vascular endothelial growth factor-121. Circulation 2002; 105:759-765.[Abstract/Free Full Text]
  6. Rim SJ, Leong-Poi H, Lindner JR, et al. Quantification of cerebral perfusion with "real-time" contrast-enhanced ultrasound. Circulation 2001; 104:2582-2587.[Abstract/Free Full Text]
  7. Seidel G, Claassen L, Meyer K, Vidal-Langwasser M. Evaluation of blood flow in the cerebral microcirculation: analysis of the refill kinetics during ultrasound contrast agent infusion. Ultrasound Med Biol 2001; 27:1059-1064.[CrossRef][Medline]
  8. Schlosser T, Pohl C, Veltmann C, et al. Feasibility of the flash-replenishment concept in renal tissue: which parameters affect the assessment of the contrast replenishment? Ultrasound Med Biol 2001; 27:937-944.[CrossRef][Medline]
  9. Schwarz KQ, Bezante GP, Chen X. When can Doppler be used in place of integrated backscatter as a measure of scattered ultrasound intensity? Ultrasound Med Biol 1995; 21:231-242.[CrossRef][Medline]
  10. Correas JM, Burns PN, Lai X, Qi X. Infusion versus bolus of an ultrasound contrast agent: in vivo dose-response measurements of BR1. Invest Radiol 2000; 35:72-79.[CrossRef][Medline]
  11. Schneider M, Arditi M, Barrau MB, et al. BR1: a new ultrasonographic contrast agent based on sulfur hexafluoride- filled microbubbles. Invest Radiol 1995; 30:451-457.[Medline]
  12. Linka AZ, Sklenar J, Wei K, Jayaweera AR, Skyba DM, Kaul S. Assessment of transmural distribution of myocardial perfusion with contrast echocardiography. Circulation 1998; 98:1912-1920.[Abstract/Free Full Text]
  13. Seidel G, Meyer K. Harmonic imaging: a new method for the sonographic assessment of cerebral perfusion. Eur J Ultrasound 2001; 14:103-113.[CrossRef][Medline]



This article has been cited by other articles:


Home page
RadiologyHome page
H. J. Paltiel, L. A. Kalish, R. A. Susaeta, F. Frauscher, P. L. O'Kane, and L. G. Freitas-Filho
Pulse-Inversion US Imaging of Testicular Ischemia: Quantitative and Qualitative Analyses in a Rabbit Model
Radiology, June 1, 2006; 239(3): 718 - 729.
[Abstract] [Full Text] [PDF]


Home page
J Ultrasound MedHome page
M. Krix, M.-A. Weber, H. Krakowski-Roosen, H. B. Huttner, S. Delorme, H.-U. Kauczor, and W. Hildebrandt
Assessment of Skeletal Muscle Perfusion Using Contrast-Enhanced Ultrasonography
J. Ultrasound Med., April 1, 2005; 24(4): 431 - 441.
[Abstract] [Full Text] [PDF]


Home page
RadiologyHome page
J. B. Kruskal
Science to Practice: Blood Flow Quantification with Contrast-enhanced US: "Entrance in the Section" Phenomenon--Phantom and Rabbit Study
Radiology, August 1, 2003; 228(2): 298 - 299.
[Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Figures Only
Right arrow Full Text (PDF)
Right arrow All Versions of this Article:
2282020699v1
228/2/473    most recent
Right arrow Submit a response
Right arrow Alert me when this article is cited
Right arrow Alert me when eLetters are posted
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Lucidarme, O.
Right arrow Articles by Berger, G.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Lucidarme, O.
Right arrow Articles by Berger, G.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
RADIOLOGY RADIOGRAPHICS RSNA JOURNALS ONLINE