1.

Introduction

The loss of bone fracture resistance as occurs in osteoporosis and with aging, is not solely due to a decrease in bone mass but likely involves deleterious changes to tissue organization, including both the mineral phase and the organic matrix. Current x-ray-based clinical diagnostics predominantly assess bone strength through analysis of mineral density^1,2 and macrostructure or microarchitecture.³ However, changes to other material properties including toughness and fracture toughness may also contribute to fracture risk. In fact, Burstein et al.⁴ found that the age-related decrease in human bone toughness is greater than the age-related decrease in material strength at the apparent-level, suggesting that brittleness, not just weakness, lowers fracture resistance with age. To complement existing clinical diagnostics, many in the photonics field have begun to explore the use of Raman Spectroscopy (RS) due to its potential to concurrently and noninvasively quantify the biochemical signature of both mineral and collagen.⁵ However, RS is traditionally limited to biochemical content without directly assessing matrix organization. Recently, the link inherent between polarization of the Raman phenomenon and structural organization of crystals and birefringent molecules^6,7 has been extended to tissue,⁸ specifically the collagen matrix.^9,10 Polarization RS may overcome the traditional structural limitations and improve the ability of RS to explore the underlying mechanisms that influence fracture resistance.

The effect of genetic disease on bone is often studied using rodent models as a means to understand the mechanisms linking cellular action to tissue biomechanics as shown in several RS studies of bone quality.^11–13 Osteogenesis imperfecta (OI) is perhaps the best characterized human disease of brittle bone; OI is known to be a disease of collagen organization;^14–16 and the RS peaks of OI mouse models change with polarization.¹⁷ This is perhaps not too surprising, given a growing body of evidence that polarized Raman intensity varies with collagen direction^10,18 and can even be used to determine collagen orientation in bone.¹⁹ Other genetic mouse models have also demonstrated brittle bone phenotypes that are not necessarily explained by composition but are related to matrix organization,^20,21 and hence may benefit from polarization RS. We recently reported that there are polarization-related changes in the RS peak ratios of bones lacking the activating transcription factor 4 (ATF4), implying an organizational component to the phenotype that includes toughness loss and decreased fracture toughness,²² but it is unclear to what extent these polarization RS changes are predictive of mechanical bone quality. This is especially interesting given that mice lacking the matrix metalloproteinase (MMP9) exhibit no changes in standard RS peak ratios even though they, too, have low bone toughness.²³

In translucent, pure media, polarization-related changes in RS intensities can be quantified relative to theory using a depolarization ratio ($I_{\text{perpendicular}}/I_{\text{parallel}}$) to assess molecular orientation;⁶ however, in turbid media like bone, this is confounded by the light scattering properties of tissue. Therefore, the effect of specific optics instrumentation, including the degree of confocality,¹⁷ on RS peaks needed to be characterized prior to implementing assessment of bone organization with polarization RS. For different objectives ($\mathrm{NA}=0.75$ and $\mathrm{NA}=0.4$), we previously showed that observable changes in Raman intensity in polarization sensitive peaks with bone orientation occur across microstructural features of human tissue samples²⁴ such that the $\nu 1$ phosphate peak is regularly out of phase with amide I. However, traditional peak ratios only report a small portion of the known spectral changes that occur with polarization as previously characterized by others.^8,18,25–27 Moreover, limiting analysis to established peak ratios may overlook crucial aspects of complex toughness phenotypes if they manifest as subtle peak widening or shifting.

RS analyses of bone are often limited to a set of strong peak intensities and validated peak ratios,²⁸ including only the characteristic frequencies and relative composition of bone. In the absence of the complexity of tissue, RS is often applied to pure crystalline and polymeric materials (i) to analyze stress and strain using Raman shift changes;^29,30 (ii) to determine orientation as a function of peak polarization;³¹ and (iii) to assess crystal quality and deformation using peak width.³² Although each would require significant empirical support to extend mechanistically to analysis of turbid bone tissue, unsupervised multivariate analysis allows for the influence of these known RS sensitivities to be related to mechanical outcomes without spurious mechanistic suppositions. Indeed, recent applications of multivariate RS analyses to bone have helped explain the fracture toughness of bone in a model of rheumatoid arthritis³³ as well as partial mechanical improvement of OI mouse bones after human stem cell transplant.³⁴

Therefore, we investigated the potential of full spectrum multivariate analysis with a “bottom-up” design built upon the principal directions of Raman variance between test groups. Using both ATF4 and MMP9 mouse models in which peak ratios alone did not fully explain the mechanical phenotype, we hypothesized that analysis of all wavenumbers within the RS “fingerprint region,” while including polarization information, would improve the ability of RS to classify knockout mice of a brittle phenotype (versus wild-type of the same strain) and subsequently improve RS correlation to mechanical measures of toughness. However, this may not be the case for mechanical properties like strength that can be attributed to composition more than organization. To this end, we scanned excised intact femurs of both Atf4−/− and Mmp9−/− mice and their corresponding controls for two different orientations at a consistent anatomical location on the anterior midshaft. Using single (less sensitivity to polarization) or paired data (more sensitivity to polarization), peak ratios and principal components were tested for the ability to explain the toughness difference between genotypes.

2.

Methods

2.4.

Raman Spectroscopy

Raman spectra were acquired from the anterior midshaft of the intact extracted femurs of each animal femur prior to mechanical testing using a standard confocal Raman microscope (Renishaw InVia Raman Microscope, Renishaw, Hoffman Estates, Illinois) set to a $35\mu \mathrm{m}$ slit opening at $1\text{-}{\mathrm{cm}}^{-1}$ spectral resolution, and equipped with a $50\times$ ($\mathrm{NA}=0.75$) objective and a mode stabilized TEM (0, 0) 785-nm laser diode source (Innovative Photonic Solutions, Monmouth Junction, New Jersey), and a $1800\text{-}\text{lines}/\mathrm{mm}$ holographic grating. Laser power was measured daily at the sample before and after measurements to ensure consistent exposure to 35-mW laser power. 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. Laser polarization was then confirmed to have an approximate extinction ratio of $1:200$ (light was not further polarized by additional optics). The spot size was approximated at $1.5\mu \mathrm{m}$ and was Gaussian 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. Daily silicon measurements before and after data collection ensured consistency of wavenumber calibration and collection arm throughput.

The third distal midshaft vessel perforation was used as a landmark to consistently select the site of analysis across bones.³⁹ Bones were thawed to room temperature, mounted on a standard microscope slide using modeling clay, then attached to a rotation stage (accurate to 5 arcsec) to ensure that the plane of the anterior midshaft would be normal to the incident laser beam. Spectra were collected at two intact bone orientations such that the long bone axis was oriented either parallel (termed the longitudinal orientation) or perpendicular (termed the orthogonal orientation) to the polarization axis of the incident laser. Colocalization of collection sites before and after rotation was manually accomplished by registering fine structural features in the bright field. Both our work²⁴ and that of Kazanci et al.^9,40 have previously shown that the RS of two orthogonal orientations of human cortical bone provides information about tissue organization. However, this does not imply that the use of a true depolarization ratio as defined by Porto et al.⁶ is straightforward for bone. Presently, using the depolarization ratio on the full spectrum from intact mouse femurs inappropriately changed data shape, where the division by “near-zero” values inflated the variance of regions with little biochemical significance (Fig. 1). Moreover, there has been recent evidence that in structured materials, the use of the depolarization ratio is not always accurate for the quantification of fiber orientation.⁴¹

Fig. 1

Depolarization ratio spectrum overlaid on average Raman spectrum of mouse bone shows how the depolarization ratio skews data shape. After the mathematical operation, data variance is located in baseline regions, instead of at peaks where changes are occurring.

Each spectrum consisted of five accumulations of 10-s integration time to yield a high signal-to-noise ratio (SNR) from 300 to $1800{\mathrm{cm}}^{-1}$ (Hydroxyproline SNR in excess of 25:1). Spectra were processed via least squares modified polynomial fit⁴² and smoothed for noise using a second-order Savitsky–Golay filter.⁴³ Custom MATLAB software (Mathworks, Natick, Massachusetts) extracted the intensity and wavenumber of the prominent spectral peaks. There was no spectral binning in this analysis. After fluorescence subtraction, a linear baseline subtraction (based on derivative zero-crossings neighboring the peak) was conducted on peaks that overlap with neighboring constituents to ensure no residual fluorescence (namely proline, hydroxyproline, $\nu 1$ phosphate, and carbonate). This data was then used to generate markers of bone composition for mineralization ($\nu 1\text{phosphate}/\text{amide}\mathrm{I}$), carbonate substitution ($\text{carbonate}/\nu 1\text{phosphate}$) and crystallinity (crystal grain size and perfection, determined by the inverse full-width at half maximum intensity of the $\nu 1$ phosphate peak).

3.

Results

3.1.

Differences in Raman Spectra between Genotypes Become Apparent upon Bone Rotation (Polarization)

Despite the strong toughness reduction in Atf4−/− mice (Fig. 2), only slight changes in RS peaks [Figs. 3(a)–3(c)] are seen at any one orientation. To investigate whether differences exist between the genotypes in ways beyond traditional peaks and peak ratios, the effect of bone rotation on the Raman spectra of bone was examined [Figs. 3(d) and 3(e)]. Upon bone rotation with respect to laser polarization, both Atf4+/+ [Fig. 3(d)] and Atf4−/− [Fig. 3(e)] RS profiles changed at previously identified polarization-sensitive regions (Fig. 4) including the shoulder of $\nu 4$ phosphate ($590{\mathrm{cm}}^{-1}$), the $\nu 3$ phosphate peak ($1045{\mathrm{cm}}^{-1}$),⁸ the amide III $\delta (\mathrm{NH})$ band ($1273{\mathrm{cm}}^{-1}$),^26,44 and amide I band,⁹ demonstrating statistically significant differences between orientations ($p<0.05$ Mann Whitney U). The spectra of Atf4−/− bone had smaller variance than Atf4+/+ bone at amide III $\delta (\mathrm{NH})$ band and the amide I band, yielding greater separation between orientations than its wild-type counterpart.

Fig. 2

Toughness differences in Atf4 and Mmp9, unlike strength changes in maturation, are not associated with changes in tissue mineral density (TMD). Atf4 demonstrates isolation of toughness loss (b) without significant changes in strength (a) or TMD (c). Mmp9 males (d to f) demonstrate similar mechanical profile of isolated toughness loss (e). MMP9 study of both genders illustrates statistically significant but less pronounced difference in toughness (h). Maturation (j to l) shows an opposite profile where strength differences (j) associate with TMD (l) but not toughness (k).

Fig. 3

Polarization changes in the Raman spectrum are more intense in brittle bone, but not apparent in any one bone orientation. Peaks from the orthogonally orientated bone spectrum show slight trends that do not significantly separate for (a) $\nu 2$ and $\nu 4$ phosphate, (b) $\nu 1$ phosphate, or (c) amide I. (d) Raman spectra ($\text{mean}\pm \text{standard}\text{deviation}$) show more pronounced differences between bone orientations for Atf4−/− (e) than Atf4 +/+ mice (d).

Fig. 4

Detail regions of Raman spectroscopy (RS) spectra ($\text{mean}\pm \mathrm{SD}$) for two orthogonal bone orientations (orthogonal in color and longitudinal in black) depict orientation (polarization)-related changes in mineral and collagen bands for ATF4-deficient bone that are more severe than those seen in wild-type bone.

3.2.

Classification of Brittle Bone Genotypes Improves with Inclusion of Full Spectrum Polarization Information

To determine whether the principal components derived from polarization-influenced Raman spectra improves the ability of RS to classify brittle bone, we performed SMLR with leave-one-out cross validation on PCA output and traditional peak ratios (Table 1). Including peak ratios from both orientations showed no improvement of classification compared to using a peak ratio from either orientation (Tables 1 and 2). The impact of increasing polarization information was further examined using models including all principal components that significantly separated Atf4−/− from Atf4+/+ (Mann-Whitney U; $p<0.05$). SMLR using PCA output from only the orthogonal bone orientation [Fig. 5(a)] classified the genotype for 70.8% of the mice (Table 1), misclassifying only 2 Atf4−/− bones. However, using spectra from both orthogonal (O) and longitudinal (L) orientations in the PCA improves SMLR-based classification capabilities as evidenced by the increase in the samples correctly classified and the improvement in class posterior probability [Fig. 5(b)]. Thus, Atf4−/− and Atf4+/+ bones have different degrees of polarization changes in RS upon bone rotation, and including multiple orientations to capture this polarization-based RS phenotype improves genotype classification (Table 1). This finding implies that the tissue organization underlies the brittleness difference between Atf4+/+ and Atf4−/− mice.

Table 1

Multivariate expression of increasing polarization information improves classification over common peak ratios.

Genotype	Orientation	Sparse multinomial logistic regression (SMLR) input	Value ($\text{mean}\pm \mathrm{SD}$) WT	Value ($\text{mean}\pm \mathrm{SD}$) KO	Genotype separation $P_{\text{value}}$a	Correctly classified	Specificity (%)	Sensitivity (%)
Atf4	–	$\mu \mathrm{CT}$ tissue mineral density (TMD)	$1303\pm 16.6$	$1313\pm 19.2$	NS	15/24	100.0	0.0
	O	Mineral-to-collagen ratio (MCR) (ν1/AmI)	$11.450\pm 2.876$	$13.200\pm 1.549$	NS	15/24	100.0	0.0
	L	MCR (ν1/AmI)	$16.330\pm 3.147$	$19.410\pm 2.253$	$p=0.017$	15/24	100.0	0.0
	O	Carb sub	$0.157\pm 0.006$	$0.151\pm 0.007$	NS	15/24	100.0	0.0
	L	Carb sub	$0.155\pm 0.006$	$0.151\pm 0.005$	NS	15/24	100.0	0.0
	O	Crystallinity	$0.0649\pm 0.0008$	$0.0646\pm 0.0009$	NS	15/24	100.0	0.0
	L	Crystallinity	$0.0646\pm 0.0009$	$0.0645\pm 0.0008$	NS	15/24	100.0	0.0
	O	PC score	$8.576\pm 23.730$	$-14.290\pm 8.013$	$p=0.015$	17/24	66.7	77.8
	L	PC score	$8.435\pm 21.420$	$-14.060\pm 10.350$	$p=0.001$	16/24	53.3	88.9
	[O, L]	PC score	$12.030\pm 26.090$	$-20.050\pm 11.230$	$p=0.0007$	19/24	73.3	88.9
Mmp9	–	$\mu \mathrm{CT}$ TMD	$1458\pm 14.8$	$1474\pm 16.1$	$p=0.004$	7/15	0.0	100.0
	O	MCR (ν1/AmI)	$14.530\pm 1.189$	$12.430\pm 1.270$	$p=0.009$	11/15	71.4	75.0
	L	MCR (ν1/AmI)	$19.930\pm 2.543$	$18.010\pm 2.040$	NS	8/15	42.9	62.5
	O	Carb sub	$0.153\pm 0.005$	$0.154\pm 0.004$	NS	7/15	100.0	0.0
	L	Carb sub	$0.152\pm 0.006$	$0.152\pm 0.005$	NS	7/15	100.0	0.0
	O	Crystallinity	$0.0645\pm 0.0005$	$0.0653\pm 0.0007$	$p=0.029$	7/15	100.0	0.0
	L	Crystallinity	$0.0650\pm 0.0005$	$0.0654\pm 0.0006$	NS	7/15	100.0	0.0
	O	PC score	$-15.330\pm 12.090$	$13.420\pm 12.120$	$p=0.001$	12/15	71.4	87.5
	L	PC score	$-11.120\pm 21.850$	$9.730\pm 11.690$	NS	10/15	57.1	75.0
	[O, L]	PC score	$-19.750\pm 20.680$	$17.280\pm 14.080$	$p=0.0006$	13/15	85.7	87.5
Mmp9 gender	–	$\mu \mathrm{CT}$ TMD	$1310\pm 19.0$	$1235\pm 328$	NS	13/26	100.0	0.0
	O	MCR (ν1/AmI)	$11.820\pm 1.901$	$11.470\pm 1.380$	NS	13/26	100.0	0.0
	L	MCR (ν1/AmI)	$18.280\pm 1.962$	$17.950\pm 3.087$	NS	13/26	100.0	0.0
	O	Carb sub	$0.152\pm 0.006$	$0.156\pm 0.004$	$p=0.034$	13/26	100.0	0.0
	L	Carb sub	$0.150\pm 0.007$	$0.152\pm 0.004$	NS	13/26	100.0	0.0
	O	Crystallinity	$0.0653\pm 0.0011$	$0.0651\pm 0.0009$	NS	13/26	100.0	0.0
	L	Crystallinity	$0.0650\pm 0.0007$	$0.0653\pm 0.0010$	NS	13/26	100.0	0.0
	O	PC score	$9.638\pm 23.610$	$-9.638\pm 15.360$	$p=0.016$	19/26	76.9	69.2
	L	PC score	$-5.882\pm 16.030$	$10.720\pm 20.970$	NS	18/26	76.9	61.5
	[O, L]	PC score	$7.300\pm 23.440$	$-13.510\pm 17.720$	$p=0.022$	21/26	76.9	84.6

^aGenotype separation is determined by Mann Whitney U test, significance at p<0.05.

Table 2

Bivariate expression of increasing polarization information from orthogonal orientations does not improve classification.

Genotype	Orientation	SMLR input	Correctly classified	Specificity (%)	Sensitivity (%)
Atf4	OL	MCR (ν1/Am1)	15/24	100.0	0.0
	OL	Carb sub	15/24	100.0	0.0
	OL	Crystallinity	15/24	100.0	0.0
Mmp9	OL	MCR (ν1/Am1)	5/15	28.6	37.5
	OL	Carb sub	7/15	100.0	0.0
	OL	Crystallinity	7/15	100.0	0.0
Mmp9 gender	OL	MCR (ν1/Am1)	18/33	0.0	100.0
	OL	Carb sub	15/18	0.0	100.0
	OL	Crystallinity	18/33	0.0	100.0

Fig. 5

Classification accuracy and posterior probability improve when including RS polarization information. (a) Sparse multinomial logistic regression (SMLR) classification for wild-type and knockout bone from the Atf4 study misclassifies two knockout and five wild-type samples when using only the orthogonal orientation. (b) SMLR classification including both orientations misclassifies one knockout and four wild-type and posterior probability increases for many samples.

Improved classification with the inclusion of paired spectra (but not paired peak ratios) from two bone orientations was consistent across studies involving genetic knock-out mice with a brittle bone phenotype [Figs. 2(e) and 2(h)]. Unlike the ATF4 study, $\mu \mathrm{CT}$-derived TMD was significantly different between genotypes for the MMP9 study involving only males, but similar to the ATF4 study, TMD had poor classification accuracy (Table 1). With the exception of $\nu 1\text{phosphate}/\text{amide}\mathrm{I}$ for the Mmp9 male data set, none of the peak ratios from either bone orientation appropriately classified knockout bones as determined by SMLR (i.e., 0% sensitivity), such that all bones were classified as wild-type. Use of the full spectrum as PCA input for either the O or L bone orientation improved SMLR classification sensitivity for all three studies (Table 1) with the two different bone orientations yielding different classification results. Including both orientations as separate variables for PCA (denoted as [O, L]) provided further improvement in overall classification accuracy in both MMP9 studies, despite the confounding factor of gender.

3.3.

Principal Components have a Stronger Correlation with Toughness than Traditional Compositional RS Measurements

Expanding upon the multivariate classification of genotype, we analyzed whether RS polarization information improves mechanical bone quality characterization using Spearman’s correlation coefficients (Table 3). Peak ratios did not correlate with toughness except for longitudinal mineral-to-collagen ratio (MCR) and the paired MCR (Spearman’s correlation coefficient of SMLR scores; $r=-0.57$) in the ATF4 study, the only study where $\mu \mathrm{CT}$-derived TMD also correlated significantly (Table 3). PC1 from either orientation correlated with toughness in the ATF4 and MMP9-gender studies, again with differing results between orientations. The inclusion of both orientations in the PCA strengthened the correlation between PC1 and toughness in all three studies. Notably, none of the peak ratios or the genotype separating PC’s from toughness models yielded a significant correlation to bending strength. Figure 6 shows linear regressions between PC1 from the [O, L] analysis of both orientations and toughness for each genotype. As one of few PCs to significantly separate genotype (Mann–Whitney U; $p<0.05$), PC1 was consistently the only PC observed to significantly explain toughness variance (Fig. 6). As such, PC1 had a significant linear relationship to toughness for each study (F test: $\text{slope}\ne 0$; $p<0.05$), although the slope of this relationship was not positive for all studies. The additional MMP9 study including both male and female femurs exhibited a significant linear regression despite a marked overlap in toughness values between the genotypes. That is, there is a notably smaller toughness difference [Fig. 2(h)] in this study than in the MMP9 male only study, but PC1 still explains the toughness variance. In summary, full spectrum polarization information improves genotype classification and correlation to toughness in all three studies when compared to standard RS peak ratios or $\mu \mathrm{CT}$-derived TMD.

Table 3

Multivariate Raman spectroscopy signatures provide consistent, improved correlation to bone toughness.

Orientation	Property	Atf4a	Mmp9a	Mmp9 gendera
		Toughness	Toughness	Toughness
–	$\mu \mathrm{CT}$ TMD	$-0.42$	NS	NS
O	MCR (ν1/AmI)	NS	NS	NS
L	MCR (ν1/AmI)	$-0.50$	NS	NS
O	Carb sub	NS	NS	NS
L	Carb sub	NS	NS	NS
O	Crystallinity	NS	NS	NS
L	Crystallinity	NS	NS	NS
O	PC1 score	0.53	NS	0.53
L	PC1 score	0.66	NS	0.57
[O, L]	PC1 score	0.62	$-0.64$	0.60

^aSpearman’s correlation coefficient, significance given by using exact test that the correlation is significantly nonzero.

Fig. 6

Linear regressions of principal component scores significantly explain a portion of bone toughness. (a) PC1 score explains toughness and separates genotype in Atf4. (b) Despite small sample size, PC1 score explains over 20% of variance in toughness of Mmp9 male samples. (c) Despite overlap due to gender, PC1 score significantly explains over 35% of Mmp9 toughness.

3.5.

Principal Components Loadings from RS of Brittle Bone Phenotypes Suggest Conserved Elements That Are Not Seen in PCs Distinguishing Bone Maturity

Examining spectral loadings of the principal components may identify elements of RS signatures that explain material property differences. Figure 7 compares principal component loadings from the [O, L] PCA that used spectra from both bone orientations and that represents the greatest explanation of variance. In the ATF4 model of toughness loss, strong negative weights of PC1, which differed between genotypes, occurred at all the mineral peaks, including $\nu 2$ phosphate ($430{\mathrm{cm}}^{-1}$), $\nu 4$ phosphate ($590{\mathrm{cm}}^{-1}$), $\nu 1$ phosphate shoulders ($960{\mathrm{cm}}^{-1}$), and carbonate ($1074{\mathrm{cm}}^{-1}$) for both orientations. Strong positive weights occurred at collagen peaks including proline ($857{\mathrm{cm}}^{-1}$), hydroxyproline ($878{\mathrm{cm}}^{-1}$), amide I ($1668{\mathrm{cm}}^{-1}$), and amide III ($1248{\mathrm{cm}}^{-1}$). For several peaks, weights differ between O and L bone orientations, including amide I and ${\mathrm{CH}}_2$ ($1454{\mathrm{cm}}^{-1}$). PC1 of the MMP9 male study, responsible for best genotype classification and correlation to toughness, shows similar weights to those seen for the ATF4 study with strong negative weights on mineral peaks and strong positive weights on collagen peaks.

Fig. 7

Principal components that explain toughness demonstrate consistent loadings among studies. Gray overlay is used to accent regions of greatest similarity between models of brittle bone. (a) PC1 of Atf4 male mice demonstrates strong positive collagen and strong negative mineral loadings. (b) Similar positive collagen weights and negative mineral weights explain toughness in PC1 of Mmp9 male mice. (c) PC3 of maturation provides best separation of class, but demonstrates strong differences in loading when compared to components from brittleness models.

Although principal components explaining toughness display some conserved spectral features, the principal components that separate maturation illustrate stark contrast from RS signatures of toughness. In Fig. 7(c), PC3 from the maturation experiment did not significantly correlate to strength (note that none of the PC’s in the [O, L] input correlated); however, the 11% of the variance explained by this PC offers the best observed separation of class. For both orientations, carbonate was negatively weighted. Surprisingly, very little to no weight was placed on the amide I peak. Most of the peaks were weighted differently between the two orientations. In summary, mouse genetic models of reduced toughness (loss of ATF4 or MMP9) have a similar RS polarization profile in the principal components that explain toughness; however, separation of maturation-based changes has strong differences in orientation that do not correlate to strength and are loaded differently than RS profiles of toughness.

4.

Discussion

In the present study, a multivariate analysis of polarization information in RS improved genotype classification and correlation to bone toughness—using bones from two genetic mouse models—over peak ratios in traditional RS analyses. As this technique is further developed, it may help establish the interplay between biochemical composition and tissue organization as a factor in the resistance of bone to fracture.

These experiments confirm that traditional peak ratio analysis is not always sufficient to explain bone mechanical quality. When classifying genotypes for which there is a difference in bone brittleness, only the polarization-sensitive $\nu 1\text{phosphate}/\text{amide}\mathrm{I}$ partially identified knockout bones (Table 1) and weakly correlated with toughness (Table 3). Overall, univariate and even bivariate analysis of RS factors were insufficient to explain genetic-related differences in brittleness, characteristics that depend on defects in matrix organization in addition to abnormal biochemical composition, even when using peak ratios from two orthogonal orientations. Full spectrum analysis of principal components derived from a single-bone orientation improved classification when compared with single-peak ratios (Table 1), suggesting that additional RS information assists in the explanation of these complex mechanical phenotypes. Notable differences in RS intensity between orthogonal bone orientations [Fig. 3(e)] occur mostly at locations previously reported by the Penel group^8,26 to be sensitive to polarization: ν3 phosphate ($1045{\mathrm{cm}}^{-1}$), amide III $\delta (\mathrm{NH})$ band ($1273{\mathrm{cm}}^{-1}$),⁴⁴ and amide I ($1668{\mathrm{cm}}^{-1}$). The most striking differences occur at peak shoulders and would, therefore, be overlooked by many standard analyses. Moreover, these orientation effects on Raman spectra were more pronounced in brittle bones than wild-type counterparts [compare Fig. 3(d) to Fig. 3(e)], suggesting that polarization RS is sensitive to the orientation and subsequently the organizational changes underlying toughness loss that occurs with the deletion of ATF4.

Although full spectrum analysis improved classification in these models, data from a single orientation alone was insufficient to produce consistent classification and correlation to toughness. This manifests as disparate classification accuracy when comparing principal components of results from longitudinal and orthogonal bone orientations (Table 1). However, use of polarization RS from both bone orientations improves overall classification in all three sets of wild-type and corresponding mutant mice (Table 1) and improved class posterior probabilities (Fig. 5). In bones lacking MMP9, this additional organization information also improved correlation of RS to toughness as well (Table 3).

PCA identified the RS signatures of brittle phenotypes, and as an unsupervised “bottom up” approach to full spectrum analysis, it lends credence to the reliability of mechanical correlation of these physiochemical signatures. In each of the ATF4 and MMP9 studies, the first principal direction of variance explained only 25% to 30% of the total variance of the Raman spectrum, yet it significantly separated genotype and correlated to toughness, explaining over 25% of this mechanical variance (Fig. 6). Utilizing unsupervised multivariate analysis potentially runs the risk of over-fitting the data, but the sensitivity of the method was validated with two different genetic mouse models and by repeating the MMP9 study using a different colony of mice. Even when the toughness measurements overlapped between genotypes in the additional MMP9 study (Fig. 2), which included both genders [Fig. 6(c)], the first principal component still significantly explained toughness.

When considering other mechanical properties and potential RS correlation, tissue organization (provided by polarization) may not play as much of a role. Using growing wild-type mice as a model for strength increase with no change in toughness, peak bending stress and $\mu \mathrm{CT}$-derived TMD increased with age [Figs. 2(j) and 2(i)]. Despite this, carbonate substitution was the only peak ratio that correlated with strength, despite a lower classification accuracy than other compositional properties (Table 4). Considering the numerous endocrine and cellular changes that occur during maturation, it is not entirely unexpected that more bone changes exist in the Raman spectrum than those that explain strength. As such, an increase in the amount of Raman wavenumbers sampled improved classification accuracy, but did not increase correlation to strength (Table 4). Moreover, the PCA of the RS data from the individual orientations of bone seem to yield consistent classifications and correlations to strength, which is strikingly different from the orientation and polarization-specific outcomes observed in models of brittleness. Since strength and toughness can be diametrically opposed in many materials including bone in certain instances,^45,46 it is perhaps not surprising that organizational information from polarization in RS specifically aids the explanation of toughness differences between genotypes but not the strength differences between age groups. Notably, since ATF4 influences collagen synthesis⁴⁷ and since MMP-9 is a gelatinase that processes denatured collagen,²³ these genes likely influence matrix organization. On the other hand, secondary mineralization accompanies skeletal maturity⁴⁸ from 8 weeks to 20 weeks of age in mice, increasing strength but not necessarily toughness.

Although this paper focuses on the use of multivariate analysis to highlight the organizational information added by multiple polarizations, multivariate analyses of RS have also been successfully used to identify compositional changes and discriminate tissue age and disease. Fu et al.⁴⁹ demonstrated how the linear discriminants of PCA separate compositional changes in the rat mandible after ovariectomy. PCA has also been used to chart developmental changes in differentiating mesenchymal stem cells⁵⁰ as well as RS signatures of cartilaginous zones during fetal femur development.⁵¹ Finally, Raghavan et al.⁵² conducted an in-depth comparison of several new multivariate techniques to help discriminate the RS signatures of tissue age as they compare to nanoindentation mechanical properties on the same length scale. These articles attest to the versatility of appropriately applied multivariate analyses of RS for skeletal tissue; however, none of these studies attempted to connect RS discrimination to the fracture resistance of bone.

Although principal components only represent one straightforward and simple application of multivariate analysis to RS, the use of a “bottom up” study design that builds upon inherent, uncorrelated dimensions of variance allows for the distinct advantage of analyzing underlying biochemical signatures for consistency between models. Different component loadings (data not shown) between the orthogonal and longitudinal orientations of the same bones within each study imply the necessity of different orientations, and subsequently organizational information, to explain toughness. There are differences in significant principal components between the two models of bone brittleness and the model of skeletal maturation (Fig. 7), noting especially heavy weightings for amide I bands ($\sim 1600$ to $1720{\mathrm{cm}}^{-1}$) in ATF4 [Fig. 7(a)] and MMP9 studies [Fig. 7(b)] but not for maturation [Fig. 7(c)]. Strong consistency between the first principal components (PC1) of ATF4 and MMP9 models includes heavy positive weights for amide III (1235 to $1280{\mathrm{cm}}^{-1}$) and amide I bands and heavy negative weights for $\nu 2$ phosphate ($430{\mathrm{cm}}^{-1}$) and $\nu 4$ phosphate ($590{\mathrm{cm}}^{-1}$). Although the mechanistic significance of these trends remains to be tested, the data imply some conserved elements in the Raman signature of toughness loss [Figs. 7(a) and 7(b)]. Note that there is no evidence that ATF4 is upstream or downstream of MMP9 activity to date. By correlating RS to bone material properties, we are effectively assuming an inherent relationship between biochemistry and biomechanics. Although the additional RS polarization information allows for a more direct link between toughness and tissue organization, the relationship between principal components and fracture resistance may not always be clear (Fig. 6). Although PC1 loading is similar for the two bone orientations between ATF4 and MMP9 male studies (Fig. 7), the spectral differences between respective genotypes are not necessarily the same across studies (MCR is higher for Atf4−/− than for Atf4+/+ mice while the opposite prevails for Mmp9 mice; Table 1).

Although polarization effects of the Raman phenomenon have been characterized since Placzek⁵³ and the application of vibrational spectroscopy to the characterization of bone has been conducted for decades (see Boskey^54,55 and Morris^5,56 for review), relatively little has been investigated at the union of these two fields until recently. Pioneers in the bone field showed that the crystalline structure was highly organized in enamel and bone;⁸ that the orientation of collagen fibers in bone could be extracted from polarization RS¹⁹ and aligned with theory;⁵⁷ and that the effect of the polarization phenomenon on bone RS could be augmented or minimized by appropriate instrumentation choices.¹⁷ Traditionally, RS probes tissue biochemistry with a high degree of molecular specificity, but polarization analysis can be sensitive to differences in tissue organization, such that RS can then be used to analyze the interaction between bone composition and organization. The additional layer of full spectrum multivariate analysis extends RS interpretation beyond the relative composition implied by peak ratios.

Since the goal of the study was to compare the performance of standard peak ratio analysis to multivariate expressions of RS polarization in their ability to explain bone mechanical quality, Raman spectra were acquired from anatomically consistent surfaces without preparation in the region of mechanical testing, and as such, this anterior midshaft site might not yield the strongest polarization difference seen in bone. The consistent, observed sensitivity of this method to bone mechanical quality is surprising considering vast under-sampling with less than $15\mu {\mathrm{m}}^3$ sampled per bone. To make an accurate comparison between standard measurements of peak ratios and full spectrum analysis that includes polarization information, the instrument was left at its inherent polarization extinction ratio. Therefore, further optimization of the Raman instrument may offer continued improvement in the explanation of bone mechanical quality. However, the current configuration makes the method readily available for any lab with a confocal RS instrument and a rotation stage.

In the present study, the subtle spectral changes that occur with polarization sensitivity in RS to matrix organization coincide with the loss of mechanical toughness and not with changes in material strength. Thus, polarization RS shows promise as a novel tool to explore the dynamic and subtle underpinnings of the mechanisms behind bone mechanical quality. However, there is indubitably more to the RS signature of bone than a complex explanation of mechanics. In addition, there are likely contributors to mechanical quality to which RS is not sensitive. Near infrared RS, used for clinical relevance, is tuned to be largely insensitive to water, and bound and pore water have a significant contribution to the fracture resistance of bone.⁵⁸ RS will not have the resolution of atomic force microscopy to map collagen d spacing,⁵⁹ nor the SNR of two photon fluorescence to examine collagen fiber orientation.⁶⁰ Nonetheless, the inherent interplay between chemical composition and tissue organization in RS polarization may prove useful in explaining changes in the fracture resistance of complex human microstructures. Ascenzi showed such complexity as mechanical properties of osteons for different loading modes (i.e., compression versus tension) were related to the primary collagen orientation relative to the directionality of the haversian canals.^61–63 Both polarized light⁶⁴ and polarization RS⁴⁰ are sensitive to osteonal lamellae. As such, if the organizational information in polarization RS continues to explain the mechanical integrity of bone as it is applied to human bone, this technique could be used to explain how microstructural heterogeneity and composition affect bone mechanical quality, ultimately producing an RS profile for healthy bone tissue.

5.

Conclusions

Multivariate analysis of Raman spectra at two bone orientations (enhancing polarization sensitivity) assisted in the explanation of a toughness loss in genetic mouse models involving two different genes (for a transcription factor and for an enzyme). Across three sets of wild-type and knockout bones, MCR, which is often used to characterize compositional differences, did not explain or poorly explained the difference in bone toughness between genotypes. This was effectively achieved with Raman spectral analysis (300 to $1800{\mathrm{cm}}^{-1}$) using principal components acquired from two orthogonal bone orientations such that the first and largest direction of variance consistently separated the brittleness phenotypes and significantly correlated with bone toughness. Using the same technique on bones from a mouse model of skeletal maturation did not improve the explanation in the age-related increase in strength when compared to a univariate approach. This implies that the multivariate analysis of Raman spectra is not simply providing more information, but rather the organizational information provided by polarization may specifically aid in explaining variance in toughness but not in strength. If the observed Raman profiles scale to the complex organization of human bone, polarization in RS may have clinical utility for bone quality assessment.