If you don't remember your password, you can reset it by entering your email address and clicking the Reset Password button. You will then receive an email that contains a secure link for resetting your password
If the address matches a valid account an email will be sent to __email__ with instructions for resetting your password
Address correspondence and reprint requests to: V.-P. Karjalainen, Research Unit of Medical Imaging, Physics and Technology, Faculty of Medicine, University of Oulu, POB 5000, FI-90014, Oulu, Finland, Tel.: 358400778107.
Lund University, Faculty of Medicine, Department of Clinical Sciences Lund, Orthopaedics, Clinical Epidemiology Unit, Lund, SwedenLund University, Faculty of Medicine, Department of Clinical Sciences Lund, Rheumatology and Molecular Skeletal Biology, Lund, Sweden
Knee osteoarthritis (OA) is associated with meniscal degeneration that may involve disorganization of the meniscal collagen fiber network. Our aims were to quantitatively analyze the microstructural organization of human meniscus samples in 3D using micro-computed tomography (μCT), and to compare the local microstructural organization between OA and donor samples.
We collected posterior horns of both medial and lateral human menisci from 10 end-stage medial compartment knee OA patients undergoing total knee replacement (medial & lateral OA) and 10 deceased donors without knee OA (medial & lateral donor). Posterior horns were dissected and fixed in formalin, dehydrated in ascending ethanol concentrations, treated with hexamethyldisilazane (HMDS), and imaged with μCT. We performed local orientation analysis of collagenous microstructure in 3D by calculating structure tensors from greyscale gradients within selected integration window to determine the polar angle for each voxel.
In donor samples, meniscus bundles were aligned circumferentially around the inner border of meniscus. In medial OA menisci, the organized structure of collagen network was lost, and main orientation was shifted away from the circumferential alignment. Quantitatively, medial OA menisci had the lowest mean orientation angle compared to all groups, −24° (95%CI -31 to −18) vs medial donor and −25° (95%CI -34 to −15) vs lateral OA.
HMDS-based μCT imaging enabled quantitative analysis of meniscal collagen fiber bundles and their orientations in 3D. In human medial OA menisci, the collagen disorganization was profound with overall lower orientation angles, suggesting collagenous microstructure disorganization as an important part of meniscus degradation.
The human meniscus is a fibrocartilaginous soft tissue that has important roles in knee joint function including load bearing, load transmission, joint lubrication, nutrition distribution, shock absorption, and the ability to withstand different mechanical forces such as compression, shear and tension
. The extracellular matrix (ECM) of a healthy meniscus consists of a thick network of collagen fibers (20–25% of wet weight) working as a principal solid component of the ECM, together with other constituents, such as water (65–75% of wet weight), proteoglycans (<1% of wet weight) and non-collagenous proteins
. The collagen fibers in the middle layer orient circumferentially around the inner border of meniscus, and with the middle layer covering most of the tissue volume in meniscus, circumferential orientation is the predominant alignment of the collagen fibers in the whole meniscus
In knee osteoarthritis (OA), it has been reported that meniscal tears and degeneration are associated with disorganization of the collagen fiber network of the meniscus, and that the posterior horn of the medial meniscus is the most susceptible region to these changes
We have previously visualized the microstructure of the posterior horn of the human meniscus ex vivo using a hexamethyldisilazane (HMDS) based sample drying technique together with micro-computed tomography (μCT)
, by dehydrating the sample and creating intrinsic contrast within the soft tissue. In our preliminary study, an HMDS-based sample processing protocol was applied to quantitatively determine the ECM orientation in AC
. However, the local orientation of the collagenous microstructure of human meniscus has not been quantitatively analyzed in 3D and compared between end-stage OA patients and healthy donors.
In this study, our objectives were to: 1) analyze the microstructural organization of ex vivo human meniscus samples in 3D using our previously implemented HMDS-based μCT imaging protocol, and 2) quantitatively compare the local microstructural organization of meniscus posterior horn between medial and lateral menisci of end-stage medial compartment knee OA patients and deceased donors without known knee OA (healthy donors). We hypothesize that there are differences in the local microstructural organization between medial OA menisci and medial donor menisci, and that the lateral OA menisci are more similar to lateral donor menisci.
Tissue sample preparation
The regional ethical review board at Lund University approved this study (Dnr 2015/39 and Dnr 2016/865). The same sample set was used as in our previous study in which the sample selection is explained in more detail
. Briefly, human meniscus samples from the knee tissue biobank MENIX at Skåne University Hospital were used for this study. We selected both the medial and lateral menisci from the right knee of 10 deceased adult donors. The donors had no known diagnosis of knee OA or rheumatoid arthritis. From the same biobank, we also retrieved both the medial and lateral menisci collected from 10 end-stage medial compartment knee OA patients (all medial samples had Outerbridge grade IV, while lateral samples had grade 0 or I) during total knee replacement (TKR) surgery. Both OA and donor groups include five men and five women. Hereafter, these medial and lateral sample groups from deceased donors and TKR patients are referred to as medial donor and lateral donor, and medial OA and lateral OA, respectively. All the menisci were frozen at −80°C within 2 h of extraction. The samples were thawed in phosphate buffered saline (PBS), and the posterior horn was separated from the meniscus body with a scalpel followed by fixation in 4% saline-buffered formaldehyde for a minimum of 11 days. A few millimeter-thick tissue sections were cut from the fixed posterior horns for μCT imaging and following analyses.
HMDS-based μCT imaging
The same μCT dataset was used as in our previous work
. Briefly, the fixed samples were dehydrated in ascending ethanol concentrations before treating them with HMDS and further leaving them to air-dry in a fume hood at room temperature overnight. For μCT image acquisition, we used a desktop device (SkyScan 1272, Bruker microCT, Kontich, Belgium) with the following settings: tube voltage 40 kV, tube current 250 μA, isotropic voxel size 2.0 μm, exposure time 1815 ms, image resolution 4032 × 2688, no additional filtering, number of projections 2400 and averaging five frames/projection. NRecon software (Bruker microCT, V 220.127.116.11) was used for image reconstruction from acquired projections. During image reconstruction, beam-hardening and ring-artifact corrections were applied
Local orientation analysis provides the local spatial 3D orientation of tissue structures. Orientation analysis was conducted with CTAn software (Bruker microCT, Local Orientation plugin, ver. 1.18) which calculates the local fiber orientation using the method of Straumit et al. (2015) to determine the principal direction of voxels using structure tensors in μCT image data (Fig. 1)
. The analysis calculates theta (Ө) angles (polar angles) for each voxel from the structure tensors to define the orientation direction of the object within the integration window.
First, a cubic 1800 μm × 1800 μm × 1800 μm sized volume-of-interest (VOI) from the middle of the meniscus posterior horn piece was selected from the reconstructed 3D μCT image stack for the local orientation analysis (Fig. 2). The selected VOI dimensions were chosen due to computational limitations and to include mostly the collagen fiber bundles in the middle layer and not in the superficial or lamellar layers, while avoiding imaging artifacts and largest calcifications that do not account for the orientation information.
Local orientation analysis was applied to the selected VOIs from all the samples with an integration window with a size of 19 voxels, i.e., the diameter that is used to calculate the direction gradient for orientation angle. In the beginning of the analysis, the μCT images are converted into 3D arrays of grey scale values from which structure tensors can be defined as shown in Fig. 1. The smallest eigenvalue and the corresponding eigenvector from decomposition of the structure tensor matrix give the direction of the local fiber orientation (polar angle Ө) for the voxel. The procedure filters falsely classified voxels by comparing alignments of neighboring voxels in the integration window
. The polar angle Ө is in the Z dimension perpendicular to the cross-section X-Y image plane. Polar angle Ө = 0 is at 90° to the X–Y plane to Z dimension. Fig. 2(D) shows how the polar angle Ө is defined between 0 and 90°.
A global background thresholding and subtraction were implemented on the VOIs after local orientation analysis to prioritize calculating orientations from intact collagen fibers over loose connective tissue. Finally, the percentage of voxels at each angle (0–90°) for each sample was calculated. Lower orientation angles indicate more disorganization or non-normal orientations for the collagen fibers. In contrast, higher orientation angles present the dominating circumferential orientation in the middle layer of intact menisci.
We present the data as the percentage of voxel per given angle (0–90°) for each sample in the figures. Further, we calculated the mean angle in each sample as a weighted average to account for the different number of voxels in each sample. We estimated unadjusted differences and also differences adjusted for age and body mass index (BMI). We did not consider adjustment for sex, as the number of males and females was the same in OA and donors. To analyze the differences in mean angle between the groups, we used a linear regression model weighted by the total number of voxels in each sample with clustered standard errors to account for the dependence of samples coming from the same individual
. The group with four levels (medial OA, lateral OA, medial donor, lateral donor) was included as independent variable, with medial OA being the reference category. Further, to analyze the shape of the distribution of voxels across angles, we used Poisson regression model. The outcome was the percentage of voxels at a given angle. The independent variables included restricted cubic splines with knots at 5th, 10th, 30th, 80th and 90th percentile of angles, group (OA vs donor) and interaction between group and the restricted cubic splines to allow for different shape of the relationship in each group. Further, to account for dependence of the observations coming from the same person we included a random intercept for individual. We fitted a separate model for medial compartment and a separate model for the lateral compartment. In a sensitivity analysis, we additionally adjusted for age and BMI. The observed and fitted data can be seen in the supplementary material (Fig. S1). From these models, we extracted predicted group specific means with 95% CIs to illustrate the group specific trajectories with uncertainty at all angles.
Descriptive statistics on the study participants, both OA patients and deceased donors in Table I. The mean age of OA patients was higher than the mean age of donors while mean height and mean weight were similar between the two groups.
Collagen fiber bundles in HMDS-treated meniscal samples were successfully depicted in 3D using μCT. In medial donor, lateral donor and lateral OA samples, meniscus bundles were aligned circumferentially around the inner border of meniscus. In the medial OA meniscus, the organized structure of collagen network was lost, and the main orientation was shifted away from the circumferential alignment. In Fig. 3, we show example images of μCT 3D volumes from OA patients and deceased donors, their selected VOIs and finally analyzed VOIs representing orientation angles of the meniscus microstructure. Furthermore, all selected and analysed 3D μCT volumes are presented in the supplementary material (Figs. S3–S6).
In quantitative local orientation analysis and statistical comparison between the medial OA group and other sample groups, we observed 24° (95% CI 18, 31), 25° (95% CI 15, 34) and 25° (95% CI 16, 35) lower orientation angles in medial OA menisci compared to medial donor, lateral OA, and lateral donor menisci respectively (Table II). Furthermore, medial donor, lateral OA and lateral donor menisci yielded similar mean angles with 68°, 66° and 68°, respectively. The orientation angles between lateral OA and lateral donor menisci were also similar with a mean difference of 2°.
Table IIThe mean angles of collagen fiber bundles within the following groups: medial OA, medial donor, lateral OA and lateral donor
Mean estimated differences vs medial OA (95% CI)
Mean estimated differences vs medial OA adjusted for age and BMI (95% CI)
In Fig. 4, we present the percentage of voxels as a function of orientation angle in the sample groups, showing diverse distribution in medial OA menisci with increase in low orientation angles and decrease in high orientation angles compared to other groups. Medial donor and lateral OA menisci had similar low and high angle distributions, while medial donors had slightly higher amount of low angle orientations of collagen fibers (Fig. 4). Lateral donor menisci had the smallest amount of low angle orientations indicating the most intact (normal) tissue. In Fig. 5(A) we show the percentage of voxels as a function of orientation angle in each sample from OA samples and in Fig. 5(B) from donor samples. The results after adjustment for age and BMI were essentially the same (Table S2).
We found that using our HMDS-based sample drying protocol together with μCT imaging, it is possible not only to visualize, but also to quantitatively analyze the local microstructural organization of collagen network in the posterior horn of human meniscus in 3D. The results of the local microstructural orientation analysis indicate substantial differences in collagen organization between osteoarthritic and donor medial posterior horn menisci. Increased disorganization of collagen fiber structure in the medial OA samples suggest a strong association between meniscus collagen (dis)organization and knee OA; the orientation of fibers in the medial OA menisci was on average 15–30° lower compared to other groups.
We found increased tissue disorganization in the medial OA menisci, quantitatively indicated by a large difference in local orientation angles of collagen fiber bundles in the medial OA menisci when compared to lateral donor, lateral OA and medial donor menisci (Table II, Fig. 4). Similarly, the estimated mean orientation angle in medial OA menisci was lower compared to all other groups (Table II). In a related mass spectrometry study using meniscus tissue plugs from same OA patients and donors as in our study, several matrix proteins such as matrix metalloproteinase three and TIMP1 were found in higher intensities in the medial OA menisci compared to medial reference menisci
. These results could be indicative of changes in the ECM, which, in our study are shown as collagen matrix disorganization. Furthermore, our results are supported by a previous study using transmission electron microscopy, where osteoarthritic meniscus tissue was reported to have degenerated ECM and disorganized collagen fiber network when compared to healthy tissue
. This affects a large proportion of people with end-stage knee OA who experience edematous swelling of the tissue, especially in the medial meniscus, which could be related to meniscal extrusion, i.e. unloading of the tissue frequently seen in end-stage medial compartment knee OA
. The disruption in collagen structure could, therefore, be caused by the swelling of these fluid compartments, which causes the collagen bundles to burst from within, as seen in our μCT images of osteoarthritic samples. This hypertrophy-induced swelling of meniscus is further supported by our previous study, in which we found that intact menisci were smaller in size than degenerated menisci
. Consequently, we believe that in OA, hypertrophy has a great impact on the disruption of large proportion of the collagen network in meniscus. However, these effects and their association with commonly observed meniscal tears and extrusion require further study.
From our quantitative μCT analyses, we observed several similarities in the microstructural local orientations between medial donor, lateral OA, and lateral donor meniscus. Additionally, the estimated mean orientation angles of medial donor, lateral OA and lateral donor were similar (68°, 66° and 68°, respectively). In our previous study, we reported similar histopathological scores for the same sample set between lateral donor, medial donor and lateral OA groups of meniscus posterior horn
. In the present study, out of these three groups, lateral donor meniscus had the lowest amount of small angle orientations suggesting higher preservation of its microstructural organization than in medial donor and lateral OA meniscus. Moreover, in Fig. 4, slight differences in the angle distribution between the lateral OA and the lateral donor menisci can be observed. In the lateral meniscus with medial compartment OA one could expect minor organizational changes compared to the lateral donor, even when their Outerbridge grades were 0 or I. However, judging from the confidence intervals, differences between donors and OA in the lateral compartment, if any, are expected to be smaller than in the medial compartment. Furthermore, the medial donor meniscus had slightly more degenerated and disorganized collagen network compared to the lateral donor meniscus. This could be expected since the weight-bearing function and biomechanical workload on the medial side of meniscus is often higher compared to the lateral side
. It is also possible that the donor menisci have some “pre-osteoarthritic” changes which would occur most likely in the medial side. Furthermore, between medial donor and lateral OA meniscus, medial donor meniscus had higher amount of small angle orientations indicating higher disorganization in its microstructure. Thus, our results suggest that the disorganization of collagen network in meniscus is more prevalent in the medial donor meniscus than in the lateral meniscus with end-stage medial compartment knee OA. In addition, this medial meniscus susceptibility to degeneration is supported by previous study showing that the medial meniscus is more vulnerable to meniscal tears and degeneration than lateral meniscus
. However, all of these techniques are limited to 5–60 μm thin tissue sections, losing the real 3D information because a single collagen fiber averages between 5 and 10 μm in diameter and collagen fiber bundles together are even larger
. We have reported that our HMDS-based μCT method enables the measurement of larger volumetric structures in human posterior horn meniscus than conventional 2D section-based histology or above-mentioned methods. Furthermore, in the beginning of degeneration the meniscal surface may remain intact while distinct changes in meniscal internal structure and composition are already observed
. For further studies, with this method it is possible to measure quantitatively the collagen fiber bundle orientations from any compartment or layer of meniscus to study the increasing disorganization of collagen network in OA.
In principle, it is possible to compare our collagen network disorganization results in meniscus with the disorganization of collagen network in AC. In fact, our HMDS protocol has been applied to quantitatively determine the orientation of ECM in AC
. However, while meniscus and AC share similarities in their structural layers, the collagen network is completely different in the circumferentially organized middle layer of meniscus compared to the deep layer of AC
. Thus, in both meniscus and AC, the disorganization of collagen network has similar features and can be quantitatively measured with our method and compared to each other with this limitation considered.
Different methods to characterize the structural orientations have been developed for different materials and tissues. Instead of structure tensors, it is possible to prepare an application of an anisotropic filter and perform a step by step rotation of a region in the image to find the most probable orientation of fibers
. In medical research, diffusion tensor image analysis is widely used in magnetic resonance imaging (MRI) to characterize structural changes by measuring orientations from diffusing water molecules in brain tissue
. In our study, we used HMDS-based protocol, in which the samples are dried to enable the μCT imaging with high resolution. Local orientation analysis was tested to be compatible with HMDS-based dried meniscus samples and μCT imaging. It has also been reported for having low uncertainty when classifying voxels
The present study has some important limitations. First, the menisci specimens were frozen after collection and then thawed before the initiation of the study. In a previous study, it was reported that freezing may cause changes in the meniscus collagen network, including shrinking in collagen fiber diameters and higher disorder in their alignment
. However, we performed the same preparation protocol to all samples, both from OA cases and donors, and therefore, the results are more likely to be comparable. Second, HMDS-based sample drying method preserves the sample to enable the μCT imaging but may not truly reflect the sample as it was in vivo. Furthermore, the increase of water content and the decrease of collagen content in the degenerated menisci compared to donor menisci introduces a different ratio of constituents between OA and donor menisci
. Thus, it is theoretically possible that HMDS-based drying would cause larger shrinking of OA samples compared to donor samples. However, this difference should be minimal in practice as the HMDS reagent is believed to cross-link proteins (collagen) and it has lower surface tension compared to water
. This should prevent the tissue from collapsing in all cases. Moreover, even if the shrinkage due to evaporation of water is minimal, it may be a possible limitation of this contrast method due to having a greater effect in the degenerated menisci with higher water content
. Third, some bony attachment of the posterior root may to various extent be left in the meniscus samples after the extraction. However, our VOI is from the posterior horn distant to the ligamentous attachment and, thus, the meniscal root is not included in the analyses. Fourth, μCT-imaging technique may introduce ring artifacts, which induce error to the orientation angle values in affected voxels. Fifth, we acknowledge that the Poisson regression and splines models work best with larger sample sizes, but considering that the splines are based on a large number of voxels, we find our modeling approach reasonable considering the exploratory nature of the study. As a final limitation, we grouped male and female specimen together due to otherwise limited sample numbers, but the number of males and females was the same in both OA and donors.
To conclude, the local orientation analysis from HMDS-treated samples imaged with conventional desktop μCT allows visualization and quantitative measurements of the collagen fiber network in 3D in the human meniscus posterior horn. We conclude a higher disorganization of the collagen network in medial OA menisci when compared to lateral donor, medial donor and lateral OA menisci.
Conception and design: VPK, IK, MF, ME and SS.
Provision of study materials and tissue preparation: ME, EF, VH, PÖ, JT, IK and MF, SS.
Micro-CT imaging: IK and MF.
Image processing/assessment: VPK and IK.
Statistical analysis: VPK and AT.
Interpretation of results: All coauthors.
Drafting of the article: VPK.
Critical revision of the article for important intellectual content: IK, MF, EF, AT, PÖ, VH, JT, ME and SS.
Final approval of the article: All coauthors.
AT works as an associate editor (statistics) in Osteoarthritis and Cartilage.
ME has received grants from European Research Council, The Swedish Research Council, the Foundation for Research in Rheumatology, the Greta and Johan Kock Foundation, the Swedish Rheumatism Association, the Österlund Foundation, the Governmental Funding of Clinical Research program within the National Health Service in Sweden, and the Faculty of Medicine, Lund University, Sweden. ME reports in 2019 serving on an advisory board for Pfizer (Tanezumab).
SS has received grants from Foundation for Research in Rheumatology and European Research Council.
Other authors (VK, IK, MF, EF, PÖ, VH, and JT) report no conflicts of interest.
Role of the funding source
The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
This project has received funding from the European Research Council funded this project under European Union's Horizon 2020 research and innovation programme and under European Union's Seventh Framework Programme (FP/2007–2013) (grant agreement No 336267). This work was also supported by the Swedish Research Council ( Dnr 2014–2348 ), the Foundation for Research in Rheumatology (FOREUM) (018EnglundPreCl), the Greta and Johan Kock Foundation , the Swedish Rheumatism Association , The Österlund Foundation, The Governmental Funding of Clinical Research program within the National Health Service (ALF) in Sweden, and the Faculty of Medicine, Lund University , Sweden.
We would like to thank the staff of Tissue Donor Bank at Skåne University Hospital, the Department of Forensic Medicine in Lund, and the MENIX clinical personnel at Trelleborg Hospital for collaborating with the sample collection.
Appendix A. Supplementary data
The following is the Supplementary data to this article: