Research Papers: Imaging

Optical spectral imaging of degeneration of articular cartilage

[+] Author Affiliations
Jussi Kinnunen

University of Eastern Finland, Department of Physics and Mathematics, P.O.B. 111, FI-80101, Joensuu, Finland

Jukka S. Jurvelin, Jaana Mäkitalo

University of Eastern Finland, Department of Physics and Mathematics, P.O.B. 1627, FI-70211, Kuopio, Finland

Markku Hauta-Kasari

University of Eastern Finland, School of Computing, P.O.B. 111, FI-80101, Joensuu, Finland

Pasi Vahimaa

University of Eastern Finland, Department of Physics and Mathematics, P.O.B. 111, FI-80101, Joensuu, Finland

Simo Saarakkala

University of Eastern Finland, Department of Physics and Mathematics, P.O.B. 1627, FI-70211, Kuopio, Finland

J. Biomed. Opt. 15(4), 046024 (August 13, 2010). doi:10.1117/1.3477190
History: Received November 03, 2009; Revised June 24, 2010; Accepted June 24, 2010; Published August 13, 2010; Online August 13, 2010
Text Size: A A A

Open Access Open Access

* Address all correspondence to: Jussi Kinnunen, University of Eastern Finland, Department of Physics and Mathematics, P.O.B. 111, FI-80101, Joensuu, Finland. Tel:358-13-25-17-919. E-mail: jussi.kinnunen@uef.fi

Osteoarthritis (OA) is a common musculoskeletal disorder often diagnosed during arthroscopy. In OA, visual color changes of the articular cartilage surface are typically observed. We demonstrate in vitro the potential of visible light spectral imaging (420to720nm) to quantificate these color changes. Intact bovine articular cartilage samples (n=26) are degraded both enzymatically using the collagenase and mechanically using the emery paper (P60 grit, 269μm particle size). Spectral images are analyzed by using standard CIELAB color coordinates and the principal component analysis (PCA). After collagenase digestion, changes in the CIELAB coordinates and projection of the spectra to PCA eigenvector are statistically significant (p<0.05). After mechanical degradation, the grinding tracks could not be visualized in the RGB presentation, i.e., in the visual appearance of the sample to the naked eye under the D65 illumination. However, after projecting to the chosen eigenvector, the grinding tracks are revealed. The tracks are also seen by using only one wavelength, i.e., 469nm, however, the contrast in the projection image is 1.6 to 2.5 times higher. Our results support the idea that the spectral imaging can be used for evaluation of the integrity of the cartilage surface.

Figures in this Article

Osteoarthritis (OA), the most common joint disease, is characterized by the degenerative changes in articular cartilage. Depletion of superficial proteoglycans,15 changes in collagen fibril orientation,34,67 and reduction in the collagen content4 are early signs of OA. Later, they lead to an increase of articular surface roughness8 and debilitate the mechanical integrity of articular cartilage by decreasing its stiffness.910 Frequently, as a complement to x-ray imaging, OA can be clinically diagnosed by arthroscopic examination.

In OA, visual color changes of the articular cartilage surface can be observed during arthroscopy, i.e., the loss of glisten and the change from white color to more yellowish. Furthermore, when evaluating the healing of the cartilage repair tissue after surgical repair procedures, similar color changes are observed.1112 These changes are visually and subjectively assessed through an arthroscope, and typically they can be clearly observed only when present in more advanced stages of OA. Consequently, there is a need for more quantitative and accurate methods for diagnostic purposes.

Optical spectral assessment was applied earlier for bovine hip joint,13 rabbit cartilage plugs after osteochondral autografting,11 rabbit cartilage after microfracture repair,12 and human articular cartilage.14 The studies have found color coordinate changes in the CIELAB colorimetric system when the cartilage integrity changes.1112,14 A model for the reflectance spectrum as a function of the cartilage thickness has also been proposed.13 Further optical properties, such as scattering and absorption coefficients and anisotropy, within the whole visible range of the spectrum have been reported for cartilage tissue.1516 Based on those studies, cartilage has been found to be highly forward-scattering tissue. However, in these earlier studies only a small area on the cartilage surface was analyzed. Thus, spatial spectroscopic assessment of the cartilage surface has not been conducted. Further, all of the earlier studies have been limited to very simple spectral analysis methods. With more sophisticated analysis methods, it may be possible to extract more detailed information on the integrity of the cartilage.

The main aim of this in vitro study was to investigate the effect of superficial degeneration on the optical spectrum of bovine articular cartilage. Controlled tissue degeneration was artificially created with enzymatic digestion using collagenase enzyme and mechanical grinding with an emery paper (P60). In addition to the color coordinate changes, the changes in the reflectance spectrum of cartilage before and after experimental procedures were analyzed with sophisticated spectral analysis methods, such as principal component analysis (PCA). Further, as a second aim we clarified the differences of the spectra between the patellar and femoral bovine cartilage.

Articular Cartilage Samples and Experimental Degradations

Bovine knee joints were obtained from the local abattoir (Atria Oyj, Kuopio, Finland). Joints were opened and cylindrical osteochondral samples (n=32, diameter=19mm) were prepared from the lateral upper quadrant of patella (PAT, n=26) and from the femoral medial condyle (FMC, n=6). All samples were visually intact except one: the patella with a spontaneous cartilage degeneration, probably associated with an earlier trauma, was also included in the study as a purpose of comparison with the artificial degradations. After preparation, the samples were kept in phosphate-buffered saline (PBS) until degradative treatments.

In enzymatic degradation, collagenase enzyme (111Uml, Collagenase Type 1a, C9891, Sigma Aldrich, St. Louis, Missouri) was used with different incubation times to obtain variable stages of degeneration. From PAT samples 8 were incubated for 6h, and 11 PAT samples and 6 FMC samples were incubated for 15h. Furthermore, 7 samples from the patella were exposed to mechanical grinding of the surface with an emery paper (P60 grit, 269μm particle size, KWH Mirka Ltd, Finland). The grinding was performed manually by sliding the cartilage surface two times against the emery paper. The slides were approximately at a 90-deg angle. The spectral images were collected for all samples before and after experimental treatments.

Acquisition of Spectral Data

The spectral images were captured by using the Nuance liquid crystal tunable filter (LCTF) spectral camera (model N-MSI-420-10, Cambridge Research & Instrumentation, Woburn, Massachusetts) (Fig. 1). The used spectral range of the camera was 420to720nm with 7-nm sampling. The measurements were conducted in standard 450 measurement geometry (45 illumination angle, normal detection angle). To simulate the arthroscopic evaluation, the samples were immersed in PBS during measurements. Illumination was by a 50-W halogen light (Accentline 50W GU5.3 12V 36D 1CT, Koninklijke Philips Electronics N.V., Amsterdam, Netherlands) driven by a mixed-mode-regulated precision laboratory power supply (EX1810R, Thurlby Thandar Instruments Ltd., Huntingdon, United Kingdom) and detection of the sample was made through a glass window of the custom-made sample container. Spectral data with the dynamic of 12bits was acquired using optimal integration times for each measured channel with the camera software (Nuance 1.6.8.2), and the data were further analyzed with a custom-made MATLAB program (MathWorks, Inc., Natick, Massachusetts). In the imaging system, the visible area was 14mm in diameter corresponding to 850pixels in the image, which resulted in a pixel resolution of the data to be 60.7pixelsmm. The spatial resolution of 14.25lpmm was assessed by using a USAF-1951 (U.S. Air Force 1951) resolution target (RES-2, Newport Corporation, Fountain Valley, California).

Graphic Jump LocationF1 :

Schematic presentation of the measurement setup as viewed on top. Measurements were conducted trough a glass window with an LCTF camera in 450 geometry while the sample was immersed in PBS.

Spectral Analysis

The acquired reflectance spectral images were analyzed using PCA. In addition, the spectral data were transformed into CIELAB uniform color space.17 The reflectance R(λ) of the sample at each pixel of the spectral image was calculated usingDisplay Formula

1R(λ)=Isample(λ)Idark(λ)Ireference(λ)Idark(λ),
where Isample is spectral power distribution of the illumination reflected from the sample, Ireference is spectral power distribution of the illumination reflected from the perfectly reflecting diffuser, and Idark is the detector dark signal. The Ireference was measured using an optical polytetrafluoroethylene (PTFE) reference material (ODM98, Gigahertz-Optik GmbH, Germany) with uniform 98% diffuse reflectance over the visible range. The reference material with a size equal to that of the tissue samples was placed in the system, and the Ireference was measured using individually optimized integration times for each wavelength channel. The reference material was replaced with the cartilage sample, and the Isample was measured by using the same integration times. The Idark, with corresponding integration times, was measured to reduce the effect of any stray light and camera bias.

The PCA for a set of spectra R from the x×y×n spectral imageDisplay Formula

2R=[R11(λ1)Rxy(λ1)R11(λn)Rxy(λn)],
where n is the number of the wavelengths, and x×y is the number of pixels, can be calculated as follows. The n×n correlation matrix C of the data is defined asDisplay Formula
3C=1x×yRRT.

Eigenvectors Φ of C are solutions of the equationDisplay Formula

4CΦ=σΦ,
where σ is the diagonal matrix of C’s eigenvalues. The eigenvectors form an orthonormal basis in the function space of measured spectra. These vectors are organized so that the one corresponding to highest eigenvalue is the first and so forth. The original image can be reconstructed accurately by a linear combination of those vectors.18 In this study, the characteristics of the data were studied by calculating projections of spectra to the eigenvectors of correlation matrix.

The CIELAB color coordinates (L*a*b*) can be calculated from the spectra by using the Commission internationale de l’eclairage (CIE) standard observer and any chosen illumination.17 The D65 standard daylight and CIE1931 standard observer were used in this study. First, the projections of the measured spectrum to the XYZ color-matching functions were calculated as follows.Display Formula

5X=kλS(λ)R(λ)x¯(λ)Δλ,
Display Formula
6Y=kλS(λ)R(λ)y¯(λ)Δλ,
Display Formula
7Z=kλS(λ)R(λ)z¯(λ)Δλ,
where S(λ) is the spectral power distribution of the chosen illumination; and x¯(λ), y¯(λ), and z¯(λ) are the color-matching functions of the standard observer. The Y has a value of 100 for the perfectly reflecting diffuser when the scaling factor k is defined asDisplay Formula
8k=100λS(λ)y¯(λ).

The nonlinear transformation from the XYZ coordinate system to the CIELAB system is defined using the following equations:17Display Formula

9L*={116(YYn)1316ifYYn>0.008856903.3(YYn)ifYYn0.008856},
Display Formula
10a*=500[(XXn)13(YYn)13],
Display Formula
11b*=200[(YYn)13(ZZn)13],
where Xn, Yn, and Zn denote the tristimulus values of the illumination. Equations 1011 are valid when X/Xn, Y/Yn, and Z/Zn are all greater than 0.008856. Otherwise the cube root in the equation in question [Eq. 10 or Eq. 11] is replaced byDisplay Formula
12(aan)13[7.787(aan)+16116],
where a is either X, Y, or Z.17

The L* denotes the brightness of the color. Coordinate a* denotes the blueness (a*>0) or the yellowness (a*<0), and the b* coordinate denotes the redness (b*>0) or the greenness (b*<0) of the color. In the CIELAB system, the color difference ΔEab* is almost uniform for all color pairs and is defined as the Euclidean distance between the coordinates:17Display Formula

13ΔEab*=[(ΔL*)2+(Δa*)2+(Δb*)2]12.

Contrast Evaluation

The contrast C between the grinding tracks and the intact areas for the mechanically degenerated cartilage was evaluated by using the Michelson contrast equation:Display Formula

14C=max(I)min(I)max(I)+min(I),
where I is the value of a pixel in a chosen pair of pixels. Totally 60 pairs of pixels, one pixel with clearly visible grinding tracks and another with intact cartilage nearby, were selected from each sample image. The pixels were selected as sparsely as possible to reduce the global variation, and a 5-×5-pixel window within the selected pixels was used for the contrast evaluation to reduce the local variations. The derived 60 values were averaged to represent the contrast of the grinding tracks for the specific sample.

Optical Microscopy

After the optical reflectance measurements, one representative sample from each sample group was processed for histological evaluation with an optical microscope (Axio Imager M2, Zeiss, Germany). Histological 3-μm-thick sections were prepared and stained with 0.5% safranin-O as described earlier.19 The sections were imaged and the integrity of the articular surface and proteoglycan content of the samples were visually evaluated. For each sample, cartilage thickness was also determined from the captured images, calibrated by using a 1-mm ruler.

Statistical Analysis

The nonparametric Wilcoxon signed rank test was applied to evaluate the statistical significance of the changes in spectral parameters (CIELAB coordinates and the projection of the spectra to eigenvector) after experimental degradations were compared to initial situation before degradations. The Kruskall-Wallis post hoc test was used to reveal statistical differences in the spectral slope between the sample groups. The level of significance was set to p<0.05. Statistical analyses were conducted using SPSS software (version 16.0, SPSS Inc., Chicago, Illinois).

Mean Reflectance Spectra

The mean reflectance spectrum at the center of each sample (200×200pixels3.3×3.3mm) was calculated before and after the degradation [Fig. 2]. After collagenase digestion, the changes in the reflectance spectra were evident in the wavelengths up to 600nm. However, in the wavelength range of 600to720nm, the changes were only minor after collagenase digestion. In contrast, after mechanical degradation, no changes could be observed in the mean reflectance spectra.

Graphic Jump LocationF2 :

(a) Mean spectra and (b) the relative change at each wavelength before and after 6h or 15h collagenase digestion or mechanical grinding using the P60 emery paper for all the patellar (PAT) and femoral (FMC) samples. The dashed lines in (b) are the linear fits to the ratios of the two average spectra, before and after the treatments.

A relative change after the experimental degradations is presented at each wavelength in Fig. 2. The mean slope (and the standard deviation) of the linear fit [Fig. 2] in each sample group is presented in Table 1. The relative change is nearly linear after the 6h collagenase digestion. After the 15-h digestion, some gaps in spectrum appear at wavelengths below 450 and around 550nm, i.e., at the same spectral bands to maximal absorption of the hemoglobin.13 The slopes after 15-h collagenase degradation are significantly different from those after mechanical degradation. After the 6-h collagenase degradation the slopes are significantly different only from the femoral samples (Table 1).

Table Grahic Jump Location
The slope of the linear fit for the relative change between the mean wavelength dependent reflectances, before and after collagenase digestion (6 or 15h) or mechanical degradation for the patellar (PAT) and fermoral (FMC) samples.
Table Footer NoteThe statistical significance (NS=nonsignificant) between the slopes is calculated by using the Kruskall-Wallis post hoc test.
CIELAB Color Coordinates

The CIELAB color coordinates were calculated for the samples at the center area (200×200pixels). After enzymatic digestions, statistically significant (p<0.05) changes, depending on the immersion time and anatomical location, were observed in CIELAB coordinates (Table 2). Similarly as for the reflectance spectra, no significant changes in coordinates could be observed after the mechanical degradation.

067004jbot2Table Grahic Jump Location
Value of the CIELAB coordinates (mean±SD) before and after mechanical degradation (mech.) and either 6- or 15-h enzymatic digestion using collagenase (coll.) for the patellar (PAT) and femoral medial condyle (FMC) samples. The RGB presentation of the colors and the ΔEab* color difference before and after treatments are also shown.
Table Footer NoteaStatistically significant difference (Wilcoxon signed rank test, p<0.05) compared to sample before degradation.

The mean color difference ΔEab* showed that the color change after 15-h collagenase digestion was two times higher for the FMC than for the patella (Table 2).

PCA

The PCA was calculated for several subsets of the optical spectral image data. The sets for PCA, classified into groups by experimental treatments, were created by aggregating the spectra at the center area (200×200pixels) from all the samples with the same treatment. The spectra before and after degradation were taken separately and the PCA was calculated. The projections of spectral images to the second vector of the femoral samples after 15-h collagenase digestion were found to give good separation between the collagenase treatments and for the grinding marks caused by the emery paper (Fig. 3).

Graphic Jump LocationF3 :

The second eigenvector after 15-h collagenase digestion of femoral samples. The projections of the cartilage spectra to the presented eigenvector were found to separate the enzymatic changes and emery paper scratches from the untreated samples. The vector was chosen to be the base vector for all the projections.

Figure 4 shows examples of the projection images. The projection shows uneven distribution of the spectral differences. The mean level of the center area (Fig. 5) in the samples treated by the emery paper showed no statistically significant difference between the samples before and after the degradation. However, for enzymatic degradations the level increases when the immersion time is extended. The increase is significant (p<0.05) after all enzyme degradations.

Graphic Jump LocationF4 :

Examples of pseudocolor projections of the spectral images to the chosen base vector the second eigenvector of the femoral samples after 15-h collagenase digestion, Fig. 3. Representative patella (PAT) and femoral (FMC) samples before and after 6- or 15-h collagenase digestion or mechanical grinding using the P60 emery paper are shown. The PAT 15-h sample showed spontaneous trauma before collagenase digestion and is presented here for the purpose of comparison. All samples show clear spatial differences including the mechanical degradation marks (“scratches”) and the spontaneous trauma in the projection image before and after the collagenase digestion.

Graphic Jump LocationF5 :

Mean values (200-×200-pixel center area) of the projections of the spectral images to the chosen base vector (the second eigenvector of the femoral samples after 15-h collagenase digestion) for samples before and after enzymatic (6-15-h collagenase) or mechanical degradations (P60 emery paper). The hollow circles indicate mean values within each sample, while the closed dots with lines represent the mean of samples within each group. The results from mechanical degradation are shown on the left side of the result of the native samples.

In the characteristic images of the mechanically degraded sample (Fig. 6) the marks are faint in the visual appearance of the surface to the naked eye under the D65 illumination. When using only one wavelength channel, the grinding is invisible at 700nm. However, specifically at 469nm the grinding can be seen. Use of one channel represents the case with the object viewed or illuminated through a narrow band filter. The projection of the sample spectra to the second eigenvector of FMC after 15-h enzyme treatment (Fig. 3) makes the grinding marks more visible in the cartilage and it was chosen to be the base vector for all the presented projections. The contrast calculation was performed for the 469-nm monochromatic and the projection image of the chosen base vector with the corresponding pixels. The contrast in projection image was found to be 1.6 to 2.5 times higher for each seven samples (Fig. 6). This kind of projection to the eigenvector indicates a combination of two images that are captured by using two separate complex optical filters.

Graphic Jump LocationF6 :

Examples of the visual appearance of one representative mechanically degraded cartilage sample in RGB, under monochromatic 700- and 469-nm light and after projection of spectra to the chosen base vector (the second eigenvector of the femoral samples after 15-h collagenase digestion). The mean contrast between the ground and intact areas in projection image is 1.6 to 2.5 timer higher for all the mechanically degraded samples when compared to image at 469nm.

Optical Microscopy

One representative sample in each group was studied under a optical microscope. Further, cartilage thickness was evaluated from the microscopic images (Fig. 7). Mechanical degradation only affected the cartilage surface but the proteoglycan content, as assessed by safranin-O staining, was uniform, compared to that of collagenase digested samples. After the 6-h collagenase digestion, the proteoglycans were depleted from the superficial cartilage. The thickness of the mechanically degraded and 6-h collagenase-digested samples was typical of normal intact bovine patellar cartilage (1901±224μm) (see, e.g., Rieppo et al. 20). In contrast, cartilage thickness was significantly reduced after the 15-h digestion with collagenase. Thickness of the one representative FMC sample was smaller (0.83mm) than those of the patellar samples.

Graphic Jump LocationF7 :

Optical microscopy images (safranin-O-stained sections) of the representative patellar samples after mechanical grinding with the P60 emery paper and after 6- and 15-h collagenase digestion. Cartilage thickness was determined from the captured images, calibrated by using the 1-mm ruler seen on the left.

In this study, the potential of visible light spectral imaging to quantify color changes in articular cartilage after experimental degradations, simulating OA, was investigated. Digestion by collagenase significantly changed the wavelength-dependent reflectance from the surface of bovine articular cartilage. The changes were more prominent at the shorter wavelengths, and statistically significant changes were observed for all the CIELAB coordinates.

The increase of the a* coordinate, indicating an increase in the tissue “redness,” was probably contributed by the red color of the subchondral bone.11,13 The characteristic absorption of the hemoglobin around 550nm and below 450nm (Ref. 13) was evident for samples after 15-h collagenase digestion (Fig. 2). A decrease in the lightness of the cartilage was related to the decrease in L* coordinate. The changes in L* and a* may be related to cartilage thickness. In the bovine knee joint, cartilage thickness at FMC is typically smaller than that in the patella.21 Interestingly, the differences in L* and a* between the FMC and PAT cartilage are consistent with the changes of coordinates after 15-h collagenase digestion, characterized by reduction in cartilage thickness. The dependence of L* and a* on cartilage thickness is also supported by the finding that the coordinate values of FMC cartilage before digestion and PAT cartilage after digestion were nearly similar, matching with the equal cartilage thickness at those states.

Because of the degenerative nature of our treatments, assessment of the sample thickness before and after treatments was not feasible. As a limitation, we could determine the cartilage thickness for only one representative sample in each sample group. However, a statistically significant loss of bovine patella articular cartilage thickness after collagenase digestion, consistent with our findings, is reported by Wang et al. 22

After 15-h collagenase digestion, the increase of the b* coordinate from negative to positive values imply the color to transform from bluish to more yellow. Inconsistently with L* and a*, this change may not be explained by the variations in tissue thickness, but the changes in the b* coordinate can be related to the cartilage structure, especially to those changes in the superficial zone of cartilage that could be revealed by the microscopy of the safranin-O stained tissue sections (Fig. 7). As cartilage is a complex tissue with multiple major components, such as collagen and proteoglycans, however, further studies should be conducted to reveal how the structural changes affect the spectrum.

In contrast, no significant changes in the mean optical reflectance spectrum were seen after mechanical degradation with the emery paper. However, by projecting the reflectance spectral image into carefully selected base vector (the second eigenvector of the femoral medial condyle data after the 15-h collagenase digestion) the visual appearance of the grinding marks after mechanical degradation was improved 1.6 to 2.5 times compared to the single channel image.

In this study, the same (native) samples before the treatments provided the control results for comparison with those after collagenase digestion and mechanical degradation. No separate control samples soaked for 6 or 15h in PBS were included. This was accepted as in our previous study,20 similar bovine patellar samples (n=12) were soaked in PBS for 44h with no significant differences in histological structure or composition, as judged by optical densitometry or polarized light microscopy. We strongly believe that the observed changes result from the collagen digestion or mechanical degradation.

Based on the studies by Ebert15 and Qu et al. 16 the value of the scattering coefficient of intact cartilage tissue is of an order of magnitude larger than the absorption coefficient. Further, the tissue is highly forward scattering. The effect of collagenase digestion on the scattering in articular cartilage, to our knowledge, has not been studied. However, the discrete particle model by Schmitt and Kumar23 indicates that the scatter-related attenuation depends on the ratio between the wavelength and the particle size in the tissue. When the collagen network is degraded by the collagenase and the proteoglycan content may also change (Fig. 7), the spectral differences (Fig. 2) may be attributed to a high concentration of collagen fragments that act as scattering centers. Thus, it seems that the spectral gaps below 450 and around 550nm are more specific to advanced cartilage degeneration, i.e., when the thickness of the tissue is significantly reduced. Consequently, for diagnosis of early changes, more advanced spectral analysis methods may be necessary.

Significant correlation between the Mankin score and the Outerbridge classification with the L* and a* values of the CIELAB coordinate system was reported for human cartilage by Ishimoto et al. 14 They found a decrease in L* coordinate and an increase in a* along the progression of cartilage degeneration. This is consistent with our findings on the effect of collagenase digestion. The increasing trend in the reflectance to the longer wavelengths in human cartilage samples differs from that found for the bovine cartilage in this study. The bovine reflectance from the intact sample showed nearly the same value at the both ends of the measured range. The differences between the reflectance spectra from human and bovine cartilage may arise from the differences in tissue thickness, as the human cartilage is known to be thicker than bovine cartilage.

While studying cartilage after microfracture repair, Hattori et al. have “hypothesized that the a* value represents the color of the fibrin clot, while the L* value indicates the newly formed cartilage”12 in the repair tissue. Based on the presented spectral reflectance and color coordinate changes, this kind of generalization, however, should be done with caution.

All the earlier spectroscopic assessments of the articular cartilage1114 have been performed by measuring only a small surface area at the time. Because the degenerative changes in cartilage are highly dependent on the location, the actual state of the cartilage might not be reliably recorded with measurements at a small area. In this study, for the first time, the optical imaging approach was utilized and tested. We demonstrated that mechanical surface degradation with emery paper could only be revealed in imaging mode after sophisticated PCA algorithms. This algorithm could be used both for evaluating the collagenase-induced changes in cartilage and the grinding marks on the surface of the cartilage, i.e., increased surface roughness. Thus, imaging-based spectral measurements may provide more information on the cartilage integrity, and consequently may be recommended for further studies.

The presented techniques yielded good results for the bovine cartilage. However, they may not be optimized for human tissue. The development of an optimized algorithm for human cartilage is among our further aims. Importantly, the reflectance measurements and spectral imaging may be carried out with a fiber optic probe during diagnostic arthroscopies. Furthermore, optimal illuminations may be designed for arthroscopy to reveal sensitive and specific diagnostic information about the cartilage.

The financial support from the Ministry of Education, Finland (University of Eastern Finland grant, Projects 5741/Kuopio and 10102/Joensuu), and from the Academy of Finland (Project 127198/Kuopio) is acknowledged.

Guilak  F., , Ratcliffe  A., , Lane  N., , Rosenwasser  M. P., , and Mow  V. C., “ Mechanical and biochemical changes in the superficial zone of articular cartilage in canine experimental osteoarthritis. ,” J. Orthop. Res..  0736-0266 12, , 474–484  ((1994)).
Arokoski  J. P., , Jurvelin  J. S., , Väätäinen  U., , and Helminen  H. J., “ Normal and pathological adaptations of articular cartilage to joint loading. ,” Scand. J. Med. Sci. Sports.  0905-7188 10, (4 ), 186–198  ((2000)).
Buckwalter  J. A., and Mankin  H. J., “ Articular cartilage, part ii: degeneration and osteoarthrosis, repair, regeneration, and transplantation. ,” J. Bone Joint Surg. Am.. 79, (4 ), 612–632  ((1997)).
Bi  X., , Yang  X., , Bostrom  M. P., , and Camacho  N. P., “ Fourier transform infrared imaging spectroscopy investigations in the pathogenesis and repair of cartilage. ,” Biochim. Biophys. Acta.  0006-3002 1758, (7 ), 934–941  ((2006)).
David-Vaudey  E., , Burghardt  A., , Keshari  K., , Brouchet  A., , Ries  M., , and Majumdar  S., “ Fourier transform infrared imaging of focal lesions in human osteoarthritic cartilage. ,” Eur. Cells Mater.  1473-2262 10, , 51–60  ((2005)).
Panula  H., , Hyttinen  M., , Arokoski  J., , Langsjo  T., , Pelttari  A., , Kiviranta  I., , and Helminen  H., “ Articular cartilage superficial zone collagen birefringence reduced and cartilage thickness increased before surface fibrillation in experimental osteoarthritis. ,” Ann. Rheum. Dis..  0003-4967 57, , 237–245  ((1998)).
Bi  X., , Li  G., , Doty  S. B., , and Camacho  N. P., “ A novel method for determination of collagen orientation in cartilage by fourier transform infrared imaging spectroscopy (FT-IRIS). ,” Osteoarthritis Cartilage.  1063-4584 13, , 1050–1058  ((2005)).
Saarakkala  S., , Laasanen  M. A., , Jurvelin  J. S., , and Töyräs  J., “ Quantitative ultrasound imaging detects degenerative changes in articular cartilage surface and subchondral bone. ,” Phys. Med. Biol..  0031-9155 51, , 5333–5346  ((2006)).
Kempson  G. E., , Spivey  C. J., , Swanson  S. A., , and Freeman  M. A., “ Patterns of cartilage stiffness on normal and degenerate human femoral heads. ,” J. Biomech..  0021-9290 4, , 597–609  ((1971)).
Armstrong  C. G., and Mow  V. C., “ Variations in the intrinsic mechanical properties of human articular cartilage with age, degeneration, and water content. ,” J. Bone Joint Surg. Am.. 64, , 88–94  ((1982)).
Hattori  K., , Uematsu  K., , Tanikake  Y., , Habata  T., , Tanaka  Y., , Yajima  H., , and Takakura  Y., “ Spectrocolorimetric assessment of cartilage plugs after autologous osteochondral grafting: correlations between color indices and histological findings in a rabbit model. ,” Arthritis Res. Ther..  1478-6354 9, , R88  ((2007)).
Hattori  K., , Uematsu  K., , Matsumori  H., , Dohi  Y., , Takakura  Y., , and Ohgushi  H., “ Spectrocolorimetric evaluation of repaired articular cartilage after a microfracture. ,” BMC Res. Notes. 1, (1 ), 87  ((2008)).
Öberg  P. Å., , Sundqvist  T., , and Johansson  A., “ Assessment of cartilage thickness utilising reflectance spectroscopy. ,” Med. Biol. Eng. Comput..  0140-0118 42, , 3–8  ((2004)).
Ishimoto  Y., , Hattori  K., , Ohgushi  H., , Uematsu  K., , Tanikake  Y., , Tanaka  Y., , and Takakura  Y., “ Spectrocolorimetric evaluation of human articular cartilage. ,” Osteoarthritis Cartilage.  1063-4584 17, (9 ), 1204–1208  ((2009)).
Ebert  D. W., “ Articular cartilage optical properties in the spectral range 300–850nm. ,” J. Biomed. Opt..  1083-3668 3, (3 ), 326–333  ((1998)).
Qu  J., , MacAulay  C., , Lam  S., , and Palcic  B., “ Optical properties of normal and carcinomatous bronchial tissue. ,” Appl. Opt..  0003-6935 33, (31 ), 7397–7405  ((1994)).
Wyszecki  G., and Stiles  W. S.,  Color Science. , 2nd ed.,  Wiley ,  New York  ((2000)).
Parkkinen  J. P. S., , Hallikainen  J., , and Jääskeläinen  T., “ Characteristic spectra of munsell colors. ,” J. Opt. Soc. Am. A.  0740-3232 6, , 318–322  ((1989)).
Kiviranta  I., , Jurvelin  J., , Tammi  M., , Säämänen  A. M., , and Helminen  H. J., “ Microspectrophotometric quantitation of glycosaminoglycans in articular cartilage sections stained with safranin o. ,” Histochemistry.  0301-5564 82, (3 ), 249–255  ((1985)).
Rieppo  J., , Töyräs  J., , Nieminen  M. T., , Kovanen  V., , Hyttinen  M. M., , Korhonen  R. K., , Jurvelin  J. S., , and Helminen  H. J., “ Structure-function relationships in enzymatically modified articular cartilage. ,” Cells Tissues Organs.  1422-6405 175, (3 ), 121–132  ((2003)).
Laasanen  M. S., , Töyräs  J., , Korhonen  R. K., , Rieppo  J., , Saarakkala  S., , Nieminen  M. T., , Hirvonen  J., , and Jurvelin  J. S., “ Biomechanical properties of knee articular cartilage. ,” Biorheology.  0006-355X 40, (1–3 ), 133–140  ((2003)).
Wang  S., , Huang  Y., , Saarakkala  S., , and Zheng  Y., “ Quantitative assessment of articular cartilage with morphologic, acoustic and mechanical properties obtained using high-frequency ultrasound. ,” Ultrasound Med. Biol..  0301-5629 36, , 512–527  (Mar. (2010)).
Schmitt  J. M., and Kumar  G., “ Optical scattering properties of soft tissue: a discrete particle model. ,” Appl. Opt..  0003-6935 37, (13 ), 2788–2797  ((1998)).

