Diffuse optical tomography (DOT) can be used to localize and characterize a tumor based on endogenous contrast. Besides, using exogenous contrast agents may provide additional information for tumor detection and characterization. For example, it may be possible to use the dynamic behavior of indocyanine green (ICG) for differential diagnosis.1, 2, 3 ICG binds to plasma proteins; thus, it is confined to the vascular compartment and behaves as an intravascular contrast agent.
To date, contrast agents have been used extensively in the magnetic resonance (MR) imaging of tumors including breast cancer.4 Clinical studies have shown that dynamic contrast enhanced–MR imaging (DCE-MRI) has a high sensitivity in breast cancer detection. Nonetheless, the major obstacle on the road to the use of DCE-MRI in breast cancer diagnosis is its low specificity. This is due to the small molecular weight of gadolinium-diethylene-triamine-pentaacetic acid that leaks out of tumor vasculature.5 As a result, some benign lesions may show malignant-type enhancement kinetics and give false positive results because of their high vascularity and high interstitial volume.
This study is motivated by the desire to use a mixture of a small and a large molecular agent to improve specificity in cancer imaging.5 Small molecular agents leak into the extravascular space in tumors, resulting in higher contrast that improves sensitivity. Meanwhile, the large molecular optical agent probes vascular permeability, resulting in improved specificity. In this letter, we report the simultaneous measurement of optical and MR contrast agent kinetics using a frequency domain DOT small animal imaging system with a time resolution of and a 4T MRI system.
We used a R3230 ac-induced small animal breast tumor model in this study. Tumor cells were injected into a Fisher rat that weighed . The tumor size was . We used a frequency domain DOT system that had already been integrated into a 4T MR scanner previously.6 The time resolution of the imaging system was a pivotal feature of the dynamic studies. Our frequency domain DOT system has eight sources and eight detectors7 with a time resolution of per time point. The time resolution of dynamic MR data acquisition was set to . Optical measurements were performed with a laser diode emitting at and modulated at a frequency of . The MR and the DOT measurements were started and performed simultaneously. A mixture of Gd-DTPA and ICG was injected intravenously into the animal after the first three weighted baseline acquisitions. Data were acquired by both systems for .
Afterward, the calibrated optical data were analyzed using the diffusion equation with the Robin boundary conditions. We used the finite element method for the numerical solution of the diffusion equation. The details of the reconstruction scheme were described in Refs. 6, 8. In the optical reconstructions, the first time point was selected as the baseline and then it was subtracted from subsequent dynamic optical data. Coregistration of MR and optical images was accomplished using agar-water–based markers attached to DOT fiber probes.9 These markers were visible on the MR images so that contact points of the fibers on the animal surface could be located. Please note that, in this preliminary study, we did not use any anatomical a priori information in the optical image reconstruction so as not to bias the inverse solution. It was important to take into account the heterogeneity of the tumor when calculating the contrast agent kinetics. The availability of both the and -weighted images made it possible to select a viable and a nonviable region in the tumor. The two main criteria that characterized a viable region were (i) high signal (bright regions) on the MR enhancement image and (ii) normal or low signal (dark regions) in the image.9 Accordingly, a nonviable region of the tumor was defined as a region that did not show signal enhancement. We divided the entire tumor into two regions: one viable and one nonviable region. These ROIs were selected from MR enhancement images corresponding to the optical imaging slice. After the DOT reconstruction was completed, the mean of reconstructed absorption in each ROI at each time point was used to obtain the fractional absorption enhancement curve that was defined as . Similarly, fractional MR signal enhancement was defined .
Figure 1a displays the enhancement image superimposed on the -weighted image of the rat cross section. The ROI that was displayed on the image covered the nonviable region. The rest of the tumor was considered to be viable. Figure 1b shows the reconstructed absorption map at , when maximum enhancement was obtained. The calculated mean was in the viable ROI for the baseline measurement, and it increased to at the peak time point. Meanwhile, mean increased to in the nonviable ROI.
Figure 2a shows the fractional enhancement for the -weighted signal obtained by MRI from the nonviable and viable ROIs. It was clear that the nonviable region exhibited low enhancement, and there was a higher signal enhancement in the viable region. In the viable ROI, Gd-DTPA kinetics showed a rapid rise. This was because small molecular contrast agent leaked out of the blood vessels in the viable region causing a rapid increase in the signal that reached the peak point immediately after the injection. The slow decay seen afterward showed that the Gd-DTPA started diffusing back into the intravascular space.
Optical absorption enhancement kinetics is displayed in Fig. 2b. Absorption kinetics showed a rapid rise in the viable region similar to MRI, and a slow decay was observed. The large molecular size of albumin-ICG made it difficult for the contrast agent to exchange between the extravascular and intravascular spaces. Meanwhile, the kinetics in the nonviable ROI showed low enhancement.
DOT imaging allowed spatially resolved measurement of absorption kinetics within the tumor. The temporal resolution of our system allowed us to measure the fast-rising slope of absorption enhancement and recover the peak value of the optical absorption. It had been shown that reduced scattering did not depend on the concentration of ICG significantly.10 Therefore, we calculated only the absorption kinetics.
The spatial resolution of DOT was low compared to the MRI. As a result, the recovered optical enhancement map did not entirely overlap with the MR enhancement map. This showed that structural a priori information obtained from MRI should be used in optical reconstruction.
Some qualitative features of both contrast agents could be observed from the kinetic curves in the viable ROI. The kinetics of Gd-DTPA displayed the characteristics of a small molecular weight agent, whereas ICG kinetics showed the features of a large molecular agent. In addition, neither agent showed a significant enhancement in the nonviable region. These results were consistent with the DCE-MRI literature.11
Using two contrast agents, one small molecular weight MR and one large molecular weight optical agent, during a single imaging session provides fully coregistered complementary information. This study was the first step toward the goal to use optical imaging in conjunction with the DCE-MRI to improve specificity in tumor characterization.
We thank Dr. Min-Ying Su and Dr. Hon J. Yu for discussions on the subject. This research was supported in part by the National Cancer Institute through Grants No. R21/33 CA-101139, No. R21/33 CA120175, No. R21 CA 121568, and No. U54 CA105480; the California Breast Cancer Research Program No. 12IB-0095; and U.S. Army BC051304.