Special Section on Advanced Biomedical Imaging and Sensing

Molecular component distribution imaging of living cells by multivariate curve resolution analysis of space-resolved Raman spectra

[+] Author Affiliations
Masahiro Ando

Waseda University, Consolidated Research Institute for Advanced Science and Medical Care, 513 Wasedatsurumaki-cho, Shinjuku-ku, Tokyo 162-0041, Japan

Hiro-o Hamaguchi

Waseda University, Consolidated Research Institute for Advanced Science and Medical Care, 513 Wasedatsurumaki-cho, Shinjuku-ku, Tokyo 162-0041, Japan

National Chiao Tung University, Institute of Molecular Science, Department of Applied Chemistry, 1001 Ta Hsueh Road, Hsinchu 300, Taiwan

J. Biomed. Opt. 19(1), 011016 (Oct 09, 2013). doi:10.1117/1.JBO.19.1.011016
History: Received May 3, 2013; Revised August 22, 2013; Accepted August 26, 2013
Text Size: A A A

Open Access Open Access

Abstract.  Label-free Raman microspectroscopy combined with a multivariate curve resolution (MCR) analysis can be a powerful tool for studying a wide range of biomedical molecular systems. The MCR with the alternating least squares (MCR-ALS) technique, which retrieves the pure component spectra from complicatedly overlapped spectra, has been successfully applied to in vivo and molecular-level analysis of living cells. The principles of the MCR-ALS analysis are reviewed with a model system of titanium oxide crystal polymorphs, followed by two examples of in vivo Raman imaging studies of living yeast cells, fission yeast, and budding yeast. Due to the non-negative matrix factorization algorithm used in the MCR-ALS analysis, the spectral information derived from this technique is just ready for physical and/or chemical interpretations. The corresponding concentration profiles provide the molecular component distribution images (MCDIs) that are vitally important for elucidating life at the molecular level, as stated by Schroedinger in his famous book, “What is life?” Without any a priori knowledge about spectral profiles, time- and space-resolved Raman measurements of a dividing fission yeast cell with the MCR-ALS elucidate the dynamic changes of major cellular components (lipids, proteins, and polysaccharides) during the cell cycle. The MCR-ALS technique also resolves broadly overlapped OH stretch Raman bands of water, clearly indicating the existence of organelle-specific water structures in a living budding yeast cell.

Figures in this Article

Raman microspectroscopy, Raman spectroscopy under a microscope, is now widely used in molecular-level investigations in various fields of bioscience and biotechnology.15 It is now well established as a strategic analytical tool in these fields. From its nature, Raman spectroscopy does not require any sample pretreatment such as dye labeling or genetic manipulation and is low invasive. It is therefore highly suitable for in vivo analysis of living cells. For example, it possesses a high potential in low-invasive screening of living cells in regenerative medicine, where the safety in the re-use of screened cells is the most crucial issue. With sub-μm spatial resolution, Raman microspectroscopy provides space-resolved information on a molecular structure and their distribution inside the cell. In the past two decades, a number of publications have shown the successful application of Raman microspectroscopy to label-free molecular-level analysis of living cells and to discrimination of cell types.623

Although Raman spectra contain rich information on molecular structure, detailed interpretation of measured spectra is often difficult because of their complexity. Each Raman spectrum obtained from space-resolved mapping measurements is usually interpreted as a superposition of several spectral components of biomolecules, as well as a background and fluorescence. In order to decompose the complicated spectra into tractable component spectra, a number of chemometric methods have been developed and applied to the analysis of Raman spectra and images. Cluster analysis (CA) is one of the widely used decomposition methods.9,2428 With the CA method, by statistically analyzing the spectral variations, distinct subsets of similar spectra are obtainable. It has been applied to diagnostic tissue discrimination24 and subcellular structure imaging.9 In the CA method, the number of clusters is important.27,28 A small cluster number may result in false allocation of some raw spectra. A large number of clusters are often needed to achieve a relevant segmentation, which complicates the interpretation. Principal component analysis (PCA) is also widely used.25,2731 This method gives an orthogonal set of dominant spectral components, called principal components (PCs), as a result of matrix factorization. Each spectrum of original data can be expressed as a linear combination of PCs. PCA method has also been proved useful for the construction of molecular images from decomposed spectral component1,31,32 as well as for the discrimination of tissue or cell types.2,13,33,34 However, physical bases of this method are rather obscure; the PCs exhibit positive and negative values, thus the physical meanings of decomposed spectral components are not clear just as they are.

Recently, the multivariate curve resolution-alternating least squares (MCR-ALS) method, also known as self-modeling curve resolution or non-negative matrix factorization,35 has been developed and applied to the spectral decomposition of Raman spectra.25,29,31,3641 It has also been exploited in many kinds of spectroscopy, high performance liquid chromatography,42,43 gas chromatograph/mass spectrometry,44 UV-VIS,45 near-infrared,46,47 FT-IR,48,49 fluorescence imaging,50 etc. In the MCR method, the experimental data is approximated by a linear combination of several spectral components.51,52 The decomposition of superposed spectral data sets is done with ALS calculation, under appropriate model constraints, such as non-negativity of spectral profiles and their concentrations. Due to these constraints, this method easily provides physically interpretable spectral components, without any a priori information on chemical components in the sample specimen such as a living cell. Using this advantage, we have successfully applied the MCR-ALS method to the analysis of molecular component distribution imaging (MCDI) in living cells, whose raw spectra contain a number of unknown spectral components and are hard to interpret without a priori information. We have developed an in-house software for the MCR-ALS of numerous Raman spectral data. In the following, we first show the capability of the MCR-ALS method by using a model system consisting of two different types of crystal polymorphs of titanium oxide (TiO2), anatase and rutile. MCDI of TiO2 polymorphs is successfully obtained. Then, we show the MCR-ALS analysis of the time- and space-resolved Raman spectra of a dividing fission yeast cell. Unexpected dynamic changes of major cellular molecular components (lipids, proteins, and polysaccharides) during the cell cycle have been elucidated. Finally, we show the results of the study of intracellular water in a living budding yeast cell. We have successfully resolved so far unknown organelle-specific water structures by the MCR-ALS method. Highly important and otherwise unobtainable MCDI information has thus been obtained for the two living cell systems in vivo.

In the MCR, the experimental data is approximated by a linear combination of several spectral components. In the matrix form, this approximation can be written as Display Formula

AWH,(1)
where A is an experimental data matrix of spectra acquired at different measuring points, written as Display Formula
A=(a11a12a1na21a22a2nam1am1amn),
where m is the number of data points per spectrum along with the wavenumber axis and n is the number of spectra in the whole data set. A is decomposed to m×k matrix W, whose columns represent pure component spectra, and k×n matrix H, whose rows represent the intensity profiles of individual corresponding spectral components. k is the number of underlying constituents, which should be set a priori by the user or estimated by singular value decomposition (SVD) or PCA. In an MCR analysis, W and H are usually obtained by iterative refinement with MCR-ALS so that the Frobenius norm ||AWH||2 is minimized, under the non-negativity constraints W0 and H0. These constraints come from the fact that, in physical terms, Raman spectra and their concentration profiles must be non-negative. Unlike other factorization methods such as SVD and PCA, the MCR-ALS does not require the orthogonality of each component but only requires their non-negativity. This results in the advantage of MCR-ALS providing solutions that are more straightforward to interpret.

In practice, the initialization method and the additional constraints should be determined in advance. Several initialization methods are proposed.25,5357 For example, the initial guess of the matrix W or H can be determined by random non-negative values, by an SVD-based manner, or SIMPLISMA. The number of components, k value, is determined by a priori information of the sample species or is estimated by the number of dominant singular values from an SVD analysis. The initialization method should be appropriately chosen by the variance of spectral data set and/or signal-to-noise ratio. In order to attain sufficient decomposition ability, in addition to non-negativity constraints, further constraints for ALS optimization can be helpful.58,59 When ill-conditioned or singular W are given, where W contains similar spectral components, then H can be easily affected by the noise of raw data A, even though the optimization of Eq. (1) is achieved. Additional constraint terms can be applied to the ALS optimization; (||AWH||2+||ΓH||2) is minimized instead of ||AWH||2. In practice, L2-norm (ridge regression) and/or L1-norm (lasso regression) penalty term can be used for this purpose. With L2-norm penalty term β, the following equations are solved: Display Formula

(WTW+β2I)HWTA,(2)
Display Formula
(HHT+β2I)WHA1,(3)
where I is a k×k identity matrix. L1-norm penalty term α can be applied as follows: Display Formula
(WTW+α2E)HWTA,(4)
Display Formula
(HHT+α2E)WHAT,(5)
where E is a k×k matrix all of whose elements are unity. These equations are iteratively solved to obtain the optimized matrices H and W, respectively. The L2-norm regularization can provide preferable solutions even though the WTW or HHT matrix is singular, whereas the L1-norm regularization can provide sparser solutions. These L2- and L1-norm regularizations are efficacious in obtaining pure spectral decomposition and sparse MCDI, especially from complex sample species and low signal-to-noise ratio Raman spectral data sets.

Example of TiO2 Crystal Polymorphs Discrimination

Here, as a model system, type discrimination of TiO2 crystal polymorphs, anatase and rutile, is performed using the MCR-ALS method. Anatase and rutile TiO2 powder were mixed and placed on quartz cover slips. Space-resolved Raman spectra were obtained by using a 785-nm excitation Raman microspectroscopic system. Raman spectra of the powder mixture were collected over a 20×20μm region with a 0.25-μm scanning interval, hence, 6400 spectra were obtained. Each Raman spectrum consisted of 1340 wavenumber points corresponding to the 1340 elements of a charge-coupled device detector. Figure 1 shows the representative Raman spectra. These spectra have a number of overlapped Raman bands, interpreted as the superposition of the two intrinsic Raman spectra of anatase and rutile TiO2 shown in Fig. 2.

Graphic Jump LocationF1 :

Space-resolved Raman spectra of TiO2 powder containing anatase and rutile forms.

Graphic Jump LocationF2 :

Raman spectra of pure anatase (solid line) and rutile (dotted line) TiO2.

In the present MCR-ALS analysis, the Raman spectra sets were combined to make a 1340×6400 matrix. The number of pure spectral components was set to k=2, and the initial guess of 1340×2 matrix W was set by random numbers. In the ALS optimization, the constraints were set as follows: (1) W0 and H0, (2) L2-norm penalty term in Eq. (2) was set to be β=0.002, (3) L1-norm penalty term in Eq. (6) was set to be α=0.002. These equations were iteratively solved with non-negative matrix factorization algorithm, and in every iteration step, column vectors of the matrix W were all normalized.60 After the 500 iteration cycle, ensuring that ||AWH||2 converged to a sufficiently small value constant value, the pure spectral components were obtained as shown in Fig. 3.

Graphic Jump LocationF3 :

(a) Multivariate curve resolution-alternating least squares (MCR-ALS) retrieved Raman spectra, (b) optical microscope image of TiO2 powder, and (c) pseudocolor distribution images of components 1 (red) and 2 (green).

The MCR-ALS based factorization successfully decomposes the raw spectral data sets into two pure component spectra, i.e., the components 1 and 2 spectra are identical to the Raman spectra of anatase and rutile, respectively [Figs. 2 and 3(a)]. Based on this discrimination, each MCDI is constructed by rearranging the row vectors of H. As shown in Fig. 3(c), a clear distribution image has been obtained, providing the qualitative and quantitative information about the TiO2 mixture sample. As shown here, the MCR-ALS technique automatically resolves, without a priori spectral information, the observed set of space-resolved Raman spectra into physically interpretable spectra of the two polymorphs.

Analysis of Living Cells

Taking advantage of the physically interpretable factorization, the MCR-ALS technique can be effectively applied to molecular-level analysis of living cells. Living cells are highly complicated molecular systems and contain a large number of spectral components with no a priori information. This is particularly the case with in vivo analysis. Furthermore, they contain many compounds that have similar molecular structures. Consequently, Raman spectra tend to show many overlapped bands. This situation makes it difficult to analyze living cell Raman spectra at a detailed molecular level. The MCR-ALS analysis has great advantages to overcome these difficulties: (1) It does not need a priori information on the spectral and concentration profiles. (2) Decomposed spectra are ready for physical and/or chemical interpretation. (3) Sparseness constraints can be effectively used to achieve high contrast MCDI.

The MCR-ALS analysis is capable of extracting dynamic information from living cells. Figure 4 shows the time-lapse MCDI of a single dividing Schizosaccharomyces pombe, fission yeast cell.61 Raman mapping measurements were performed at 600 to 800 points (depending on the image size) at an interval of 0.5 μm and at nine different times (1, 2, 4, 6, 6.5, 10, 14, 18, and 22 h after inoculation of yeast cells into medium) in the cell cycle. The resultant 6885 Raman spectra were assembled to construct one A matrix; two spatial and one temporal dimensions were combined to a single dimension. The ALS optimization was conducted with an SVD-based initialization (six SVD spectral components were used as the initial guess of W matrix) and L1-norm regularization, yielding sparse solutions. The resulting six components are denoted 1 to 6 as given in Figs. 4(a) and 4(b). It should be noted that the MCR-ALS optimization started with random initialization was in vain for decomposition in this complicated cell system. In order to avoid falling into a false local minimum rather than the global minimum, the initialization of W or H was a key step.

Graphic Jump LocationF4 :

Time-lapse Raman imaging of single dividing fission yeast cell. (a) Six spectral components (1 to 6) derived from MCR-ALS and (b) distribution images of components 1 to 6 together with the optical microscope images. Scale bar=2μm.

By the MCR-ALS method, the separation of background signals is easily carried out; the component 1 is interpreted as due to the background because it shows a featureless spectral profile [Fig. 4(a) 1] and high intensity distribution patterns outside the cells [Fig. 4(b) 1]. The distributions of polysaccharides are clearly observed from the component 2, whose spectrum [Fig. 4(a) 2] very much resembles the spectra of glucan and mannan.7 It localizes at the cell wall and septum at all measurement times [Fig. 4(b) 2]. The component 3 [Fig. 4(a) 3] shows typical lipid spectrum with unsaturated lipids7 and ergosterol.20 From the time dependent changes of the corresponding spatial distribution [Fig. 4(a) 3], we know that the concentration of lipids shows a drastic decrease just after the cell division and gradually comes back as time goes on. The spectral components 4 and 5 [Fig. 4(a) 4 and 5] contain well-known protein bands. The concentrations and distributions of proteins [Fig. 4(b) 4 and 5] show time dependence that is totally different from that of lipids. The fact that we have two spectral components for proteins means that we have two groups of protein molecules with different structures (and hence different spectra) that show distinct time- and space-dependence during the cell cycle. We need additional information on those protein groupings in order to resolve this set of Raman spectra into more physically meaningful protein spectra. The origin of the component 6 is still unclear, although the MCR-ALS analysis ends up with much less clear results without this component. In this way, the intrinsic spectra and MCDI obtained from MCR-ALS has elucidated unknown and unexpected molecular-level dynamics taking place during the process of cell division.

Another study with the MCR-ALS shows the existence of organelle-specific water structures in a living budding yeast cell.62 Water molecules inside a cell are believed to play critical roles in physiological processes, creating distinct structural and chemical properties as compared to bulk water.63 A detailed intracellular water structural information in living yeast cells (diploid Saccharomyces cerevisiae) has been obtained using Raman microspectroscopy, in which the OH stretch Raman band of water is sensitive to the changes in the hydrogen-bonding networks. In the following MCR-ALS analysis, 247 mapping measured Raman spectra ranging from 3100 to 3800cm1 (572 wavenumber points) were used as an input. From an SVD analysis, the number of components was determined to be five. For the initial guess of W, the bulk water spectrum was used as the only “fixed” component and the other spectral profiles were randomly set. The ALS iteration was performed with L2-norm regularization, which is known to be effective for MCR-ALS analysis of data sets including like component spectra.

The resultant five components are shown in Fig. 5. The five pseudocolor MCDIs are almost mutually exclusive and the five spectral components show varying relative intensities of the OH stretch bands. The component 1 corresponds to bulk water whose spectrum is fixed. Its spatial distribution shows high value outside the cell. The component 5 has a significantly lower intensity in the 3200cm1 OH stretch region as compared to the intensity in the 3400cm1 region, indicating a lower proportion of stable hydrogen-bonding network in component 5 than in bulk water. From the analysis of the average spectra in the fingerprint region (data not shown here), the components 2 to 5 have been indicated to originate from the cell wall, cytoplasm, nuclear, and lipid bodies, respectively. Thus, organelle-specific water structures in living yeast cells are successfully retrieved and elucidated by using the MCR-ALS method. It is well recognized that water is the most difficult molecule to study in living cells. Only the combination of Raman microspectroscopy and MCR-ALS can provide this unique information of organelle-specific water structures in a cell in vivo.

Graphic Jump LocationF5 :

MCR-ALS analysis of the OH stretch Raman band of intracellular water in a living budding yeast cell. Raman spectral components 1 to 5 and their distribution images. The top image is the optical microscope image.

As Raman microspectroscopy has been exploited in biomedical analysis, thorough interpretations of complicated spectra have often been beset with difficulties. Even a single cell has complicated subcellular structure containing a large number of components. Various multivariate methods like CA and PCA have thus been attempted to resolve the observed spectra into a pure spectral component, but with only limited success. Due to appropriate constraints, non-negativity and L1-norm regularization for fission yeast and non-negativity and L2-norm regularization for water, the MCR-ALS method can provide physically sound spectra and high-contrast MCDI of living cells as has been shown by the two examples described here.

The authors acknowledge the work of Mr. Hirokazu Minami who developed the basis software for the MCR-ALS analysis described in the present article. They also thank the coauthors of the two original papers discussed here, Dr. C.-K. Huang, professor S. Shigeto, and Dr. S. Tiwari.

Gendrin  C., Roggo  Y., Collet  C., “Pharmaceutical applications of vibrational chemical imaging and chemometrics: a review,” J. Pharm. Biomed. Anal.. 48, (3 ), 533 –553 (2008). 0731-7085 CrossRef
Kendall  C. et al., “Vibrational spectroscopy: a clinical tool for cancer diagnostics,” Analyst. 134, (6 ), 1029 –1045 (2009). 0365-4885 CrossRef
Wu  H. W. et al., “In vivo lipidomics using single-cell Raman spectroscopy,” Proc. Natl. Acad. Sci. U. S. A.. 108, (9 ), 3809 –3814 (2011). 0027-8424 CrossRef
Lopes  P. C. et al., “Discriminating adenocarcinoma from normal colonic mucosa through deconvolution of Raman spectra,” J. Biomed. Opt.. 16, (12 ), 127001  (2011). 1083-3668 CrossRef
Krafft  C. et al., “Raman and coherent anti-Stokes Raman scattering microspectroscopy for biomedical applications,” J. Biomed. Opt.. 17, (4 ), 040801  (2012). 1083-3668 CrossRef
Huang  Y. S. et al., “Raman spectroscopic signature of life in a living yeast cell,” J. Raman Spectrosc.. 35, (7 ), 525 –526 (2004). 0377-0486 CrossRef
Huang  Y. S. et al., “Molecular-level investigation of the structure, transformation, and bioactivity of single living fission yeast cells by time- and space-resolved Raman spectroscopy,” Biochemistry. 44, (30 ), 10009 –10019 (2005). 0144-0578 CrossRef
Naito  Y., Toh-e  A., Hamaguchi  H., “In vivo time-resolved Raman imaging of a spontaneous death process of a single budding yeast cell,” J. Raman Spectrosc.. 36, (9 ), 837 –839 (2005). 0377-0486 CrossRef
van Manen  H. J. et al., “Single-cell Raman and fluorescence microscopy reveal the association of lipid bodies with phagosomes in leukocytes,” Proc. Natl. Acad. Sci. U. S. A.. 102, (29 ), 10159 –10164 (2005). 0027-8424 CrossRef
Notingher  I., Hench  L. L., “Raman microspectroscopy: a noninvasive tool for studies of individual living cells in vitro,” Expert Rev. Med. Devices. 3, (2 ), 215 –234 (2006). 1743-4440 CrossRef
Onogi  C., Hamaguchi  H., “Photobleaching of the “Raman spectroscopic signature of life” and mitochondrial activity in Rho(−) budding yeast cells,” J. Phys. Chem. B. 113, (31 ), 10942 –10945 (2009). 1520-6106 CrossRef
Onogi  C., Torii  H., Hamaguchi  H., “Raman spectra of isotope-substituted mitochondria of living budding yeast cells: possible origin of the “Raman spectroscopic signature of life,” Chem. Lett.. 38, (9 ), 898 –899 (2009). 0366-7022 CrossRef
Neugebauer  U. et al., “Identification and differentiation of single cells from peripheral blood by Raman spectroscopic imaging,” J. Biophotonics. 3, (8–9 ), 579 –587 (2010). 1864-063X CrossRef
Okuno  M., Hamaguchi  H., “Multifocus confocal Raman microspectroscopy for fast multimode vibrational imaging of living cells,” Opt. Lett.. 35, (24 ), 4096 –4098 (2010). 0146-9592 CrossRef
Okuno  M. et al., “Quantitative CARS molecular fingerprinting of single living cells with the use of the maximum entropy method,” Angew. Chem. Int. Ed.. 49, (38 ), 6773 –6777 (2010). 1433-7851 CrossRef
Onogi  C., Hamaguchi  H., “In vivo resonance Raman detection of ferrous cytochrome c from mitochondria of single living yeast cells,” Chem. Lett.. 39, (3 ), 270 –271 (2010). 0366-7022 CrossRef
Ando  M. et al., “1064 nm deep near-infrared (NIR) excited Raman microspectroscopy for studying photolabile organisms,” Appl. Spectrosc.. 65, (5 ), 488 –492 (2011). 0003-7028 CrossRef
Huang  C.-K., Hamaguchi  H., Shigeto  S., “In vivo multimode Raman imaging reveals concerted molecular composition and distribution changes during yeast cell cycle,” Chem. Commun.. 47, (33 ), 9423 –9425 (2011). 1364-548X CrossRef
Noothalapati Venkata  H. N., Nomura  N., Shigeto  S., “Leucine pools in Escherichia coli biofilm discovered by Raman imaging,” J. Raman Spectrosc.. 42, (11 ), 1913 –1915 (2011). 0377-0486 CrossRef
Chiu  L. et al., “On the origin of the 1602  cm−1 Raman band of yeasts; contribution of ergosterol,” J. Biophotonics. 5, (10 ), 724 –728 (2012). 1864-063X CrossRef
Kakita  M., Kaliaperumal  V., Hamaguchi  H., “Resonance Raman quantification of the redox state of cytochromes b and c in-vivo and in-vitro,” J. Biophotonics. 5, (1 ), 20 –24 (2012). 1864-063X CrossRef
Kaliaperumal  V., Hamaguchi  H., “Casting new physicochemical light on the fundamental biological processes in single living cells by using Raman microspectroscopy,” Chem. Rec.. 12, (6 ), 567 –580 (2012). 1527-8999 CrossRef
Kakita  M., Okuno  M., Hamaguchi  H., “Quantitative analysis of the redox states of cytochromes in a living L929 (NCTC) cell by resonance Raman microspectroscopy,” J. Biophotonics. 6, (3 ), 256 –259 (2013). 1864-063X CrossRef
Nijssen  A. et al., “Discriminating basal cell carcinoma from its surrounding tissue by Raman spectroscopy,” J. Invest. Dermatol.. 119, (1 ), 64 –69 (2002). 0022-202X CrossRef
Zhang  L., Henson  M. J., Sekulic  S. S., “Multivariate data analysis for Raman imaging of a model pharmaceutical tablet,” Anal. Chim. Acta. 545, (2 ), 262 –278 (2005). 0003-2670 CrossRef
van Manen  H.-J., Lenferink  A., Otto  C., “Noninvasive imaging of protein metabolic labeling in single human cells using stable isotopes and Raman microscopy,” Anal. Chem.. 80, (24 ), 9576 –9582 (2008). 0003-2700 CrossRef
Miljković  M. et al., “Label-free imaging of human cells: algorithms for image reconstruction of Raman hyperspectral datasets,” Analyst. 135, (8 ), 2002 –2013 (2010). 0365-4885 CrossRef
Hedegaard  M. et al., “Spectral unmixing and clustering algorithms for assessment of single cells by Raman microscopic imaging,” Theor. Chem. Acc.. 130, (4–6 ), 1249 –1260 (2011). 1432-881X CrossRef
Shafer-Peltier  K. E. et al., “Model-based biological Raman spectral imaging,” J. Cell. Biochem.. 87, (S39 ), 125 –137 (2002). 0730-2312 CrossRef
Huang  W. E. et al., “Raman microscopic analysis of single microbial cells,” Anal. Chem.. 76, (15 ), 4452 –4458 (2004). 0003-2700 CrossRef
Shinzawa  H. et al., “Multivariate data analysis for Raman spectroscopic imaging,” J. Raman Spectrosc.. 40, (12 ), 1720 –1725 (2009). 0377-0486 CrossRef
Drumm  C. A., Morris  M. D., “Microscopic Raman line-imaging with principal component analysis,” Appl. Spectrosc.. 49, (9 ), 1331 –1337 (1995). 0003-7028 CrossRef
Chan  J. W. et al., “Micro-Raman spectroscopy detects individual neoplastic and normal hematopoietic cells,” Biophys. J.. 90, (2 ), 648 –656 (2006). 0006-3495 CrossRef
Zhou  Y. et al., “Human brain cancer studied by resonance Raman spectroscopy,” J. Biomed. Opt.. 17, (11 ), 116021  (2012). 1083-3668 CrossRef
Lee  D. D., Seung  H. S., “Learning the parts of objects by non-negative matrix factorization,” Nature. 401, (6755 ), 788 –791 (1999). 0028-0836 CrossRef
Andrew  J. J., Hancewicz  T. M., “Rapid analysis of Raman image data using two-way multivariate curve resolution,” Appl. Spectrosc.. 52, (6 ), 797 –807 (1998). 0003-7028 CrossRef
de Juan  A., Tauler  R., “Chemometrics applied to unravel multicomponent processes and mixtures—revisiting latest trends in multivariate resolution,” Anal. Chim. Acta. 500, (1–2 ), 195 –210 (2003). 0003-2670 CrossRef
Perera  P. et al., “Solute-induced perturbations of solvent-shell molecules observed using multivariate Raman curve resolution,” J. Am. Chem. Soc.. 130, (14 ), 4576 –4577  (2008). 0002-7863 CrossRef
Simpson  J. V., Balakrishnan  G., Jiji  R. D., “MCR-ALS analysis of two-way UV resonance Raman spectra to resolve discrete protein secondary structural motifs,” Analyst. 134, (1 ), 138 –147 (2009). 0365-4885 CrossRef
Patel  I. I. et al., “High contrast images of uterine tissue derived using Raman microspectroscopy with the empty modelling approach of multivariate curve resolution-alternating least squares,” Analyst. 136, (23 ), 4950 –495 9 (2011). 0365-4885 CrossRef
Kwok  K., Taylor  L. S., “Analysis of counterfeit Cialis tablets using Raman microscopy and multivariate curve resolution,” J. Pharm. Biomed. Anal.. 66, , 126 –135 (2012). 0731-7085 CrossRef
Tauler  R., “Multivariate curve resolution applied to second order data,” Chemom. Intell. Lab. Syst.. 30, (1 ), 133 –146 (1995). 0169-7439 CrossRef
Tauler  R., Smilde  A., Kowalski  B., “Selectivity, local rank, three-way data analysis and ambiguity in multivariate curve resolution,” J. Chemom.. 9, (1 ), 31 –58 (1995). 0886-9383 CrossRef
Jalali-Heravi  M., Parastar  H., “Recent trends in application of multivariate curve resolution approaches for improving gas chromatography-mass spectrometry analysis of essential oils,” Talanta. 85, (2 ), 835 –849 (2011). 0039-9140 CrossRef
Bijlsma  S., Smilde  A. K., “Application of curve resolution based methods to kinetic data,” Anal. Chim. Acta. 396, (2–3 ), 231 –240 (1999). 0003-2670 CrossRef
Awa  K. et al., “Self-modeling curve resolution (SMCR) analysis of near-infrared (NIR) imaging data of pharmaceutical tablets,” Anal. Chim. Acta. 619, (1 ), 81 –86 (2008). 0003-2670 CrossRef
Amigo  J. M., Ravn  C., “Direct quantification and distribution assessment of major and minor components in pharmaceutical tablets by NIR-chemical imaging,” Eur. J. Pharm. Sci.. 37, (2 ), 76 –82 (2009). 0928-0987 CrossRef
Budevska  B. O., Sum  S. T., Jones  T. J., “Application of multivariate curve resolution for analysis of FT-IR microspectroscopic images of in situ plant tissue,” Appl. Spectrosc.. 57, (2 ), 124 –131 (2003). 0003-7028 CrossRef
Garrido  M., Rius  F. X., Larrechi  M. S., “Multivariate curve resolution-alternating least squares (MCR-ALS) applied to spectroscopic data from monitoring chemical reactions processes,” Anal. Bioanal. Chem.. 390, (8 ), 2059 –2066 (2008). 1618-2642 CrossRef
Vermaas  W. F. J. et al., “In vivo hyperspectral confocal fluorescence imaging to determine pigment localization and distribution in cyanobacterial cells,” Proc. Natl. Acad. Sci. U. S. A.. 105, (10 ), 4050 –4055 (2008). 0027-8424 CrossRef
Windig  W., “Self-modeling mixture analysis of spectral data with continuous concentration profiles,” Chemom. Intell. Lab. Syst.. 16, (1 ), 1 –16 (1992). 0169-7439 CrossRef
Jiang  J.-H., Ozaki  Y., “Self-modeling curve resolution (SMCR): principles, techniques, and applications,” Appl. Spectrosc. Rev.. 37, (3 ), 321 –345 (2002). 0570-4928 CrossRef
Gampp  H. et al., “Calculation of equilibrium constants from multiwavelength spectroscopic data—III: model-free analysis of spectrophotometric and ESR titrations,” Talanta. 32, (12 ), 1133 –1139 (1985). 0039-9140 CrossRef
Windig  W., Guilment  J., “Interactive self-modeling mixture analysis,” Anal. Chem.. 63, (14 ), 1425 –1432 (1991). 0003-2700 CrossRef
Wild  S., Curry  J., Dougherty  A., “Improving non-negative matrix factorizations through structured initialization,” Pattern Recognit.. 37, (11 ), 2217 –2232 (2004). 0031-3203 CrossRef
Berry  M. W. et al., “Algorithms and applications for approximate nonnegative matrix factorization,” Comput. Stat. Data Anal.. 52, (1 ), 155 –173 (2007). 0167-9473 CrossRef
Boutsidis  C., Gallopoulos  E., “SVD based initialization: a head start for nonnegative matrix factorization,” Pattern Recognit.. 41, (4 ), 1350 –1362 (2008). 0031-3203 CrossRef
Paatero  P., Tapper  U., “Positive matrix factorization: a non-negative factor model with optimal utilization of error estimates of data values,” Environmetrics. 5, (2 ), 111 –126 (1994). 1180-4009 CrossRef
Hoyer  P. O., “Non-negative matrix factorization with sparseness constraints,” J. Mach. Learn. Res.. 5, , 1457 –1469 (2004). 1532-4435 
Paatero  P., Tapper  U., “Positive matrix factorization—a nonnegative factor model with optimal utilization of error estimates of data values,” Environmetrics. 5, (2 ), 111 –126 (1994). 1180-4009 CrossRef
Huang  C.-K. et al., “Disentangling dynamic changes of multiple cellular components during the yeast cell cycle by in vivo multivariate Raman imaging,” Anal. Chem.. 84, (13 ), 5661 –5668 (2012). 0003-2700 CrossRef
Tiwari  S., Ando  M., Hamaguchi  H., “Investigation of organelle-specific intracellular water structures with Raman microspectroscopy,” J. Raman Spectrosc.. 44, (2 ), 167 –169 (2013). 0377-0486 CrossRef
Ball  P., “Water as an active constituent in cell biology,” Chem. Rev.. 108, (1 ), 74 –108 (2008). 0009-2665 CrossRef

© The Authors. Published by SPIE under a Creative Commons Attribution 3.0 Unported License. Distribution or reproduction of this work in whole or in part requires full attribution of the original publication, including its DOI.

Citation

Masahiro Ando and Hiro-o Hamaguchi
"Molecular component distribution imaging of living cells by multivariate curve resolution analysis of space-resolved Raman spectra", J. Biomed. Opt. 19(1), 011016 (Oct 09, 2013). ; http://dx.doi.org/10.1117/1.JBO.19.1.011016


Figures

Graphic Jump LocationF1 :

Space-resolved Raman spectra of TiO2 powder containing anatase and rutile forms.

Graphic Jump LocationF2 :

Raman spectra of pure anatase (solid line) and rutile (dotted line) TiO2.

Graphic Jump LocationF3 :

(a) Multivariate curve resolution-alternating least squares (MCR-ALS) retrieved Raman spectra, (b) optical microscope image of TiO2 powder, and (c) pseudocolor distribution images of components 1 (red) and 2 (green).

Graphic Jump LocationF4 :

Time-lapse Raman imaging of single dividing fission yeast cell. (a) Six spectral components (1 to 6) derived from MCR-ALS and (b) distribution images of components 1 to 6 together with the optical microscope images. Scale bar=2μm.

Graphic Jump LocationF5 :

