1 March 2007 Intrinsic tumor biomarkers revealed by novel double-differential spectroscopic analysis of near-infrared spectra
Author Affiliations +
J. of Biomedical Optics, 12(2), 020509 (2007). doi:10.1117/1.2709701
We develop a double-differential spectroscopic analysis method for broadband near-infrared (NIR, 650 to 1000 nm) absorption spectra. Application of this method to spectra of tumor-containing breast tissue reveals specific cancer biomarkers. In this method, patient-specific variations in molecular composition are removed by using the normal tissue as an internal control. The effects of concentration differences of the four major tissue absorbers (oxyhemoglobin, deoxyhemoglobin, water, and bulk lipid) between the tumor and normal tissue are accounted for to reveal small spectral components unique to cancer. From a pilot study of 15 cancer patients, we find these spectral components to be characterized by specific NIR absorption bands. Based on the spectral regions of absorption at about 760, 930, and 980 nm, we identify these biomarkers with changes in state or addition of lipid and/or water. To quantify spectral variation in the absorption bands, we construct the specific tumor component (STC) index. The STC index identifies regions of the breast with tumors.
Kukreti, Cerussi, Tromberg, and Gratton: Intrinsic tumor biomarkers revealed by novel double-differential spectroscopic analysis of near-infrared spectra

The basis for near-infrared (NIR) spectroscopy and imaging in breast cancer has been that tumor alterations of tissue vascularization/angiogenesis and oxygen consumption can be measured through hemoglobin concentration and oxygenation state, respectively. Several investigators have developed diffuse optical imaging (DOI) and diffuse optical spectroscopy (DOS) instruments at discrete wavelengths in the 650- to 850-nm range to detect and characterize breast tumors due to the absorption of both oxy- and deoxyhemoglobin. NIR tissue absorption spectra are typically fitted with hemoglobin extinction spectra obtained in vitro to quantify tissue hemoglobin concentrations.

Increased levels of hemoglobin are often found in malignant breast tissues.1 However, hemoglobin is not a specific marker for tumors. Several groups have increased the wavelength range and number of absorbing components to which the absorption data are fitted to include water and bulk lipid.1, 2, 3 These four components form the basis spectra for breast tissues. However, like hemoglobin, water and bulk lipid are present in both diseased and normal tissues, and are not specific to cancer.

In our experience, fitting broadband absorption spectra to the known basis spectra for breast chromophores produces small mismatch residuals. There are several possible sources of these mismatches. One possibility is that the four component basis spectra, which have been obtained in vitro, do not correctly quantify absorption in vivo. Published hemoglobin extinction spectra have small but clear differences between their shapes.4 There is debate whether or not the spectrum of water found in tissues is the same as that of bulk water. Ideally a patient-specific basis spectra set should be used for analysis.

Another possibility for the mismatches is that the current basis spectra do not fully account for all NIR tissue absorption. For example, tumors alter protein expression, modify lipid and water states, and generate hemoglobin breakdown products.5 Multivariate spectral analysis methods have addressed this shortcoming of currently used basis spectra, but finding the meaning behind cancer-related coordinates can be difficult.

To account for all possible compositional differences between tumor and normal tissues, we hypothesize that a unique NIR tumor component spectrum exists, which we refer to as the specific tumor component (STC) spectrum. The STC spectrum would be highly specific to the tumor-containing tissue, and should not be found in normal tissues. However, since mismatches between absorption spectra and the basis spectra fits exist for both tumor and normal tissues, we need a method that can account for patient-specific shifts in basis spectra, as well as individual biochemical variations in normal tissue. To detect a potential STC spectrum, we have developed a double-differential method. We begin with the difference spectrum between tumor and normal tissue, fit this spectrum to the four component basis spectra, then focus on the residuals to this fit.

It is imperative that we obtain a broadband tissue absorption spectrum that is corrected for the effects of multiple scattering. We used a DOS instrument to measure broadband NIR absorption and scattering spectra from 650 to 1000nm .1 Using a handheld probe absorption spectra were obtained at several spatial locations over both tumor-containing and normal tissue for each subject, with the normal region serving as a patient-specific control. Fifteen subjects were studied, and all lesions were confirmed by standard pathology. Data were acquired in compliance with an institutionally approved human subjects research protocol [University of California, Irvine (UCI) 02-2306 and 95-563].

Here we present the algorithm for the double-differential spectroscopic method for computing the residual due to components that are unaccounted for by the basis spectra. Tissue absorption and scattering spectra were calculated using custom software in MATLAB .1 The resulting absorption spectra were analyzed by in-house Elantest software to calculate the STC spectra.

We begin by assuming that the absorption spectrum at a particular location can be expressed as a linear combination of the basis spectra plus the unaccounted STC:


