Biopolymers - 2006 - Chowdary
Biopolymers - 2006 - Chowdary
Biopolymers - 2006 - Chowdary
Chowdary1,2
K. Kalyan Kumar1,2
Jacob Kurien2
Discrimination of Normal,
Stanley Mathew3 Benign, and Malignant Breast
C. Murali Krishna1
1
Tissues by Raman Spectroscopy
Center for Laser
Spectroscopy, Manipal
Academy of Higher Education,
Manipal-576 104,
Karnataka, India
2
Department of Surgical
Oncology, Shirdi Sai Baba
Cancer Hospital and Research
Center, Manipal Academy of
Higher Education, Manipal-576
104, Karnataka, India
3
Department of General
Surgery, Kasturba Medical
College, Manipal Academy of
Higher Education, Manipal-576
104, Karnataka, India
Received 3 May 2006;
revised 15 July 2006;
accepted 2 August 2006
Published online 8 August 2006 in Wiley InterScience (www.interscience.wiley.com). DOI 10.1002/bip.20586
Abstract: Breast cancers are the leading cancers among females. Diagnosis by fine needle aspira-
tion cytology (FNAC) is the gold standard. The widely practiced screening method, mammography,
suffers from high false positive results and repeated exposure to harmful ionizing radiation. As with
all other cancers survival rates are shown to heavily depend on stage of the cancers (Stage 0, 95%;
Stage IV, 75%). Hence development of more reliable screening and diagnosis methodology is of
considerable interest in breast cancer management. Raman spectra of normal, benign, and malig-
nant breast tissue show significant differences. Spectral differences between normal and diseased
breast tissues are more pronounced than between the two pathological conditions, malignant and
benign tissues. Based on spectral profiles, the presence of lipids (1078, 1267, 1301, 1440, 1654,
1746 cm–1) is indicated in normal tissue and proteins (stronger amide I, red shifted DCH2, broad
and strong amide III, 1002, 1033, 1530, 1556 cm–1) are found in benign and malignant tissues. The
major differences between benign and malignant tissue spectra are malignant tissues seem to have
an excess of lipids (1082, 1301, 1440 cm–1) and presence of excess proteins (amide I, amide III, red
shifted DCH2, 1033, 1002 cm–1) is indicated in benign spectra. The multivariate statistical tool,
principal components analysis (PCA) is employed for developing discrimination methods. A score
of factor 1 provided a reasonable classification of all three tissue types. The analysis is further fine-
tuned by employing Mahalanobis distance and spectral residuals as discriminating parameters.
This approach is tested both retrospectively and prospectively. The limit test, which provides the
556
10970282, 2006, 5, Downloaded from https://onlinelibrary.wiley.com/doi/10.1002/bip.20586 by Agricultural University of Faisalabad, Wiley Online Library on [16/09/2023]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
Breast Tissue Discrimination Using Raman Spectroscopy 557
most unambiguous discrimination, is also considered and this approach clearly discriminated all
three tissue types. These results further support the efficacy of Raman spectroscopic methods in dis-
criminating normal and diseased breast tissues. # 2006 Wiley Periodicals, Inc. Biopolymers 83:
556–569, 2006
This article was originally published online as an accepted preprint. The ‘‘Published Online’’ date
corresponds to the preprint version. You can request a copy of the preprint by emailing the
Biopolymers editorial office at [email protected]
Table I Details of Samples Used in the Study logical conditions, benign and malignant. The spectral
features of normal tissue, bands at around 1078, 1267,
Standard Samples 1301, 1440, 1654, and 1746 cm–1, were assigned to
vibrational modes of lipids (Figure 1a). The spectral
Spectrum Number Nature Histopathology
profiles of benign and malignant tissues, strong amide
I, broad amide III, and peaks at 1002, 1083, 1530, and
1–105 Normal Normal
106–206 Malignant Malignant 1556 cm–1 were attributed to vibrational modes of
207–258 Fibroadenoma Fibroadenoma proteins (Figure 1b and 1c). The major differences
between benign spectra with respect to malignant
spectra are: red shifted and broad DCH2 peak, rela-
recorded at six or more different sites on each tissue sam- tively strong amide I, and stronger and broader amide
ple. Only good spectra from clear normal, malignant, and III (Figure 1c). All of the observed spectral features
benign tissues were employed for analysis. A total of 258
may provide vital clues in understanding the differen-
sites/spectra (105 normal, 101 malignant, 52 benign) from
29 normal, 24 malignant, and 7 benign subjects were
ces in the biochemical composition of tissues.
recruited in the study. Details of the samples used in the To bring out the differences in spectral profiles
study are shown in the Table I. more clearly, difference spectra were computed by
subtracting mean normal spectrum from mean malig-
nant and benign spectra, respectively (Figure 2a and
Laser Raman Spectroscopy
b). All of the negative peaks, 1066, 1301, 1437, and
Raman spectra were recorded using the set up assembled by 1744 cm–1, were due to normal spectrum and these
us.19–21 This setup consisted of a diode laser (785 nm, SDL spectral features can be assigned to lipids. The posi-
diode laser 8530, 150 mW) for excitation. Raman signals tive peaks in Figure 2a and b were due to malignant
were detected by combination of spectrometer (Spex Triax
and benign spectra, respectively. Most of the positive
320, 600 g/mm, blazed at 900 nm) and Spectrum one liquid
nitrogen cooled CCD (Jobin Yvon-Spex, Instruments S.A.,
bands seen around 857, 940, 1002, 1033, 1169, 1244,
Inc.). A holographic filter (HLBF-785.0, Kaiser Optics, 1281, 1341, 1381, 1417, 1465, 1530, 1556, 1587,
Ann Arbor, MI, USA) was used to reject unwanted lines 1637, and 1666 cm–1 can be assigned to different
from the excitation source. Rayleigh scattering was
removed by a notch filter (HSPF-5812, Kaiser Optics, Ann
Arbor, MI, USA). Other experimental conditions were inte-
gration time – 30 s and number of accumulations – 20.
These settings were kept constant for all the measurements.
The recorded spectra were calibrated with a cubic fit to
known frequencies of Tylenol (4-acetamidophenol).
Data Analysis
Baseline correction, smoothing, calibration, and normaliza-
tion with respect to DCH2 band were carried out using algo-
rithms of Grams 32 (Galactic Industries Corporation,
Salem, NH, USA). PCA was carried out using Grams PLS
Plus (Galactic Industries Corporation). PCA was carried out
under different conditions: different number of factors,
entire range (800–1800 cm–1), several selected regions, de-
rivative spectra to explore the best discrimination. In our
analysis, the 1400–1750 cm–1 spectral region with nine fac-
tors provided the best discrimination. Further analysis was
carried out under these conditions.
Table II PCA of Raman Spectra of Normal, The classification based on scores of factor though
Malignant, and Benign Breast Tissues: a widely used approach could be subjective as dis-
Eigen values, Total % Variance criminating components are selected. Generally algo-
rithms of discriminant analysis are used to identify
Factor Number Eigen Value Total % Variance the most discriminating factors or principal compo-
nents.18 However, this may not provide a very good
1 1047.349 93.67731
2 27.7394 96.15838
classification under all circumstances.19–21 In view of
3 22.31309 98.15412 this we have considered another approach. In this sec-
4 7.245638 98.80218 ond approach, standard sets of each class, normal,
5 4.583775 99.21217 malignant, and benign, using pathologically certified
6 3.097776 99.48924 samples were developed. Spectral residuals (sum of
7 2.279846 99.69315 the squares of the difference between observed and
8 1.869763 99.86039 simulated spectra) and Mahalanobis distance (meas-
9 1.560891 100 ured in terms standard deviation from mean of the
training set) were computed for the spectra against
ues versus factors are shown in Figure 5. The first three the standard sets in both retrospective as well as pro-
factors accounted for 98% variation in the spectra (Ta- spective modality. And these parameters have been
ble II). Profiles of factors, ‘‘factor loads,’’ of the first explored to achieve better discrimination of tissue
three factors and the last two factors are shown in Fig- type. This approach has been successfully employed
ure 4. In our analysis, scores of factor 1 provided a rea- by us for discrimination of oral and cervical tis-
sonable classification among three groups (Figure 4). sues.19–21 The mathematical basis of the Mahalanobis
As can be seen from the figure, the contribution of fac- distance calculation is well known24–26 and widely
tor 1 is negative for normal spectra and largely positive used for spectral discrimination. It is very sensitive to
for benign and malignant spectra. Cluster of malignant intervariable changes in the calibration data and it is
spectra occupy between clusters corresponding to nor- measured in terms of standard deviations from the
mal and benign tissues. The mean and standard devia- mean of the training samples. In addition to providing
tions of score of factor 1 for normal, malignant, and spectral discrimination, it also gives a statistical mea-
benign spectra are 0.068 6 0.006, 0.028 6 0.023, sure of how well the unknown sample spectrum
and 0.083 6 0.007, respectively. Thus a reasonable matches or does not match. Thus it is a statistical
classification of three clusters can be seen up to mean measure of proximity of two spectra. In the present
6 1 standard deviation, which indicates 75% sensitiv- analysis, Mahalanobis distance is coupled with PCA,
ity/specificity. which is a very effective data reduction technique for
FIGURE 5 PCA of Raman spectra of normal, malignant, and benign breast tissues. Normal;
* Malignant; ~ Benign.
spectroscopic data and thus allows full spectral cover- nign standard sets. Randomly picked spectra based on
age for all samples. In our analysis, Mahalanobis dis- scores of factor 1 and pathological certifications are
tance is computed using the expression employed to develop standard calibration sets for nor-
mal (36 spectra), malignant (35 spectra), and benign
D2 ¼ ðStest Þ M1 ðStest Þ0 (21 spectra) conditions. The standard sets were veri-
Where Stest is the vector of the scores and sum fied retrospectively where spectra from the standard
of squared spectral residuals for a given test sample, sets are matched to compute Mahalanobis distance
and M is an f by f Mahalanobis matrix given by and spectral residuals. It is expected that if spectra and
0
M ¼ ðnSS1Þ, where S contains the corresponding para- standard sets are of same class, values for Mahalano-
meters for the calibration set (n standards). bis distance and spectral residuals are low and vice
The advantages of using D2 as a discriminating pa- versa. As an example, retrospective analysis against
rameter are described elsewhere.23–25 Mahalanobis normal standard is shown in Figure 6. As can be seen
distance and spectral residuals are computed by from the figure, very good classification of tissue type
matching spectra against normal, malignant, and be- is achieved. Mahalanobis distance for all normal spec-
FIGURE 6 PCA of Raman spectra of breast tissue—retrospective study against normal standard
set: (a) plot of spectral residual versus Mahalanobis distance; (b) plot of sample number versus
Mahalanobis distance. Normal; * Malignant; ~ Benign.
tra are around 0.881 6 0.44, and values for malignant þ 1.24 and 10.52 þ 1.75, respectively. Spectral resid-
and benign spectra are around 104.42 6 28.64 and ual values for malignant, benign and normal spectra
259.85 6 28.43, respectively. Spectral residual values were around 0.65 þ 0.40, 2.60 þ 0.56 and 5.23 þ
for normal spectra are around 0.09 6 0.05 and for ma- 1.24, respectively.
lignant and benign spectra are around 8.40 6 2.27 and The above approach was further extended to com-
20.76 6 2.26, respectively. The approach was eval- pute match/mismatch tables or a ‘‘Limit test,’’ which
uated prospectively, matching all normal, benign, and is considered to obtain unambiguous discrimination
malignant spectra that were not included in the stand- of tissue types. This methodology is successfully ap-
ard sets, against standard sets. Typical results obtained plied to discriminate normal, inflammatory, premalig-
in prospective analysis against malignant standard sets nant, and malignant oral tissues and normal and ma-
are shown in Figure 7. The clear clustering, barring lignant conditions of the cervix.19–21 In this method-
minor overlap among the clusters, which could not be ology, spectra are matched against standard sets with
assigned pathological status, obtained in the analysis inclusions/exclusions set for three discrimination pa-
validates standard sets. Mean Mahalanobis distance rameters: scores of factors, Mahalanobis distance,
for all malignant spectra were around 1.19 þ 0.64, and spectral residuals. If the values of the given spec-
and values for benign and normal spectra around 5.53 trum fall within set limits, spectra show YES or PASS,
FIGURE 7 PCA of Raman spectra of breast tissues— prospective study against malignant stand-
ard set: (a) plot of spectral residual versus Mahalanobis distance; (b) plot of sample number versus
Mahalanobis distance. Normal; * Malignant; ~ Benign.
Table III Multiparametric Match/Mismatch Table Compared Against Benign Standard Set
otherwise NO or FAIL. Thus any given spectrum is and malignant, show no match (NO/FAIL). Thus the
matched against all the available standard sets and spec- results obtained in the study further support the efficacy
tra to be classified as a particular class should show a of Raman spectroscopic methods in discrimination of
match (YES/PASS) against that particular standard set normal, malignant, and benign breast tissues.
and no match (NO/FAIL) against all other sets. Thus,
this methodology provides the most unambiguous dis- CONCLUSION
crimination. A typical match/mismatch table compared
against benign standard sets is given in Table III. All Raman spectra obtained from normal and pathologi-
spectra belonging to the benign condition show match cal tissues show significant differences. The spectral
(YES/PASS) and all other classes of spectra, normal profile of normal tissues is indicative of higher levels
Biopolymers DOI 10.1002/bip
10970282, 2006, 5, Downloaded from https://onlinelibrary.wiley.com/doi/10.1002/bip.20586 by Agricultural University of Faisalabad, Wiley Online Library on [16/09/2023]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
Breast Tissue Discrimination Using Raman Spectroscopy 569
of lipids. In comparison, the spectral profile of patho- 7. Majumder, S. K.; Gupta, P. K.; Jain, B.; Uppal, A.
logical tissues, both benign and malignant, indicates Lasers Life Sci 1998, 8, 249–264.
the presence of more proteins and fewer lipids. Further, 8. Yang, Y.; Katz, A.; Celmer, E. J.; Zurawska-Szczepa-
among the pathological tissues, malignant tissue con- niak, M.; Alfano, R. R.; Photochem Photobiol 1997,
66, 518–522.
tains relatively more lipids in comparison to benign tis-
9. Bigio, I. J.; Bown, S. G.; Briggs, G.; Kelley, C.;
sue. These findings are in corroboration with our basic
Lakhani, S.; Pickard, D.; Ripley, P. M.; Rose, I. G.;
understanding that neoplastic processes are fundamen- Saunders, C. J. Biomed Opt 2000, 5, 221–228.
tally more cellular and hence would biochemically be 10. Haka, A. S.; Volynskaya, Z.; Gardecki, J. A.; Nazemi,
observed to have higher levels of protein. Also, benign J.; Hicks, D.; Fitzumaurice, M.; Dasari, R. R.; Crowe,
tissue having more stromal components relative to the J. P.; Feld, M. S. Cancer Res 2006, 66, 3317–3322.
more cellular components of malignant tissue would 11. Abigail, S. H.; Shaefer-Peltier, K. E.; Maryann, F.; Jo-
explain our observation of relatively more lipid in ma- seph, C.; Dasari, R. R.; Feld, M. S. Proc Natl Acad Sci
lignant tissue in comparison to benign tissue. PCA is U S A 2005, 102, 12371–12376.
employed for developing discrimination methods. 12. Abigail, S. H.; Shaefer-Peltier, K. E.; Maryann, F.;
Scores of factor 1 provided a reasonable discrimination. Crowe, J.; Dasari, R. R.; Feld, M. S. Cancer Res 2002,
The analysis is further finetuned by employing Mahala- 62, 5375–5380.
nobis distance and spectral residuals as discriminating 13. Shaefer-Peltier, K. E.; Haka, A. S.; Fitzumaurice, M.;
Crowe, J.; Myles, J.; Dasari, R. R.; Feld, M. S. J Raman
parameters. Limit tests or computation of match–mis-
Spectrosc 2002, 33, 552–563.
match tables showed most unambiguous discrimina-
14. Shaefer-Peltier, K. E.; Haka, A. S.; Motz, J. T.; Fitzu-
tion. Thus the results obtained in the study further sup- maurice, M.; Dasari, R. R.; Feld, M. S. J Cell Biochem
ports the efficacy of Raman spectroscopic methods in 2002, 87, S125–S137.
discrimination of normal, malignant, and benign breast 15. Kneipp, J.; Tom, B. S.; Kliffen, M.; Marian, M. P.;
tissues. Perceptively, these results are significant as Puppels, G. Vib Spectrosc 2003, 32, 67–74.
with improvement of technique (development of fiber- 16. Manoharan, R.; Shafer, K.; Perelman, I.; Wu, J.; Chen,
optic probes), reproducibility, and data analysis, develop- K.; Deinum, G.; Fitzumaurice, M.; Myles, J.; Crowe,
ment of objective, rapid, and minimally invasive Raman J.; Dasari, R. R.; Feld, M. S. Photochem Photobiol
spectroscopic technique for breast cancer screening and 1998, 67, 15–22.
diagnosis seems quite feasible. 17. Hanlon, E. B.; Manoharan, R.; Koo, T. W.; Shafer, K.
E.; Motz, J. T.; Fitzumaurice, M.; Kramer, J. R.; Itzkan,
This work is carried out under the project entitled: ‘‘Devel- I.; Dasari, R. R.; Feld, M. S. Phys Med Biol 2000, 45,
opment of laser spectroscopy techniques for early detection R1–R59.
and follow up of therapy in breast malignancy’’ No. 5/13/ 18. Murali Krishna, C.; Sockalingam, G. D.; Kurien, J.;
23/2003-NCD-III, Indian Council for Medical Research, Lakshmi, R.; Venteo, L.; Pluot, M.; Manfait, M.; Kar-
Government of India. The authors (M.V.P.C. and K.K.K.) tha, V. B. Appl Spectrosc 2004, 58, 107–114.
are thankful to ICMR for senior research fellowships. Ms. 19. Murali Krishna, C.; Prathima, N. B.; Malini, R.; Vad-
Keerthi and Mr. Chethan N. Anand are acknowledged for hiraja, B. M.; Rani, Bhatt, A.; Donald, J. F.; Pralhad,
their technical assistance in the sample handling/preserva- K.; Vidyasagar, M. S.; Kartha, V. B. Vib Spectrosc
tion and data acquisition. 2006, 41, 136–141.
20. Malini, R.; Venkatakrishna, K.; Kurien, J.; Keethilatha
Pai, M.; Kartha, V. B.; Murali Krishna, C. Biopolymers
REFERENCES 2006, 81, 179–193.
21. Manjunath, B. K.; Kurien, J.; Rao, L.; Murali Krishna,
C.; Chidananda, M. S.; Venkatakrishna, K.; Kartha, V.
1. Max, P. D.; Freddie, B.; Ferlay, J.; Poala, P. A Cancer J
B. J Photochem Photobiol B: Biology, 2004, 73, 49–58.
Clin 2005, 55, 74–108.
2. Ciatto, S.; Rosselli, D. T. M.; Zappa, M. Br J Cancer 22. Clark, R. J. H.; Hester, R. E. In Advances in Spectros-
1995, 71, 337–339. copy; New York: Wiley, 1993.
3. Esserman, L.; Cowley, H.; Eberie, C.; Eberle, C.; Kirkpa- 23. Parker, E. S. Applications of Infrared and Raman and
trick, A.; Chang, S.; Berbaum, K.; Gale, A. J Natl Cancer Resonance Raman Spectroscopy in Biochemistry;
Inst 2002, 369–375. New York: Plenum Press, 1983.
4. Destounis, S. V.; DiNitto, P.; Logan-Young, W.; 24. Mahalanobis, P. C. Proc Natl Inst Sci India, 1936, 2, 49.
Bonaccio, E.; Zuley, M. L.; Willison, K. M. Radiology 25. Mark, H. L. Anal Chem 1987, 59, 790.
2004, 232, 578–584. 26. PLSplus/IQ User’s Guide; Galactic Industries Corpora-
5. Liberman, L.; Dershaw, D. D.; Rosen, P. P.; Abramson, A. F.; tion: Salem, New Hampshire, USA, 1999.
Deutch, B. M.; Hann, L. E. Radiology 1994, 192, 793–795.
6. Gupta, P. K.; Majumder, S. K.; Uppal, A. Lasers Surg
Med 1997, 21, 417–422. Reviewing Editor: Laurence Nafie