Methodology development for three-dimensional MR-guided near infrared spectroscopy of breast tumors

Multiwavelength photo-magnetic imaging algorithm improved for direct chromophore concentration recovery using spectral constraints

Multiwavelength photo-magnetic imaging (PMI) is a novel combination of diffuse optics and magnetic resonance imaging, to the best of our knowledge, that yields tissue chromophore concentration maps with high resolution and quantitative accuracy. Here, we present the first experimental results, to the best of our knowledge, obtained using a spectrally constrained PMI image reconstruction method, where chromophore concentration maps are directly recovered, unlike the conventional two-step approach that requires an intermediate step of reconstructing wavelength-dependent absorption coefficient maps. The imposition of the prior spectral information into the PMI inverse problem improves the reconstructed image quality and allows recovery of highly quantitative concentration maps, which are crucial for effective cancer detection and characterization. The obtained results demonstrate the higher performance of the direct reconstruction method. Indeed, the reconstructed concentration maps are not only of higher quality but also obtained approximately 2 times faster than the conventional method.

Performance assessment of MRI guided continuous wave near-infrared spectral tomography for breast imaging

Integration of magnetic resonance imaging (MRI) and near-infrared spectral tomography (NIRST) has yielded promising diagnostic performance for breast imaging in the past. This study focused on whether MRI-guided NIRST can quantify hemoglobin concentration using only continuous wave (CW) measurements. Patients were classified into four breast density groups based on their MRIs. Optical scattering properties were assigned based on average values obtained from these density groups, and MRI-guided NIRST images were reconstructed from calibrated CW data. Total hemoglobin (HbT) contrast between suspected lesions and surrounding normal tissue was used as an indicator of the malignancy. Results obtained from simulations and twenty-four patient cases indicate that the diagnostic power when using only CW data to differentiate malignant from benign abnormalities is similar to that obtained from combined frequency domain (FD) and CW data. These findings suggest that eliminating FD detection to reduce the cost and complexity of MRI-guided NIRST is possible.

Hybrid time-domain and continuous-wave diffuse optical tomography instrument with concurrent, clinical magnetic resonance imaging for breast cancer imaging

Diffuse optical tomography has demonstrated significant potential for clinical utility in the diagnosis and prognosis of breast cancer, and its use in combination with other structural imaging modalities improves lesion localization and the quantification of functional tissue properties. Here, we introduce a hybrid diffuse optical imaging system that operates concurrently with magnetic resonance imaging (MRI) in the imaging suite, utilizing commercially available MR surface coils. The instrument acquires both continuous-wave and time-domain diffuse optical data in the parallel-plate geometry, permitting both absolute assignment of tissue optical properties and three-dimensional tomography; moreover, the instrument is designed to incorporate diffuse correlation spectroscopic measurements for probing tissue blood flow. The instrument is described in detail here. Image reconstructions of a tissue phantom are presented as an initial indicator of the system's ability to accurately reconstruct optical properties and the concrete benefits of the spatial constraints provided by concurrent MRI. Last, we briefly discuss how various data combinations that the instrument could facilitate, including tissue perfusion, can enable more comprehensive assessment of lesion physiology.

Direct Regularization From Co-Registered Contrast MRI Improves Image Quality of MRI-Guided Near-Infrared Spectral Tomography of Breast Lesions

An approach using direct regularization from co-registered dynamic contrast enhanced magnetic reson- ance images was used to reconstruct near-infrared spectral tomography patient images, which does not need image segmentation. 20 patients with mammography/ultrasound confirmed breast abnormalities were involved in this paper, and the resulting images indicated that tumor total hemoglobin concentration contrast differentiated malignant from benign cases (p-value = 0.021). The approach prod- uced reconstructed images, which significantly reduced surface artifacts near the source-detector locations (p-value = 4.16e-6).

Near-Infrared Diffuse Optical Imaging for Early Prediction of Breast Cancer Response to Neoadjuvant Chemotherapy: A Comparative Study Using 18F-FDG PET/CT

Diffuse optical spectroscopic imaging (DOSI) is used as an indicator of tumor blood volume quantified by tissue hemoglobin concentrations. We aimed to determine whether early changes in tumor total hemoglobin (tHb) concentration can predict a pathologic complete response (pCR) to neoadjuvant chemotherapy in patients with operable breast cancer, and we compared the predictive value of pCR between DOSI and (18)F-FDG PET combined with CT.

METHODS

Of the 100 patients enrolled, 84 patients were prospectively evaluated for primary objective analysis. Sixty-four of the patients underwent both sequential DOSI scans at baseline after their first and second chemotherapy courses and (18)F-FDG PET/CT at baseline and after their second chemotherapy course. The mean tHb (tHbmean) concentration and SUVmax of the lesion were measured using DOSI and (18)F-FDG PET/CT, respectively, and the percentage change in tHbmean (∆tHbmean) and change in SUVmax (∆SUVmax) were calculated. We compared the diagnostic performances of DOSI and (18)F-FDG PET/CT for predicting pCR via the analysis of the receiver-operating-characteristic curves.

RESULTS

pCR was achieved in 16 patients, and neoadjuvant chemotherapy caused a significant reduction of ∆tHbmean in pCR compared with non-pCR after the 2 chemotherapy courses. When the tentative ∆tHbmean cutoff values after the first and second courses were used, the ability to predict pCR was as follows: 81.2% sensitivity/47.0% specificity and 93.7% sensitivity/47.7% specificity, respectively. Comparison of the diagnostic performances of DOSI and (18)F-FDG PET/CT revealed areas under the curve of 0.69 and 0.75 of ∆tHbmean after the first and second courses, respectively, which were lower than those of ∆SUVmax (0.90).

CONCLUSION

DOSI predicted pCR in patients with breast cancer with moderate accuracy. The diagnostic performance of DOSI was inferior to that of the early metabolic response as monitored by (18)F-FDG PET/CT.

Direct regularization from co-registered anatomical images for MRI-guided near-infrared spectral tomographic image reconstruction

Combining anatomical information from high resolution imaging modalities to guide near-infrared spectral tomography (NIRST) is an efficient strategy for improving the quality of the reconstructed spectral images. A new approach for incorporating image information directly into the inversion matrix regularization was examined using Direct Regularization from Images (DRI), which encodes the gray-scale data into the NIRST image reconstruction problem. This process has the benefit of eliminating user intervention such as image segmentation of distinct regions. Specifically, the Dynamic Contrast Enhanced Magnetic Resonance (DCE-MR) image intensity value differences within the anatomical image were used to implement an exponentially-weighted regularization function between the image pixels. The algorithm was validated using simulated reconstructions with noise, and the results showed that spatial resolution and robustness of the reconstructed images were significantly improved by appropriate choice of the regularization weight parameters. The proposed approach was also tested on in vivo breast data acquired in a recent clinical trial combining NIRST / MRI for cancer tumor characterization. Relative to the standard "no priors" diffuse recovery, the contrast of the tumor to the normal surrounding tissue increased from 2.4 to 3.6, and the difference between the tumor size segmented from DCE-MR images and reconstructed optical images decreased from 18% to 6%, while there was an overall decrease in surface artifacts.

MR-Guided Near-Infrared Spectral Tomography Increases Diagnostic Performance of Breast MRI

PURPOSE

The purpose of this study was to determine the diagnostically most important molecular biomarkers quantified by magnetic resonance-guided (MR) near-infrared spectral tomography (NIRST) that distinguish malignant breast lesions from benign abnormalities when combined with outcomes from clinical breast MRI.

EXPERIMENTAL DESIGN

The study was HIPAA compliant and approved by the Dartmouth Institutional Review Board, the NIH, the United States State Department, and Xijing Hospital. MR-guided NIRST evaluated hemoglobin, water, and lipid content in regions of interest defined by concurrent dynamic contrast-enhanced MRI (DCE-MRI) in the breast. MRI plus NIRST was performed in 44 subjects (median age, 46, age range, 20-81 years), 28 of whom had subsequent malignant pathologic diagnoses, and 16 had benign conditions. A subset of 30 subject examinations yielded optical data that met minimum sensitivity requirements to the suspicious lesion and were included in the analyses of diagnostic performance.

RESULTS

In the subset of 30 subject examinations meeting minimum optical data sensitivity criterion, the MR-guided NIRST separated malignant from benign lesions using total hemoglobin (HbT; P < 0.01) and tissue optical index (TOI; P < 0.001). Combined MRI plus TOI data caused one false positive and 1 false negative, and produced the best diagnostic performance, yielding an AUC of 0.95, sensitivity of 95%, specificity of 89%, positive predictive value of 95%, and negative predictive value of 89%, respectively.

CONCLUSIONS

MRI plus NIRST results correlated well with histopathologic diagnoses and could provide additional information to reduce the number of MRI-directed biopsies.

Sensitivity of MRI-guided near-infrared spectroscopy clinical breast exam data and its impact on diagnostic performance

In this study, data from breast MRI-guided near infrared spectroscopy (NIRS) exams delivered to 44 patients scheduled for surgical resection (ending in 16 benign and 28 malignant diagnoses) were analyzed using a spatial sensitivity metric to quantify the adequacy of the optical measurements for interrogating the tumor region of interest, as derived from the concurrent MRI scan. Along with positional sensitivity, the incorporation of spectral priors and the selection of an appropriate regularization parameter in the image reconstruction were considered, and found to influence the diagnostic accuracy of the recovered images. Once optimized, the MRI/NIRS data was able to differentiate the malignant from benign lesions through both total hemoglobin (p = 0.0037) and tissue optical index (p = 0.00019), but required the relative spatial sensitivity of the optical measurement data to each lesion to be above 1%. Spectral constraints implemented during the reconstruction were required to obtain statistically significant diagnostic information from images of H2O, lipids, and Tissue Optical Index (TOI). These results confirm the need for optical systems that have homogenous spatial coverage of the breast while still being able to accommodate the normal range of breast sizes.

Adaptable near-infrared spectroscopy fiber array for improved coupling to different breast sizes during clinical MRI

RATIONALE AND OBJECTIVES

Near-infrared spectroscopy (NIRS) of breast can provide functional information on the vascular and structural compartments of tissues in regions identified during simultaneous magnetic resonance imaging (MRI). NIRS can be acquired during dynamic contrast-enhanced MRI (DCE-MRI) to accomplish image-guided spectroscopy of the enhancing regions, potentially increasing the diagnostic specificity of the examination and reducing the number of biopsies performed as a result of inconclusive MRI breast imaging studies.

MATERIALS AND METHODS

We combine synergistic attributes of concurrent DCE-MRI and NIRS with a new design of the clinical NIRS breast interface that couples to a standard MR breast coil and allows imaging of variable breast sizes. Spectral information from healthy volunteers and cancer patients is recovered, providing molecular information in regions defined by the segmented MR image volume.

RESULTS

The new coupling system significantly improves examination utility by allowing improved coupling of the NIR fibers to breasts of all cup sizes and lesion locations. This improvement is demonstrated over a range of breast sizes (cup size A through D) and normal tissue heterogeneity using a group of eight healthy volunteers and two cancer patients. Lesions located in the axillary region and medial-posterior breast are now accessible to NIRS optodes. Reconstructed images were found to have biologically plausible hemoglobin content, oxygen saturation, and water and lipid fractions.

CONCLUSIONS

In summary, a new NIRS/MRI breast interface was developed to accommodate the variation in breast sizes and lesion locations that can be expected in clinical practice. DCE-MRI-guided NIRS quantifies total hemoglobin, oxygenation, and scattering in MR-enhancing regions, increasing the diagnostic information acquired from MR examinations.

Hybrid photomultiplier tube and photodiode parallel detection array for wideband optical spectroscopy of the breast guided by magnetic resonance imaging

A new optical parallel detection system of hybrid frequency and continuous-wave domains was developed to improve the data quality and accuracy in recovery of all breast optical properties. This new system was deployed in a previously existing system for magnetic resonance imaging (MRI)-guided spectroscopy, and allows incorporation of additional near-infrared wavelengths beyond 850 nm, with interlaced channels of photomultiplier tubes (PMTs) and silicon photodiodes (PDs). The acquisition time for obtaining frequency-domain data at six wavelengths (660, 735, 785, 808, 826, and 849 nm) and continuous-wave data at three wavelengths (903, 912, and 948 nm) is 12 min. The dynamic ranges of the detected signal are 105 and 106 for PMT and PD detectors, respectively. Compared to the previous detection system, the SNR ratio of frequency-domain detection was improved by nearly 103 through the addition of an RF amplifier and the utilization of programmable gain. The current system is being utilized in a clinical trial imaging suspected breast cancer tumors as detected by contrast MRI scans.

Towards non-invasive characterization of breast cancer and cancer metabolism with diffuse optics

We review recent developments in diffuse optical imaging and monitoring of breast cancer, i.e. optical mammography. Optical mammography permits non-invasive, safe and frequent measurement of tissue hemodynamics oxygen metabolism and components (lipids, water, etc.), the development of new compound indices indicative of the risk and malignancy, and holds potential for frequent non-invasive longitudinal monitoring of therapy progression.

Combined hemoglobin and fluorescence diffuse optical tomography for breast tumor diagnosis: a pilot study on time-domain methodology

A combined time-domain fluorescence and hemoglobin diffuse optical tomography (DOT) system and the image reconstruction methods are proposed for enhancing the reliability of breast-dedicated optical measurement. The system equipped with two pulsed laser diodes at wavelengths of 780 nm and 830 nm that are specific to the peak excitation and emission of the FDA-approved ICG agent, and works with a 4-channel time-correlated single photon counting device to acquire the time-resolved distributions of the light re-emissions at 32 boundary sites of tissues in a tandem serial-to-parallel mode. The simultaneous reconstruction of the two optical (absorption and scattering) and two fluorescent (yield and lifetime) properties are achieved with the respective featured-data algorithms based on the generalized pulse spectrum technique. The performances of the methodology are experimentally assessed on breast-mimicking phantoms for hemoglobin- and fluorescence-DOT alone, as well as for fluorescence-guided hemoglobin-DOT. The results demonstrate the efficacy of improving the accuracy of hemoglobin-DOT based on a priori fluorescence localization.

Automatic and robust calibration of optical detector arrays for biomedical diffuse optical spectroscopy

The design and testing of a new, fully automated, calibration approach is described. The process was used to calibrate an image-guided diffuse optical spectroscopy system with 16 photomultiplier tubes (PMTs), but can be extended to any large array of optical detectors and associated imaging geometry. The design goals were accomplished by developing a routine for robust automated calibration of the multi-detector array within 45 minutes. Our process was able to characterize individual detectors to a median norm of the residuals of 0.03 V for amplitude and 4.4 degrees in phase and achieved less than 5% variation between all the detectors at the 95% confidence interval for equivalent measurements. Repeatability of the calibrated data from the imaging system was found to be within 0.05 V for amplitude and 0.2 degrees for phase, and was used to evaluate tissue-simulating phantoms in two separate imaging geometries. Spectroscopic imaging of total hemoglobin concentration was recovered to within 5% of the true value in both cases. Future work will focus on streamlining the technology for use in a clinical setting with expectations of achieving accurate quantification of suspicious lesions in the breast.

NIRViz: 3D visualization software for multimodality optical imaging using visualization toolkit (VTK) and insight segmentation toolkit (ITK)

Optical imaging using near-infrared light is used for noninvasive probing of tissues to recover vascular and molecular status of healthy and diseased tissues using hemoglobin contrast arising due to absorption of light. While multimodality optical techniques exist, visualization techniques in this area are limited. Addressing this issue, we present a simple framework for image overlay of optical and magnetic resonance (MRI) or computerized tomographic images which is intuitive and easily usable, called NIRViz. NIRViz is a multimodality software platform for the display and navigation of Digital Imaging and Communications in Medicine (DICOM) MRI datasets and 3D optical image solutions geared toward visualization and coregistration of optical contrast in diseased tissues such as cancer. We present the design decisions undertaken during the design of the software, the libraries used in the implementation, and other implementation details as well as preliminary results from the software package. Our implementation uses the Visualization Toolkit library to do most of the work, with a Qt graphical user interface for the front end. Challenges encountered include reslicing DICOM image data and coregistration of image space and mesh space. The resulting software provides a simple and customized platform to display surface and volume meshes with optical parameters such as hemoglobin concentration, overlay them on magnetic resonance images, allow the user to interactively change transparency of different image sets, rotate geometries, clip through the resulting datasets, obtain mesh and optical solution information, and successfully interact with both functional and structural medical image information.

Remote positioning optical breast magnetic resonance coil for slice-selection during image-guided near-infrared spectroscopy of breast cancer

