| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
ORIGINAL ARTICLES |
From the Departments of Surgery (TMB), Anatomic Pathology (FX, KWG), Biomedical Engineering (GMP, NR), and Electrical and Computer Engineering (CZ), University of Wisconsin-Madison, Madison, Wisconsin.
Correspondence: Address correspondence and reprint requests to: Tara M. Breslin, MD, Department of Surgery, University of Wisconsin, 600 Highland Ave., H4/744 CSC, Madison, WI 53792-7375; Fax: 608-263-7652; E-mail: breslin{at}surgery.wisc.edu
| ABSTRACT |
|---|
|
|
|---|
Methods: Optical measurements were performed on 56 samples of tumor or benign breast tissue. Autofluorescence spectra were measured at excitation wavelengths ranging from 300 to 460 nm, and diffuse reflectance was measured between 300 and 600 nm. Principal component analysis to dimensionally reduce the spectral data and a Wilcoxon ranked sum test were used to determine which wavelengths showed statistically significant differences. A support vector machine algorithm compared classification results with the histological diagnosis (gold standard).
Results: Several excitation wavelengths and diffuse reflectance spectra showed significant differences between tumor and benign tissues. By using the support vector machine algorithm to incorporate relevant spectral differences, a sensitivity of 70.0% and specificity of 91.7% were achieved.
Conclusions: A statistically significant difference was demonstrated in the diffuse reflectance and fluorescence emission spectra of benign and malignant breast tissue. These differences could be exploited in the development of adjuncts to diagnostic and surgical procedures.
Key Words: Detection Breast cancer Fluorescence spectroscopy Imaging
| INTRODUCTION |
|---|
|
|
|---|
Fluorescence spectroscopy has been successfully used to provide fast and minimally invasive detection of cancers and precancers in a variety of organ systems in vivo, including the cervix, colon, bronchus, bladder, and oral mucosa.38 These studies involved modification of existing endoscopic equipment to include fluorescence spectroscopy capability. Several groups have evaluated fluorescence spectroscopy for breast cancer detection in ex vivo studies. Alfano et al.9 were the first to measure spectra of normal and malignant breast tissues from two patients at 488 and 457.9 nm excitation. Subsequently, Yang et al.1012 reported 93% sensitivity and 95% specificity rates for discrimination between malignant and benign tissue by using 300 nm excitation spectra.
Diffuse reflectance spectroscopy for breast cancer detection has also been studied. Using wavelengths between 330 and 750 nm, Bigio et al.13 measured diffuse reflectance spectra through a core biopsy needle and during breast cancer surgery and showed that this technique can differentiate tumor from normal tissue with a sensitivity of 60% to 70% and a specificity of 85% to 95%.
Despite the promising results of previous work, significant gaps in knowledge remain. The main limitations of the previous studies are that fluorescence spectra were obtained only at one excitation wavelength or several excitation wavelengths, and the utility of combining fluorescence and diffuse reflectance spectroscopy has not been evaluated. In addition, the challenge exists of adapting this technology, which is ideally suited to evaluating epithelial surfaces, to a three-dimensional organ system. The primary goal of this study was to characterize the multiexcitation fluorescence spectra (at nine excitation wavelengths in the UV/VIS range) and UV/VIS diffuse reflectance spectra of benign and malignant breast tissue and to identify the optimal spectral features for breast cancer diagnosis. We hypothesized that this technique would be able to demonstrate a statistically significant difference between benign and malignant breast tissue on the basis of their unique fluorescence and reflectance properties. In a previous report,14 we described the development of the optical system and a novel nonparametric algorithm for statistical analysis and sample classification. This report focuses on the clinical application of this exciting technology to breast cancer diagnosis and surgery.
| PATIENTS AND METHODS |
|---|
|
|
|---|
1 cm away from the grossly visible tumor margin to minimize the potential for measuring adjacent ductal carcinoma-in-situ. A Skinscan spectrofluorometer (J. Y. Horiba, Edison, NJ) was used for all measurements. This instrument consists of a fiberoptic probe with a central collection region and an outer ring of excitation fibers. Fluorescence emission spectra were recorded in 20-nm increments at excitation wavelengths of 300 to 460 nm. At each excitation wavelength, fluorescence emission was recorded in 5-nm increments, beginning at a wavelength 10 nm longer than the excitation wavelength, up to 600 nm (e.g., 310 to 600 nm for a 300-nm excitation). Fluorescence emission spectra were thus obtained at nine excitation wavelengths. Diffuse reflectance was measured by illuminating and collecting at the same wavelength ranging from 300 to 600 nm in 5-nm increments.
After each measurement, the probe position on each tissue sample was inked (TMD-BK; Triangle Biomedical Sciences, Durham, NC), and the specimen was formalin-fixed and processed for routine histopathology. Microscopic evaluation of each histological section was performed (F.X. and K.W.G.) and a consensus diagnosis was reached. When a sample exhibited a heterogeneous diagnosis at the site of measurement, the worst-case diagnosis was used (e.g., for samples that contained both normal glandular tissue and malignant tissue, the diagnosis was coded as malignant). In cases in which normal adipose and fibrous/glandular tissues were present, the histology was determined by the predominant tissue type at the measurement site.
Data analysis was performed with the Matlab (Mathworks Inc., Natick, MA) software package. Principal component analysis (PCA) was used as a data-reduction technique. PCA characterizes a majority of the variance while greatly reducing the input data set into a few orthogonal variables. The principal components (PCs) are extracted such that the first PC (PC1) accounts for the largest amount of the total variance of the input data. The second PC (PC2) accounts for the second largest amount of the variance while being orthogonal to PC1, and so on. There are two advantages to this transformation: (1) the input data can be represented by a few subsets of PCs with minimal mean square error, which reduces the dimensionality of the data set; and (2) the projection onto the PC subspace maximizes the separation of data clusters.15,16 The primary drawback of this technique is that it is not trivial to determine which parts of the data input (i.e., spectra) are diagnostically useful. To balance these concerns, the PCA was performed individually on all fluorescence emission spectra, one excitation wavelength at a time. A similar process was performed on the diffuse reflectance data, thus yielding a set of PCs to characterize each spectrum.
To determine which of the PCs represented variance due to malignancy (rather than normal variability), a Wilcoxon ranked sum test (P < .0005) was used to test which PCs demonstrated significant differences between malignant and nonmalignant samples. Because PCA was performed on each excitation wavelength separately, the wavelength responsible for each PC could be determined, and the most diagnostically useful wavelengths were deduced.
A support vector machine (SVM) algorithm was used for classification.17,18 SVM is a classification algorithm based on statistical learning theory. The principal idea of an SVM is to determine an optimal separating hyperplane that maximizes the margin between two classes in a multidimensional data space. With the largest separation of the two data clusters, the SVM classifier gives a lower expected risk, which means that future error can be minimized if more data are added to the sample pool. The final step in data analysis involved comparing the sensitivity and specificity of those PCs found to be statistically significant with the gold standard histopathological diagnosis. For comparison with the classification scheme, the samples were labeled as either tumor or benign (fibrous/glandular and adipose). The rationale for grouping adipose tissue with benign (fibrous/glandular) tissue despite the observation that they are histologically distinct, with different optical properties, was based on our previous work. In developing the statistical classification scheme, we found that the algorithm is not yet able to discriminate between adipose and nonadipose tissue with 100% accuracy. As a result, the overall classification was more accurate with a one-step process as opposed to a two-step process.14
A full cross-validation was performed on the classification scheme by sequential removal of individual samples from the training datasets used in the PCA, Wilcoxon, and SVM algorithms. The accuracy of the resulting classification scheme was then tested on the deleted sample. This method allows for unbiased evaluation of the diagnostic algorithm when the sample size is too small for division into separate training and testing datasets. This procedure was performed for all 56 samples; each case had 55 training samples and 1 testing sample.
| RESULTS |
|---|
|
|
|---|
The mean age of all participants was 48.4 years and was 51.5 years for those with cancer. Twenty-seven subjects had a diagnosis of infiltrating carcinoma, and five underwent reduction mammoplasty for benign conditions. For each sample, the histopathologic diagnosis was classified as tumor, fibrous/glandular, or adipose. Table 1 summarizes the histopathologic diagnosis of the tissue samples used for spectroscopic analysis.
|
|
|
After determining which PCs demonstrated statistically significant differences between tumor and benign tissue, the spectra were analyzed by using the SVM algorithm to classify samples as tumor or benign. These results were compared with the histopathologic diagnosis (gold standard) to determine the sensitivity and specificity of the technique. Table 2 depicts the diagnostic data in 2 x 2 format and demonstrates a sensitivity of 70.0% and specificity of 91.7% for fluorescence spectroscopy alone. Six tumor samples and three benign samples were misclassified, yielding a positive predictive value of 82.3% and a negative predictive value of 84.6%.
|
| DISCUSSION |
|---|
|
|
|---|
Diffuse reflectance exploits some of the same optical features as fluorescence spectroscopy. The most notable similarities are biological chromophores that result in absorption, such as oxygenated and deoxygenated hemoglobin and ß-carotene, and tissue scattering properties. ß-Carotene is abundant in adipose tissue and may provide a unique indicator for this tissue type. The primary advantage to using diffuse reflectance is that it can be performed at a fraction of the cost of fluorescence spectroscopy. However, with a similar classification technique, diffuse reflectance spectroscopy alone was not capable of discriminating between tumor and nontumor tissue with the same accuracy as the combined technique.14 Therefore, there is a potential advantage to using fluorescence spectroscopy instead of diffuse reflectance. This advantage is likely due to additional chemical specificity arising from characterization of the many intrinsic fluorophores present in human tissue. Future studies will focus on performing measurements with only those fluorescence and reflectance spectra found to be statistically significant by PCA.
A nontrivial issue with using endogenous fluorophores as used in this study is that is difficult to determine the true biologic basis of the spectroscopic differences observed between tumor and nontumor tissue. There are, however, known biologic fluorophores that contribute to the fluorescence spectra at the optimum excitation wavelengths as determined by this study. The excitation wavelengths used, ranging from 300 to 460 nm, allow for characterization of a number of biologic fluorophores, including tryptophan, nicotinamide adenine dinucleotide, flavoproteins, and collagen, all of which are present in tissue systems.19 The fluorescence excitation-emission wavelengths identified by the PCA/SVM algorithm suggest that tryptophan, nicotinamide adenine dinucleotide, and flavoproteins are important for breast cancer diagnosis. Hemoglobin is another important absorber that may contribute to differences between tumor and nontumor tissue.
The statistical algorithm used for this study was developed for analysis of this dataset, taking into account both the small sample size and the large number of individual data points. Previously published studies using cervical biopsies had a much larger sample size, which allowed for splitting the data into training and testing subsets. In one study, Ramanujam et al.20 divided the samples into training and testing datasets and used PCA followed by logistic discrimination to classify tissue types and calculate the posterior probability of a correct diagnosis. This technique allows for a more robust evaluation of the classification algorithm. As our dataset enlarges, we plan to develop independent training and testing datasets and explore alternate classification schemes by using neural networks.
These exciting preliminary data derived from measurements performed on ex vivo breast tissues form an essential background for performing future in vivo experiments. The spectral differences demonstrated to exist between tumor and nontumor tissue could be exploited in the development of diagnostic adjuncts to core biopsy, fiberoptic ductoscopy, or evaluation of surgical margins. Future challenges include issues such as the development of a sterilizable probe appropriate for intraparenchymal or intraductal use, determining the optimum measurement depth, and elucidating the biologic mechanisms behind these observed optical differences.
| ACKNOWLEDGMENTS |
|---|
Supported by the Whitaker Foundation and the Department of Defense Breast Cancer Research Program.
| FOOTNOTES |
|---|
Received for publication March 7, 2003. Accepted for publication August 25, 2003.
| REFERENCES |
|---|
|
|
|---|
srg/publications/pdf/SVM.pdf. Southampton, UK: University of Southampton, Department of Electronics and Computer Science: 1998.
This article has been cited by other articles:
![]() |
M. P. L. Bard, A. Amelink, M. Skurichina, V. Noordhoek Hegt, R. P. W. Duin, H. J. C. M. Sterenborg, H. C. Hoogsteden, and J. G. J. V. Aerts Optical spectroscopy for the classification of malignant lesions of the bronchial tree. Chest, April 1, 2006; 129(4): 995 - 1001. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. B. Wolbarst and W. R. Hendee Evolving and Experimental Technologies in Medical Imaging Radiology, January 1, 2006; 238(1): 16 - 39. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |