|
1.IntroductionApplications of spectroscopy techniques to biological and biomedical problems have been rapidly advancing in recent years.1–3 Based on inelastic scattering, Raman spectroscopy, as a type of vibrational spectroscopy, provides extremely rich molecular information about multiple chemical species present in a sample/specimen simultaneously.4–13 The intensity of a Raman signal bears a linear relationship to the analyte concentrations, therefore, Raman spectroscopy can be used as a quantitative tool in concentration measurements as well.1,4,5,14–17 An ultimate goal in this field is to develop Raman spectroscopy-based techniques for biomedical applications through instrumentation,18–22 plasmonic substrates,23–27 devices,28,29 assays,30,31 and techniques.32,33 Raman spectroscopic measurements, like other optical techniques, pose minimal danger from exposure to ionizing radiation due to the low-energy optical radiation exposure. One additional advantage of the NIR source used in Raman is that the tissue sampling region is much deeper than those provided by other optical approaches due to the reduced tissue scattering and reduced water and chromophore absorption at those wavelengths. As a result, NIR Raman spectroscopy satisfies two critical prerequisites for a truly noninvasive technique to transcutaneously monitor clinically important chemicals in vivo.14,21,34–36 Nevertheless, noninvasive techniques of this kind will be valuable in a wide variety of clinical settings and laboratory tests. We have employed multivariate calibration techniques such as partial least squares (PLS) and constrained regularization to noninvasive glucose sensing using Raman spectroscopy.1,14,16,17 Glucose is a convenient analyte to study because its concentration can be conveniently altered and monitored in humans and other living animals such as dogs. The dog model provides several advantages such as a similar physiological glucose response as humans, no motion artifacts owing to the anesthesia that can be administered, and the flexibility to perform glucose clamping studies. The dog study described in this paper was performed on a beagle anesthetized for , during which its blood glucose concentration was clamped at several different levels. Raman spectra were continuously acquired from the ear and reference blood glucose measurements were taken using plasma glucose obtained from venous blood draws. Results from PLS analyses demonstrate that the calibration model can predict samples that were not included in the calibration set. The PLS analysis included a leave-one-level-out analysis, and thus represents a step in the direction of a completely prospective analysis. 2.Methods2.1.Raman SystemThe Raman system for this study was modified from the previously published version for human volunteers.14 The geometry of the light delivery path was modified to allow the excitation laser to have normal incidence from beneath the dog ear through a hole in the paraboloidal mirror that subsequently collects and collimates the backscattered Raman signal. In other words, a mirror was placed underneath the paraboloidal mirror to redirect the originally in-plane optical path to be perpendicular to the optical plane. Figure 1 shows the equivalent optical path after folding the perpendicular path in-plane. The dog ear was placed in contact with a sapphire window embedded within an aluminum platform, with the backside of the window serving as a reference plane. The optimal distance between the reference plane and the paraboloidal mirror was determined by maximizing the intralipid Raman signal from a tissue phantom contained in a sample holder simulating the dog ear geometry, i.e., a cylindrical tissue phantom solution with optical properties and a thickness close to the dog ear. Figure 2 shows an aluminum sample stage where a dog subject can lie on its stomach with its ear positioned over the sapphire window aperture. 2.2.Experimental ProtocolThe dog was maintained under an isoflourane inhaled anesthetic for the duration of the study. The blood glucose levels were clamped at eight levels for a period of 45 min at each level. These glucose levels were achieved and maintained by infusing 20% dextrose and insulin into ear veins. Blood samples were drawn every 5 min and were analyzed using a glucose analyzer. The temperature of the ear was constantly maintained with a closed-loop thermoelectric cooling temperature control of the plate containing the sapphire window. The dog’s glucose concentration was clamped at eight different levels within the range 5.6 to 25.6 mM (100 to ). Each clamping level lasted for approximately 35 min. During the course of the experiment, Raman spectra were continuously collected with 1.8 s per frame and a 1.6 s data transfer time, giving a frame every 3.4 s. (The duty cycle was limited by data transfer.) The laser was not shuttered during the file transfer. Each frame has pixel dimensions as hardware binning of every five vertical pixels was chosen. After data collection, the curvature correction algorithm was applied to all frames before vertical binning.37 Since various frame-averaging schemes were adopted, the individual spectra are referred to as “frames” though they are one-dimensional, and the subsequent averaged spectra are referred to as “sample spectra.” Figure 3 shows the examples of the 33-frame averaged sample spectra with between successive spectra. Apparent sapphire Raman peaks and a broadband decreasing background are observed. To better accentuate Raman peaks from the dog ear, a fifth-order polynomial background subtraction routine was employed and the background removed spectrum is also shown in Fig. 3 with prominent sapphire peaks. 3.Results and Discussions3.1.Minimum Detection Error AnalysisFor a nearly shot-noise limited spectrum measurement, the minimum detection error can be estimated using experimental parameters such as signal-to-noise ratio (SNR) and an overlap factor.15 To estimate spectral random noise, we calculated the variance of each pixel among 10 adjacent frames (frame 6485 to 6494). The rational for selecting these frames is to minimize the apparent variance owing to the background decay. The decay was observed to diminish with time. The calculated two-dimensional variance map was then processed by the curvature correction algorithm previously described and a single spectrum of variance was obtained. The estimated noise value, 360, was obtained from the average across the square root of the variance spectrum. The Raman spectrum of glucose was obtained from a 50-mM glucose water solution contained in the dog-ear-like sapphire sample holder. The norm of the glucose signal was calculated to be using either a pixel range 240 to 1040 or 200 to 1200. The overlap factor for the experiment was estimated to be to 1.4 using the nine-component model described earlier. Using a previously developed theory, the minimum detection error based on these experimental parameters is to 9 mM (using raw frames).15 If frame averaging is performed, is 1.46 to 1.57 mM and 1.04 to 1.11 mM for 33- and 65- frame averaging, respectively. Note that the formalism considers only random noise in the predicted spectra, not the calibration spectra nor the reference concentration, i.e., an absolutely correct model. 3.2.PreprocessingAdditional preprocessing steps were implemented besides the background removal mentioned earlier. Among the 6498 frames, we observed that the laser intensity fluctuated at two fixed frequencies, causing fluctuations at the same frequencies in the collected frames. Fourier filtering was employed to effectively remove the slowly varying laser intensity fluctuations. Owing to the high SNR, the charge-coupled device fixed pattern noise was very significant. We first heavily smoothed the sample spectrum using a 101-point Savitzky–Golay filter and then subtracted the smoothed spectrum from the original sample spectra to identify the fixed pattern noise. The fixed pattern noise in individual frames was subsequently removed according to intensity levels. 3.3.Partial Least Squares Analysis with Cross ValidationVarious datasets were formed for PLS analysis using leave-one-out cross validation with differences in the following aspects: number of frames averaged; with or without 25-pt Savitzky–Golay smoothing; and spectral range selection. The results give us a general evaluation of the performance of our technique. Figures 4(a) to 4(d) show example results from one analysis with 33-frame averaging, 25-pt smoothing, and the plasma glucose as the reference concentration. Figure 4(a) shows the calculated root-mean-square error of the cross validation (RMSECV) versus number of PLS factors. The observed minimum indicates that the optimal calibration model contains eight factors. The Clark’s error grid is plotted in Fig. 4(b) using the predicted concentrations obtained from the cross-validation procedure. This type of grid analysis is used by physicians to evaluate the performance of glucose analysis technologies. Predictions falling in zones A and B are clinically considered acceptable. Figure 4(c) compares the reference to the predicted glucose concentration over time ( between two samples). The regression vector and glucose Raman spectrum are plotted in Fig. 4(d). Distinctive similarities are observed between the two, indicating that glucose was indeed measured because there was no prior glucose spectral information supplied to the PLS model. Table 1 lists all results from the cross-validation analyses with various calibration set formations. Table 1Cross-validation analysis with various preprocessing and parameters.
3.4.Partial Least Squares Analysis with Cross Validation and PredictionWe then picked one set of parameters, i.e., 33-frame averaging and 25-pt smoothing, to perform further analysis with level splitting. Because the 65-frame averaging scheme did not give much improved RMSECV and results in fewer samples, analyses here were done using a 33-frame averaging. All the samples collected at the clamping levels were divided into a calibration set and a prediction set. Building calibration models solely based on the leveled data avoids additional confounding factors during the glucose rise and fall phases. PLS was performed on the calibration set to calculate RMSECV and the vector, which was subsequently used to predict on the prediction set with the root-mean-square error of prediction (RMSEP). The vector was then used to predict on all samples except the calibration samples, excluding the samples during glucose rise and fall phases, to calculate , and including the samples during glucose rise or fall phases to calculate and . Table 2 lists all results from the level-splitting analysis. The higher values observed in suggest that the divergence between ISF and plasma glucose is significant during the glucose rise and fall phases, which could not be corrected for using the calibration models based on leveled regions. These RMSEP values agree with the minimum detection error estimated earlier. Table 2Level-splitting analysis with various preprocessing and parameters.
The next analysis was to form the calibration set with one level entirely left out, and then predict on the left-out level () and all samples not included in the calibration set (). Results are summarized in Table 3. It is observed that the for level 1 is much higher than for other levels. This is because fluorescence photobleaching was most significant during that time and also the instrument and experimental subject needed a warm up time. Table 3Leave-one-level-out analysis with various preprocessing and parameters.
Finally, two randomized concentration profiles were used to demonstrate that the previous calibration models are indeed predictive. In the first case, random concentrations in the experimental range were paired with measured spectra. In the second case, the order of the reference concentration measurements was randomly scrambled. As shown in Table 4, result from these tests suggests RMSEP with a model that lacks prediction capability. Therefore, results from previous calibration models are predictive for glucose concentration. Table 4Randomized concentration analysis.
4.ConclusionsThis paper describes an in vivo survival dog study performed on a beagle anesthetized for , during which its blood glucose concentration was clamped at several different levels. A glucose clamping study allows better disentangling of systematic effects from real glucose changes. Raman spectra were continuously acquired from the ear and reference blood glucose measurements were taken via venous blood draw. Using only the level data, RMSEP on the order of 1.5 to 2 mM (10% of the average concentration) was obtained, agreeing with the minimum detection error analysis. This RMSEP is higher than needed for diabetic patients. However, the average glucose level in the dog study () was significantly higher than the fasting concentration in human subjects. Since it is critically important to accurately determine glucose concentration near hypoglycemia, future experimental designs will include lower clamping concentrations. Distinctive similarities were observed between the resulting vector and the glucose Raman spectrum measured in water, indicating that glucose was indeed measured. We have identified the photobleaching of tissue autofluorescence to be a key error source from the results of the leave-level 1-out analysis. A potential strategy for a mitigating scheme was suggested.38 Results from this study demonstrate the feasibility of detecting glucose in vivo using transcutaneous Raman spectroscopy. In addition, the analyses and results provide valuable insights for improving our technique for future studies. AcknowledgmentsWCS acknowledges funding from the National Science Foundation (NSF) CAREER Award (No. CBET-1151154), the National Aeronautics and Space Administration (NASA) Early Career Faculty Grant (No. NNX12AQ44G), Gulf of Mexico Research Initiative (GoMRI-030), and Cullen College of Engineering at the University of Houston. This work was initially performed at the MIT Laser Biomedical Research Center supported by the NIH National Center for Research Resources, Grant No. P41-RR02594. We thank Professor Michael S. Feld (1940 to 2010) and Dr. Ramachandra R. Dasari for their mentorship and guidance. ReferencesW. C. Shih, K. L. Bechtel and M. S. Feld, Handbook Of Optical Sensing Of Glucose In Biological Fluids And Tissues, CRC Press, New York
(2008). Google Scholar
V. Tuchin, Tissue Optics: Light Scattering Methods and Instruments for Medical Diagnosis, 2nd ed.SPIE Press Monograph, Bellingham, Washington
(2007). Google Scholar
M. G. Ghosn, V. V. Tuchin and K. V. Larin,
“Depth-resolved monitoring of glucose diffusion in tissues by using optical coherence tomography,”
Opt. Lett., 31
(15), 2314
–2316
(2006). http://dx.doi.org/10.1364/OL.31.002314 OPLEDP 0146-9592 Google Scholar
K. L. Bechtel, W. C. Shih and M. S. Feld,
“Intrinsic Raman spectroscopy for quantitative biological spectroscopy Part II: experimental applications,”
Opt. Express, 16
(17), 12737
–12745
(2008). http://dx.doi.org/10.1364/OE.16.012737 OPEXFF 1094-4087 Google Scholar
W. C. Shih, K. L. Bechtel and M. S. Feld,
“Intrinsic Raman spectroscopy for quantitative biological spectroscopy Part I: theory and simulations,”
Opt. Express, 16
(17), 12726
–12736
(2008). http://dx.doi.org/10.1364/OE.16.012726 OPEXFF 1094-4087 Google Scholar
J. Zeng et al.,
“Analysis of ethyl and methyl centralite vibrational spectra for mapping organic gunshot residues,”
Analyst, 139
(17), 4270
–4278
(2014). http://dx.doi.org/10.1039/C4AN00657G ANLYAG 0365-4885 Google Scholar
G. J. Puppels et al.,
“Studying single living cells and chromosomes by confocal Raman microspectroscopy,”
Nature, 347
(6290), 301
–303
(1990). http://dx.doi.org/10.1038/347301a0 NATUAS 0028-0836 Google Scholar
A. Mahadevan-Jansen et al.,
“Near-infrared Raman spectroscopy for in vitro detection of cervical precancers,”
Photochem. Photobiol., 68
(1), 123
–132
(1998). http://dx.doi.org/10.1111/php.1998.68.issue-1 PHCBAP 0031-8655 Google Scholar
A. J. Berger et al.,
“Multicomponent blood analysis by near-infrared Raman spectroscopy,”
Appl. Opt., 38
(13), 2916
–2926
(1999). http://dx.doi.org/10.1364/AO.38.002916 APOPAI 0003-6935 Google Scholar
K. Maquelin et al.,
“Raman spectroscopic method for identification of clinically relevant microorganisms growing on solid culture medium,”
Anal. Chem., 72
(1), 12
–19
(2000). http://dx.doi.org/10.1021/ac991011h ANCHAM 0003-2700 Google Scholar
H. Ding et al.,
“Development of Raman spectral markers to assess metastatic bone in breast cancer,”
J. Biomed. Opt., 19
(11), 111606
(2014). http://dx.doi.org/10.1117/1.JBO.19.11.111606 JBOPFO 1083-3668 Google Scholar
N. Sudheendran et al.,
“Line-scan Raman microscopy complements optical coherence tomography for tumor boundary detection,”
Laser Phys. Lett., 11
(10), 105602
(2014). http://dx.doi.org/10.1088/1612-2011/11/10/105602 1612-2011 Google Scholar
Y. Li et al.,
“Micro-Raman spectroscopy study of cancerous and normal nasopharyngeal tissues,”
J. Biomed. Opt., 18
(2), 027003
(2013). http://dx.doi.org/10.1117/1.JBO.18.2.027003 JBOPFO 1083-3668 Google Scholar
A. M. K. Enejder et al.,
“Raman spectroscopy for noninvasive glucose measurements,”
J. Biomed. Opt., 10
(3), 031114
(2005). http://dx.doi.org/10.1117/1.1920212 JBOPFO 1083-3668 Google Scholar
O. R. Scepanovic et al.,
“Determination of uncertainty in parameters extracted from single spectroscopic measurements,”
J. Biomed. Opt., 12
(6), 064012
(2007). http://dx.doi.org/10.1117/1.2815692 JBOPFO 1083-3668 Google Scholar
W. C. Shih, K. L. Bechtel and M. S. Feld,
“Constrained regularization: hybrid method for multivariate calibration,”
Anal. Chem., 79
(1), 234
–239
(2007). http://dx.doi.org/10.1021/ac060732v ANCHAM 0003-2700 Google Scholar
W.-C. Shih, K. Bechtel and M. Feld, In Vivo Glucose Measurements, John Wiley & Sons, Inc.(2009). Google Scholar
J. Qi and W.-C. Shih,
“Parallel Raman microspectroscopy using programmable multi-point illumination,”
Opt. Lett., 37
(8), 1289
–1291
(2012). http://dx.doi.org/10.1364/OL.37.001289 OPLEDP 0146-9592 Google Scholar
J. Qi, J. Li and W.-C. Shih,
“High-speed hyperspectral Raman imaging for label-free compositional microanalysis,”
Biomed. Opt. Express, 4
(11), 2376
–2382
(2013). http://dx.doi.org/10.1364/BOE.4.002376 BOEICL 2156-7085 Google Scholar
J. Qi and W. C. Shih,
“Performance of line-scan Raman microscopy (LSRM) for high-throughput chemical imaging of cell population,”
Appl. Opt., 53
(13), 2881
–2885
(2014). http://dx.doi.org/10.1364/AO.53.002881 APOPAI 0003-6935 Google Scholar
Z. W. Huang et al.,
“Rapid near-infrared Raman spectroscopy system for real-time in vivo skin measurements,”
Opt. Lett., 26
(22), 1782
–1784
(2001). http://dx.doi.org/10.1364/OL.26.001782 OPLEDP 0146-9592 Google Scholar
M. S. Bergholt, W. Zheng and Z. Huang,
“Development of a multiplexing fingerprint and high wavenumber Raman spectroscopy technique for real-time in vivo tissue Raman measurements at endoscopy,”
J. Biomed. Opt., 18
(3), 030502
(2013). http://dx.doi.org/10.1117/1.JBO.18.3.030502 JBOPFO 1083-3668 Google Scholar
M. M. P. Arnob et al.,
“Laser rapid thermal annealing enables tunable plasmonics in nanoporous gold nanoparticles,”
Nanoscale, 6
(21), 12470
–12475
(2014). http://dx.doi.org/10.1039/C4NR03672G 1556-276X Google Scholar
F. Zhao et al.,
“Monolithic NPG nanoparticles with large surface area, tunable plasmonics, and high-density internal hot-spots,”
Nanoscale, 6 8199
–8207
(2014). http://dx.doi.org/10.1039/c4nr01645a 1556-276X Google Scholar
J. Zneg et al.,
“Internal and external morphology-dependent plasmonic resonance in monolithic nanoporous gold nanoparticles,”
RSC Adv., 4 36682
–36688
(2014). http://dx.doi.org/10.1039/C4RA04922E RSCACL 2046-2069 Google Scholar
J. Qi et al.,
“Surface-enhanced Raman spectroscopy with monolithic nanoporous gold disk substrates,”
Nanoscale, 5
(10), 4105
–4109
(2013). http://dx.doi.org/10.1039/c2nr33242f 1556-276X Google Scholar
H. Wang et al.,
“Quantitative analysis of creatinine in urine by metalized nanostructured parylene,”
J. Biomed. Opt., 15
(2), 027004
(2010). http://dx.doi.org/10.1117/1.3369002 JBOPFO 1083-3668 Google Scholar
M. Li et al.,
“Microfluidic surface-enhanced Raman scattering sensor with monolithically integrated nanoporous gold disk arrays for rapid and label-free biomolecular detection,”
J. Biomed. Opt., 19
(11), 111611
(2014). http://dx.doi.org/10.1117/1.JBO.19.11.111611 JBOPFO 1083-3668 Google Scholar
K. W. Kho et al.,
“Polymer-based microfluidics with surface-enhanced Raman-spectroscopy-active periodic metal nanostructures for biofluid analysis,”
J. Biomed. Opt., 13
(5), 054026
(2008). http://dx.doi.org/10.1117/1.2976140 JBOPFO 1083-3668 Google Scholar
J. Qi et al.,
“Label-free, in situ SERS monitoring of individual DNA hybridization in microfluidics,”
Nanoscale, 6
(15), 8521
–8526
(2014). http://dx.doi.org/10.1039/C4NR01951B 1556-276X Google Scholar
P. Z. McVeigh et al.,
“Widefield quantitative multiplex surface enhanced Raman scattering imaging in vivo,”
J. Biomed. Opt., 18
(4), 046011
(2013). http://dx.doi.org/10.1117/1.JBO.18.4.046011 JBOPFO 1083-3668 Google Scholar
M. Li et al.,
“Stamping surface-enhanced Raman spectroscopy for label-free, multiplexed, molecular sensing and imaging,”
J. Biomed. Opt., 19
(5), 050501
(2014). http://dx.doi.org/10.1117/1.JBO.19.5.050501 JBOPFO 1083-3668 Google Scholar
S. M. Stranahan and K. A. Willets,
“Super-resolution optical imaging of single-molecule SERS hot spots,”
Nano Lett., 10
(9), 3777
–3784
(2010). http://dx.doi.org/10.1021/nl102559d NALEFD 1530-6984 Google Scholar
J. L. Lambert et al.,
“Measurement of aqueous glucose in a model anterior chamber using Raman spectroscopy,”
J. Raman Spectrosc., 33
(7), 524
–529
(2002). http://dx.doi.org/10.1002/(ISSN)1097-4555 JRSPAF 0377-0486 Google Scholar
P. J. Caspers, G. W. Lucassen and G. J. Puppels,
“Combined in vivo confocal Raman spectroscopy and confocal microscopy of human skin,”
Biophys. J., 85
(1), 572
–580
(2003). http://dx.doi.org/10.1016/S0006-3495(03)74501-9 BIOJAU 0006-3495 Google Scholar
C. A. Lieber et al.,
“Raman microspectroscopy for skin cancer detection in vitro,”
J. Biomed. Opt., 13
(2), 024013
(2008). http://dx.doi.org/10.1117/1.2899155 JBOPFO 1083-3668 Google Scholar
J. Qi, K. L. Bechtel and W.-C. Shih,
“Automated image curvature assessment and correction for high-throughput Raman spectroscopy and microscopy,”
Biomed. Spectrosc. Imaging, 3
(4), 359
–368
(2014). http://dx.doi.org/10.3233/BSI-140092 BSIIAX 2212-8794 Google Scholar
H. Wang et al.,
“Improving skin Raman spectral quality by fluorescence photobleaching,”
Photodiagn. Photodyn. Ther., 9
(4), 299
–302
(2012). http://dx.doi.org/10.1016/j.pdpdt.2012.02.001 PPTHBF 1572-1000 Google Scholar
|