The design and testing of a pneumatic optical positioning interface produced with the goal of improving fiber positioning in magnetic resonance (MR)-guided diffuse spectral imaging of breast cancer is presented. The system was created for vertical positioning of optical fibers inside the MR bore during a patient exam to target suspicious lesions with MR scans for reference and collect multiple planes of optical data. The interface includes new fiber plates for mechanical and optical coupling to the breast, and was tested in phantoms and human imaging. Reconstructions with data taken in the new interface show acceptable linearity over different absorber concentrations (residual norm = 0.067), and exhibit good contrast recovery at different imaging planes, which is consistent with previous work. An example of human breast imaging through the new interface is shown and a discussion of how it compares to other patient interfaces for breast imaging is presented. Design goals of increasing the available degrees of freedom for fiber positioning while maintaining good patient-fiber contact and comfort were accomplished. This interface allows improved volumetric imaging with interactive and accurate slice selection to quantify targeted suspicious lesions.

Microwave dielectric contrast imaging in a magnetic resonant environment and the effect of using magnetic resonant spatial information in image reconstruction

Microwave Tomography (MT) can determine the permittivity and conductivity of a volume of interest; it has been shown that a contrast exists between these electrical properties in healthy and malignant tissues, and MT can be used to discern the dielectric contrast image of these tissues by recovering their electrical property values. Simulation and phantom experiments of objects with known spatial locations have shown that using boundary information derived from internal structures in the imaged volume greatly increases the accuracy of the recovered property values. In practice this spatial information, which will be used for reconstructing the tissue's electrical property images, must be determined with high enough resolution to segment boundary regions and internal structures of interest. This experiment investigates the use of Magnetic Resonant Imaging (MRI) in obtaining the desired spatial information used in mesh generation for image reconstruction and provides microwave image results comparing electrical properties recovered with and without this prior spatial information.

MR water quantitative priors improves the accuracy of optical breast imaging

Magnetic resonance (MR) guided optical breast imaging is a promising modality to improve the specificity of breast imaging, because it provides high-resolution quantitative maps of total hemoglobin, oxygen saturation, water content, and optical scattering. These properties have been shown to distinguish malignant from benign lesions. However, the optical detection hardware required for deep tissue imaging has poor spectral sensitivity which limits accurate water quantification; this reduces the accuracy of hemoglobin quantification. We present a methodology to improve optical quantification by utilizing the ability of Dixon MR imaging to quantitatively estimate water and fat; this technique effectively reduces optical crosstalk between water and oxyhemoglobin. The techniques described in this paper reduce hemoglobin quantification error by as much as 38%, as shown in a numerical phantom, and an experimental phantom. Error is reduced by as much 20% when imperfect MR water quantification is given. These techniques may also increase contrast between diseased and normal tissue, as shown in breast tissue in vivo. It is also shown that using these techniques may permit fewer wavelengths to be used with similar quantitative accuracy, enabling higher temporal resolution. In addition, it is shown that these techniques can improve the ability of MRI to quantify water in the presence of bias in the Dixon water/fat separation.

Image guided near-infrared spectroscopy of breast tissue in vivo using boundary element method

We demonstrate quantitative functional imaging using image-guided near-infrared spectroscopy (IG-NIRS) implemented with the boundary element method (BEM) for reconstructing 3-D optical property estimates in breast tissue in vivo. A multimodality MRI-NIR system was used to collect measurements of light reflectance from breast tissue. The BEM was used to model light propagation in 3-D based only on surface discretization in order to reconstruct quantitative values of total hemoglobin (HbT), oxygen saturation, water, and scatter. The technique was validated in experimental measurements from heterogeneous breast-shaped phantoms with known values and applied to a total of seven subjects comprising six healthy individuals and one participant with cancer imaged at two time points during neoadjuvant chemotherapy. Using experimental measurements from a heterogeneous breast phantom, BEM for IG-NIRS produced accurate values for HbT in the inclusion with a <3% error. Healthy breast tissues showed higher HbT and water in fibroglandular tissue than in adipose tissue. In a subject with cancer, the tumor showed higher HbT compared to the background. HbT in the tumor was reduced by 9 μM during treatment. We conclude that 3-D MRI-NIRS with BEM provides quantitative and functional characterization of breast tissue in vivo through measurement of hemoglobin content. The method provides potentially complementary information to DCE-MRI for tumor characterization.

A coupled finite element-boundary element method for modeling Diffusion equation in 3D multi-modality optical imaging

Three dimensional image reconstruction for multi-modality optical spectroscopy systems needs computationally efficient forward solvers with minimum meshing complexity, while allowing the flexibility to apply spatial constraints. Existing models based on the finite element method (FEM) require full 3D volume meshing to incorporate constraints related to anatomical structure via techniques such as regularization. Alternate approaches such as the boundary element method (BEM) require only surface discretization but assume homogeneous or piece-wise constant domains that can be limiting. Here, a coupled finite element-boundary element method (coupled FE-BEM) approach is demonstrated for modeling light diffusion in 3D, which uses surfaces to model exterior tissues with BEM and a small number of volume nodes to model interior tissues with FEM. Such a coupled FE-BEM technique combines strengths of FEM and BEM by assuming homogeneous outer tissue regions and heterogeneous inner tissue regions. Results with FE-BEM show agreement with existing numerical models, having RMS differences of less than 0.5 for the logarithm of intensity and 2.5 degrees for phase of frequency domain boundary data. The coupled FE-BEM approach can model heterogeneity using a fraction of the volume nodes (4-22%) required by conventional FEM techniques. Comparisons of computational times showed that the coupled FE-BEM was faster than stand-alone FEM when the ratio of the number of surface to volume nodes in the mesh (N(s)/N(v)) was less than 20% and was comparable to stand-alone BEM ( ± 10%).

Compositional-prior-guided image reconstruction algorithm for multi-modality imaging

The development of effective multi-modality imaging methods typically requires an efficient information fusion model, particularly when combining structural images with a complementary imaging modality that provides functional information. We propose a composition-based image segmentation method for X-ray digital breast tomosynthesis (DBT) and a structural-prior-guided image reconstruction for a combined DBT and diffuse optical tomography (DOT) breast imaging system. Using the 3D DBT images from 31 clinically measured healthy breasts, we create an empirical relationship between the X-ray intensities for adipose and fibroglandular tissue. We use this relationship to then segment another 58 healthy breast DBT images from 29 subjects into compositional maps of different tissue types. For each breast, we build a weighted-graph in the compositional space and construct a regularization matrix to incorporate the structural priors into a finite-element-based DOT image reconstruction. Use of the compositional priors enables us to fuse tissue anatomy into optical images with less restriction than when using a binary segmentation. This allows us to recover the image contrast captured by DOT but not by DBT. We show that it is possible to fine-tune the strength of the structural priors by changing a single regularization parameter. By estimating the optical properties for adipose and fibroglandular tissue using the proposed algorithm, we found the results are comparable or superior to those estimated with expert-segmentations, but does not involve the time-consuming manual selection of regions-of-interest.

Diffuse Optics for Tissue Monitoring and Tomography

This review describes the diffusion model for light transport in tissues and the medical applications of diffuse light. Diffuse optics is particularly useful for measurement of tissue hemodynamics, wherein quantitative assessment of oxy- and deoxy-hemoglobin concentrations and blood flow are desired. The theoretical basis for near-infrared or diffuse optical spectroscopy (NIRS or DOS, respectively) is developed, and the basic elements of diffuse optical tomography (DOT) are outlined. We also discuss diffuse correlation spectroscopy (DCS), a technique whereby temporal correlation functions of diffusing light are transported through tissue and are used to measure blood flow. Essential instrumentation is described, and representative brain and breast functional imaging and monitoring results illustrate the workings of these new tissue diagnostics.

Sensitivity of hemoglobin concentration on optical probe positioning in image-guided near infrared spectroscopy

Multi-modality image-guided near infrared spectroscopy provides volume-based quantification of total hemoglobin, oxygen saturation, water and scatter in various tissue types in-vivo. The accuracy of these parameters depends on the location of the imaging probe and its distance from the tumor. In a numerical study, we have performed simulation to analyze this effect in a breast-specific imaging domain. Results show that the accuracy of total hemoglobin decreases by 25% for every centimeter away from the tumor center, for a typical size of 35mm cancer in the breast. Image guidance is necessary for accurate positioning, and multi-plane acquisition can improve this accuracy.