Coregistration of diffuse optical spectroscopy and magnetic resonance imaging in a rat tumor model

Quantitative In Vivo Monitoring of Hypoxia and Vascularization of Patient-Derived Murine Xenografts of Mantle Cell Lymphoma Using Photoacoustic and Ultrasound Imaging

Tumor oxygenation and vascularization are important parameters that determine the aggressiveness of the tumor and its resistance to cancer therapies. We introduce dual-modality ultrasound and photoacoustic imaging (US-PAI) for the direct, non-invasive real-time in vivo evaluation of oxygenation and vascularization of patient-derived xenografts (PDXs) of B-cell mantle cell lymphomas. The different optical properties of oxyhemoglobin and deoxyhemoglobin make it possible to determine oxygen saturation (sO₂) in tissues using PAI. High-frequency color Doppler imaging enables the visualization of blood flow with high resolution. Tumor oxygenation and vascularization were studied in vivo during the growth of three different subcutaneously implanted patient-derived xenograft (PDX) lymphomas (VFN-M1, VFN-M2 and VFN-M5 R1). Similar values of sO₂ (sO₂ Vital), determined from US-PAI volumetric analysis, were obtained in small and large VFN-M1 tumors ranging from 37.9 ± 2.2 to 40.5 ± 6.0 sO₂ Vital (%) and 37.5 ± 4.0 to 35.7 ± 4.6 sO₂ Vital (%) for small and large VFN-M2 PDXs. In contrast, the higher sO₂ Vital values ranging from 57.1 ± 4.8 to 40.8 ± 5.7 sO₂ Vital (%) (small to large) of VFN-M5 R1 tumors corresponds with the higher aggressiveness of that PDX model. The different tumor percentage vascularization (assessed as micro-vessel areas) of VFN-M1, VFN-M2 and VFN-M5 R1 obtained by color Doppler (2.8 ± 0.1%, 3.8 ± 0.8% and 10.3 ± 2.7%) in large-stage tumors clearly corresponds with their diverse growth and aggressiveness. The data obtained by color Doppler were validated by histology. In conclusion, US-PAI rapidly and accurately provided relevant and reproducible information on tissue oxygenation in PDX tumors in real time without the need for a contrast agent.

Improving vascular imaging with co-planar mutually guided photoacoustic and diffuse optical tomography: a simulation study

Diffuse optical tomography (DOT) and photoacoustic tomography (PAT) are functional imaging modalities that provide absorption coefficient maps of the tissue. Spatial resolution of DOT is relatively low due to light scattering characteristics of the tissue. On the other hand, although PAT can resolve regions of different absorptions with a high spatial resolution, measuring the absolute value of optical absorptions using PAT is challenging due to unknown light fluence distribution in the tissue. Development of image guidance techniques using a priori information of imaging target structure has been shown to increase the accuracy of DOT. PAT is one such method that can be used as a complementary modality to serve as a guide for DOT image reconstruction. On the other hand, estimated fluence map provided by DOT can be used to quantitatively correct PAT images. In this study we introduce a mutually-guided imaging system for fast and simultaneous optical and photoacoustic measurements of tissue absorption map, where DOT is guided by the PAT image and vice versa. Using the obtained absorption map of the tissue, we then estimate the tissue scattering map. We conducted this study using a series of simulations on digital phantoms and demonstrated the effectiveness of the proposed method.

Combined diffuse optical tomography and photoacoustic tomography for enhanced functional imaging of small animals: a methodological study on phantoms