MCR-ALS analysis of the OH stretch Raman band of intracellular water in a living budding yeast cell. Raman spectral components 1 to 5 and their distribution images. The top image is the optical microscope image.

Tables

References

Gendrin  C., Roggo  Y., Collet  C., “Pharmaceutical applications of vibrational chemical imaging and chemometrics: a review,” J. Pharm. Biomed. Anal.. 48, (3 ), 533 –553 (2008). 0731-7085 CrossRef
Kendall  C. et al., “Vibrational spectroscopy: a clinical tool for cancer diagnostics,” Analyst. 134, (6 ), 1029 –1045 (2009). 0365-4885 CrossRef
Wu  H. W. et al., “In vivo lipidomics using single-cell Raman spectroscopy,” Proc. Natl. Acad. Sci. U. S. A.. 108, (9 ), 3809 –3814 (2011). 0027-8424 CrossRef
Lopes  P. C. et al., “Discriminating adenocarcinoma from normal colonic mucosa through deconvolution of Raman spectra,” J. Biomed. Opt.. 16, (12 ), 127001  (2011). 1083-3668 CrossRef
Krafft  C. et al., “Raman and coherent anti-Stokes Raman scattering microspectroscopy for biomedical applications,” J. Biomed. Opt.. 17, (4 ), 040801  (2012). 1083-3668 CrossRef
Huang  Y. S. et al., “Raman spectroscopic signature of life in a living yeast cell,” J. Raman Spectrosc.. 35, (7 ), 525 –526 (2004). 0377-0486 CrossRef
Huang  Y. S. et al., “Molecular-level investigation of the structure, transformation, and bioactivity of single living fission yeast cells by time- and space-resolved Raman spectroscopy,” Biochemistry. 44, (30 ), 10009 –10019 (2005). 0144-0578 CrossRef
Naito  Y., Toh-e  A., Hamaguchi  H., “In vivo time-resolved Raman imaging of a spontaneous death process of a single budding yeast cell,” J. Raman Spectrosc.. 36, (9 ), 837 –839 (2005). 0377-0486 CrossRef
van Manen  H. J. et al., “Single-cell Raman and fluorescence microscopy reveal the association of lipid bodies with phagosomes in leukocytes,” Proc. Natl. Acad. Sci. U. S. A.. 102, (29 ), 10159 –10164 (2005). 0027-8424 CrossRef
Notingher  I., Hench  L. L., “Raman microspectroscopy: a noninvasive tool for studies of individual living cells in vitro,” Expert Rev. Med. Devices. 3, (2 ), 215 –234 (2006). 1743-4440 CrossRef
Onogi  C., Hamaguchi  H., “Photobleaching of the “Raman spectroscopic signature of life” and mitochondrial activity in Rho(−) budding yeast cells,” J. Phys. Chem. B. 113, (31 ), 10942 –10945 (2009). 1520-6106 CrossRef
Onogi  C., Torii  H., Hamaguchi  H., “Raman spectra of isotope-substituted mitochondria of living budding yeast cells: possible origin of the “Raman spectroscopic signature of life,” Chem. Lett.. 38, (9 ), 898 –899 (2009). 0366-7022 CrossRef
Neugebauer  U. et al., “Identification and differentiation of single cells from peripheral blood by Raman spectroscopic imaging,” J. Biophotonics. 3, (8–9 ), 579 –587 (2010). 1864-063X CrossRef
Okuno  M., Hamaguchi  H., “Multifocus confocal Raman microspectroscopy for fast multimode vibrational imaging of living cells,” Opt. Lett.. 35, (24 ), 4096 –4098 (2010). 0146-9592 CrossRef
Okuno  M. et al., “Quantitative CARS molecular fingerprinting of single living cells with the use of the maximum entropy method,” Angew. Chem. Int. Ed.. 49, (38 ), 6773 –6777 (2010). 1433-7851 CrossRef
Onogi  C., Hamaguchi  H., “In vivo resonance Raman detection of ferrous cytochrome c from mitochondria of single living yeast cells,” Chem. Lett.. 39, (3 ), 270 –271 (2010). 0366-7022 CrossRef
Ando  M. et al., “1064 nm deep near-infrared (NIR) excited Raman microspectroscopy for studying photolabile organisms,” Appl. Spectrosc.. 65, (5 ), 488 –492 (2011). 0003-7028 CrossRef
Huang  C.-K., Hamaguchi  H., Shigeto  S., “In vivo multimode Raman imaging reveals concerted molecular composition and distribution changes during yeast cell cycle,” Chem. Commun.. 47, (33 ), 9423 –9425 (2011). 1364-548X CrossRef
Noothalapati Venkata  H. N., Nomura  N., Shigeto  S., “Leucine pools in Escherichia coli biofilm discovered by Raman imaging,” J. Raman Spectrosc.. 42, (11 ), 1913 –1915 (2011). 0377-0486 CrossRef
Chiu  L. et al., “On the origin of the 1602  cm−1 Raman band of yeasts; contribution of ergosterol,” J. Biophotonics. 5, (10 ), 724 –728 (2012). 1864-063X CrossRef
Kakita  M., Kaliaperumal  V., Hamaguchi  H., “Resonance Raman quantification of the redox state of cytochromes b and c in-vivo and in-vitro,” J. Biophotonics. 5, (1 ), 20 –24 (2012). 1864-063X CrossRef
Kaliaperumal  V., Hamaguchi  H., “Casting new physicochemical light on the fundamental biological processes in single living cells by using Raman microspectroscopy,” Chem. Rec.. 12, (6 ), 567 –580 (2012). 1527-8999 CrossRef
Kakita  M., Okuno  M., Hamaguchi  H., “Quantitative analysis of the redox states of cytochromes in a living L929 (NCTC) cell by resonance Raman microspectroscopy,” J. Biophotonics. 6, (3 ), 256 –259 (2013). 1864-063X CrossRef
Nijssen  A. et al., “Discriminating basal cell carcinoma from its surrounding tissue by Raman spectroscopy,” J. Invest. Dermatol.. 119, (1 ), 64 –69 (2002). 0022-202X CrossRef
Zhang  L., Henson  M. J., Sekulic  S. S., “Multivariate data analysis for Raman imaging of a model pharmaceutical tablet,” Anal. Chim. Acta. 545, (2 ), 262 –278 (2005). 0003-2670 CrossRef
van Manen  H.-J., Lenferink  A., Otto  C., “Noninvasive imaging of protein metabolic labeling in single human cells using stable isotopes and Raman microscopy,” Anal. Chem.. 80, (24 ), 9576 –9582 (2008). 0003-2700 CrossRef
Miljković  M. et al., “Label-free imaging of human cells: algorithms for image reconstruction of Raman hyperspectral datasets,” Analyst. 135, (8 ), 2002 –2013 (2010). 0365-4885 CrossRef
Hedegaard  M. et al., “Spectral unmixing and clustering algorithms for assessment of single cells by Raman microscopic imaging,” Theor. Chem. Acc.. 130, (4–6 ), 1249 –1260 (2011). 1432-881X CrossRef
Shafer-Peltier  K. E. et al., “Model-based biological Raman spectral imaging,” J. Cell. Biochem.. 87, (S39 ), 125 –137 (2002). 0730-2312 CrossRef
Huang  W. E. et al., “Raman microscopic analysis of single microbial cells,” Anal. Chem.. 76, (15 ), 4452 –4458 (2004). 0003-2700 CrossRef
Shinzawa  H. et al., “Multivariate data analysis for Raman spectroscopic imaging,” J. Raman Spectrosc.. 40, (12 ), 1720 –1725 (2009). 0377-0486 CrossRef
Drumm  C. A., Morris  M. D., “Microscopic Raman line-imaging with principal component analysis,” Appl. Spectrosc.. 49, (9 ), 1331 –1337 (1995). 0003-7028 CrossRef
Chan  J. W. et al., “Micro-Raman spectroscopy detects individual neoplastic and normal hematopoietic cells,” Biophys. J.. 90, (2 ), 648 –656 (2006). 0006-3495 CrossRef
Zhou  Y. et al., “Human brain cancer studied by resonance Raman spectroscopy,” J. Biomed. Opt.. 17, (11 ), 116021  (2012). 1083-3668 CrossRef
Lee  D. D., Seung  H. S., “Learning the parts of objects by non-negative matrix factorization,” Nature. 401, (6755 ), 788 –791 (1999). 0028-0836 CrossRef
Andrew  J. J., Hancewicz  T. M., “Rapid analysis of Raman image data using two-way multivariate curve resolution,” Appl. Spectrosc.. 52, (6 ), 797 –807 (1998). 0003-7028 CrossRef
de Juan  A., Tauler  R., “Chemometrics applied to unravel multicomponent processes and mixtures—revisiting latest trends in multivariate resolution,” Anal. Chim. Acta. 500, (1–2 ), 195 –210 (2003). 0003-2670 CrossRef
Perera  P. et al., “Solute-induced perturbations of solvent-shell molecules observed using multivariate Raman curve resolution,” J. Am. Chem. Soc.. 130, (14 ), 4576 –4577  (2008). 0002-7863 CrossRef
Simpson  J. V., Balakrishnan  G., Jiji  R. D., “MCR-ALS analysis of two-way UV resonance Raman spectra to resolve discrete protein secondary structural motifs,” Analyst. 134, (1 ), 138 –147 (2009). 0365-4885 CrossRef
Patel  I. I. et al., “High contrast images of uterine tissue derived using Raman microspectroscopy with the empty modelling approach of multivariate curve resolution-alternating least squares,” Analyst. 136, (23 ), 4950 –495 9 (2011). 0365-4885 CrossRef
Kwok  K., Taylor  L. S., “Analysis of counterfeit Cialis tablets using Raman microscopy and multivariate curve resolution,” J. Pharm. Biomed. Anal.. 66, , 126 –135 (2012). 0731-7085 CrossRef
Tauler  R., “Multivariate curve resolution applied to second order data,” Chemom. Intell. Lab. Syst.. 30, (1 ), 133 –146 (1995). 0169-7439 CrossRef
Tauler  R., Smilde  A., Kowalski  B., “Selectivity, local rank, three-way data analysis and ambiguity in multivariate curve resolution,” J. Chemom.. 9, (1 ), 31 –58 (1995). 0886-9383 CrossRef
Jalali-Heravi  M., Parastar  H., “Recent trends in application of multivariate curve resolution approaches for improving gas chromatography-mass spectrometry analysis of essential oils,” Talanta. 85, (2 ), 835 –849 (2011). 0039-9140 CrossRef
Bijlsma  S., Smilde  A. K., “Application of curve resolution based methods to kinetic data,” Anal. Chim. Acta. 396, (2–3 ), 231 –240 (1999). 0003-2670 CrossRef
Awa  K. et al., “Self-modeling curve resolution (SMCR) analysis of near-infrared (NIR) imaging data of pharmaceutical tablets,” Anal. Chim. Acta. 619, (1 ), 81 –86 (2008). 0003-2670 CrossRef
Amigo  J. M., Ravn  C., “Direct quantification and distribution assessment of major and minor components in pharmaceutical tablets by NIR-chemical imaging,” Eur. J. Pharm. Sci.. 37, (2 ), 76 –82 (2009). 0928-0987 CrossRef
Budevska  B. O., Sum  S. T., Jones  T. J., “Application of multivariate curve resolution for analysis of FT-IR microspectroscopic images of in situ plant tissue,” Appl. Spectrosc.. 57, (2 ), 124 –131 (2003). 0003-7028 CrossRef
Garrido  M., Rius  F. X., Larrechi  M. S., “Multivariate curve resolution-alternating least squares (MCR-ALS) applied to spectroscopic data from monitoring chemical reactions processes,” Anal. Bioanal. Chem.. 390, (8 ), 2059 –2066 (2008). 1618-2642 CrossRef
Vermaas  W. F. J. et al., “In vivo hyperspectral confocal fluorescence imaging to determine pigment localization and distribution in cyanobacterial cells,” Proc. Natl. Acad. Sci. U. S. A.. 105, (10 ), 4050 –4055 (2008). 0027-8424 CrossRef
Windig  W., “Self-modeling mixture analysis of spectral data with continuous concentration profiles,” Chemom. Intell. Lab. Syst.. 16, (1 ), 1 –16 (1992). 0169-7439 CrossRef
Jiang  J.-H., Ozaki  Y., “Self-modeling curve resolution (SMCR): principles, techniques, and applications,” Appl. Spectrosc. Rev.. 37, (3 ), 321 –345 (2002). 0570-4928 CrossRef
Gampp  H. et al., “Calculation of equilibrium constants from multiwavelength spectroscopic data—III: model-free analysis of spectrophotometric and ESR titrations,” Talanta. 32, (12 ), 1133 –1139 (1985). 0039-9140 CrossRef
Windig  W., Guilment  J., “Interactive self-modeling mixture analysis,” Anal. Chem.. 63, (14 ), 1425 –1432 (1991). 0003-2700 CrossRef
Wild  S., Curry  J., Dougherty  A., “Improving non-negative matrix factorizations through structured initialization,” Pattern Recognit.. 37, (11 ), 2217 –2232 (2004). 0031-3203 CrossRef
Berry  M. W. et al., “Algorithms and applications for approximate nonnegative matrix factorization,” Comput. Stat. Data Anal.. 52, (1 ), 155 –173 (2007). 0167-9473 CrossRef
Boutsidis  C., Gallopoulos  E., “SVD based initialization: a head start for nonnegative matrix factorization,” Pattern Recognit.. 41, (4 ), 1350 –1362 (2008). 0031-3203 CrossRef
Paatero  P., Tapper  U., “Positive matrix factorization: a non-negative factor model with optimal utilization of error estimates of data values,” Environmetrics. 5, (2 ), 111 –126 (1994). 1180-4009 CrossRef
Hoyer  P. O., “Non-negative matrix factorization with sparseness constraints,” J. Mach. Learn. Res.. 5, , 1457 –1469 (2004). 1532-4435 
Paatero  P., Tapper  U., “Positive matrix factorization—a nonnegative factor model with optimal utilization of error estimates of data values,” Environmetrics. 5, (2 ), 111 –126 (1994). 1180-4009 CrossRef
Huang  C.-K. et al., “Disentangling dynamic changes of multiple cellular components during the yeast cell cycle by in vivo multivariate Raman imaging,” Anal. Chem.. 84, (13 ), 5661 –5668 (2012). 0003-2700 CrossRef
Tiwari  S., Ando  M., Hamaguchi  H., “Investigation of organelle-specific intracellular water structures with Raman microspectroscopy,” J. Raman Spectrosc.. 44, (2 ), 167 –169 (2013). 0377-0486 CrossRef
Ball  P., “Water as an active constituent in cell biology,” Chem. Rev.. 108, (1 ), 74 –108 (2008). 0009-2665 CrossRef

Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging & repositioning the boxes below.

Related Book Chapters

Topic Collections

Advertisement
  • Don't have an account?
  • Subscribe to the SPIE Digital Library
  • Create a FREE account to sign up for Digital Library content alerts and gain access to institutional subscriptions remotely.
Access This Article
Sign in or Create a personal account to Buy this article ($20 for members, $25 for non-members).
Access This Proceeding
Sign in or Create a personal account to Buy this article ($15 for members, $18 for non-members).
Access This Chapter

Access to SPIE eBooks is limited to subscribing institutions and is not available as part of a personal subscription. Print or electronic versions of individual SPIE books may be purchased via SPIE.org.