JBO Letters

Transcutaneous fiber optic Raman spectroscopy of bone using annular illumination and a circular array of collection fibers

[+] Author Affiliations
Matthew V. Schulmerich, Kathryn A. Dooley, Michael D. Morris

University of Michigan, Department of Chemistry, Ann Arbor, Michigan 48109-1055

Thomas M. Vanasse, Steven A. Goldstein

University of Michigan Medical School, Orthopaedic Research Laboratories, Ann Arbor, MI 48109-2200

J. Biomed. Opt. 11(6), 060502 (December 13, 2006). doi:10.1117/1.2400233
History: Received August 01, 2006; Revised September 06, 2006; Accepted September 18, 2006; Published December 13, 2006; Online December 13, 2006
Text Size: A A A

Open Access Open Access

* Tel: +1-734-764-7360; E-mail: mdmorris@umich.edu

Transcutaneous bone Raman spectroscopy with an exciting annulus of 785-nm laser light surrounding the field of view of a circular array of collection fibers is demonstrated. The configuration provides distributed laser light. The annulus is located 2to3mm beyond the edge of the field of view of the collection fibers to reject contributions from skin and other overlying tissues. Data are presented for rat and chicken tissue. For rat tibia, the carbonate/phosphate ratio measured at a depth of 1mm below the skin is in error by 2.3% at an integration time of 120s and within 10% at a 30-s integration time. For chicken tibia 4mm below the skin surface, the error is less than 8% with a 120-s integration time.

Figures in this Article

Dual energy x-ray absorptiometry (DEXA) has revolutionized the diagnosis and treatment of osteoporosis, a disorder characterized by loss of bone mass and quality resulting in more than 1 million fractures per year in the United States, causing disability and even death. DEXA measures only bone mineral density and not bone quality, a term that encompasses bone architecture, material properties, and remodeling dynamics. It thus predicts only 60 to 70% of fracture risk. There is significant overlap in measurements of bone mineral density between individuals with and without fractures.13 DEXA also grossly underpredicts the influence of drugs used to treat osteoporosis.4

The limitations of DEXA have fueled the search for new diagnostic and monitoring technologies. We used Raman microscopic imaging of trabecular bone (density and age matched) from the proximal femur to distinguish between tissue from females who had suffered osteoporotic fracture and those who died from causes unrelated to bone tissue problems.5 Iliac crest biopsies from females with and without the clinical diagnosis of osteoporosis were also evaluated. The best predictor for osteoporosis was the carbonate/phosphate ratio at the periosteal surface of the biopsy specimens. Matrix bands were less predictive because the archived specimens were fixed and embedded in polymethyl methacrylate which can partially denature proteins.67

The role of diffusive light transport in highly scattering systems is widely known in biomedical optics,89 but was not fully appreciated by the Raman community until recently.1012 As in other spectroscopies, fiber optic Raman probes can be used to gather subsurface information from highly scattering specimens if the detecting fibers are located some distance from the point(s) of introduction of the exciting laser light.13

We have used a probe with distributed laser light (250to270mW) and an array of collection fibers14 to recover bone tissue spectra in cadaveric human tissue15 at depths of 3mm. We have used glycerol optical clearing16 to improve transcutaneous signal recovery.17 Using a probe with two concentric rings of collection fibers and a single illumination fiber Matousek et al. demonstrated in vivo transcutaneous Raman spectroscopy of phalanges of human volunteers.18

Here we demonstrate transcutaneous Raman spectroscopy using a fiber optic probe with an axicon/lens doublet19 that distributes laser power in an annulus surrounding the field of view of a circular array of 50 collection fibers. In our previously reported global illumination configuration, which enabled subsurface mapping, coincidence of the illumination region and the field of view emphasized Raman scatter arising from the surface or close to the surface of the specimen. A completely offset collection region results in a higher ratio of subsurface to surface spectra,13,20 enhancing contrast at the cost of losing subsurface spatial information.

Figure 1 is a schematic of the probe. A 785-nm diode laser (Invictus, Kaiser Optical Systems, Ann Arbor, Michigan) is focused into a 62.5-μm core-fused silica fiber, which is collimated and focused into a fused silica axicon (cone angle 179.5deg, Delmar Photonics, San Diego, California). A planoconvex and negative-focal-length fused silica lenses focus the beam to an annulus. A dichroic mirror aligns the annulus to the field of view of the collection fibers (PhAT probe, Kaiser Optical Systems). The assembly is mounted on the distal end of the PhAT probe. The PhAT probe output is presented to an imaging spectrograph (HoloSpec, Kaiser Optical Systems). The dispersed Raman scatter is imaged onto a 1024×256 back-illuminated deep-depletion CCD (Andor Classic, Andor Technologies, Belfast, United Kingdom) and analyzed by band target entropy minimization (BTEM).14 Band height ratios were calculated in GRAMS/AI (Thermo Galactic, Madison, Wisconsin). At the probe head the laser beam was about 110mW (chicken) and 170mW (rat).

