Bruce J. Tromberg

BOUNDARY-CONDITIONS FOR THE DIFFUSION EQUATION IN RADIATIVE-TRANSFER

Using the method of images, we examine the three boundary conditions commonly applied to the surface of a semi-infinite turbid medium. We find that the image-charge configurations of the partial-current and extrapolated-boundary conditions have the same dipole and quadrupole moments and that the two corresponding solutions to the diffusion equation are approximately equal. In the application of diffusion theory to frequency-domain photon-migration (FDPM) data, these two approaches yield values for the scattering and absorption coefficients that are equal to within 3%. Moreover, the two boundary conditions can be combined to yield a remarkably simple, accurate, and computationally fast method for extracting values for optical parameters from FDPM data. FDPM data were taken both at the surface and deep inside tissue phantoms, and the difference in data between the two geometries is striking. If one analyzes the surface data without accounting for the boundary, values deduced for the optical coefficients are in error by 50% or more. As expected, when aluminum foil was placed on the surface of a tissue phantom, phase and modulation data were closer to the results for an infinite-medium geometry. Raising the reflectivity of a tissue surface can, in principle, eliminate the effect of the boundary. However, we find that phase and modulation data are highly sensitive to the reflectivity in the range of 80-100%, and a minimum value of 98% is needed to mimic an infinite-medium geometry reliably. We conclude that noninvasive measurements of optically thick tissue require a rigorous treatment of the tissue boundary, and we suggest a unified partial-current--extrapolated boundary approach.

Imaging cells and extracellular matrix in vivo by using second-harmonic generation and two-photon excited fluorescence

Multiphoton microscopy relies on nonlinear light–matter interactions to provide contrast and optical sectioning capability for high-resolution imaging. Most multiphoton microscopy studies in biological systems have relied on two-photon excited fluorescence (TPEF) to produce images. With increasing applications of multiphoton microscopy to thick-tissue “intravital” imaging, second-harmonic generation (SHG) from structural proteins has emerged as a potentially important new contrast mechanism. However, SHG is typically detected in transmission mode, thus limiting TPEF/SHG coregistration and its practical utility for in vivo thick-tissue applications. In this study, we use a broad range of excitation wavelengths (730–880 nm) to demonstrate that TPEF/SHG coregistration can easily be achieved in unstained tissues by using a simple backscattering geometry. The combined TPEF/SHG technique was applied to imaging a three-dimensional organotypic tissue model (RAFT). The structural and molecular origin of the image-forming signal from the various tissue constituents was determined by simultaneous spectroscopic measurements and confirming immunofluorescence staining. Our results show that at shorter excitation wavelengths (<800 nm), the signal emitted from the extracellular matrix (ECM) is a combination of SHG and TPEF from collagen, whereas at longer excitation wavelengths the ECM signal is exclusively due to SHG. Endogenous cellular signals are consistent with TPEF spectra of cofactors NAD(P)H and FAD at all excitation wavelengths. The reflected SHG intensity follows a quadratic dependence on the excitation power, decays exponentially with depth, and exhibits a spectral dependence in accordance with previous theoretical studies. The use of SHG and TPEF in combination provides complementary information that allows noninvasive, spatially localized in vivo characterization of cell–ECM interactions in unstained thick tissues.

In vivo absorption, scattering, and physiologic properties of 58 malignant breast tumors determined by broadband diffuse optical spectroscopy

Diffuse optical imaging (DOI) may be a beneficial diagnostic method for women with mammographically dense breast tissue. In order to evaluate the utility of DOI, we are developing broadband diffuse optical spectroscopy (DOS) to characterize the functional origins of optical signals in breast cancer patients. Broadband DOS combines multifrequency intensity-modulated and continuous-wave near-infrared light to quantify tissue absorption and scattering spectra from 650 to 1000 nm. Values of intrinsic physiological properties (oxy- and deoxy-hemoglobin, water, lipid, and scatter power) derived from absorption and scattering spectra provide detailed information on breast physiology. We present the results of clinical studies of 58 stage II/III malignant breast tumors using a noninvasive, handheld, broadband DOS probe. On average, eight positions were scanned over tumor and contralateral normal breast for each subject. Intrinsic physiological properties were statistically significantly different for malignant vs. normal tissues for all subjects, without patient age or tumor size/type stratification. Breast tissues containing malignant tumors displayed reduced lipid content ( approximately 20%) and increased water, deoxy-, and oxy-hemoglobin (>50% each) compared to normal breast tissues. Functional perturbations by the tumor were significantly larger than functional variations in normal tissues. A tissue optical index (TOI) derived from intrinsic physiological properties yielded an average two-fold contrast difference between malignant tumors and intrinsic tissue properties. Our results demonstrate that intrinsic optical signals can be influenced by functional perturbations characteristic of malignant transformation; cellular metabolism, extracellular matrix composition, and angiogenesis. Our findings further underscore the importance of broadband measurements and patient age stratification in breast cancer DOI.

In vivo local determination of tissue optical properties: applications to human brain

Local and superficial near-infrared (NIR) optical-property characterization of turbid biological tissues can be achieved by measurement of spatially resolved diffuse reflectance at small source-detector separations (<1.4 mm). However, in these conditions the inverse problem, i.e., calculation of localized absorption and the reduced scattering coefficients, is necessarily sensitive to the scattering phase function. This effect can be minimized if a new parameter of the phase function gamma, which depends on the first and the second moments of the phase function, is known. If gamma is unknown, an estimation of this parameter can be obtained by the measurement, but the uncertainty of the absorption coefficient is increased. A spatially resolved reflectance probe employing multiple detector fibers (0.3-1.4 mm from the source) is described. Monte Carlo simulations are used to determine gamma, the reduced scattering and absorption coefficients from reflectance data. Probe performance is assessed by measurements on phantoms, the optical properties of which were measured by other techniques [frequency domain photon migration (FDPM) and spatially resolved transmittance]. Our results show that changes in the absorption coefficient, the reduced scattering coefficient, and gamma can be measured to within +/-0.005 mm(-1), +/-0.05 mm(-1), and +/-0.2, respectively. In vivo measurements performed intraoperatively on a human skull and brain are reported for four NIR wavelengths (674, 811, 849, 956 nm) when the spatially resolved probe and FDPM are used. The spatially resolved probe shows optimum measurement sensitivity in the measurement volume immediately beneath the probe (typically 1 mm(3) in tissues), whereas FDPM typically samples larger regions of tissues. Optical-property values for human skull, white matter, scar tissue, optic nerve, and tumors are reported that show distinct absorption and scattering differences between structures and a dependence on the phase-function parameter gamma.

Noninvasive functional optical spectroscopy of human breast tissue

Near infrared diffuse optical spectroscopy and diffuse optical imaging are promising methods that eventually may enhance or replace existing technologies for breast cancer screening and diagnosis. These techniques are based on highly sensitive, quantitative measurements of optical and functional contrast between healthy and diseased tissue. In this study, we examine whether changes in breast physiology caused by exogenous hormones, aging, and fluctuations during the menstrual cycle result in significant alterations in breast tissue optical contrast. A noninvasive quantitative diffuse optical spectroscopy technique, frequency-domain photon migration, was used. Measurements were performed on 14 volunteer subjects by using a hand-held probe. Intrinsic tissue absorption and reduced scattering parameters were calculated from frequency-domain photon migration data. Wavelength-dependent absorption (at 674, 803, 849, and 956 nm) was used to determine tissue concentration of oxyhemoglobin, deoxyhemoglobin, total hemoglobin, tissue hemoglobin oxygen saturation, and bulk water content. Results show significant and dramatic differences in optical properties between menopausal states. Average premenopausal intrinsic tissue absorption and reduced scattering values at each wavelength are 2.5- to 3-fold higher and 16–28% greater, respectively, than absorption and scattering for postmenopausal subjects. Absorption and scattering properties for women using hormone replacement therapy are intermediate between premenopausal and postmenopausal populations. Physiological properties show differences in mean total hemoglobin (7.0 μM, 11.8 μM, and 19.2 μM) and water concentration relative to pure water (10.9%, 15.3%, and 27.3%) for postmenopausal, hormone replacement therapy, and premenopausal subjects, respectively. Because of their unique, quantitative information content, diffuse optical methods may play an important role in breast diagnostics and improving our understanding of breast disease.

Modulated imaging: quantitative analysis and tomography of turbid media in the spatial-frequency domain

Experiments performed on turbid phantoms demonstrate that spatially modulated illumination facilitates quantitative wide-field optical property mapping and tomographic imaging in turbid media.

Assessing the future of diffuse optical imaging technologies for breast cancer management

Diffuse optical imaging (DOI) is a noninvasive optical technique that employs near-infrared (NIR) light to quantitatively characterize the optical properties of thick tissues. Although NIR methods were first applied to breast transillumination (also called diaphanography) nearly 80 years ago, quantitative DOI methods employing time- or frequency-domain photon migration technologies have only recently been used for breast imaging (i.e., since the mid-1990s). In this review, the state of the art in DOI for breast cancer is outlined and a multi-institutional Network for Translational Research in Optical Imaging (NTROI) is described, which has been formed by the National Cancer Institute to advance diffuse optical spectroscopy and imaging (DOSI) for the purpose of improving breast cancer detection and clinical management. DOSI employs broadband technology both in near-infrared spectral and temporal signal domains in order to separate absorption from scattering and quantify uptake of multiple molecular probes based on absorption or fluorescence contrast. Additional dimensionality in the data is provided by integrating and co-registering the functional information of DOSI with x-ray mammography and magnetic resonance imaging (MRI), which provide structural information or vascular flow information, respectively. Factors affecting DOSI performance, such as intrinsic and extrinsic contrast mechanisms, quantitation of biochemical components, image formation/visualization, and multimodality co-registration are under investigation in the ongoing research NTROI sites. One of the goals is to develop standardized DOSI platforms that can be used as stand-alone devices or in conjunction with MRI, mammography, or ultrasound. This broad-based, multidisciplinary effort is expected to provide new insight regarding the origins of breast disease and practical approaches for addressing several key challenges in breast cancer, including: Detecting disease in mammographically dense tissue, distinguishing between malignant and benign lesions, and understanding the impact of neoadjuvant chemotherapies.

Predicting response to breast cancer neoadjuvant chemotherapy using diffuse optical spectroscopy

Diffuse optical spectroscopy (DOS) and imaging are emerging diagnostic techniques that quantitatively measure the concentration of deoxy-hemoglobin (ctHHb), oxy-hemoglobin (ctO2Hb), water (ctH2O), and lipid in cm-thick tissues. In early-stage clinical studies, diffuse optical imaging and DOS have been used to characterize breast tumor biochemical composition and monitor therapeutic response in stage II/III neoadjuvant chemotherapy patients. We investigated whether DOS measurements obtained before and 1 week into a 3-month adriamycin/cytoxan neoadjuvant chemotherapy regimen can predict final, postsurgical pathological response. Baseline DOS measurements of 11 patients before therapy revealed significant increases in tumor ctHHb, ctO2Hb, ctH2O, and spectral scattering slope, and decreases in bulk lipids, relative to normal breast tissue. Tumor concentrations of ctHHb, ctO2Hb, and ctH2O dropped 27 ± 15%, 33 ± 7%, and 11 ± 15%, respectively, within 1 week (6.5 ± 1.4 days) of the first treatment for pathology-confirmed responders (n = 6), whereas nonresponders (n = 5) and normal side controls showed no significant changes in these parameters. The best single predictor of therapeutic response 1 week posttreatment was ctHHb (83% sensitivity, 100% specificity), while discrimination analysis based on combined ctHHb and ctH2O changes classified responders vs. nonresponders with 100% sensitivity and specificity. In addition, the pretreatment tumor-to-normal ctO2Hb ratio was significantly higher in responders (2.82 ± 0.44) vs. nonresponders (1.82 ± 0.49). These results highlight DOS sensitivity to tumor cellular metabolism and biochemical composition and demonstrate its potential for predicting and monitoring an individuals response to treatment.

Frequency-domain photon migration measurements of normal and malignant tissue optical properties in a human subject

A 1-GHz multifrequency, multiwavelength frequency-domain photon migration instrument is used to measure quantitatively the optical absorption (mu(a)) and effective optical scattering (mu(s) ?) of normal and malignant tissues in a human subject. Large ellipsoidal (~10-cm major axis, ~6-cm minor axes) subcutaneous malignant lesions were compared with adjacent normal sites in the abdomen and back. Absorption coefficients recorded at 674, 811, 849, and 956 nm were used to calculate tissue hemoglobin concentration (oxyhemoglobin, deoxyhemoglobin, and total), water concentration, hemoglobin oxygen saturation, and blood volume fraction in vivo. Our results show that the normal and the malignant tissues measured in the patient have clearly resolvable optical and physiological property differences that may be broadly useful in identifying and characterizing tumors.

Phase measurement of light absorption and scatter in human tissue

Analog and digital technologies are presented for precise measurement of propagation delay of photons from source and detector placed on portions of the human body. The goal of the apparatus design is to quantify absorption (μa) and scattering (μs′) induced by biological pigments and biological structures, respectively. Body tissues are highly scattering with a mean distance between scatterers of less than a mm (at 700–850 nm). Significant absorption is mainly due to 5%–10% of the tissue volume occupied by blood. Measurement of μa and μs′ is done by both time and frequency domain equipment. This article focuses upon frequency domain equipment because of its simplicity, reduced noise bandwidth, versatility, and the strong analogy to very high frequency/ultrahigh frequency communication devices, particularly those using phase modulation. Comparisons are made of homodyne and heterodyne systems together with evaluation of single and multiple side band systems, with particular emphasis on methods for multiplex...