3D analysis and grading of calcifications from ex vivo human meniscus.

OBJECTIVE: Meniscal calcifications are associated with the pathogenesis of knee osteoarthritis (OA). We propose a micro-computed tomography (µCT) based 3D analysis of meniscal calcifications ex vivo, including a new grading system. METHOD: Human medial and lateral menisci were obtained from 10 patients having total knee replacement for medial compartment OA and 10 deceased donors without knee OA (healthy references). The samples were fixed; one subsection was imaged with µCT, and the adjacent tissue was processed for histological evaluation. Calcifications were examined from the reconstructed 3D µCT images, and a new grading system was developed. To validate the grading system, meniscal calcification volumes (CVM) were quantitatively analyzed and compared between the calcification grades. Furthermore, we estimated the relationship between histopathological degeneration and the calcification severity. RESULTS: 3D µCT images depict calcifications in every sample, including diminutive calcifications that are not visible in histology. In the new grading system, starting from grade 2, each grade results in a CVM that is 20.3 times higher (95% CI 13.3-30.5) than in the previous grade. However, there was no apparent difference in CVM between grades 1 and 2. The calcification grades appear to increase with the increasing histopathological degeneration, although histopathological degeneration is also observed with small calcification grades. CONCLUSIONS: 3D µCT grading of meniscal calcifications is feasible. Interestingly, it seems that there are two patterns of degeneration in the menisci of our sample set: 1) with diminutive calcifications (calcification grades 1-2), and 2) with large to widespread calcifications (calcification grades 3-5).

The CMS19 disease model specifies a pivotal role for collagen XIII in bone homeostasis.

Mutations in the COL13A1 gene result in congenital myasthenic syndrome type 19 (CMS19), a disease of neuromuscular synapses and including various skeletal manifestations, particularly facial dysmorphisms. The phenotypic consequences in Col13a1 null mice (Col13a1-/-) recapitulate the muscle findings of the CMS19 patients. Collagen XIII (ColXIII) is exists as two forms, a transmembrane protein and a soluble molecule. While the Col13a1-/- mice have poorly formed neuromuscular junctions, the prevention of shedding of the ColXIII ectodomain in the Col13a1tm/tm mice results in acetylcholine receptor clusters of increased size and complexity. In view of the bone abnormalities in CMS19, we here studied the tubular and calvarial bone morphology of the Col13a1-/- mice. We discovered several craniofacial malformations, albeit less pronounced ones than in the human disease, and a reduction of cortical bone mass in aged mice. In the Col13a1tm/tm mice, where ColXIII is synthesized but the ectodomain shedding is prevented due to a mutation in a protease recognition sequence, the cortical bone mass decreased as well with age and the cephalometric analyses revealed significant craniofacial abnormalities but no clear phenotypical pattern. To conclude, our data indicates an intrinsic role for ColXIII, particularly the soluble form, in the upkeep of bone with aging and suggests the possibility of previously undiscovered bone pathologies in patients with CMS19.

Quantitative three-dimensional collagen orientation analysis of human meniscus posterior horn in health and osteoarthritis using micro-computed tomography.

OBJECTIVE: 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. METHOD: 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. RESULTS: 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. CONCLUSIONS: 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.

Design and development of poly-L/D-lactide copolymer and barium titanate nanoparticle 3D composite scaffolds using breath figure method for tissue engineering applications.

In tissue engineering, the scaffold topography influences the adhesion, proliferation, and function of cells. Specifically, the interconnected porosity is crucial for cell migration and nutrient delivery in 3D scaffolds. The objective of this study was to develop a 3D porous composite scaffold for musculoskeletal tissue engineering applications by incorporating barium titanate nanoparticles (BTNPs) into a poly-L/D-lactide copolymer (PLDLA) scaffold using the breath figure method. The porous scaffold fabrication utilised 96/04 PLDLA, dioleoyl phosphatidylethanolamine (DOPE), and different types of BTNPs, including uncoated BTNPs, Al2O3-coated BTNPs, and SiO2-coated BTNPs. The BTNPs were incorporated into the polymer scaffold, which was subsequently analysed using field emission scanning electron microscopy (FE-SEM). The biocompatibility of each scaffold was tested using ovine bone marrow stromal stem cells. The cell morphology, viability, and proliferation were evaluated using FE-SEM, LIVE/DEAD staining, and Prestoblue assay. Porous 3D composite scaffolds were successfully produced, and it was observed that the incorporation of uncoated BTNPs increased the average pore size from 1.6 µm (PLDLA) to 16.2 µm (PLDLA/BTNP). The increased pore size in the PLDLA/BTNP scaffolds provided a suitable porosity for the cells to migrate inside the scaffold, while in the pure PLDLA scaffolds with their much smaller pore size, cells elongated on the surface. To conclude, the breath figure method was successfully used to develop a PLDLA/BTNP scaffold. The use of uncoated BTNPs resulted in a composite scaffold with an optimal pore size while maintaining the honeycomb-like structure. The composite scaffolds were biocompatible and yielded promising structures for future tissue engineering applications.

Comparison of material models for anterior cruciate ligament in tension: from poroelastic to a novel fibril-reinforced nonlinear composite model.

Computational models of the knee joint are useful for evaluating stresses and strains within the joint tissues. However, the outcome of those models is sensitive to the material model and material properties chosen for ligaments, the collagen reinforced tissues connecting bone to bone. The purpose of this study was to investigate different compositionally motivated material models and further to develop a model that can accurately reproduce experimentally measured stress-relaxation data of bovine anterior cruciate ligament (ACL). Tensile testing samples were extracted from ACLs of bovine knee joints (N = 10) and subjected to a three-step stress-relaxation test at the toe region. Data from the experiments was averaged and one average finite element model was generated to replicate the experiment. Poroelastic and different fibril-reinforced poro(visco)elastic material models were applied, and their material parameters were optimized to reproduce the experimental force-time response. Material models with only fluid flow mediated relaxation were not able to capture the stress-relaxation behavior (R2 = 0.806, 0.803 and 0.938). The inclusion of the viscoelasticity of the fibrillar network improved the model prediction (R2 = 0.978 and 0.976), but the complex stress-relaxation behavior was best captured by a poroelastic model with a nonlinear two-relaxation-time strain-recruited viscoelastic fibrillar network (R2 = 0.997). The results suggest that in order to replicate the multi-step stress-relaxation behavior of ACL in tension, the fibrillar network formulation should include the complex nonlinear viscoelastic phenomena.

Automating three-dimensional osteoarthritis histopathological grading of human osteochondral tissue using machine learning on contrast-enhanced micro-computed tomography.

OBJECTIVE: To develop and validate a machine learning (ML) approach for automatic three-dimensional (3D) histopathological grading of osteochondral samples imaged with contrast-enhanced micro-computed tomography (CEµCT). DESIGN: A total of 79 osteochondral cores from 24 total knee arthroplasty patients and two asymptomatic donors were imaged using CEµCT with phosphotungstic acid -staining. Volumes-of-interest (VOI) in surface (SZ), deep (DZ) and calcified (CZ) zones were extracted depth-wise and subjected to dimensionally reduced Local Binary Pattern -textural feature analysis. Regularized linear and logistic regression (LR) models were trained zone-wise against the manually assessed semi-quantitative histopathological CEµCT grades (diameter = 2 mm samples). Models were validated using nested leave-one-out cross-validation and an independent test set (4 mm samples). The performance was primarily assessed using Mean Squared Error (MSE) and Average Precision (AP, confidence intervals are given in square brackets). RESULTS: Highest performance on cross-validation was observed for SZ, both on linear regression (MSE = 0.49, 0.69 and 0.71 for SZ, DZ and CZ, respectively) and LR (AP = 0.9 [0.77-0.99], 0.46 [0.28-0.67] and 0.65 [0.41-0.85] for SZ, DZ and CZ, respectively). The test set evaluations yielded increased MSE on all zones. For LR, the performance was also best for the SZ (AP = 0.85 [0.73-0.93], 0.82 [0.70-0.92] and 0.8 [0.67-0.9], for SZ, DZ and CZ, respectively). CONCLUSION: We present the first ML-based automatic 3D histopathological osteoarthritis (OA) grading method which also adequately perform on grading unseen data, especially in SZ. After further development, the method could potentially be applied by OA researchers since the grading software and all source codes are publicly available.

Comparison of water, hydroxyproline, uronic acid and elastin contents of bovine knee ligaments and patellar tendon and their relationships with biomechanical properties.

Mechanical material properties of ligaments originate from their biochemical composition and structural organization. However, it is not yet fully elucidated how biochemical contents vary between knee ligaments and patellar tendon (PT) and how they relate with mechanical properties. The purpose of this study was to compare water, collagen, proteoglycan and elastin contents between bovine knee ligaments and PT and correlate them with tensile material properties. Hydroxyproline (collagen), uronic acid (proteoglycan) and elastin contents per wet and dry weights were measured using colorimetric biochemical methods for bovine knee ligament and PT samples (n = 10 knees). Direct comparison and correlation with multiple linear regression were performed against biomechanical properties measured in our earlier study. Anterior cruciate ligament (ACL) and PT exhibited lower hydroxyproline content per wet weight compared with other ligaments (p < 0.05). Cruciate ligaments had higher uronic acid content per dry weight compared with collateral ligaments (p < 0.05). Posterior cruciate ligament had higher elastin content than ACL (p < 0.05). Higher hydroxyproline content per wet weight implied higher Young's modulus, strength and toughness. Quantitatively, higher elastin content per wet weight predicted higher toe region nonlinearity and Young's modulus whereas higher uronic acid content per dry weight predicted lower Young's modulus, yield stress and toughness. Differences between ligaments in biochemical composition highlight differences in their physiological function and loading regimes. As expected, collagen content showed similar trend with stiffness and strength. The predictive role of proteoglycan and elastin contents on the mechanical properties might indicate their important functional role in ligaments.

Quantifying Subresolution 3D Morphology of Bone with Clinical Computed Tomography.

The aim of this study was to quantify sub-resolution trabecular bone morphometrics, which are also related to osteoarthritis (OA), from clinical resolution cone beam computed tomography (CBCT). Samples (n = 53) were harvested from human tibiae (N = 4) and femora (N = 7). Grey-level co-occurrence matrix (GLCM) texture and histogram-based parameters were calculated from CBCT imaged trabecular bone data, and compared with the morphometric parameters quantified from micro-computed tomography. As a reference for OA severity, histological sections were subjected to OARSI histopathological grading. GLCM and histogram parameters were correlated to bone morphometrics and OARSI individually. Furthermore, a statistical model of combined GLCM/histogram parameters was generated to estimate the bone morphometrics. Several individual histogram and GLCM parameters had strong associations with various bone morphometrics (|r| > 0.7). The most prominent correlation was observed between the histogram mean and bone volume fraction (r = 0.907). The statistical model combining GLCM and histogram-parameters resulted in even better association with bone volume fraction determined from CBCT data (adjusted R2 change = 0.047). Histopathology showed mainly moderate associations with bone morphometrics (|r| > 0.4). In conclusion, we demonstrated that GLCM- and histogram-based parameters from CBCT imaged trabecular bone (ex vivo) are associated with sub-resolution morphometrics. Our results suggest that sub-resolution morphometrics can be estimated from clinical CBCT images, associations becoming even stronger when combining histogram and GLCM-based parameters.

Three-dimensional microstructure of human meniscus posterior horn in health and osteoarthritis.

OBJECTIVE: To develop and perform ex vivo 3D imaging of meniscus posterior horn microstructure using micro-computed tomography (µCT), and to compare specimens from healthy references against end-stage osteoarthritis (OA) using conventional section-based histology and qualitative µCT. DESIGN: We retrieved human medial and lateral menisci from 10 deceased donors without knee OA (healthy references) and medial and lateral menisci from 10 patients having total knee replacement for medial compartment OA. Meniscal posterior horns were dissected and fixed in formalin. One subsection underwent hexamethyldisilazane (HMDS) treatment and µCT imaging. Pauli's histopathological scoring was performed for 3 other subsections. The differences in histopathological scores were estimated using mixed linear regression, resulting in fixed effects estimates for within-knee comparisons and adjusted for age and body mass index for between-subjects comparisons. RESULTS: 3D visualization with µCT qualitatively revealed similar microstructural changes in the posterior horns as conventional histology. The mean histopathological score was higher for medial menisci from OA knees vs both medial reference menisci (mean difference [95% CI], 3.9 [2.6,5.3]), and lateral menisci from OA knees (3.9 [2.9,5.0]). The scores were similar between lateral menisci from OA knees and lateral reference menisci (0.8 [-0.6,2.2]), and between medial and lateral reference menisci (0.8 [-0.3,1.9]). CONCLUSIONS: HMDS-based µCT protocol allows unique 3D visualization of meniscus microstructures. Posterior horns of medial menisci from medial compartment OA knees had higher histopathological scores than both the lateral posterior horns from the same OA knees and medial reference menisci, suggesting a strong association between meniscus degradation and unicompartmental knee OA.

3D morphometric analysis of calcified cartilage properties using micro-computed tomography.

OBJECTIVE: Our aim is to establish methods for quantifying morphometric properties of calcified cartilage (CC) from micro-computed tomography (µCT). Furthermore, we evaluated the feasibility of these methods in investigating relationships between osteoarthritis (OA), tidemark surface morphology and open subchondral channels (OSCCs). METHOD: Samples (n = 15) used in this study were harvested from human lateral tibial plateau (n = 8). Conventional roughness and parameters assessing local 3-dimensional (3D) surface variations were used to quantify the surface morphology of the CC. Subchondral channel properties (percentage, density, size) were also calculated. As a reference, histological sections were evaluated using Histopathological osteoarthritis grading (OARSI) and thickness of CC and subchondral bone (SCB) was quantified. RESULTS: OARSI grade correlated with a decrease in local 3D variations of the tidemark surface (amount of different surface patterns (rs = -0.600, P = 0.018), entropy of patterns (EP) (rs = -0.648, P = 0.018), homogeneity index (HI) (rs = 0.555, P = 0.032)) and tidemark roughness (TMR) (rs = -0.579, P = 0.024). Amount of different patterns (ADP) and EP associated with channel area fraction (CAF) (rp = 0.876, P < 0.0001; rp = 0.665, P = 0.007, respectively) and channel density (CD) (rp = 0.680, P = 0.011; rp = 0.582, P = 0.023, respectively). TMR was associated with CAF (rp = 0.926, P < 0.0001) and average channel size (rp = 0.574, P = 0.025). CC topography differed statistically significantly in early OA vs healthy samples. CONCLUSION: We introduced a µ-CT image method to quantify 3D CC topography and perforations through CC. CC topography was associated with OARSI grade and OSCC properties; this suggests that the established methods can detect topographical changes in tidemark and CC perforations associated with OA.

Composition, structure and tensile biomechanical properties of equine articular cartilage during growth and maturation.

Articular cartilage undergoes structural and biochemical changes during maturation, but the knowledge on how these changes relate to articular cartilage function at different stages of maturation is lacking. Equine articular cartilage samples of four different maturation levels (newborn, 5-month-old, 11-month-old and adult) were collected (N = 25). Biomechanical tensile testing, Fourier transform infrared microspectroscopy (FTIR-MS) and polarized light microscopy were used to study the tensile, biochemical and structural properties of articular cartilage, respectively. The tensile modulus was highest and the breaking energy lowest in the newborn group. The collagen and the proteoglycan contents increased with age. The collagen orientation developed with age into an arcade-like orientation. The collagen content, proteoglycan content, and collagen orientation were important predictors of the tensile modulus (p < 0.05 in multivariable regression) and correlated significantly also with the breaking energy (p < 0.05 in multivariable regression). Partial least squares regression analysis of FTIR-MS data provided accurate predictions for the tensile modulus (r = 0.79) and the breaking energy (r = 0.65). To conclude, the composition and structure of equine articular cartilage undergoes changes with depth that alter functional properties during maturation, with the typical properties of mature tissue reached at the age of 5-11 months.

In vitro method for 3D morphometry of human articular cartilage chondrons based on micro-computed tomography.

OBJECTIVE: The aims of this study were: to 1) develop a novel sample processing protocol to visualize human articular cartilage (AC) chondrons using micro-computed tomography (µCT), 2) develop and validate an algorithm to quantify the chondron morphology in 3D, and 3) compare the differences in chondron morphology between intact and osteoarthritic AC. METHOD: The developed protocol is based on the dehydration of samples with hexamethyldisilazane (HMDS), followed by imaging with a desktop µCT. Chondron density and depth, as well as volume and sphericity, were calculated in 3D with a custom-made and validated algorithm employing semi-automatic chondron selection and segmentation. The quantitative parameters were analyzed at three AC depth zones (zone 1: 0-10%; zone 2: 10-40%; zone 3: 40-100%) and grouped by the OARSI histological grades (OARSI grades 0-1.0, n = 6; OARSI grades 3.0-3.5, n = 6). RESULTS: After semi-automatic chondron selection and segmentation, 1510 chondrons were approved for 3D morphometric analyses. The chondrons especially in the deeper tissue (zones 2 and 3) were significantly larger (P < 0.001) and less spherical (P < 0.001), respectively, in the OARSI grade 3-3.5 group compared to the OARSI grade 0-1.0 group. No statistically significant difference in chondron density between the OARSI grade groups was observed at different depths. CONCLUSION: We have developed a novel sample processing protocol for chondron imaging in 3D, as well as a high-throughput algorithm to semi-automatically quantify chondron/chondrocyte 3D morphology in AC. Our results also suggest that 3D chondron morphology is affected by the progression of osteoarthritis (OA).

Imaging of Osteoarthritic Human Articular Cartilage using Fourier Transform Infrared Microspectroscopy Combined with Multivariate and Univariate Analysis.

The changes in chemical composition of human articular cartilage (AC) caused by osteoarthritis (OA) were investigated using Fourier transform infrared microspectroscopy (FTIR-MS). We demonstrate the sensitivity of FTIR-MS for monitoring compositional changes that occur with OA progression. Twenty-eight AC samples from tibial plateaus were imaged with FTIR-MS. Hyperspectral images of all samples were combined for K-means clustering. Partial least squares regression (PLSR) analysis was used to compare the spectra with the OARSI grade (histopathological grading of OA). Furthermore, the amide I and the carbohydrate regions were used to estimate collagen and proteoglycan contents, respectively. Spectral peak at 1338 cm(-1) was used to estimate the integrity of the collagen network. The layered structure of AC was revealed using the carbohydrate region for clustering. Statistically significant correlation was observed between the OARSI grade and the collagen integrity in the superficial (r = -0.55) and the deep (r = -0.41) zones. Furthermore, PLSR models predicted the OARSI grade from the superficial (r = 0.94) and the deep (r = 0.77) regions of the AC with high accuracy. Obtained results suggest that quantitative and qualitative changes occur in the AC composition during OA progression, and these can be monitored by the use of FTIR-MS.

Effects of developmental exposure to perfluorooctanoic acid (PFOA) on long bone morphology and bone cell differentiation.

Perfluorooctanoic acid (PFOA) is a ubiquitous and persistent environmental chemical, which has been used extensively due to its stability and surface tension-lowering properties. Toxicological effects include induction of neonatal mortality and reproductive toxicity. In this study, pregnant C57BL/6 mice were exposed orally to 0.3mg PFOA/kg/day throughout pregnancy, and female offspring were studied at the age of 13 or 17months. Morphometrical and biomechanical properties of femurs and tibias were analyzed with micro-computed tomography and 3-point bending, and bone PFOA concentrations were determined by mass spectrometry. The effects of PFOA on bone cell differentiation were studied in osteoclasts from C57BL/6 mice and in the MC3T3 pre-osteoblast cell line. PFOA exposed mice showed increased femoral periosteal area as well as decreased mineral density of tibias. Biomechanical properties of these bones were not affected. Bone PFOA concentrations were clearly elevated even at the age of 17months. In osteoblasts, low concentrations of PFOA increased osteocalcin (OCN) expression and calcium secretion, but at PFOA concentrations of 100µM and above osteocalcin (OCN) expression and calcium secretion were decreased. The number of osteoclasts was increased at all PFOA concentrations tested and resorption activity dose-dependently increased from 0.1-1.0µM, but decreased at higher concentrations. The results show that PFOA accumulates in bone and is present in bones until the old age. PFOA has the potential to influence bone turnover over a long period of time. Therefore bone is a target tissue for PFOA, and altered bone geometry and mineral density seem to persist throughout the life of the animal.