Despite recent advances in the ability to assess fracture risk,1,2 definitive metrics do not yet exist to identify individuals in need of an intervention that lowers fracture risk. Complementary to established x-ray-based diagnostics of bone, Raman spectroscopy (RS) is an emerging technology that offers nondestructive34.5.6.–7 measures of the biochemical nature of tissue. As an indication of its potential to assess fracture risk, RS detected differences in carbonate concentration relative to phosphate between bone samples from nonfracture patients and bones from osteoporotic fracture cases.8 In addition to quantifying the amount of carbonate in calcified tissue,9 RS is sensitive to local changes in mineral accumulation through mineral to collagen peak ratios, as well as local changes in mineral maturation through measurements of crystallinity.10 These properties become less heterogeneous with aging,11 change in response to tissue damage,12 and correlate to mechanical strength of rodent bones,13,14 as well as human cortical bone.10
While these attributes make RS a candidate for clinical diagnosis of bone quality and disease states, unresolved issues regarding instrument polarization and its impact on analysis hamper unambiguous derivation of quantities reflecting the biochemical properties of bone tissue, referred to henceforth as biomarkers. For example, probe-based instruments have been developed to acquire Raman spectra from bone through the overlying tissue.1516.–17 However, most RS studies assessing bone use laser confocal microscopes18 in which the laser is polarized.19,20 Since most fiber optic instruments do not preserve polarization, there is potential for significant discrepancy between the relevant biomarkers of fracture resistance as obtained from RS microscopes and those obtained from existing clinically relevant instruments.
Earlier studies in using Fourier transform infrared spectroscopy (FTIR) identified that vibrational spectroscopy markers for mineralization, crystallinity, carbonation, and collagen cross-linking all significantly associate with fragility fracture,21 and while the influence of polarization on FTIR was previously characterized,22,23 the polarization state for this instrument was not reported. FTIR and RS are linked as vibrational spectroscopy methods, but fundamental differences give RS an apparent clinical advantage that has fueled cross-correlation and validation in RS.
Recent correlation studies linking RS to the fracture resistance of bone have reported different biomarker sensitivity of RS to the biomechanical properties of bone, possibly due to instrument polarization differences. In one RS study involving a commercial confocal system, differences in Phosphate/-wag between trabecular and cortical bone were related to differences in nanomechanical properties between the tissue types.24In another study analyzing cortical bone from genetic mouse models involving matrix metalloproteinases (Mmp2-/- and Mmp9-/-) and using a similar commercial instrument, correlations among nanoindentation modulus, bending strength, and RS were reported for several mineral-to-collagen ratios, but not Phosphate/Amide I.25 On the other hand, strength of bone from vehicle- and glucocorticoid-treated mice was correlated to various RS peaks when normalized to Amide I in a fiber optic system.26 Despite differences in modes of biomechanical testing among these studies, the polarization state of the instrument likely influences which RS biomarkers are sensitive to experimental groups.
Even though bone is a birefringent material, only a few investigations have intentionally examined the effect of polarization on RS peaks of bone.2728.29.–30 The vast majority of polarization RS studies of mineralized tissues utilized isolator and analyzer polarizers to “fully polarize” both the input light and collection arm of the Raman system (Table 1).To the best of our knowledge, the extent to which polarization may affect various peak ratios used to assess composition of bone is not well understood for RS instruments without these added optics. Addressing this is important because Raman scattering bands are inherently and differentially affected by polarization due to vibrational modes that give rise to the Raman effect.31,32 Even though Legendre polynomials can be used to extract the distribution of collagen and mineral orientation by modulating polarization,30 inherent polarization within the instrument can affect bands even in the absence of molecular organization (e.g., analyzing carbon tetrachloride).33 Because sample volume, molecular organization, scattering anisotropy, and tissue turbidity influence the scattered light,30 it is difficult to predict how polarization affects RS spectra of bone tissue (which is both turbid and organized), especially when the input laser light is unaltered (Table 1). In the context of RS instrumentation, we refer to changes in Raman peak intensities due to bone rotation (relative to light polarization) as a “polarization bias.”
Previous polarized Raman studies used altered instrumentation to analyze mineralized tissues and collagen.
|Tissue||Instrumentation||Polarization||Laser (nm)||Factor(s)||Raman property||Reported trends||Ref.|
|Human tooth enamel||Spectrophotometer coupled to microscope with 100x and 50x objectives||Fully polarized||632.8, 514||Crystalline orientation||ν1PO4||Carbonation preserves the orientation of enamel crystals. Vibrational modes of enamel rods are reconciled to theory.||34|
|Human tooth enamel||Confocal microscope, 10x objective||Fully polarized||830||Caries and lesions||ν1PO4||Depolarization ratio and anisotropy of||35|
|Sheep tendon, fixed human femur||Confocal microscope, 100x objective||Not reported||633||Collagen Fiber distribution along osteon radius, and bone rotation||ν3PO4/ν1PO4||Polarization-sensitive peak ratios varied with distance from haversian canal.||29|
|Rat tail tendon collagen fibrils||Confocal microscope, 60x||Fully polarized||785||Mechanical stress||Amide I, Amide III||Stress changes orientation distribution but not vibration frequencies of collagen.||36|
|Whole and powdered mouse tibia||Custom confocal Raman microprobe, various objectives||Fully polarized||785||Osteogenesis Imperfecta (oim)||ν1PO4||Orientation distribution of mineral and collagen differ in oim mice.||30|
|Collagen from bovine achilles tendon||Confocal microscope 60x objective||Collection arm polarized||633||Collagen fiber orientation||Amide I, Amide III, and Amide III ratio (||Amide III doublet is sensitive to orientation.||33|
|Turkey leg tendon, embedded mouse femurs||Confocal microscope, 60x||Fully polarized||785||Orientation of cut, depolarization||ν1PO4/Amide I||Tissue age affects Raman polarization. Collagen orientation obtained from polarized mapping.||19|
|Embedded human femur osteons||Confocal microscope, 100x||Not reported||532||Tissue type and orientation||ν1PO4/Amide I||Tissue types have differential polarization sensitivity.||28|
|Femoral head spongy bone||Confocal microscope, 50x objective||Fully polarized||785||Specimens from osteoporotic fracture patients.||ν1PO4/Amide I||Polarized Raman helps determine the location of collagen fibers within trabeculae.||38|
|Human femur||Confocal microscope, 50x objective||Inherent polarization||785||Tissue type and orientation||Various||Polarization sensitivity inherent in peaks and Raman instrumentation can be minimized by optimizing peak phase.||Current Study|
Fully polarized indicates the use of polarizing optics in the input and collection arm.
Rather than intentionally polarizing Raman collection, the present study compares RS biomarkers of human cortical bone with respect to inherent system polarization using commercial instrumentation. The hypothesis of this study was that the intensity of polarization-sensitive Raman bands would oscillate relative to input polarization. Thus, optimization of the relative phase of these Raman bands’ oscillations could yield biomarkers that are better suited for the study of either composition (phase matched and less polarization sensitive) or structure (phase mismatched and polarization sensitive). This study uses polarization theory to quantify the phase and amplitude for a number of Raman peaks of bone, and in doing so, identifies phase matching as the source of polarization insensitivity in known and newly characterized peak ratios of bone composition. Whereas other studies have observed the intensity change of specific peak ratios (Table 1), the present work establishes phase profiles for many of the prominent peaks arising from bone with and without added polarization optics.
Transverse human cortical bone specimens from the lateral femur midshaft were prepared as per previously published methods.18 Briefly, bone samples were mounted to slides using cyanoacrylate and ground on silicon carbide papers of sequential grit, then polished with 0.05 μm alumina beads in solution to an ultimate surface area of and thickness of approximately 4 mm. One sample from each of 6 donors was used (4 males ages 48, 80, 82, and 94 and 2 females ages 86 and 95). To generate a control sample, a human molar was embedded in polymethylmethacrylate; a thick section was cut in the longitudinal direction; and the surfaced polished as previously described.18
To fully examine the influence of instrument polarization and bone structure on collected Raman spectra, we conducted several experiments, each with a different collection protocol or degree of polarization. Raman spectra were acquired from the polished surface of the bone tissue in air using a standard confocal Raman microscope (Ramanscope Mark III and InVia Raman Microscope, Renishaw, Hoffman Estates, Illinois) equipped with Renishaw EasyConfocal, a 35 μm slit opening, and a spectral resolution of , equipped with a 785 nm laser diode source with a polarization extinction ratio (PER) of (Innovative Photonic Solutions, Monmouth Junction, New Jersey). To eliminate grating bias according to Renishaw specifications, the polarization was aligned upright within the instrument (left-right when operator faces stage), confirmed with known polarizers and silicon standard intensity. Placing a mirror in the sample plane, the PER was also measured as after the dichroic and after the grating. Additional optics increased polarization of the Raman microscope, such that the system operated in a highly polarized regime. An isolator (NIR linear polarizer, extinction ratio, Thorlabs, Newton, New Jersey) was used to isolate a polarization angle of input laser light prior to sample incidence. An analyzer (additional linear polarizer, same specifications) isolates a particular polarization angle of light reflected off the sample. A quartz wedge depolarizer (AR coated achromatic depolarizer DPU-25-B, Thorlabs, Newton, New Jersey) effectively scrambles the polarization state of light in space prior to the spectral grating to prevent instrumentation bias by transmitting a pseudo-random polarized beam. Removal of the analyzer ( extinctionratio) decreased system polarization sensitivity, but retaining the input polarizer provided an “input polarization regime.” In this regime, the bone sample is rotated to examine bias and the depolarizer remains in the system to minimize instrumentation bias of the grating. Without added optics, the system retains a degree of inherent polarization sensitivity, henceforth referred to as an “unaltered polarization regime.”
To preserve system throughput across experiments despite differences in added optics, spectral acquisition exposure times were scaled to ensure 480 mW·s apparent exposure at the sample. This provided a signal to noise ratio (SNR) for the low intensity Proline peak in excess of in highly polarized experiments, translating to at least in unaltered experiments. Unless otherwise stated, spectra were obtained with 3 accumulations after 5 s photobleaching. Spectra were then binned to a resolution of , and processed via least squares modified polynomial fit39 and smoothed for noise using an 2nd order Savitsky-Golay filter.26,40 After fluorescence subtraction, a linear baseline subtraction (based on derivative zero-crossings neighboring the peak) was conducted on peaks that overlap with other constituents to ensure no residual fluorescence, namely Proline, Hydroxyproline, Phosphate, and Carbonate. Spatial resolution for each objective used was approximated via edge detection on a polished silicon standard. System Raman shift calibration was accomplished using a neon lamp and a silicon standard with Renishaw software to account for grating motion. Silicon measurements before and after each beam path change and at system “startup” ensured wavenumber calibration consistency.
Since dentin has less heterogeneity in collagen fibril orientation than bone, we collected Raman spectra from the same site as a human tooth rotated from 0 deg to 180 deg in 20 deg increments in order to characterize the polarization sensitivity of our RS instrument without additional polarization optics. In these dentin measurements, known polarization sensitive peaks oscillated through rotation with percent changes in mean normalized intensity of 6.6% and 22.6% for Phosphate and Amide I, respectively.
Highly polarized analyzer rotation
Polarization analysis used known bias from previous work19,28,37 to confirm the ability of Malus’s law to model phase and amplitude of Raman peaks. In effect, our first experiment was designed to evaluate phase oscillation for sensitive RS peaks. To account for within sample variation, five osteons and neighboring interstitial sites were selected from a single bone sample.18 In brief, selected osteons were spaced evenly over the surface and distributed by osteon size and pore size. Using upright input polarization through our 50×, numerical aperture objective (lateral resolution 3 to 4 μm, as measured by edge detection) and stationary bone orientation, the analyzer was rotated at each site from 0 deg to 180 deg in 20 deg increments. This study used an adaptation of Porto’s notations based on microscope translation stage directions since Porto’s notations traditionally depend on sample crystallographic axis,41,42 which varies within cortical bone. In this adaptation, the polarization regime is Z(XB)-Z, where B denotes analyzer rotation relative to instrument input X (always left-to-right as viewed by operator). Intensity for each of the prominent peaks in the bone spectrum was then modeled as a function of polarization angle B to compare the degree of oscillation between quantities.
Bone rotation for two polarization regimes
Next, spectra were collected as a function of bone orientation to evaluate peak and peak ratio sensitivity in less polarization sensitive systems. We analyzed a single osteon and neighboring interstitial site from each of three bones under both input polarized (added isolator and depolarizer) or unaltered (no added optics) polarization regimes. As was done with tooth, the bone sample was rotated around the optical axis using a custom stage to preserve collection location while obtaining spectra (50×, objective) from 0 deg to 180 deg rotation in 20 deg increments. The polarization regime is Z(X)-Z, where denotes bone rotation around Z relative to instrument input X (left-to-right as viewed by operator).
Spectral mapping of bone tissue rotation
Using the unaltered polarization regime, we acquired confocal Raman maps of spatial heterogeneity to demonstrate the effects of phase-matching on compositional discrimination of known osteonal and interstitial tissue differences. Phase-matching of peak ratios is defined as minimizing the phase difference of the ratio components, effectively choosing peaks that have the most similar rotation angle of maximum intensity, subsequently reducing the impact of rotation angle upon the observed ratio intensity. One osteon and the neighboring interstitial area (20×, objective, lateral resolution of 12 μm) from each of three bones was mapped using unaltered instrumentation at a pixel size of for 0 deg, 45 deg, and 90 deg rotations of the bone sample about the optical axis. To analyze discrimination of osteonal from interstitial tissue, intensity maps were generated for selected peak ratios applying a uniform scale based upon full intensity range, such that a polarization insensitive spectral constituent will show the same intensity image in all three acquisitions. Instrument polarization in direction X is denoted with X-Y stage directions in each figure panel. For one bone, the mapping process was repeated using the 50× objective for an osteonal-interstitial border within the original 20× map to demonstrate Raman maps of polarization bias with a smaller sample volume.
Data Modeling and Statistics
Data modeling and statistics were performed on peak heights extracted from each processed spectrum [Fig. 1(a)]. Peak intensities were modeled to Malus’ Law43,44 (intensity varies with polarization angle as a function of cosine squared) for phase and amplitude of oscillation [Fig. 1(b)]. The custom algorithm employed a least squares fit for amplitude nested inside a mean squared error driven optimization (Matlab implementation of Nelder-Mead simplex,45 Mathworks, Natick, Massachusetts), outputting peak phase, amplitude, and mean intensity as illustrated in Fig. 1(b).The degree of orientation sensitivity across the three generated polarization regimes was quantified for each prominent peak as a function of oscillation amplitude normalized to mean peak intensity. For less sensitive peaks, individual sample oscillations could become noisy or undetectable, such that data fails the underlying assumptions of the Malus’ law model. Modeled data were excluded from quantitative analysis if the model fit was not significant () via analysis of variance (ANOVA) regression (all fits shown in figures and tables are significant). For each peak, the number of samples with statistically significant models and the observable percent change in intensity measurements of the same quantity were recorded.
Phase Differences in Raman Peaks of Bone Under Highly Polarized Light
Acquired under a highly polarized regime, RS biomarker peaks exhibited differential polarization behavior in both degree and relative phase of intensity oscillation. For the most part, relative phase varied insignificantly between osteonal and interstitial tissue types for any given peak. However, phase oscillation varied distinctly between different peaks representing the same bone compositional element (i.e., Amide I at versus Amide III at , both biomarkers of collagen in Fig. 2).In reference to the ν1 Phosphate peak (), the strongest spectral signal for bone mineral, Proline () was found to have the best phase match for the generation of a mineral to collagen ratio, a metric commonly used as an indicator of bone quality.
Building upon our previous findings,18 the observed difference between osteonal and interstitial tissue composition (Fig. 3) was small (2% to 30% difference) relative to intensity change as a function of polarization angle (100% to 300% difference). However, under the traditional calculation of a mineral to collagen ratio using Phosphate (mineral) and Amide I (collagen) as biomarkers, different quantities would be observed at different polarization angles [e.g., 60 deg versus 140 deg in Fig. 3(a)]. As an alternative mineral to collagen ratio that still utilizes the signal strength of Phosphate, the phase-matched Proline peak can be used to represent collagen [Fig. 3(b)]. Also, as indicated by Kazanci et al.,28 other RS mineral quantities can be substituted for Phosphate [Fig. 3(c) and 3(d)]. The distinct phase mismatch between Phosphate and Amide I [Fig. 3(c)] was reversed by using Amide III for collagen [Fig. 3(d)].
Susceptibility of Certain Raman Peaks to Polarization Bias
When defined as the model amplitude normalized to mean peak intensity, the peak sensitivity to polarization decreased from the highly polarized regime to the input polarized and unaltered polarized regimes (Fig. 4). Hydroxyproline () and Amide I, the two most sensitive peaks in the highly polarized regime, remained polarization sensitive in the input polarized regime. Despite the fact that spectra were acquired from different bone samples, the oscillation sensitivity trends among most peaks remain consistent between the highly polarized and input polarized regimes. Comparing input polarized and unaltered polarization regimes (paired measurements of the same sample locations at the same rotation increments), polarization sensitivity dropped off markedly for some peaks like Carbonate and Hydroxyproline (Fig. 4). However, for other peaks like Amide I and Phosphate, degree of oscillation amplitude remains relatively unchanged. In the unaltered polarization regime, less sensitive peaks like Amide III fell into the noise floor, as evidenced by decrease in number of significant model fits by ANOVA regression (Table 2).
Peak sensitivity ranking as a percent change in intensity during bone rotation shows that some peaks still oscillate with unaltered system polarization.
|Highly polarizeda||Input polarized||Unaltered polarization|
|Raman peakb||% Changec||Sig. modelsd||% Change||Sig. models||% Change||Sig. models|
Derived from separate samples.
Rank is based on the degree of mean normalized intensity fluctuation for unaltered polarization.
% Change is defined as the range of observed intensities normalized to the mean intensity observed during specimen rotation (unaltered and input polarized) or during analyzer rotation (highly polarized) for significant (sig.) model fits. When more than one model fit is significant, the mean is given.
Number of significant models. Significance of model fit is defined as p<0.05 of the ANOVA regression (goodness of fit).
N/A: No significant models.
An RS surface plot for a single osteon acquired with unaltered polarization [Fig. 5(a)] illustrates that Phosphate peak intensity fluctuations [Fig. 5(b)] were out of phase with Amide I intensity fluctuations [Fig. 5(c)] but matched to the fluctuations of Proline [Fig. 5(d)]. Although noise has a significant impact on model fit in the unaltered regime, the trends of polarization phase between mineral and collagen peaks (Fig. 5) remained consistent with trends observed when the analyzer was rotated with the bone sample stationary (Fig. 3). Phase mismatch trends of RS biomarkers from highly polarized data persisted in unaltered polarization.
Performance of Phase-Matched Ratios for Compositional Differences
RS maps demonstrate how phase mismatch in RS peak ratios confounds the consistent measurement of spatial heterogeneity, even in an unaltered polarization regime (Fig. 6). Expected differences in mineral to collagen ratio between an osteon and surrounding interstitial tissue is not maintained throughout bone rotation for polarization sensitive Phosphate/Amide I [Fig. 6(b)]; whereas, Phosphate/Amide III [Fig 6(c)] shows consistent overall intensity differences between the tissue types despite rotation. Yet, this latter image is noisier than the former image due to significantly lower SNR of the Phosphate and Amide III peaks, relative to Phosphate. Maps of Phosphate/Proline [Fig. 6(d)] illustrate a relatively consistent image of compositional heterogeneity throughout rotation, differentiating the osteonal tissue from the more mineralized interstitial tissue. This peak ratio map is independent of bone rotation because of the low phase difference between Phosphate and Proline (Table 3).
Paired phase difference between selected Raman peak ratios and overall variance of peak ratios were estimated for bone and tooth rotation using the unaltered polarization instrument.
|Bone samples||Dentin control|
|Peak ratio||Representation||Phase difference (deg)a||Coefficient of variationb||Phase difference (deg)||Coefficient of variation|
|Carbonate/ν1 Phosphate||carbonate substitution||0.1||0.020||3.1||0.009|
|Carbonate/ν2 Phosphate||carbonate substitution||N/A||0.025||6.5||0.016|
|Carbonate/ν4 Phosphate||carbonate substitution||31.6||0.022||9. 3||0.022|
To obtain phase difference, the difference in phase between paired peaks (same spectra for the same site of rotation) was averaged among the multiple sites of data collection in which model fit was significant for both peaks.
To obtain coefficient of variation, the standard deviation per mean of the peak ratio for a given rotation was averaged across multiple sites of data collection.
N/A: No phase angle could be generated from significant models within the same sample. Significance of model fit is defined as p<0.05 of the ANOVA regression (goodness of fit).
Figure 7 shows how NA and subsequent differences in sample volume averaging affect the apparent sensitivity of mineral to matrix calculations to tissue type. The calculations compare ν1 Phosphate/Amide I sensitivity to the phase-matched Phosphate/Proline using a map at 20× magnification and a 50× map of a portion of the same area [Fig. 7(a)]. The polarization sensitive Amide I ratio produced a distinct intensity change at 50× magnification in the 45 deg map [Fig. 7(c)] that was less pronounced but arguably still apparent at 20× [Fig. 7(b)]. The mineral to matrix ratio with Proline [Fig. 7(d) and 7(e)] was relatively consistent throughout.
Development of RS methods towards bone diagnostics requires a firm understanding of which RS measures are sensitive to bone tissue composition and which are sensitive to bone tissue organization. While the potential for highly polarized RS to discriminate bone organization is known (Table 1), the present study provides a thorough characterization of the contribution of polarization bias in a standard Raman microscope to RS measurements of bone. When the goal of RS analysis is to assess compositional differences in bone, polarization bias adds uncertainty to the measurements. The addition of polarization optics, even to reduce polarization bias, leads to increased data collection time. Therefore, this study characterized peak and peak ratio specificity without altering instrumentation and found that polarization bias exists in a standard microscope and needs to be addressed.
Consistent with findings from the highly polarized RS analysis of mouse bone by Raghavan et al.,30 polarization bias can persist for a low NA objective () (Fig. 7) that effectively averages the signal over larger spatial volumes than a 50× objective with a NA of 0.75. Nonetheless, matching polarization-orientation phase, effectively matching the organizational component of RS peaks, allows for consistent measures of composition [Fig. 3(b) and 3(d), as well as Fig. 7(c) and 7(d)]. This study used spectral maps of spatial heterogeneity within bone (Fig. 6) to establish that peak ratio sensitivity identified by polarized RS studies (Table 1) remains in unaltered polarization regimes (i.e., standard confocal RS instruments). In effect, polarization phase can be exploited to distinguish compositional heterogeneity from organizational heterogeneity. Finally, throughput costs of added optics can be avoided by using less sensitive peaks and phase-optimized peak ratios.
Spectral fingerprints of disease may rely upon organization of collagen and mineral as much as the relative composition of these constituents. This is of particular importance because bone disorders and disease states including osteoporosis and osteogenesis imperfecta can involve deleterious changes in bone organization.30 Moreover, the organization of the constituents of bone tissue influences the fracture resistance of bone.4647.48.49.–50 The observations of phase difference (Fig. 2) from our modeling of spectral data as a function of analyzer rotation is consistent with findings from previous highly polarized studies,19,28,37 as well as recent theoretical models of collagen orientation within osteons.51 The phase differences between RS biomarker peaks arise from the organization of mineral crystals and collagen fibrils, such that information encoded within polarization phase may provide new insight for future disease diagnosis or fracture risk assessments. However, clear associations of phase and bone’s fracture resistance remain to be established. Nonetheless, as RS matures toward clinical use, consistency in discriminating composition from organization may contribute to accurate assessment of fracture risk.
Alternatively, RS could prove especially effective in the diagnosis of bone diseases that are pathologically based upon bone composition, including diabetes, chronic kidney disease, and the discrimination of grades of osteogenesis imperfecta. Recent investigations employing spatial sampling regimens to average out polarization bias concur that alternative mineral to collagen ratios (Fig. 6) exhibit increased sensitivity to known osteonal and interstitial differences in composition.18,52 Specifically, improvement in mineral to collagen ratio variance and microstructure delineation when utilizing Proline instead of Amide I as a collagen component can be explained by a decrease in underlying polarization phase difference between ν1 Phosphate (mineral constituent) and Proline (collagen constituent) (Table 3).
RS acquisition of bone spectra often spans a range of 300 to to cover prominent peaks, though it is expanded in some studies to to capture a CH peak. Collecting Proline, Phosphate, and Carbonate would require spanning only , further reducing necessary instrumentation, data processing, and collection time. Despite the availability of commercial systems, RS instrumentation differs largely between research groups and studies, implying that polarization bias from orientation sensitive peaks may vary between studies.
Changes in Raman spectra as the bone or tooth rotates relative to the incident laser light reveal the persistence of a polarization-orientation bias for sensitive peaks like Amide I even within less sensitive unaltered RS setups. Polarization sensitivity trends (Fig. 4) are conserved among prominent peaks in both the highly polarized regime and in the input polarized regime, despite the fact that spectra were acquired from different bone samples. Phase trends seen in highly polarized acquisition of bone appear to persist in unaltered instrumentation as well (Fig. 5). Phase-matched peak ratios of mineral to collagen ratio demonstrate lower coefficients of variation and therefore greater consistency (Table 3). When determined from the same site of bone rotation and then averaged across the various sites, the phase difference between Phosphate and Proline was 2.2 deg for bone and 9 deg for more highly organized dentin (Table 3), suggesting that these trends in phase difference may be conserved between tissues and anatomical locations. Despite the low intensity of Proline, the high intensity of the Phosphate peak may make Phosphate/Proline a more practical compositional metric than Phosphate/Amide III, which also has a low paired phased difference (Table 3). In addition, the use of peak phase difference confirmed the compositional nature of carbonate substitution (Table 3). While results suggest optimal metrics for bone composition and caution against possible inconsistent use of other metrics, the polarization-orientation information of RS biomarkers may have greater implications for future clinical bone diagnostics.
Consistent use of less polarization sensitive peaks or phase-matched ratios may allow for clearer comparisons between instruments and studies. RS sensitivity to glucocorticoid-treatment in rheumatoid arthritis bone shows compositional difference despite normalization to Amide I when using a fiber optic (polarization insensitive) system.26 These biomechanical correlations are likely separate and distinct from RS correlations to collagen tension changes seen in formal polarization analysis.36 Given low instrument polarization and results from less polarization sensitive carbonate and Amide III bands, analysis of bone from osteoarthritic patients on different load bearing surfaces can be interpreted as a largely compositional effect.53 Phase mismatch of Phosphate/ [see Figs. 2 and 6(e)] may have contributed to biomechanical correlation due to use of a commercial confocal system,24 thereby indicating a predominantly organizational phenotype. Interpretation of results from these and future studies in light of instrument polarization may help to define consistent Raman signatures for compositional and organizational disease.
Polarization-orientation information in bone biomarkers, as seen in highly polarized studies involving RS, persists in unaltered commercial systems with lower inherent sensitivity to polarization. Modeling this consistent bias shows that matched phase information between peaks yields biomarker ratios that are less sensitive to polarization-orientation, without the loss of throughput necessitated by additional optics. Bias in compositional measures can be minimized by phase matching; specifically, findings support using Phosphate/Proline for mineral to collagen and Phosphate for carbonate substitution. In the diagnosis of organizational phenotypes, polarization-orientation can be maximized by phase mismatch (i.e., Phosphate/Amide I) without necessarily including polarization optics. Optimizing polarization in the instrument and in biomarkers should help to increase discrimination and consistency in future studies of bone.
This material is based upon work supported by the Department of Veterans Affairs, Veterans Health Administration, Office of Research and Development, Biomedical Laboratory Research and Development. The authors would also like to acknowledge financial support from NSF Grant 1068988.
S. Amin et al., “Association of hip strength estimates by finite-element analysis with fractures in women and men,” J. Bone Miner. Res. 26(7), 1593–1600 (2011).JBMREJ0884-0431http://dx.doi.org/10.1002/jbmr.347Google Scholar
P. Matousek et al., “Subsurface probing in diffusely scattering media using spatially offset Raman spectroscopy,” Appl. Spectrosc. 59(4), 393–400 (2005).APSPA40003-7028http://dx.doi.org/10.1366/0003702053641450Google Scholar
M. V. Schulmerich et al., “Transcutaneous fiber optic Raman spectroscopy of bone using annular illumination and a circular array of collection fibers,” J. Biomed. Opt. 11(6), 060502 (2006).JBOPFO1083-3668http://dx.doi.org/10.1117/1.2400233Google Scholar
B. R. McCreadie et al., “Bone tissue compositional differences in women with and without osteoporotic fracture,” Bone 39(6), 1190–1195 (2006).BONEDL8756-3282http://dx.doi.org/10.1016/j.bone.2006.06.008Google Scholar
A. Awonusi, M. D. Morris and M. M. Tecklenburg, “Carbonate assignment and calibration in the Raman spectrum of apatite,” Calcif. Tissue Int. 81(1), 46–52 (2007).CTINDZ0171-967Xhttp://dx.doi.org/10.1007/s00223-007-9034-0Google Scholar
J. S. Yerramshetty and O. Akkus, “The associations between mineral crystallinity and the mechanical properties of human cortical bone,” Bone 42(3), 476–482 (2008).BONEDL8756-3282http://dx.doi.org/10.1016/j.bone.2007.12.001Google Scholar
J. S. Yerramshetty, C. Lind and O. Akkus, “The compositional and physicochemical homogeneity of male femoral cortex increases after the sixth decade,” Bone 39(6), 1236–1243 (2006).BONEDL8756-3282http://dx.doi.org/10.1016/j.bone.2006.06.002Google Scholar
J. A. Timlin et al., “Raman spectroscopic imaging markers for fatigue-related microdamage in bovine bone,” Anal. Chem. 72(10), 2229–2236 (2000).ANCHAM0003-2700http://dx.doi.org/10.1021/ac9913560Google Scholar
J. G. Ramasamy and O. Akkus, “Local variations in the micromechanical properties of mouse femur: the involvement of collagen fiber orientation and mineralization,” J. Biomech. 40(4), 910–918 (2007).JBMCB50021-9290http://dx.doi.org/10.1016/j.jbiomech.2006.03.002Google Scholar
O. Akkus, F. Adar and M. B. Schaffler, “Age-related changes in physicochemical properties of mineral crystals are related to impaired mechanical function of cortical bone,” Bone 34(3), 443–453 (2004).BONEDL8756-3282http://dx.doi.org/10.1016/j.bone.2003.11.003Google Scholar
E. R. Draper et al., “Novel assessment of bone using time-resolved transcutaneous Raman spectroscopy,” J. Bone Miner. Res. 20(11), 1968–1972 (2005).JBMREJ0884-0431http://dx.doi.org/10.1359/JBMR.050710Google Scholar
C. A. Patil et al., “Integrated system for combined Raman spectroscopy-spectral domain optical coherence tomography,” J. Biomed. Opt. 16(1), 011007 (2011).JBOPFO1083-3668http://dx.doi.org/10.1117/1.3520132Google Scholar
J. S. Nyman et al., “Measuring differences in compositional properties of bone tissue by confocal Raman spectroscopy,” Calcif. Tissue Int. 89(2), 111–122 (2011).CTINDZ0171-967Xhttp://dx.doi.org/10.1007/s00223-011-9497-xGoogle Scholar
S. Gamsjaeger et al., “Cortical bone composition and orientation as a function of animal and tissue age in mice by Raman spectroscopy,” Bone 47(2), 392–399 (2010).BONEDL8756-3282http://dx.doi.org/10.1016/j.bone.2010.04.608Google Scholar
C. B. Juang, L. Finzi and C. J. Bustamante, “Design and application of a computer-controlled confocal scanning differential polarization microscope,” Rev. Sci. Instrum. 59(11), 2399–2408 (1988).RSINAK0034-6748http://dx.doi.org/10.1063/1.1139918Google Scholar
S. Gourion-Arsiquaud et al., “Use of FTIR spectroscopic imaging to identify parameters associated with fragility fracture,” J. Bone Miner. Res. 24(9), 1565–1571 (2009).JBMREJ0884-0431http://dx.doi.org/10.1359/jbmr.090414Google Scholar
S. J. Gadaleta et al., “Fourier transform infrared microscopy of calcified turkey leg tendon,” Calcif. Tissue Int. 58(1), 17–23 (1996).CTINDZ0171-967Xhttp://dx.doi.org/10.1007/BF02509541Google Scholar
M. E. Szabo et al., “Similar damage initiation but different failure behavior in trabecular and cortical bone tissue,” J. Mech. Behav. Biomed. Mater. 4(8), 1787–1796 (2011).JMBBCP1751-6161http://dx.doi.org/10.1016/j.jmbbm.2011.05.036Google Scholar
J. S. Nyman et al., “Differential effects between the loss of MMP-2 and MMP-9 on structural and tissue-level properties of bone,” J. Bone Miner. Res. 26(6), 1252–1260 (2011).JBMREJ0884-0431http://dx.doi.org/10.1002/jbmr.326Google Scholar
J. R. Maher et al., “Raman spectroscopy detects deterioration in biomechanical properties of bone in a glucocorticoid-treated mouse model of rheumatoid arthritis,” J. Biomed. Opt. 16(8), 087012 (2011).JBOPFO1083-3668http://dx.doi.org/10.1117/1.3613933Google Scholar
B. Z. Gevorkian, N. E. Arnotskaia and E. N. Fedorova, “Study of bone tissue structure using polarized Raman spectra,” Biofizika 29(6), 1046–1052 (1984).BIOFAI0006-3029Google Scholar
M. Kazanci et al., “Bone osteonal tissues by Raman spectral mapping: orientation-composition,” J. Struct. Biol. 156(3), 489–496 (2006).JSBIEM1047-8477http://dx.doi.org/10.1016/j.jsb.2006.06.011Google Scholar
G. Falgayrac et al., “New method for Raman investigation of the orientation of collagen fibrils and crystallites in the Haversian system of bone,” Appl. Spectrosc. 64(7), 775–780 (2010).APSPA40003-7028http://dx.doi.org/10.1366/000370210791666255Google Scholar
G. Placzek, “Rayleigh-streuung und Raman-effekt,” in Handbuch der Radiologie, E. Marx, Ed., pp. 205–374, Akademische Verlagsgesellschaft, Leipzig, Germany (1934).Google Scholar
A. Bonifacio and V. Sergo, “Effects of sample orientation in Raman microspectroscopy of collagen fibers and their impact on the interpretation of the amide III band,” Vib. Spec. 53(2), 314–317 (2010).VISPEK0924-2031http://dx.doi.org/10.1016/j.vibspec.2010.04.004Google Scholar
A. Masic et al., “Observations of multiscale, stress-induced changes of collagen orientation in tendon by polarized Raman spectroscopy,” Biomacromolecules 12(11), 3989–3996 (2011).BOMAF61525-7797http://dx.doi.org/10.1021/bm201008bGoogle Scholar
M. Kozielski et al., “Determination of composition and structure of spongy bone tissue in human head of femur by Raman spectral mapping,” J. Mater. Sci. Mater. Med. 22(7), 1653–1661 (2011).JSMMEL0957-4530http://dx.doi.org/10.1007/s10856-011-4353-0Google Scholar
C. A. Lieber and A. Mahadevan-Jansen, “Automated method for subtraction of fluorescence from biological Raman spectra,” Appl. Spectrosc. 57(11), 1363–1367 (2003).APSPA40003-7028http://dx.doi.org/10.1366/000370203322554518Google Scholar
A. Robichaux-Viehoever et al., “Characterization of Raman spectra measured in vivo for the detection of cervical dysplasia,” Appl. Spectrosc. 61(9), 986–993 (2007).APSPA40003-7028http://dx.doi.org/10.1366/000370207781746053Google Scholar
E. Hecht and A. Zajac, Optics, Addison-Wesley Pub. Co., Reading, Massachusetts (1987).Google Scholar
E. Malus, “Mémoires de physique et de chimie de la Société d’Arcueil,” (1809).Google Scholar
M. G. Ascenzi and A. Lomovtsev, “Collagen orientation patterns in human secondary osteons, quantified in the radial direction by confocal microscopy,” J. Struct. Biol. 153(1), 14–30 (2006).JSBIEM1047-8477http://dx.doi.org/10.1016/j.jsb.2005.08.007Google Scholar
R. B. Martin and D. L. Boardman, “The effects of collagen fiber orientation, porosity, density, and mineralization on bovine cortical bone bending properties,” J. Biomech. 26(9), 1047–1054 (1993).JBMCB50021-9290http://dx.doi.org/10.1016/S0021-9290(05)80004-1Google Scholar
R. B. Martin and J. Ishida, “The relative effects of collagen fiber orientation, porosity, density, and mineralization on bone strength,” J. Biomech. 22(5), 419–426 (1989).JBMCB50021-9290http://dx.doi.org/10.1016/0021-9290(89)90202-9Google Scholar
R. B. Martin et al., “Collagen fiber organization is related to mechanical properties and remodeling in equine bone. A comparison of two methods,” J. Biomech. 29(12), 1515–1521 (1996).JBMCB50021-9290http://dx.doi.org/10.1016/S0021-9290(96)80002-9Google Scholar
X. Bi et al., “Raman and mechanical properties correlate at whole bone- and tissue-levels in a genetic mouse model,” J. Biomech. 44(2), 297–303 (2011).JBMCB50021-9290http://dx.doi.org/10.1016/j.jbiomech.2010.10.009Google Scholar