Grahic Jump LocationF1 :

Schematic of fiber optic Raman system. The components are as marked in the figure. Lenses of different focal lengths were used to adjust dimensions of illuminating annulus as needed. The drawing is not to scale.

Sprague-Dawley rat tibiae and femora were used to study integration time dependence. Measurements were made at an overlying skin thickness of approximately 1 mm using integration times between 5 and 120 s. Depilatory lotion was used to remove hair from the rat skin. Chicken tibiae were used to study penetration at a depth of 2 and 4 mm. For all specimens, tissues were treated with an optical clearing agent, glycerol, for approximately 1 min before the start of each experiment, and after completing a sequence of measurements, the overlying tissue was removed and Raman spectra of the exposed bone were collected with the PhAT probe.

Figure 2 shows recovered rat tibia spectra as a function of integration time. The spectra are presented as recovered by BTEM without low-pass filtering or smoothing. The mineral spectrum recovery is good above about 30s. The intense phosphate ν1 band is visible, though distorted, even at integration times as low as 5s. As Table 1 shows, it is possible to obtain reasonable agreement of carbonate/phosphate ratios between measurements made on exposed bone and through 1mm of skin at integration times of 30s or longer. BTEM did not completely separate the skin spectrum from the bone spectrum. Amide III and the CH2 bending mode at ca. 1446cm1 are about four times larger in the recovered spectra than in the spectrum of exposed bone. Incomplete separation results from oversampling the skin spectrum. In these experiments the annulus outer and inner diameters were about 7 and 6mm and the collection fiber field of view about 3mm. It was not possible with the available optics to generate a collection diameter below 3mm. Increasing the annulus i.d. to larger than 7mm to reduce contribution from skin caused a large fraction of laser light to miss the tibia entirely.

Table Grahic Jump Location
Carbonate/phosphate ratios measured transcutaneously and on exposed bone.
Table Footer NoteCalculated as height ratio of 1070to958cm1 bands.
Table Footer NoteRatio calculated using data from all collection fibers.
Table Footer NoteRatio calculated without using data from the outer ring of fibers.

Grahic Jump LocationF2 :

Transcutaneous Raman spectra of rat tibia. The spectrum integration times are shown in the figure. The bands are phosphate, 958cm1; carbonate, 1070cm1; amide III, 1240to1270cm1; matrix CH2 deformation, 1446cm1.

We collected spectra of chicken tibiae using the same excitation annulus. Figure 3 shows the recovered and corresponding exposed bone spectra at depth of 4 mm (120-s integration time), and also the average spectrum (dotted line) collected by the probe. The spectra have not been smoothed or filtered. Figure 3 shows the effect of using all 50 collection fibers for BTEM calculations. Comparison of the recovered spectrum with the exposed bone spectrum shows incomplete separation of the bone factor from the overlying tissue spectrum. This is most easily seen as the erroneously high intensity of the CH21446-cm1 band. As Table 1 shows, at a 4-mm depth, the error in measurement of carbonate/phosphate is 11.8%. If the outer ring of 18 collection fibers is excluded from the calculation, the recovery of the bone factor is improved [Fig. 3] and the error is reduced to 7.6% because the operation is equivalent to increasing the distance between excitation annulus and collection fibers. SNR is not improved because the number of spectra included in the calculation is reduced.

Grahic Jump LocationF3 :

Measurements made through 4mm of overlying tissue on a chicken tibia at the mid diaphysis. Transcutaneous (dotted), recovered bone factor (gray), exposed bone (black): (a) recovered bone factor using data from all 50 collection fibers and (b) recovered bone factor using data from the 32 innermost collection fibers.

These preliminary results can be substantially improved by optimization of the design based on realistic modeling of the tissue optics of these complex skin, bone, and connective tissue systems. While much remains to be done, Raman spectroscopy shows promise as a complement to DEXA. Periosteal surfaces may be probed noninvasively up to a few millimeters below the skin surface (e.g., distal radius, vertebrae), or minimally invasive if deeper (e.g., iliac crest). Similarly, Raman spectroscopy may be useful in diagnosis and monitoring of patients with other bone diseases or genetic defects, most of which affect bone quality.

Note that after this manuscript was submitted we learned that a similar annular illumination approach has been under independent development by P. Matousek with a manuscript reporting results in press.21

Acknowledgments

This research supported in part by the University of Michigan Musculoskeletal Core Research Center through National Institutes of Health (NIH) Grant No. P30 AR46024.

Buchanan  J. R., , Myers  C., , Greer  R. B., , Lloyd  T., , and Varano  L. A., “ Assessment of the risk of vertebral fracture in menopausal women. ,” J. Bone Jt. Surg., Am. Vol..  0021-9355 69, , 212–218  ((1987)).
Riggs  B. L., , Wahner  H. W., , Dunn  W. L., , Mazess  R. B., , Offord  K., , and Melton  L. J., “ Differential changes in bone mineral density of the appendicular and axial skeleton with aging: Relationship to spinal osteoporosis. ,” J. Clin. Invest..  0021-9738 67, , 328–335  ((1981)).
Riggs  B. L., , Wahner  H. W., , Seeman  E., , Offord  K. P., , Dunn  W. L., , Mazess  R. B., , Johnson  K. A., , and Melton  L. J. d., “ Changes in bone mineral density of the proximal femur and spine with aging: differences between the posteneopausal and senine osteoporosis syndromes. ,” J. Clin. Invest..  0021-9738 67, , 328–335  ((1982)).
Delmas  P. D., and Seeman  E., “ Changes in bone mineral density explain little of the reduction in vertebral or nonvertebral fracture risk with anti-resorptive therapy. ,” Bone (N.Y.).  8756-3282 34, , 599–604  ((2004)).
McCreadie  B. R., , Morris  M. D., , Chen  T.-C., , Rao  D. S., , Finney  W. F., , Widjaja  E., , and Goldstein  S. A., “ Bone tissue compositional differences in women with and without osteoporotic fracture. ,” Bone (N.Y.).  8756-3282  (in press).
Pleshko  N. L., , Boskey  A. L., , and Mendelsohn  R., “ An FT-IR microscopic investigation of the effects of tissue preservation on bone. ,” Calcif. Tissue Int..  0171-967X 51, , 72–77  ((1992)).
Yeni  Y., , Yerramshetty  J., , Akkus  O., , Pechey  C., , and Les  C., “ Effect of fixation and embedding on Raman spectroscopic analysis of bone tissue. ,” Calcif. Tissue Int..  0171-967X 78, , 363–371  ((2006)).
Tuchin  V.,  Tissue Optics. Light Scattering Methods and Instruments for Medical Diagnosis. ,  SPIE Press ,  Bellingham, WA  ((2000)).
Vo-Dinh  T.,  Biomedical Photonics Handbook. ,  CRC Press ,  Boca Raton, FL  ((2003)).
Everall  N., , Hahn  T., , Matousek  P., , and Parker  A. W., “ Photon migration in Raman spectroscopy. ,” Appl. Spectrosc..  0003-7028 58, , 591–597  ((2004)).
Everall  N., , Hahn  T., , Matousek  P., , Parker  A. W., , and Towrie  M., “ Picosecond time-resolved Raman spectroscopy of solids: capabilities and limitations for fluorescence rejection and the influence of diffuse reflectance. ,” Appl. Spectrosc..  0003-7028 55, , 1701–1708  ((2001)).
Matousek  P., , Everall  N., , Towrie  M., , and Parker  A. W., “ Depth profiling in diffusely scattering media using Raman spectroscopy and picosecond Kerr gating. ,” Appl. Spectrosc..  0003-7028 59, , 200–205  ((2005)).
Matousek  P., , Clark  I. P., , Draper  E., , Morris  M. D., , Goodship  A., , Everall  N., , Towrie  M., , Finney  W. F., , and Parker  A. W., “ Subsurface probing of diffusely scattering media using spatially offset Raman spectroscopy (SORS). ,” Appl. Spectrosc..  0003-7028 59, , 393–400  ((2005)).
Schulmerich  M. V., , Finney  W. F., , Fredericks  R. A., , and Morris  M. D., “ Subsurface Raman spectroscopy and mapping using a globally illuminated non-confocal fiber-optic array probe in the presence of Raman photon migration. ,” Appl. Spectrosc..  0003-7028 60, , 109–114  ((2006)).
Schulmerich  M. V., , Finney  W. F., , Popescu  V., , Morris  M. D., , Vanasse  T. M., , and Goldstein  S. A., “ Transcutaneous Raman spectroscopy of bone tissue using a nonconfocal fiber optic array probe. ,” Proc. SPIE.  0277-786X 6093, , 60930O  ((2006)).
Tuchin  V. V.,  Optical Clearing of Tissues and Blood. ,  SPIE Press ,  Bellingham, WA  ((2006)).
Schulmerich  M. V., , Morris  M. D., , Vanasse  T. M., , and Goldstein  S. A., “ Optical clearing in transcutaneous Raman spectroscopy of bone tissue. ” (submitted for publication).
Matousek  P., , Draper  E. R. C., , Goodship  A. E., , Clark  I. P., , Ronayne  K. L., , and Parker  A. W., “ Noninvasive Raman spectroscopy of human tissue in vivo. ” Appl. Spectrosc..  0003-7028 60, , 758–763  ((2006)).
Bélanger  P.-A., and Rioux  M., “ Ring pattern of a lens-axicon doublet illuminated by a Gaussian beam. ,” Appl. Opt..  0003-6935 17, , 1080–1086  ((1978)).
Matousek  P., , Morris  M. D., , Everall  N., , Clark  I. P., , Towrie  M., , Draper  E., , Goodship  A., , and Parker  A. W., “ Numerical simulations of subsurface probing in diffusely scattering media using spatially offset Raman spectroscopy. ,” Appl. Spectrosc..  0003-7028 59, , 1485–1492  ((2005)).
Matousek  P., “ Inverse spatially offset Raman spectroscopy for deep non-invasive probing of turbid media. ” Appl. Spectrosc.. (in press).
© 2006 Society of Photo-Optical Instrumentation Engineers

Citation

Matthew V. Schulmerich ; Kathryn A. Dooley ; Michael D. Morris ; Thomas M. Vanasse and Steven A. Goldstein
"Transcutaneous fiber optic Raman spectroscopy of bone using annular illumination and a circular array of collection fibers", J. Biomed. Opt. 11(6), 060502 (December 13, 2006). ; http://dx.doi.org/10.1117/1.2400233


Figures

Grahic Jump LocationF1 :

Schematic of fiber optic Raman system. The components are as marked in the figure. Lenses of different focal lengths were used to adjust dimensions of illuminating annulus as needed. The drawing is not to scale.

Grahic Jump LocationF2 :

Transcutaneous Raman spectra of rat tibia. The spectrum integration times are shown in the figure. The bands are phosphate, 958cm1; carbonate, 1070cm1; amide III, 1240to1270cm1; matrix CH2 deformation, 1446cm1.

Grahic Jump LocationF3 :

Measurements made through 4mm of overlying tissue on a chicken tibia at the mid diaphysis. Transcutaneous (dotted), recovered bone factor (gray), exposed bone (black): (a) recovered bone factor using data from all 50 collection fibers and (b) recovered bone factor using data from the 32 innermost collection fibers.

Tables

Table Grahic Jump Location
Carbonate/phosphate ratios measured transcutaneously and on exposed bone.
Table Footer NoteCalculated as height ratio of 1070to958cm1 bands.
Table Footer NoteRatio calculated using data from all collection fibers.
Table Footer NoteRatio calculated without using data from the outer ring of fibers.

References

Buchanan  J. R., , Myers  C., , Greer  R. B., , Lloyd  T., , and Varano  L. A., “ Assessment of the risk of vertebral fracture in menopausal women. ,” J. Bone Jt. Surg., Am. Vol..  0021-9355 69, , 212–218  ((1987)).
Riggs  B. L., , Wahner  H. W., , Dunn  W. L., , Mazess  R. B., , Offord  K., , and Melton  L. J., “ Differential changes in bone mineral density of the appendicular and axial skeleton with aging: Relationship to spinal osteoporosis. ,” J. Clin. Invest..  0021-9738 67, , 328–335  ((1981)).
Riggs  B. L., , Wahner  H. W., , Seeman  E., , Offord  K. P., , Dunn  W. L., , Mazess  R. B., , Johnson  K. A., , and Melton  L. J. d., “ Changes in bone mineral density of the proximal femur and spine with aging: differences between the posteneopausal and senine osteoporosis syndromes. ,” J. Clin. Invest..  0021-9738 67, , 328–335  ((1982)).
Delmas  P. D., and Seeman  E., “ Changes in bone mineral density explain little of the reduction in vertebral or nonvertebral fracture risk with anti-resorptive therapy. ,” Bone (N.Y.).  8756-3282 34, , 599–604  ((2004)).
McCreadie  B. R., , Morris  M. D., , Chen  T.-C., , Rao  D. S., , Finney  W. F., , Widjaja  E., , and Goldstein  S. A., “ Bone tissue compositional differences in women with and without osteoporotic fracture. ,” Bone (N.Y.).  8756-3282  (in press).
Pleshko  N. L., , Boskey  A. L., , and Mendelsohn  R., “ An FT-IR microscopic investigation of the effects of tissue preservation on bone. ,” Calcif. Tissue Int..  0171-967X 51, , 72–77  ((1992)).
Yeni  Y., , Yerramshetty  J., , Akkus  O., , Pechey  C., , and Les  C., “ Effect of fixation and embedding on Raman spectroscopic analysis of bone tissue. ,” Calcif. Tissue Int..  0171-967X 78, , 363–371  ((2006)).
Tuchin  V.,  Tissue Optics. Light Scattering Methods and Instruments for Medical Diagnosis. ,  SPIE Press ,  Bellingham, WA  ((2000)).
Vo-Dinh  T.,  Biomedical Photonics Handbook. ,  CRC Press ,  Boca Raton, FL  ((2003)).
Everall  N., , Hahn  T., , Matousek  P., , and Parker  A. W., “ Photon migration in Raman spectroscopy. ,” Appl. Spectrosc..  0003-7028 58, , 591–597  ((2004)).
Everall  N., , Hahn  T., , Matousek  P., , Parker  A. W., , and Towrie  M., “ Picosecond time-resolved Raman spectroscopy of solids: capabilities and limitations for fluorescence rejection and the influence of diffuse reflectance. ,” Appl. Spectrosc..  0003-7028 55, , 1701–1708  ((2001)).
Matousek  P., , Everall  N., , Towrie  M., , and Parker  A. W., “ Depth profiling in diffusely scattering media using Raman spectroscopy and picosecond Kerr gating. ,” Appl. Spectrosc..  0003-7028 59, , 200–205  ((2005)).
Matousek  P., , Clark  I. P., , Draper  E., , Morris  M. D., , Goodship  A., , Everall  N., , Towrie  M., , Finney  W. F., , and Parker  A. W., “ Subsurface probing of diffusely scattering media using spatially offset Raman spectroscopy (SORS). ,” Appl. Spectrosc..  0003-7028 59, , 393–400  ((2005)).
Schulmerich  M. V., , Finney  W. F., , Fredericks  R. A., , and Morris  M. D., “ Subsurface Raman spectroscopy and mapping using a globally illuminated non-confocal fiber-optic array probe in the presence of Raman photon migration. ,” Appl. Spectrosc..  0003-7028 60, , 109–114  ((2006)).
Schulmerich  M. V., , Finney  W. F., , Popescu  V., , Morris  M. D., , Vanasse  T. M., , and Goldstein  S. A., “ Transcutaneous Raman spectroscopy of bone tissue using a nonconfocal fiber optic array probe. ,” Proc. SPIE.  0277-786X 6093, , 60930O  ((2006)).
Tuchin  V. V.,  Optical Clearing of Tissues and Blood. ,  SPIE Press ,  Bellingham, WA  ((2006)).
Schulmerich  M. V., , Morris  M. D., , Vanasse  T. M., , and Goldstein  S. A., “ Optical clearing in transcutaneous Raman spectroscopy of bone tissue. ” (submitted for publication).
Matousek  P., , Draper  E. R. C., , Goodship  A. E., , Clark  I. P., , Ronayne  K. L., , and Parker  A. W., “ Noninvasive Raman spectroscopy of human tissue in vivo. ” Appl. Spectrosc..  0003-7028 60, , 758–763  ((2006)).
Bélanger  P.-A., and Rioux  M., “ Ring pattern of a lens-axicon doublet illuminated by a Gaussian beam. ,” Appl. Opt..  0003-6935 17, , 1080–1086  ((1978)).
Matousek  P., , Morris  M. D., , Everall  N., , Clark  I. P., , Towrie  M., , Draper  E., , Goodship  A., , and Parker  A. W., “ Numerical simulations of subsurface probing in diffusely scattering media using spatially offset Raman spectroscopy. ,” Appl. Spectrosc..  0003-7028 59, , 1485–1492  ((2005)).
Matousek  P., “ Inverse spatially offset Raman spectroscopy for deep non-invasive probing of turbid media. ” Appl. Spectrosc.. (in press).

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

PubMed Articles
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.