† This is equal to the RGB presentation of the images.

© 2010 Society of Photo-Optical Instrumentation Engineers

Citation

Jussi Kinnunen ; Jukka S. Jurvelin ; Jaana Mäkitalo ; Markku Hauta-Kasari ; Pasi Vahimaa, et al.
"Optical spectral imaging of degeneration of articular cartilage", J. Biomed. Opt. 15(4), 046024 (August 13, 2010). ; http://dx.doi.org/10.1117/1.3477190


Figures

Graphic Jump LocationF1 :

Schematic presentation of the measurement setup as viewed on top. Measurements were conducted trough a glass window with an LCTF camera in 450 geometry while the sample was immersed in PBS.

Graphic Jump LocationF2 :

(a) Mean spectra and (b) the relative change at each wavelength before and after 6h or 15h collagenase digestion or mechanical grinding using the P60 emery paper for all the patellar (PAT) and femoral (FMC) samples. The dashed lines in (b) are the linear fits to the ratios of the two average spectra, before and after the treatments.

Graphic Jump LocationF7 :

Optical microscopy images (safranin-O-stained sections) of the representative patellar samples after mechanical grinding with the P60 emery paper and after 6- and 15-h collagenase digestion. Cartilage thickness was determined from the captured images, calibrated by using the 1-mm ruler seen on the left.

Graphic Jump LocationF6 :

Examples of the visual appearance of one representative mechanically degraded cartilage sample in RGB, under monochromatic 700- and 469-nm light and after projection of spectra to the chosen base vector (the second eigenvector of the femoral samples after 15-h collagenase digestion). The mean contrast between the ground and intact areas in projection image is 1.6 to 2.5 timer higher for all the mechanically degraded samples when compared to image at 469nm.

Graphic Jump LocationF5 :

Mean values (200-×200-pixel center area) of the projections of the spectral images to the chosen base vector (the second eigenvector of the femoral samples after 15-h collagenase digestion) for samples before and after enzymatic (6-15-h collagenase) or mechanical degradations (P60 emery paper). The hollow circles indicate mean values within each sample, while the closed dots with lines represent the mean of samples within each group. The results from mechanical degradation are shown on the left side of the result of the native samples.

Graphic Jump LocationF4 :

Examples of pseudocolor projections of the spectral images to the chosen base vector the second eigenvector of the femoral samples after 15-h collagenase digestion, Fig. 3. Representative patella (PAT) and femoral (FMC) samples before and after 6- or 15-h collagenase digestion or mechanical grinding using the P60 emery paper are shown. The PAT 15-h sample showed spontaneous trauma before collagenase digestion and is presented here for the purpose of comparison. All samples show clear spatial differences including the mechanical degradation marks (“scratches”) and the spontaneous trauma in the projection image before and after the collagenase digestion.

Graphic Jump LocationF3 :

The second eigenvector after 15-h collagenase digestion of femoral samples. The projections of the cartilage spectra to the presented eigenvector were found to separate the enzymatic changes and emery paper scratches from the untreated samples. The vector was chosen to be the base vector for all the projections.

Tables

Table Grahic Jump Location
The slope of the linear fit for the relative change between the mean wavelength dependent reflectances, before and after collagenase digestion (6 or 15h) or mechanical degradation for the patellar (PAT) and fermoral (FMC) samples.
Table Footer NoteThe statistical significance (NS=nonsignificant) between the slopes is calculated by using the Kruskall-Wallis post hoc test.
067004jbot2Table Grahic Jump Location
Value of the CIELAB coordinates (mean±SD) before and after mechanical degradation (mech.) and either 6- or 15-h enzymatic digestion using collagenase (coll.) for the patellar (PAT) and femoral medial condyle (FMC) samples. The RGB presentation of the colors and the ΔEab* color difference before and after treatments are also shown.
Table Footer NoteaStatistically significant difference (Wilcoxon signed rank test, p<0.05) compared to sample before degradation.

References

Guilak  F., , Ratcliffe  A., , Lane  N., , Rosenwasser  M. P., , and Mow  V. C., “ Mechanical and biochemical changes in the superficial zone of articular cartilage in canine experimental osteoarthritis. ,” J. Orthop. Res..  0736-0266 12, , 474–484  ((1994)).
Arokoski  J. P., , Jurvelin  J. S., , Väätäinen  U., , and Helminen  H. J., “ Normal and pathological adaptations of articular cartilage to joint loading. ,” Scand. J. Med. Sci. Sports.  0905-7188 10, (4 ), 186–198  ((2000)).
Buckwalter  J. A., and Mankin  H. J., “ Articular cartilage, part ii: degeneration and osteoarthrosis, repair, regeneration, and transplantation. ,” J. Bone Joint Surg. Am.. 79, (4 ), 612–632  ((1997)).
Bi  X., , Yang  X., , Bostrom  M. P., , and Camacho  N. P., “ Fourier transform infrared imaging spectroscopy investigations in the pathogenesis and repair of cartilage. ,” Biochim. Biophys. Acta.  0006-3002 1758, (7 ), 934–941  ((2006)).
David-Vaudey  E., , Burghardt  A., , Keshari  K., , Brouchet  A., , Ries  M., , and Majumdar  S., “ Fourier transform infrared imaging of focal lesions in human osteoarthritic cartilage. ,” Eur. Cells Mater.  1473-2262 10, , 51–60  ((2005)).
Panula  H., , Hyttinen  M., , Arokoski  J., , Langsjo  T., , Pelttari  A., , Kiviranta  I., , and Helminen  H., “ Articular cartilage superficial zone collagen birefringence reduced and cartilage thickness increased before surface fibrillation in experimental osteoarthritis. ,” Ann. Rheum. Dis..  0003-4967 57, , 237–245  ((1998)).
Bi  X., , Li  G., , Doty  S. B., , and Camacho  N. P., “ A novel method for determination of collagen orientation in cartilage by fourier transform infrared imaging spectroscopy (FT-IRIS). ,” Osteoarthritis Cartilage.  1063-4584 13, , 1050–1058  ((2005)).
Saarakkala  S., , Laasanen  M. A., , Jurvelin  J. S., , and Töyräs  J., “ Quantitative ultrasound imaging detects degenerative changes in articular cartilage surface and subchondral bone. ,” Phys. Med. Biol..  0031-9155 51, , 5333–5346  ((2006)).
Kempson  G. E., , Spivey  C. J., , Swanson  S. A., , and Freeman  M. A., “ Patterns of cartilage stiffness on normal and degenerate human femoral heads. ,” J. Biomech..  0021-9290 4, , 597–609  ((1971)).
Armstrong  C. G., and Mow  V. C., “ Variations in the intrinsic mechanical properties of human articular cartilage with age, degeneration, and water content. ,” J. Bone Joint Surg. Am.. 64, , 88–94  ((1982)).
Hattori  K., , Uematsu  K., , Tanikake  Y., , Habata  T., , Tanaka  Y., , Yajima  H., , and Takakura  Y., “ Spectrocolorimetric assessment of cartilage plugs after autologous osteochondral grafting: correlations between color indices and histological findings in a rabbit model. ,” Arthritis Res. Ther..  1478-6354 9, , R88  ((2007)).
Hattori  K., , Uematsu  K., , Matsumori  H., , Dohi  Y., , Takakura  Y., , and Ohgushi  H., “ Spectrocolorimetric evaluation of repaired articular cartilage after a microfracture. ,” BMC Res. Notes. 1, (1 ), 87  ((2008)).
Öberg  P. Å., , Sundqvist  T., , and Johansson  A., “ Assessment of cartilage thickness utilising reflectance spectroscopy. ,” Med. Biol. Eng. Comput..  0140-0118 42, , 3–8  ((2004)).
Ishimoto  Y., , Hattori  K., , Ohgushi  H., , Uematsu  K., , Tanikake  Y., , Tanaka  Y., , and Takakura  Y., “ Spectrocolorimetric evaluation of human articular cartilage. ,” Osteoarthritis Cartilage.  1063-4584 17, (9 ), 1204–1208  ((2009)).
Ebert  D. W., “ Articular cartilage optical properties in the spectral range 300–850nm. ,” J. Biomed. Opt..  1083-3668 3, (3 ), 326–333  ((1998)).
Qu  J., , MacAulay  C., , Lam  S., , and Palcic  B., “ Optical properties of normal and carcinomatous bronchial tissue. ,” Appl. Opt..  0003-6935 33, (31 ), 7397–7405  ((1994)).
Wyszecki  G., and Stiles  W. S.,  Color Science. , 2nd ed.,  Wiley ,  New York  ((2000)).
Parkkinen  J. P. S., , Hallikainen  J., , and Jääskeläinen  T., “ Characteristic spectra of munsell colors. ,” J. Opt. Soc. Am. A.  0740-3232 6, , 318–322  ((1989)).
Kiviranta  I., , Jurvelin  J., , Tammi  M., , Säämänen  A. M., , and Helminen  H. J., “ Microspectrophotometric quantitation of glycosaminoglycans in articular cartilage sections stained with safranin o. ,” Histochemistry.  0301-5564 82, (3 ), 249–255  ((1985)).
Rieppo  J., , Töyräs  J., , Nieminen  M. T., , Kovanen  V., , Hyttinen  M. M., , Korhonen  R. K., , Jurvelin  J. S., , and Helminen  H. J., “ Structure-function relationships in enzymatically modified articular cartilage. ,” Cells Tissues Organs.  1422-6405 175, (3 ), 121–132  ((2003)).
Laasanen  M. S., , Töyräs  J., , Korhonen  R. K., , Rieppo  J., , Saarakkala  S., , Nieminen  M. T., , Hirvonen  J., , and Jurvelin  J. S., “ Biomechanical properties of knee articular cartilage. ,” Biorheology.  0006-355X 40, (1–3 ), 133–140  ((2003)).
Wang  S., , Huang  Y., , Saarakkala  S., , and Zheng  Y., “ Quantitative assessment of articular cartilage with morphologic, acoustic and mechanical properties obtained using high-frequency ultrasound. ,” Ultrasound Med. Biol..  0301-5629 36, , 512–527  (Mar. (2010)).
Schmitt  J. M., and Kumar  G., “ Optical scattering properties of soft tissue: a discrete particle model. ,” Appl. Opt..  0003-6935 37, (13 ), 2788–2797  ((1998)).

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

Related Content

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

Related Book Chapters

Topic Collections

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

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