Hybrid imaging methods combining diffuse optical tomography (DOT) and other anatomical or nonoptical functional modalities have been widely investigated to improve imaging performance degraded by the strong optical scattering of biological tissues, through constraining the reconstruction process by prior structures. However, these modalities with different contrast mechanisms may be ineffective in revealing early-staged lesions with high optical contrast but no morphological changes. Photoacoustic tomography (PAT) is particularly useful for visualizing light-absorbing structures embedded in soft tissues with high spatial resolution. Although it is still challenging for PAT to quantitatively disclose the absorption distribution, the modality does provide reliable and specific a priori information differentiating light-absorbing structures of soft tissues and might be more appropriate to guide DOT in lesion diagnosis, as compared with other anatomical or nonoptical functional modalities. In this study, a PAT-guided DOT approach is introduced with both soft- and hard-prior regularizations. The methodology is experimentally validated on small-animal-sized phantoms using a computed-tomography-analogous (CT-analogous) PAT/DOT dual-modality system, focusing on future whole-body applications. The results show that the proposed scheme is capable of effectively improving the quantitative accuracy and spatial resolution of DOT reconstruction.

Combined Diffuse Optical Tomography (DOT) and MRI System for Cancer Imaging in Small Animals

Recently, there has been a great amount of interest in developing multi-modality imaging techniques for oncologic research and clinical studies with the aim of obtaining complementary information and, thus, improving the detection and characterization of tumors. In this present work, the details of a combined MR-diffuse optical imaging system for dual-modality imaging of small animals are given. As a part of this effort, a multi-spectral frequency domain diffuse optical tomography system is integrated with an MRI system. Here, a network analyzer provides the rf modulation signal for the laser diodes and measures the amplitude and the phase of the detected signals. Photomultiplier tubes are utilized to measure low-level signals. The integration of this optical imaging system with the 4T MRI system is realized by incorporating a fiber adaptive interface inside the MR magnet. Coregistration is achieved by a special probe design utilizing fiducial markers. A finite element algorithm is used to solve the diffusion equation and an inverse solver based on this forward solver is implemented to calculate the absorption and scattering maps from the acquired data. The MR a priori information is used to guide the optical reconstruction algorithm. Phantom studies show that the absorption coefficient of a 7 mm inclusion in an irregular object located in 64 mm phantom is recovered with 11% error when MR a priori information is used. ENU induced tumor model is used to test the performance of the system in vivo.

Coregistration of Dynamic Contrast Enhanced MRI and Broadband Diffuse Optical Spectroscopy for Characterizing Breast Cancer

A handheld scanning probe based on broadband Diffuse Optical Spectroscopy (DOS) was used in combination with dynamic contrast enhanced MRI (DCE-MRI) to quantitatively characterize locally-advanced breast cancers in six patients. Measurements were performed sequentially using external fiducial markers for co-registration. Tumor patterns were categorized according to MRI morphological data, and 3D DCE-MRI slices were converted into a volumetric matrix with isotropic voxels to generate views that coincided with the DOS scanning plane. Tumor volume and depth at each DOS measurement site were determined, and a tissue optical index (TOI) that reflects both angiogenic and stromal characteristics was derived from broadband DOS data.

In all six cases, optical scans showed significant TOI contrast corresponding to MRI morphological information. Sharp TOI peaks were recovered for well-circumscribed masses. A reduction in TOI was found inside a tumor with a necrotic center. A broadened peak was observed for a diffuse tumor pattern, and an inflammatory septal case provided two TOI peaks that correlated qualitatively with MRI enhancement. These results provide qualitative confirmation of the common signal origin and complementary information content that can be achieved by combining optical and MR imaging for breast cancer detection and clinical management.

Validity of the semi-infinite tumor model in diffuse reflectance spectroscopy for epithelial cancer diagnosis: a Monte Carlo study

The accurate understanding of optical properties of human tissues plays an important role in the optical diagnosis of early epithelial cancer. Many inverse models used to determine the optical properties of a tumor have assumed that the tumor was semi-infinite, which infers infinite width and length but finite thickness. However, this simplified assumption could lead to large errors for small tumor, especially at the early stages. We used a modified Monte Carlo code, which is able to simulate light transport in a layered tissue model with buried tumor-like targets, to investigate the validity of the semi-infinite tumor assumption in two common epithelial tissue models: a squamous cell carcinoma (SCC) tissue model and a basal cell carcinoma (BCC) tissue model. The SCC tissue model consisted of three layers, i.e. the top epithelium, the middle tumor and the bottom stroma. The BCC tissue model also consisted of three layers, i.e. the top epidermis, the middle tumor and the bottom dermis. Diffuse reflectance was simulated for two common fiber-optic probes. In one probe, both source and detector fibers were perpendicular to the tissue surface; while in the other, both fibers were tilted at 45 degrees relative to the normal axis of the tissue surface. It was demonstrated that the validity of the semi-infinite tumor model depends on both the fiber-optic probe configuration and the tumor dimensions. Two look-up tables, which relate the validity of the semi-infinite tumor model to the tumor width in terms of the source-detector separation, were derived to guide the selection of appropriate tumor models and fiber optic probe configuration for the optical diagnosis of early epithelial cancers.

Dual-contrast dynamic MRI-DOT for small animal imaging

In this paper we present first-of-its-kind spatially resolved enhancement kinetics of optical and magnetic resonance (MR) agents obtained by a combined MR and Diffuse Optical Tomography (MR-DOT) animal imaging system. A unique MR compatible fiber optic interface allows co-registration of MR and DOT data in space and time. High temporal resolution of the hybrid system permits acquisition of data in dynamic mode. Rats bearing a R3230 AC breast cancer tumor model are used for in vivo studies. Thirty-two optical and thirty MR images are acquired during a single imaging session that lasts nearly ten minutes. Both optical, indocyanine green (ICG), and MR contrast agents, gadolinium-DTPA (Gd-DTPA), are injected simultaneously after the acquisition of several baseline frames. Contrast enhancement time curves obtained by MR and DOT systems both indicate higher average enhancement in tumor regions, up to ten-fold for MRI and 3-fold for DOT, compared to close by non-tumor regions. This feasibility study is the first step towards clinical translation of this hybrid imaging platform. The ultimate aim is to use the enhancement kinetics of the optical agent ICG, which binds to plasma proteins, as complementary information to the kinetics of the MR agent Gd-DTPA, a small molecular agent that does not bind to plasma proteins, to better differentiate benign and malignant lesions.

Simultaneous in vivo dynamic magnetic resonance-diffuse optical tomography for small animal imaging

We present simultaneous measurement of enhancement kinetics of an optical and a magnetic resonance (MR) contrast agent in a small animal breast tumor model (R3230 ac) using a combined MR-diffuse optical tomographic (MR-DOT) imaging system. A mixture of a small molecular-weight MR contrast agent gadolinium-diethylene-triamine-pentaacetic acid (Gd-DTPA) and a large molecular-weight optical contrast agent indocyanine green (ICG) was administered intravenously for multimodal dynamic imaging. Coregistration of optical and MR images was accomplished using agar-water-based markers. Using T(2) and dynamic T(1) weighted MR images, we divided the entire tumor into two regions of interest (ROI): a viable and a nonviable region. The absorption enhancements in the ROIs were calculated. An enhancement of the ICG was observed in the viable region. On the contrary, there was a lower enhancement in the nonviable region.

A non-invasive, in vivo technique for monitoring vascular status of glioblastoma during angiogenesis

The growth of solid tumors dependent on the process of angiogenesis in which growth factors secreted by tumor and stromal cells promote endothelial cell proliferation, migration, and maturation. This process generates a tumor-specific vascular supply and enables small or dormant tumors to grow rapidly with exponential increases in tumor volume. Determination of tumor oxygenation at the microvascular level will provide important insight into tumor growth, angiogenesis, necrosis, and therapeutic response, and will facilitate to develop protocols for studying tumor behavior. A non-invasive multi-modality approach based on near infrared spectroscopy (NIRS) technique, namely: Steady State Diffuse Optical Spectroscopy (SSDOS) along with Magnetic Resonance Imaging (MRI) is applied for monitoring the concentration of oxyhemoglobin, deoxyhemoglobin and water within tumor region and for studying the vascular status of tumor and the patho-physiological changes that occur during angiogenesis. Since, the growth of solid tumors depends on the formation of new blood vessels, an association between intramural microvessel density (MVD) and tumor oxygenation is also investigated. The relative decrease in oxygenation value with tumor growth indicates that though blood vessels infiltrate and proliferate the tumor region, a hypoxic trend is clearly present.

A noninvasive multimodal technique to monitor brain tumor vascularization

Determination of tumor oxygenation at the microvascular level will provide important insight into tumor growth, angiogenesis, necrosis and therapeutic response and will facilitate to develop protocols for studying tumor behavior. The non-ionizing near infrared spectroscopy (NIRS) technique has the potential to differentiate lesion and hemoglobin dynamics; however, it has a limited spatial resolution. On the other hand, magnetic resonance imaging (MRI) has achieved high spatial resolution with excellent tissue discrimination but is more susceptible to limited ability to monitor the hemoglobin dynamics. In the present work, the vascular status and the pathophysiological changes that occur during tumor vascularization are studied in an orthotopic brain tumor model. A noninvasive multimodal approach based on the NIRS technique, namely steady state diffuse optical spectroscopy (SSDOS) along with MRI, is applied for monitoring the concentrations of oxyhemoglobin, deoxyhemoglobin and water within tumor region. The concentrations of oxyhemoglobin, deoxyhemoglobin and water within tumor vasculature are extracted at 15 discrete wavelengths in a spectral window of 675-780 nm. We found a direct correlation between tumor size, intratumoral microvessel density and tumor oxygenation. The relative decrease in tumor oxygenation with growth indicates that though blood vessels infiltrate and proliferate the tumor region, a hypoxic trend is clearly present.

Scaling method for fast Monte Carlo simulation of diffuse reflectance spectra from multilayered turbid media

A scaling Monte Carlo method has been developed to calculate diffuse reflectance from multilayered media with a wide range of optical properties in the ultraviolet-visible wavelength range. This multilayered scaling method employs the photon trajectory information generated from a single baseline Monte Carlo simulation of a homogeneous medium to scale the exit distance and exit weight of photons for a new set of optical properties in the multilayered medium. The scaling method is particularly suited to simulating diffuse reflectance spectra or creating a Monte Carlo database to extract optical properties of layered media, both of which are demonstrated in this paper. Particularly, it was found that the root-mean-square error (RMSE) between scaled diffuse reflectance, for which the anisotropy factor and refractive index in the baseline simulation were, respectively, 0.9 and 1.338, and independently simulated diffuse reflectance was less than or equal to 5% for source-detector separations from 200 to 1500 microm when the anisotropy factor of the top layer in a two-layered epithelial tissue model was varied from 0.8 to 0.99; in contrast, the RMSE was always less than 5% for all separations (from 0 to 1500 microm) when the anisotropy factor of the bottom layer was varied from 0.7 to 0.99. When the refractive index of either layer in the two-layered tissue model was varied from 1.3 to 1.4, the RMSE was less than 10%. The scaling method can reduce computation time by more than 2 orders of magnitude compared with independent Monte Carlo simulations.

Extended Kalman Filtering for the Modeling and Analysis of ICG Pharmacokinetics in Cancerous Tumors Using NIR Optical Methods

Compartmental modeling of indocyanine green (ICG) pharmacokinetics, as measured by near infrared (NIR) techniques, has the potential to provide diagnostic information for tumor differentiation. In this paper, we present three different compartmental models to model the pharmacokinetics of ICG in cancerous tumors. We introduce a systematic and robust approach to model and analyze ICG pharmacokinetics based on the extended Kalman filtering (EKF) framework. The proposed EKF framework effectively models multiple-compartment and multiple-measurement systems in the presence of measurement noise and uncertainties in model dynamics. It provides simultaneous estimation of pharmacokinetic parameters and ICG concentrations in each compartment. Moreover, the recursive nature of the Kalman filter estimator potentially allows real-time monitoring of time varying pharmacokinetic rates and concentration changes in different compartments. Additionally, we introduce an information theoretic criteria for the best compartmental model order selection, and residual analysis for the statistical validation of the estimates. We tested our approach using the ICG concentration data acquired from four Fischer rats carrying adenocarcinoma tumor cells. Our study indicates that, in addition to the pharmacokinetic rates, the EKF model may provide parameters that may be useful for tumor differentiation.

Spectral and temporal near-infrared imaging of ex vivo cancerous and normal human breast tissues

Cancerous and normal ex vivo human breast tissues were investigated using spectroscopic and time-sliced two-dimensional (2-D) transillumination imaging methods in order to demonstrate the importance and potential of spectral and temporal measurements in breast cancer detection and diagnosis. The experimental arrangement for time-sliced optical imaging used 120 fs, 1 kHz repetition-rate, 800 nm light pulses from a Ti:sapphire laser system for sample illumination, and a 80 ps resolution ultrafast gated intensified camera system for recording 2-D time-sliced images. The spectroscopic imaging arrangement used 1225-1300 nm tunable output of a Cr: forsterite laser for sample illumination, a Fourier space gate to discriminate against multiple-scattered light, and a near-infrared area camera to record 2-D images. Images recorded with earlier temporal slices of transmitted light highlighted tumors, while those recorded with later slices accentuated normal tissues. When light was tuned closer to the 1203 nm absorption resonance of adipose tissues, a marked enhancement in contrast between the images of adipose and fibrous tissues was observed. A similar wavelength-dependent difference between normal and cancerous tissues was observed. These results correlate well with pathology and nuclear magnetic resonance based analyses of the samples.

Combined diffuse optical spectroscopy and contrast-enhanced magnetic resonance imaging for monitoring breast cancer neoadjuvant chemotherapy: a case study

Monitoring tumor response to therapy can enable assessment of treatment efficacy, maximizing patient outcome and survival. We employ a noninvasive, handheld laser breast scanner (LBS) based on broadband diffuse optical spectroscopy (DOS) in conjunction with contrast-enhanced magnetic resonance imaging (cMRI) to assess tumor response to presurgical neoadjuvant chemotherapy. DOS and cMRI scans are performed after the first and fourth cycles of a doxorubicin/cyclophosphamide regimen in a patient with invasive ductal carcinoma. DOS measurements are used to quantify bulk tissue optical and physiological parameters, which are mapped to T2- and T1-weighted cMRI images. Initial DOS measurements show high tumor/normal contrast in total hemoglobin concentration (THC, 56+/-7 versus 27+/-4 microM) and water fraction (81.4+/-1% versus 24+/-3%) colocalized with regions of strongly enhancing T2-weighted and cMRI signals. After the fourth cycle of chemotherapy, we observe decreases in peak MRI contrast-enhancement values (37.6%) and apparent lesion volume (21.9 versus 13.7 cm3), which corresponds to physiological changes measured by DOS, including a 20 to 25% reduction in the spatial extent of the tumor and a 38.7% drop in mean total hemoglobin content (THC, 41.6 versus 23.4 microM). These data provide in vivo validation of the accuracy of broadband DOS and the sensitivity of optical methods to changes in tumor physiology.

Non-PET functional imaging techniques: optical

The strong and steady development of diffuse optical spectroscopy and tomography as new biomedical optics technologies promises to bring these optical techniques into clinical practice. This article provides a brief review of the light-tissue interaction, the instrumentation, and the theory relevant to this field. This is followed by a survey of the three main applications: brain imaging, muscle imaging, and breast imaging. Lastly, the future outlook of the technology is presented, highlighting the new promises based on recent breakthroughs.

Comparison of water and lipid content measurements using diffuse optical spectroscopy and MRI in emulsion phantoms

We present a quantitative comparison of lipid and water signals obtained from broadband Diffuse Optical Spectroscopy (DOS) and Magnetic Resonance Imaging (MRI). DOS and MRI measurements were performed on an identical set of emulsion phantoms that were composed of different water/soybean oil fractions. Absolute concentrations of water and lipid ranging from 35-94% and 63-6%, respectively were calculated from quantitative broadband near-infrared (NIR) absorption spectra (650-1000 nm). MR images of fat and water were separated using the three-point Dixon technique. DOS and MRI measured water and lipid were highly correlated (R(2) = 0.98 and R(2) = 0.99, respectively) suggesting that these techniques are complementary over a broad range of physiologically relevant water and lipid values. In addition, comparison of DOS derived concentrations to the MRI "gold standard" technique validates our quantitation approach and permits estimation of DOS accuracy and sensitivity in vivo.

Spatial variations in optical and physiological properties of healthy breast tissue

Near-infrared (NIR) diffuse optical spectroscopy (DOS) and diffuse optical imaging (DOI) show promise as noninvasive clinical techniques for breast cancer screening and diagnosis. Since NIR methods are based on optical contrast between healthy and diseased tissue, it is essential to characterize the sources of endogenous contrast in normal subjects. We report intra- and inter-subject variation and bilateral asymmetry of the optical and physiological parameters of 31 women using a seven-wavelength NIR frequency-domain photon migration (FDPM) instrument. Wavelength-dependent absorption and reduced scattering parameters (micro(a) and micro(s'), respectively) were measured in four major quadrants and the areolar regions of left and right breasts. These values were used to determine tissue concentrations of oxy-(HbO(2)) and deoxy-(Hb-R) hemoglobin, lipid content, water concentration, and tissue "scatter power." Mean total hemoglobin for premenopausal (PRE) women (20 to 30 microM) is approximately two-fold higher than for postmenopausal (POST) subjects at all positions. POST women have approximately 50% higher lipid content (50 to 60%) than PRE at all positions. Water concentration on average is 1.8-fold higher for PRE subjects (30 to 40%) than POST. These differences are most pronounced when comparing the areolar complex to the other regions of the breast. In premenopausal women, the areolar regions have 40 to 45% increased total hemoglobin concentration (THC), 20 to 25% lower lipid content, and 30 to 60% higher scatter power versus the quadrants. Small-scale (3 cm) changes in optical properties are negligible compared to large-scale variations over all quadrants, where the intrinsic spatial heterogeneity of healthy breast tissue is 20 to 40% for micro(a) and 5 to 12% for micro(s'). Although no consistent right-left differences are observed in the study population, relative differences between symmetric positions ranged from 18 to 30% for THC, 10 to 40% for adipose, 10 to 25% for water, and 4 to 9% for scattering (674 nm) within an individual.

In vivo quantification of optical contrast agent dynamics in rat tumors by use of diffuse optical spectroscopy with magnetic resonance imaging coregistration

We present a study of the dynamics of optical contrast agents indocyanine green (ICG) and methylene blue (MB) in an adenocarcinoma rat tumor model. Measurements are conducted with a combined frequency-domain and steady-state optical technique that facilitates rapid measurement of tissue absorption in the 650-1000-nm spectral region. Tumors were also imaged by use of contrast-enhanced magnetic resonance imaging (MRI) and coregistered with the location of the optical probe. The absolute concentrations of contrast agent, oxyhemoglobin, deoxyhemoglobin, and water are measured simultaneously each second for approximately 10 min. The differing tissue uptake kinetics of ICG and MB in these late-stage tumors arise from differences in their effective molecular weights. ICG, because of its binding to plasma proteins, behaves as a macromolecular contrast agent with a low vascular permeability. A compartmental model describing ICG dynamics is used to quantify physiologic parameters related to capillary permeability. In contrast, MB behaves as a small-molecular-weight contrast agent that leaks rapidly from the vasculature into the extravascular, extracellular space, and is sensitive to blood flow and the arterial input function.