| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Regular Articles |
Department of Biological Regulation, Weizmann Institute of Science, Rehovot, Israel
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
Overall, the aforementioned vascular and tissue abnormalities lead to the inadequate delivery of diagnostic and therapeutic agents to solid tumors (3) . Therefore, quantifying these abnormalities is of major importance for predicting the efficacy of chemotherapeutic drug delivery, in particular, in breast cancer because neoadjuvant, preoperative treatment is increasingly being used (5 , 6) . Quantification of transport-related parameters in vivo during tumor growth and progression, as well as in response to therapy requires the use of high-resolution noninvasive imaging methods. Although in vivo microscopy provides the desired spatial resolution required to monitor events at the cellular level, this method is depth limited. Dynamic contrast-enhanced magnetic resonance imaging (MRI) effectively provides means to measure the physiological parameters of perfusion in a noninvasive manner (Refs. 7, 8, 9, 10, 11 and references cited therein). This method has been extensively used to diagnose breast and other tumors in preclinical (12 , 13) and clinical (9 , 14) studies. Dynamic acquisition of images following injection of a contrast agent enables the tracking of tracer uptake and clearance over time, as well as providing information about tissue perfusion. In a recent study, the uptake kinetics of Gadolinium-diethylene-triamino-pentaacetic acid (Gd-DTPA) as measured by dynamic contrast-enhanced MRI were shown to correlate with the delivery of the anticancer agent phenyl acetate (15) .
The uptake of diffusible tracers such as Gd-DTPA has been extensively studied, using model-based equations to describe physiological parameters (16, 17, 18) . Here, we used the generalized scheme proposed by Tofts (18) that uses a compartmental model to describe the transfer of the contrast agent from the intravascular to the interstitial space and back, thus yielding the transcapillary influx and outflux transfer constants, respectively. In most previous studies, it was assumed that the transfer constant is the same in both directions; however, disparate influx and outflux transfer constants more accurately describes the dispersal of the contrast agent into and out of the interstitium, particularly when there are differences in diffusion or pressure on either side of the capillary walls (18) . This disparity is more likely to be present in tumors characterized by high interstitial pressure.
In this article, we present data from MRI studies of Gd-DTPA perfusion and diffusion in MDA-MB-231 breast tumors orthotopically inoculated in severe combined immunodeficient (SCID) mice. These tumors exhibited rapid growth accompanied with impaired drainage and formation of cysts and revealed invasive features which includes the formation of metastases. Analyses of dynamic contrast-enhanced MR images of the tumors revealed two distinct patterns of the contrast agent uptake: transcapillary exchange in the cellular regions and diffusion in the cyst. Each process was, therefore, analyzed according to a related model-based equation to yield the appropriate physiological parameters. Monitoring changes in these parameters enabled us to follow alterations in tumor perfusion during growth.
| MATERIALS AND METHODS |
|---|
|
|
|---|
During the experiments, mice were anesthetized by inhalation of 1% Isoflurane (Medeva Pharmaceuticals PA, Inc., Bethlehem, PA) in an O2:N2O (3:7) mixture, applied through a nose cone. After euthanasia, tumors were surgically removed, fixed with 4% formaldehyde, and sectioned in a plane parallel to that of the MRI. Tumors were then embedded in paraffin, sectioned to obtain 4-µm slices, and stained with H&E. All experimental protocols were reviewed and approved by the Institutional Animal Care and Use Committee of the Weizmann Institute of Science.
Contrast Agent Pharmacokinetics.
Gd-DTPA (Schering, Berlin, Germany) was injected as a bolus into the tail vein of the mice, at a dose of 0.5 mmol/kg body weight. Changes in the blood concentration of the contrast agent after injection can be monitored in vivo (19)
or ex vivo (20)
. We used the latter method in a separate group of SCID mice by drawing blood samples from the anesthetized mice by retro-orbital sinus puncture. Measurements of the T1-relaxation times of the plasma samples and calculations of the pharmacokinetic parameters were performed as described previously (12)
. Table 1
summarizes the values of the pharmacokinetic parameters obtained for SCID mice.
|
T1-relaxation times of MDA-MB-231 tumors were measured in vivo by applying a spin-echo imaging sequence with a constant echo time (TE) of 15 ms and a series of variable repetition time (TR) values in a range of 150-5000 ms. T2-relaxation times were measured using a constant TR of 3200 ms and TEs ranging from 18146 ms. Average proton relaxation times in cellular regions of the tumors are summarized in Table 1
.
The protocol of each dynamic study included an initial multislice rapid acquisition with relaxation enhancement (RARE) T2-weighted sequence, with TE/TR of 80/3193 ms; 8 averages; a matrix of 256 x 256; a slice thickness of 1 mm; an interslice thickness of 0.2 mm; and a field of view of 3 x 3 or 5 x 5 cm depending on tumor size. Subsequently, precontrast, T1-weighted, three-dimensional gradient echo images were recorded at the following parameter settings: TE/TR of 4.3/18.3 ms; one acquisition; flip angle 30°; a 256 x 128 x 16 matrix reconstructed to a 256 x 256 x 16 matrix; slice thickness, 1.2 mm; acquisition time of 37 s, and the same field of view as in the T2-weighted images. The precontrast image acquisition was followed by a bolus injection of Gd-DTPA; a series of 3040 postcontrast images were then collected, using the aforementioned three-dimensional T1-weighted protocol.
Image Analysis Using a Kinetic Model of Gd-DTPA Perfusion.
Image analysis of the dynamic contrast-enhanced data were based on the generalized kinetic model developed by Tofts (17
, 18)
. This model predicts the changes in concentration of the contrast agent in the tissue due to its transfer from the plasma into the interstitial space and back, according to Ref. 18
:
![]() |
(kin for short) is the influx volume transfer constant, koutPS
(kout for short) is the outflux volume transfer constant, and
e and
p are the interstitial and intravascular volume fractions, respectively. The term kout/ve = kep is defined as the efflux rate constant from the interstitial space into the plasma. The other parameters are associated with the pharmacokinetic parameters, obtained as described above. From equation A
, it is clear that kout cannot be determined independently; hence, the free parameters are kin, kep (= kout/ ve), and
p. Analysis of the data yielded low intravascular volume fractions,
p of
1% (02%), in most of the pixels in agreement with results obtained in this type of tumor using a macromolecular blood pool contrast agent (21)
. Several assumptions were made in Tofts model to relate the MRI signal enhancement to the changes in the tissue concentration of the contrast agent (18) . One of them concern the effect of water exchange (intravascular to extravascular and extracellular to intracellular) on the determination of contrast agent concentration (22, 23, 24, 25, 26, 27) . We have analyzed the enhancement curves assuming fast exchange for both processes. Additional assessment of the validity of both assumptions was performed by simulations that took into account the exchange rates, using the appropriate equations for a two compartment exchange process (27, 28, 29) , as well as the MRI experimental parameters and typical enhancement ranges (see Supplementary Material). It was found that possible deviations due to the usage of the fast exchange approximation were negligible.
The time course of signal enhancement during the 20 min after bolus injection was analyzed at pixel resolution, using a nonlinear Levenberg-Marquardt least-square fitting algorithm (12)
. The quality of the fit was assessed by calculating a correlation coefficient value R2 for each pixel (12)
. The output provided parametric maps of the influx transfer constant kin and the efflux rate constant kep (= kout/ve). Additional analysis yielded estimated maps of the outflux transfer constant kout using a lower limit for the interstitial volume fraction
e of 0.2 (30)
.
Diffusion Studies.
The diffusion coefficient of Gd-DTPA in the cysts, which formed in the tumors, was estimated by adapting the "null point" method and tracking the displacement of a specified concentration of Gd-DTPA as it diffused through the cyst fluid (31)
. In short, the diffusion of Gd-DTPA in a homogeneous medium is given by the following solution to Ficks diffusion equation (32)
using predetermined initial conditions and boundaries (31)
:
![]() |
, where D is the diffusion coefficient, x is the distance from the start of the tracked route that at time t reached a constant value of concentration Cnull, C0 is an initial steady concentration of Gd-DTPA at x = 0, and erf-1 designates the inverse error function. C0 was estimated as the average concentration of Gd-DTPA at the boundaries of the cyst 5 min after injection. Different routes taken by Gd-DTPA from the edge of the cyst toward its center were tracked, and the time to reach 20% signal enhancement (corresponding to Cnull = 0.03 mM Gd-DTPA) in each pixel along these routes was determined. The distance x of each pixel from the edge of the cyst versus the square root of this time yielded linear plots. For each tumor, 610 plots were created, one for each extracted route. The diffusion coefficient was then determined from the corresponding slope, using equation B . The mean value of the diffusion coefficient was calculated by averaging values obtained for tumors of similar age.
| RESULTS |
|---|
|
|
|---|
3-fold to 0.91 ± 0.14 cm3. Inoculation of a reduced amount of cells (2 x 106) resulted in a
5-fold slower rate of tumor growth. All of the fast-growing tumors exhibited a high intensity central region in the early T2-weighted images (811 days after inoculation), suggesting accumulation of fluid in a cyst-like structure (Fig. 1A)
|
45 days after inoculation of 107 tumor cells into the mammary fat pad of the mice.
Dynamic MRI of Transcapillary Transport.
The transcapillary transfer rates of the contrast agent Gd-DTPA were measured in the tumors over the course of their growth by dynamic contrast-enhanced MRI. Fig. 2
demonstrates the time evolution of Gd-DTPA uptake as seen by signal enhancement post contrast-agent injection in MR images of a fast-growing MDA-MB-231 breast tumor, 8 days after inoculation. Temporal differences in uptake and clearance of the contrast agent were observed throughout the entire tumor and were particularly distinct between the cyst and the viable cellular regions (Fig. 2, C and D)
. The cyst region is clearly delineated in the corresponding T2-weighted image (Fig. 2A)
. The viable areas demonstrated rapid enhancement by the contrast agent followed by a slow decay, whereas the enhancement of the contrast agent seen within the cyst was delayed and far more gradual.
|
0.5 in the viable regions and failed to fit at cyst or nonperfused regions. A typical example of the change in the kin and kout transfer rate constants during tumor growth is illustrated by color-coded parametric images of these constants (Fig. 3)
|
0.5) in the tumor demonstrated a heterogeneous, skewed distribution; hence, statistical analysis was based on calculating the 25th, median, and 75th percentiles of each parameter (Table 2)
0.05) in the values of kin and kout at all percentiles was observed (Table 2
0.05; Table 2
|
|
|
|
t, yielded a linear increase that was in agreement with Ficks Second Law of Diffusion (32)
. Fig. 7C
|
|
| DISCUSSION |
|---|
|
|
|---|
In general, the transcapillary exchange of a tracer molecule between the intravascular and the extracellular spaces occurs by diffusion and convection (3) . To separate the effects of these mechanisms and to quantify the parameters underlying each, it is necessary to track, over time, the changes in the tracer concentration gradient, as well as in the pressure gradient across the capillary walls (33) . Alternatively, tracking the global changes in contrast agent distribution, discriminating inward and outward fluxes, enabled us to monitor changes in the transcapillary delivery of the contrast agent during tumor growth, as well as to discover the disparity of the transcapillary exchange.
The analysis of the data using Tofts model (18)
directly yielded the influx transfer constant of the contrast agent kin, as well as the efflux rate, kep, which is the quotient of the outflux transfer constant, kout, and the extravascular extracellular volume fraction,
e (kep = kout/
e). In most previous studies of dynamic contrast-enhanced MRI, it was generally assumed that the influx and outflux transfer constants were similar (kin
kout); consequently, the two parameters kin and
e could be readily obtained by fitting the data to the proper equations (12
, 17)
. Analysis of our data using this assumption yielded a large fraction of
e < 0.1. These excessively low
e values indicated that the transcapillary rate constants might be different, kin
kout. Microscopic examination and analysis of the histological slices of the tumors revealed high cell density in the viable regions in the MDA-MB-231 tumors; accordingly, we estimated a
e of
0.2, which is the lower limit for this parameter in cancer tissues (30)
. This value was used in our analysis for estimation of the lower limit for kout from the obtained kep maps. Our dynamic contrast-enhanced study, therefore, provided absolute transfer constants for the influx transfer constant, lower limits for the outflux transfer constant, and the degree of disparity among the fluxes. Clearly, to obtain accurate values for kout and
e, it is necessary to use some other method to independently estimate
e at the same resolution as that of the dynamic data.
The significant reduction that we observed in the transcapillary transfer rates in the course of tumor growth can be explained by decreased vascular permeability. Simultaneous measurements for both flow- and permeability-limited tracer revealed that the perfusion of Gd-DTPA is mainly permeability limited (34) . Newly formed microcapillaries in tumors are highly permeable due to primary organization, but remodeling and thickening of the vessels walls leads to less permeable vessels (1) . This progressive process suggests that the duration of tumor growth, rather than tumor size per se, determines the vascular perfusion parameters. The reduction in both the influx and outflux transfer rates as the tumor develops could, therefore, reflect decreased capillary permeability because of vessel maturation.
Interestingly, the decrease in the transfer rates of the contrast agent was not identical for both directions. Rather, the decrease in the influx transfer rate was higher than that of the outflux rate. Consequently, the disparity in the transfer rates increased with tumor growth. A similar disparity was observed in the slow-growing tumors, although they reached smaller size. This suggests that the increase in the disparity is attributed to physiological processes that occur during tumor development rather than tumor size.
We hypothesize that the increased disparity in transfer rates as the tumor grows is because of increased interstitial hypertension as high interstitial pressure will force fluid to reenter the blood vessels thus increasing outflux to influx ratio. It was previously shown that the difference in hydrostatic pressure across the capillary walls is the main mechanism underlying interstitial hypertension (35)
. Moreover, IFP has been shown to increase with tumor size and more specifically to be related to the number of days after tumor implantation (35, 36, 37)
. Additionally, the spatial distribution of the disparity between the constants shown in Fig. 6
is also in accordance with measurements of IFP, showing elevated levels throughout the tumor that drops in the tumor periphery (38)
.
We speculate that staging of the tumors as to their disparity between outflux to influx transcapillary transfer may correlate with their net physiological barrier to transcapillary transfer and hence drug delivery. Additional investigations, however, are necessary to prove the correlation between the disparity in the transcapillary transfer constants and IFP.
It has been shown that a combination of rapid fluid extravasation and poor drainage of fluid within the tumor leads to the formation of cysts and to the development of IFP (30)
. Hence, the formation of cysts is also in agreement with the hypothesis that these tumors may develop high interstitial pressure. The dynamics of the contrast agent uptake in the cyst matched a diffusive behavior, thus enabled direct measurement of the diffusion coefficient of Gd-DTPA in vivo, by tracking the displacement of the contrast agent within the cyst fluid. This in vivo measurement yielded a diffusion coefficient of
1 x 10-5 cm2/s at body temperature. To the best of our knowledge, the Gd-DTPA diffusion coefficient was previously measured only in polyvinyl chloride hydrogel, which serves as a tissue mimic, yielding a diffusion coefficient of 2.6 x 10-6 cm2/s (31)
. The differing diffusion coefficients seen in the cystic fluid and in the hydrogel are in accord with previous studies in which the interstitial diffusion coefficients of various molecules were shown to be about one-third of those in solution (39)
. Hence, the estimated time that it took for Gd-DTPA to diffuse within a plane in the interstitium across a distance L of
200 µm (
pixel dimension) was
32 s (t = L2/4D). This supports our interpretation of the dynamic data in the cellular regions of the tumors (obtained at a temporal resolution of 37 s), which assumed a fast and nearly even distribution of the contrast agent in the interstitial space of each pixel during the acquisition time.
In summary, our results showed that, by the use of dynamic contrast-enhanced MRI, it is possible to distinguish between the influx and outflux transfer constants and that these constants decrease at different rates as the tumor develops, leading to an increase in disparity between the corresponding transfer constants, with the outflux rate exceeding the influx rate. Extending this approach may help predict the efficacy of delivering molecules with properties similar to those of the contrast agents.
| ACKNOWLEDGMENTS |
|---|
| FOOTNOTES |
|---|
The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
Note: Supplementary data for this article can be found at Cancer Research Online (http://cancerres.aacrjournals.org).
Requests for reprints: Hadassa Degani, Department of Biological Regulation, Weizmann Institute of Science, Rehovot 76100, Israel. Phone: 972-8-934-3920; Fax: 972-8-934-4186; E-mail: hadassa.degani{at}weizmann.ac.il
Received 8/26/03. Revised 2/ 5/04. Accepted 2/19/04.
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
M. Dadiani, V. Kalchenko, A. Yosepovich, R. Margalit, Y. Hassid, H. Degani, and D. Seger Real-time Imaging of Lymphogenic Metastasis in Orthotopic Human Breast Cancer Cancer Res., August 15, 2006; 66(16): 8037 - 8041. [Abstract] [Full Text] [PDF] |
||||
![]() |
N. Papo, D. Seger, A. Makovitzki, V. Kalchenko, Z. Eshhar, H. Degani, and Y. Shai Inhibition of tumor growth and elimination of multiple metastases in human prostate and breast xenografts by systemic inoculation of a host defense-like lytic Peptide. Cancer Res., May 15, 2006; 66(10): 5371 - 5378. [Abstract] [Full Text] [PDF] |
||||
![]() |
Y. Hassid, E. Furman-Haran, R. Margalit, R. Eilam, and H. Degani Noninvasive magnetic resonance imaging of transport and interstitial fluid pressure in ectopic human lung tumors. Cancer Res., April 15, 2006; 66(8): 4159 - 4166. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Cancer Research | Clinical Cancer Research |
| Cancer Epidemiology Biomarkers & Prevention | Molecular Cancer Therapeutics |
| Molecular Cancer Research | Cancer Prevention Research |
| Cancer Prevention Journals Portal | Cancer Reviews Online |
| Annual Meeting Education Book | Meeting Abstracts Online |