where AT is the absorption spectrum of tumor-containing tissue as a function of position on the breast given by x and y coordinates, and wavelength λ . Following, ai is the fractional contribution to the overall absorption. Si* contains the four component basis spectra. The * in this notation indicates that the spectra are patient specific (i.e., not known). In other words, we account for individual molecular differences of hemoglobin, water, and lipid by defining basis spectra for each subject. The goal of our analysis is to recover the STC.

In Fig. 1 we show a schematic of the algorithm. We begin by using a representative AT (650 to 990nm ) at a single spatial location over the tumor [Fig. 1a]. This algorithm should also be performed for normal tissue. We then calculate the average absorption spectrum of the normal breast AN(λ) . This is obtained by averaging the spectra taken at the various positions on the unaffected normal breast tissue. By taking an average we acquire a baseline spectrum of “normal” tissue for the patient [Fig. 1b]:


where ci is the fractional contribution of the patient-specific basis spectra to the overall absorption in normal tissue.

Fig. 1

Algorithm for the double-differential method: (a) Scatter-corrected absorption spectrum at a single spatial location over tumor-containing breast tissue. (b) Representation of average normal absorption spectra. (c) Difference spectrum obtained by subtracting (b) from (a). (d) STC spectrum obtained from subtracting the fit of the difference spectrum using the four basis components from the difference spectrum (c).


The difference spectrum D [Fig. 1c], at each spatial location is obtained by subtracting the AN [Fig. 1b] from the AT [Fig. 1a]:


Here, the symbol Δi(x,y)=(aici) indicates differences between the basis component concentrations at different spatial locations. In calculating STC, the spectra are referenced to a patient-specific baseline of normal tissue. In other words, we correct the tumor spectra for any intrasubject variation to reveal only differences between tumor and normal tissue. Solving Eq. 3 for the STC spectrum leaves:


In Eq. 4, the last term is the fit of D to the standard basis spectra. By fitting the difference spectrum D as opposed to AT , the variation due to differences between the S and S* basis spectra are minimized; thus the STC component can be recovered with relatively high precision. If we had fitted Eq. 1 instead of Eq. 3 using the standard basis spectra, the coefficients of the fit ai would have been large and the differences between S* and S (the patient-specific basis spectra and the standard basis spectra) would overwhelm the subtle differences due to the STC component. Thus the single differential residual mainly reflects intersubject differences, thereby masking the subtle differences between tumor and normal tissue for a given patient. In the double-differential method, we subtract out the common residual for tumor and normal from fitting to the standard basis spectra; thus the selection of the basis spectra has little impact in revealing spectral differences due only to tumors.

For the normal tissue, the difference spectrum (at each location) has been completely fitted by using the four components, resulting in residual spectra which are essentially a flat, featureless line with values close to zero. This means that the spectral differences between the normal and reference spectrum should be accounted for by the “natural” compositional differences, thus no new component should be needed.

Two crucial internal controls must be performed for each patient. If we calculate the residual from the double-differential method for different locations of the unaffected breast, we expect the residual spectra to be that of normal tissue. If we repeat this second differential analysis for the breast with tumor, we expect to observe a spatially dependent residual. The outcome for the normal surrounding tissue should be the same as that for the normal tissue spectra of the unaffected side. However, at the tumor regions we expect a residual different from normal tissue. We performed a phantom experiment to verify that the double-differential method can recover an unknown spectrum that is not part of the basis set. Unaccounted for components were recovered using the double-differential method.

In Fig. 2 , we show a comparison of the average of 15 tumor spectra along with the average of normal spectra obtained from the contralateral normal side. Patients ranged in age from 32 to 57 years, nine pre- and six postmenopausal, with pathologically confirmed diagnosis of invasive ductal carcinoma, with one patient diagnosed with adenocarcinoma with lobular features.

Fig. 2

Average STC spectra (double-differential residual) for 15 patients at tumor locations exhibiting the largest variation in comparison to surrounding normal tissue on tumor-containing breast, as well as normal tissue at the equivalent position on normal contralateral breast. Standard error of mean bars calculated every 20nm indicate STC spectral variation in population. Regions are defined in the text.


From inspection of the average STC spectra, we find systematic differences in roughly five spectral regions: region 1: 650 to 665nm ; region 2: 730 to 800nm ; region 3: 875 to 930nm ; region 4: 930 to 960nm ; and region 5: 980 to 990nm (Fig. 1). In an effort to quantify the amount of STC spectrum, or the amount of “tumor biomarker,” we calculate the local residual variance for each spectral region, defined by:


The local variance Lk is a function of position on the breast given by x and y coordinates. The index k indicates a given spectral region and Nk indicates the total number of wavelengths in the spectral region. STCi(λi,x,y) is the value of the STC spectra at a given wavelength. The STC index is the sum of all local variance Lk . Across all patients, the STC index displays the maximum value over regions of tumor-containing tissue, as opposed to normal tissue surrounding the lesion or normal contralateral breast tissue. Note that the STC index is our first step toward quantifying the unique spectrum and using it to identify lesions. Improved methods of spectral content analysis must be further developed.

With regards to tissue heterogeneity, we make two observations. First, we used reference spectra from different spatial locations and found that tissue heterogeneity does not affect the recovered STC spectra. Second, we believe that the amount of the STC component could depend on the particular geometry of the lesion. However, the overall spectra shape of the STC should be less affected.

The investigation of the biochemical/physical origin of the STC spectrum will be the subject of a future manuscript. Changes in region 1 (650 to 665nm ) may be due to other types of hemoglobins or breakdown products in tumors. In region 2 (730 to 800nm ), we observe a distinct negative peak in the STC spectra. As per the definition of STC, a negative peak indicates absence of a component or a broadening of the NIR band. Spectral changes in this region may be suggestive of different lipid composition in the tumor. With regards to region 3 (875 to 930nm ) and region 4 (930 to 960nm ), we see a negative peak immediately followed by a positive peak. This may be indicative of a spectral shift toward longer wavelengths. This is also the spectral region in which lipids have characteristic absorption spectra due to absorption of C-H bonds at 920 to 930nm . Thus if there are differences in lipid composition in the tumor-containing tissue with respect to the normal tissue, then this spectral region is likely to be affected. In region 5 ( 980990nm ) we observe spectral changes that may be attributed to differences in the O-H overtone region. Two possible candidates for these changes are the O-H of water and lipid oxidation products.

In conclusion, we present an algorithm for a novel analysis method of NIR absorption spectra. We show an application of method to breast tissue absorption spectra. Through this method we retrieve the spectral differences between normal and diseased tissue for a patient by fitting the difference spectra using the four basis component spectra and then analyzing the residuals of this fit. This differential approach requires regions of normal and tumor breast tissue within a patient. With this method, we show the presence of spatially localized intrinsic specific spectroscopic biomarkers of breast tumors.


This work was made possible by the NIH NCRR Biomedical Technology Center, Laser Microbeam and Medical Program (LAMMP, P41-RR01192), NCI Network for Translational Research in Optical Imaging (U54-CA105480-01), California Breast Cancer Research Program, and the Chao Family Comprehensive Cancer Center (P30-CA62203). The research was also supported by NIH RO1 (EB00559), and NIH (P41-RR03155). Beckman Laser Institute programmatic support is provided by the AFOSR MFEL and the Beckman Foundation. Finally, we thank our staff and patients who generously assisted with these studies.



A. Cerussi, N. Shah, D. Hsiang, A. Durkin, J. Butler, and B. J. Tromberg, “In vivo absorption, scattering and physiologic properties of 58 malignant breast tumors determined by broadband diffuse optical spectroscopy,” J. Biomed. Opt.1083-366810.1117/1.2337546 11(4), 044005 (2006).Google Scholar


B. W. Pogue, S. Jiang, H. Dehghani, C. Kogel, S. Soho, S. Srinivasan, X. Song, T. D. Tosteson, S. P. Poplack, and K. Paulsen, “Characterization of hemoglobin, water, and NIR scattering in breast tissue: analysis of intersubject variability and menstrual cycle changes,” J. Biomed. Opt.1083-366810.1117/1.1691028 9(3), 541–52 (2004).Google Scholar


P. Taroni, G. Danesini, A. Torricelli, A. Pifferi, L. Spinelli, and R. Cubeddu, “Clinical trial of time-resolved scanning optical mammography at 4 wavelengths between 683 and 975nm,” J. Biomed. Opt.1083-366810.1117/1.1695561 9(3), 464–73 (2004).Google Scholar


J. G. Kim, M. N. Xia, and H. L. Liu, “Extinction coefficients of hemoglobin for near-infrared spectroscopy of tissue,” IEEE Eng. Med. Biol. Mag.0739-517510.1109/MEMB.2005.1411359 24(2), 118–121 (2005).Google Scholar


R. Katz-Brull, P. T. Lavin, and R. E. Lenkinski, “Clinical utility of proton magnetic resonance spectroscopy in characterizing breast lesions,” J. Natl. Cancer Inst.0027-8874 94(16), 1197–1203 (2002).Google Scholar

Shwayta Kukreti, Albert E. Cerussi, Bruce Jason Tromberg, Enrico Gratton, "Intrinsic tumor biomarkers revealed by novel double-differential spectroscopic analysis of near-infrared spectra," Journal of Biomedical Optics 12(2), 020509 (1 March 2007). http://dx.doi.org/10.1117/1.2709701





Near infrared



Back to Top