Clinical Super-Resolution Computed Tomography of Bone Microstructure: Application in Musculoskeletal and Dental Imaging

PURPOSE

Clinical cone-beam computed tomography (CBCT) devices are limited to imaging features of half a millimeter in size and cannot quantify the tissue microstructure. We demonstrate a robust deep-learning method for enhancing clinical CT images, only requiring a limited set of easy-to-acquire training data.

METHODS

Knee tissue from five cadavers and six total knee replacement patients, and 14 teeth from eight patients were scanned using laboratory CT as training data for the developed super-resolution (SR) technique. The method was benchmarked against ex vivo test set, 52 osteochondral samples are imaged with clinical and laboratory CT. A quality assurance phantom was imaged with clinical CT to quantify the technical image quality. To visually assess the clinical image quality, musculoskeletal and maxillofacial CBCT studies were enhanced with SR and contrasted to interpolated images. A dental radiologist and surgeon reviewed the maxillofacial images.

RESULTS

The SR models predicted the bone morphological parameters on the ex vivo test set more accurately than conventional image processing. The phantom analysis confirmed higher spatial resolution on the SR images than interpolation, but image grayscales were modified. Musculoskeletal and maxillofacial CBCT images showed more details on SR than interpolation; however, artifacts were observed near the crown of the teeth. The readers assessed mediocre overall scores for both SR and interpolation. The source code and pretrained networks are publicly available.

CONCLUSION

Model training with laboratory modalities could push the resolution limit beyond state-of-the-art clinical musculoskeletal and dental CBCT. A larger maxillofacial training dataset is recommended for dental applications.

Site-specific elastic and viscoelastic biomechanical properties of healthy and osteoarthritic human knee joint articular cartilage

Articular cartilage exhibits site-specific biomechanical properties. However, no study has comprehensively characterized site-specific cartilage properties from the same knee joints at different stages of osteoarthritis (OA). Cylindrical osteochondral explants (n = 381) were harvested from donor-matched lateral and medial tibia, lateral and medial femur, patella, and trochlea of cadaveric knees (N = 17). Indentation test was used to measure the elastic and viscoelastic mechanical properties of the samples, and Osteoarthritis Research Society International (OARSI) grading system was used to categorize the samples into normal (OARSI 0-1), early OA (OARSI 2-3), and advanced OA (OARSI 4-5) groups. OA-related changes in cartilage mechanical properties were site-specific. In the lateral and medial tibia and trochlea sites, equilibrium, instantaneous and dynamic moduli were higher (p < 0.001) in normal tissue than in early and advanced OA tissue. In lateral and medial femur, equilibrium, instantaneous and dynamic moduli were smaller in advanced OA, but not in early OA, than in normal tissue. The phase difference (0.1-0.25 Hz) between stress and strain was significantly smaller (p < 0.05) in advanced OA than in normal tissue across all sites except medial tibia. Our results indicated that in contrast to femoral and patellar cartilage, equilibrium, instantaneous and dynamic moduli of the tibia and trochlear cartilage decreased in early OA. These may suggest that the tibia and trochlear cartilage degrades faster than the femoral and patellar cartilage. The information is relevant for developing site-specific computational models and engineered cartilage constructs.

Sensitivity analysis of paediatric knee kinematics to the graft surgical parameters during anterior cruciate ligament reconstruction: A sequentially linked neuromusculoskeletal-finite element analysis

BACKGROUND AND OBJECTIVE

Incidence of paediatric anterior cruciate ligament (ACL) rupture has increased substantially over recent decades. Following ACL rupture, ACL reconstruction (ACLR) surgery is typically performed to restore passive knee stability. This surgery involves replacing the failed ACL with a graft, however, surgeons must select from range of surgical parameters (e.g., type, size, insertion, and pre-tension) with no robust evidence guiding these decisions. This study presents a systemmatic computational approach to study effects of surgical parameter variation on kinematics of paediatric knees.

METHODS

This study used sequentially-linked neuromusculoskeletal (NMSK) finite element (FE) models of three paediatric knees to estimate the: (i) sensitivity of post-operative knee kinematics to four surgical parameters (type, size, insertion, and pre-tension) through multi-input multi-output sensitivity analysis; (ii) influence of motion and loading conditions throughout stance phase of walking gait on sensitivity indices; and (iii) influence of subject-specific anatomy (i.e., knee size) on sensitivivty indices. A previously validated FE model of the intact knee for each subject served as a reference against which ACLR knee kinematics were compared.

RESULTS

Sensitivity analyses revealed significant influences of surgical parameters on ACLR knee kinematics, albeit without discernible trend favouring any one parameter. Graft size and pre-tension were primary drivers of variation in knee translations and rotations, however, their effects fluctuated across stance indicating motion and loading conditions affect system sensitivity to surgical parameters. Importantly, the sensitivity of knee kinematics to surgical parameter varied across subjects, indicating geometry (i.e., knee size) influenced system sensitivity. Notably, alterations in graft parameters yielded substantial effects on kinematics (normalized root-mean-square-error > 10 %) compared to intact knee models, indicating surgical parameters vary post-operative knee kinematics.

CONCLUSIONS

Overall, this initial study highlights the importance of surgical parameter selection on post-operative kinematics in the paediatric ACLR knee, and provides evidence of the need for personalized surgical planning to ultimately enhance patient outcomes.

Muscle Morphology Does Not Solely Determine Knee Flexion Weakness After Anterior Cruciate Ligament Reconstruction with a Semitendinosus Tendon Graft: A Combined Experimental and Computational Modeling Study

The distal semitendinosus tendon is commonly harvested for anterior cruciate ligament reconstruction, inducing substantial morbidity at the knee. The aim of this study was to probe how morphological changes of the semitendinosus muscle after harvest of its distal tendon for anterior cruciate ligament reconstruction affects knee flexion strength and whether the knee flexor synergists can compensate for the knee flexion weakness. Ten participants 8-18 months after anterior cruciate ligament reconstruction with an ipsilateral distal semitendinosus tendon autograft performed isometric knee flexion strength testing (15°, 45°, 60°, and 90°; 0° = knee extension) positioned prone on an isokinetic dynamometer. Morphological parameters extracted from magnetic resonance images were used to inform a musculoskeletal model. Knee flexion moments estimated by the model were then compared with those measured experimentally at each knee angle position. A statistically significant between-leg difference in experimentally-measured maximal isometric strength was found at 60° and 90°, but not 15° or 45°, of knee flexion. The musculoskeletal model matched the between-leg differences observed in experimental knee flexion moments at 15° and 45° but did not well estimate between-leg differences with a more flexed knee, particularly at 90°. Further, the knee flexor synergists could not physiologically compensate for weakness in deep knee flexion. These results suggest additional factors other than knee flexor muscle morphology play a role in knee flexion weakness following anterior cruciate ligament reconstruction with a distal semitendinosus tendon graft and thus more work at neural and microscopic levels is required for informing treatment and rehabilitation in this demographic.

Independent and combined effects of obesity and traumatic joint injury to the structure and composition of rat knee cartilage

Osteoarthritis (OA) is a multifactorial joint disease characterized by articular cartilage degradation. Risk factors for OA include joint trauma, obesity, and inflammation, each of which can affect joint health independently, but their interaction and the associated consequences of such interaction were largely unexplored. Here, we studied compositional and structural alterations in knee joint cartilages of Sprague-Dawley rats exposed to two OA risk factors: joint injury and diet-induced obesity. Joint injury was imposed by surgical transection of anterior cruciate ligaments (ACLx), and obesity was induced by a high fat/high sucrose diet. Depth-dependent proteoglycan (PG) content and collagen structural network of cartilage were measured from histological sections collected previously in Collins et al.. (2015). We found that ACLx primarily affected the superficial cartilages. Compositionally, ACLx led to reduced PG content in lean animals, but increased PG content in obese rats. Structurally, ACLx caused disorganization of collagenous network in both lean and obese animals through increased collagen orientation in the superficial tissues and a change in the degree of fibrous alignment. However, the cartilage degradation attributed to joint injury and obesity was not necessarily additive when the two risk factors were present simultaneously, particularly for PG content and collagen orientation in the superficial tissues. Interestingly, sham surgeries caused a through-thickness disorganization of collagen network in lean and obese animals. We conclude that the interactions of multiple OA risk factors are complex and their combined effects cannot be understood by superposition principle. Further research is required to elucidate the interactive mechanism between OA subtypes.

Deep learning based detection of osteophytes in radiographs and magnetic resonance imagings of the knee using 2D and 3D morphology

In this study, we investigated the discriminative capacity of knee morphology in automatic detection of osteophytes defined by the Osteoarthritis Research Society International atlas, using X-ray and magnetic resonance imaging (MRI) data. For the X-ray analysis, we developed a deep learning (DL) based model to segment femur and tibia. In case of MRIs, we utilized previously validated segmentations of femur, tibia, corresponding cartilage tissues, and menisci. Osteophyte detection was performed using DL models in four compartments: medial femur (FM), lateral femur (FL), medial tibia (TM), and lateral tibia (TL). To analyze the confounding effects of soft tissues, we investigated their morphology in combination with bones, including bones+cartilage, bones+menisci, and all the tissues. From X-ray-based 2D morphology, the models yielded balanced accuracy of 0.73, 0.69, 0.74, and 0.74 for FM, FL, TM, TL, respectively. Using 3D bone morphology from MRI, balanced accuracy was 0.80, 0.77, 0.71, and 0.76, respectively. The performance was higher than in 2D for all the compartments except for TM, with significant improvements observed for femoral compartments. Adding menisci or cartilage morphology consistently improved balanced accuracy in TM, with the greatest improvement seen for small osteophyte. Otherwise, the models performed similarly to bones-only. Our experiments demonstrated that MRI-based models show higher detection capability than X-ray based models for identifying knee osteophytes. This study highlighted the feasibility of automated osteophyte detection from X-ray and MRI data and suggested further need for development of osteophyte assessment criteria in addition to OARSI, particularly, for early osteophytic changes.

Triple contrast computed tomography reveals site-specific biomechanical differences in the human knee joint-A proof of concept study

Cartilage and synovial fluid are challenging to observe separately in native computed tomography (CT). We report the use of triple contrast agent (bismuth nanoparticles [BiNPs], CA4+, and gadoteridol) to image and segment cartilage in cadaveric knee joints with a clinical CT scanner. We hypothesize that BiNPs will remain in synovial fluid while the CA4+ and gadoteridol will diffuse into cartilage, allowing (1) segmentation of cartilage, and (2) evaluation of cartilage biomechanical properties based on contrast agent concentrations. To investigate these hypotheses, triple contrast agent was injected into both knee joints of a cadaver (N = 1), imaged with a clinical CT at multiple timepoints during the contrast agent diffusion. Knee joints were extracted, imaged with micro-CT (µCT), and biomechanical properties of the cartilage surface were determined by stress-relaxation mapping. Cartilage was segmented and contrast agent concentrations (CA4+ and gadoteridol) were compared with the biomechanical properties at multiple locations (n = 185). Spearman's correlation between cartilage thickness from clinical CT and reference µCT images verifies successful and reliable segmentation. CA4+ concentration is significantly higher in femoral than in tibial cartilage at 60 min and further timepoints, which corresponds to the higher Young's modulus observed in femoral cartilage. In this pilot study, we show that (1) large BiNPs do not diffuse into cartilage, facilitating straightforward segmentation of human knee joint cartilage in a clinical setting, and (2) CA4+ concentration in cartilage reflects the biomechanical differences between femoral and tibial cartilage. Thus, the triple contrast agent CT shows potential in cartilage morphology and condition estimation in clinical CT.

Effects of gait modifications on tissue-level knee mechanics in individuals with medial tibiofemoral osteoarthritis: A proof-of-concept study towards personalized interventions

Gait modification is a common nonsurgical approach to alter the mediolateral distribution of knee contact forces, intending to decelerate or postpone the progression of mechanically induced knee osteoarthritis (KOA). Nevertheless, the success rate of these approaches is controversial, with no studies conducted to assess alterations in tissue-level knee mechanics governing cartilage degradation response in KOA patients undertaking gait modifications. Thus, here we investigated the effect of different conventional gait conditions and modifications on tissue-level knee mechanics previously suggested as indicators of collagen network damage, cell death, and loss of proteoglycans in knee cartilage. Five participants with medial KOA were recruited and musculoskeletal finite element analyses were conducted to estimate subject-specific tissue mechanics of knee cartilages during two gait conditions (i.e., barefoot and shod) and six gait modifications (i.e., 0°, 5°, and 10° lateral wedge insoles, toe-in, toe-out, and wide stance). Based on our results, the optimal gait modification varied across the participants. Overall, toe-in, toe-out, and wide stance showed the greatest reduction in tissue mechanics within medial tibial and femoral cartilages. Gait modifications could effectually alter maximum principal stress (~20 ± 7%) and shear strain (~9 ± 4%) within the medial tibial cartilage. Nevertheless, lateral wedge insoles did not reduce joint- and tissue-level mechanics considerably. Significance: This proof-of-concept study emphasizes the importance of the personalized design of gait modifications to account for biomechanical risk factors associated with cartilage degradation.

Microscale investigation of the anisotropic swelling of cartilage tissue and cells in response to hypo-osmotic challenges

Tissue swelling represents an early sign of osteoarthritis, reflecting osmolarity changes from iso- to hypo-osmotic in the diseased joints. Increased tissue hydration may drive cell swelling. The opposing cartilages in a joint may swell differently, thereby predisposing the more swollen cartilage and cells to mechanical injuries. However, our understanding of the tissue-cell interdependence in osmotically loaded joints is limited as tissue and cell swellings have been studied separately. Here, we measured tissue and cell responses of opposing patellar (PAT) and femoral groove (FG) cartilages in lapine knees exposed to an extreme hypo-osmotic challenge. We found that the tissue matrix and most cells swelled during the hypo-osmotic challenge, but to a different extent (tissue: <3%, cells: 11%-15%). Swelling-induced tissue strains were anisotropic, showing 2%-4% stretch and 1%-2% compression along the first and third principal directions, respectively. These strains were amplified by 5-8 times in the cells. Interestingly, the first principal strains of tissue and cells occurred in different directions (60-61° for tissue vs. 8-13° for cells), suggesting different mechanisms causing volume expansion in the tissue and the cells. Instead of the continuous swelling observed in the tissue matrix, >88% of cells underwent regulatory volume decrease to return to their pre-osmotic challenge volumes. Cell shapes changed in the early phase of swelling but stayed constant thereafter. Kinematic changes to tissue and cells were larger for PAT cartilage than for FG cartilage. We conclude that the swelling-induced deformation of tissue and cells is anisotropic. Cells actively restored volume independent of the surrounding tissues and seemed to prioritize volume restoration over shape restoration. Our findings shed light on tissue-cell interdependence in changing osmotic environments that is crucial for cell mechano-transduction in swollen/diseased tissues.

3D T1 relaxation time measurements in an equine model of subtle post-traumatic osteoarthritis using MB-SWIFT

The aim of this study is to assess whether articular cartilage changes in an equine model of post-traumatic osteoarthritis (PTOA), induced by surgical creation of standard (blunt) grooves, and very subtle sharp grooves, could be detected with ex vivo T1 relaxation time mapping utilizing three-dimensional (3D) readout sequence with zero echo time. Grooves were made on the articular surfaces of the middle carpal and radiocarpal joints of nine mature Shetland ponies and osteochondral samples were harvested at 39 weeks after being euthanized under respective ethical permissions. T1 relaxation times of the samples (n = 8 + 8 for experimental and n = 12 for contralateral controls) were measured with a variable flip angle 3D multiband-sweep imaging with Fourier transform sequence. Equilibrium and instantaneous Young's moduli and proteoglycan (PG) content from OD of Safranin-O-stained histological sections were measured and utilized as reference parameters for the T1 relaxation times. T1 relaxation time was significantly (p < 0.05) increased in both groove areas, particularly in the blunt grooves, compared with control samples, with the largest changes observed in the superficial half of the cartilage. T1 relaxation times correlated weakly (Rs  ≈ 0.33) with equilibrium modulus and PG content (Rs  ≈ 0.21). T1 relaxation time in the superficial articular cartilage is sensitive to changes induced by the blunt grooves but not to the much subtler sharp grooves, at the 39-week timepoint post-injury. These findings support that T1 relaxation time has potential in detection of mild PTOA, albeit the most subtle changes could not be detected.

Prediction of Knee Joint Compartmental Loading Maxima Utilizing Simple Subject Characteristics and Neural Networks

Joint loading may affect the development of osteoarthritis, but patient-specific load estimation requires cumbersome motion laboratory equipment. This reliance could be eliminated using artificial neural networks (ANNs) to predict loading from simple input predictors. We used subject-specific musculoskeletal simulations to estimate knee joint contact forces for 290 subjects during over 5000 stance phases of walking and then extracted compartmental and total joint loading maxima from the first and second peaks of the stance phase. We then trained ANN models to predict the loading maxima from predictors that can be measured without motion laboratory equipment (subject mass, height, age, gender, knee abduction-adduction angle, and walking speed). When compared to the target data, our trained models had NRMSEs (RMSEs normalized to the mean of the response variable) between 0.14 and 0.42 and Pearson correlation coefficients between 0.42 and 0.84. The loading maxima were predicted most accurately using the models trained with all predictors. We demonstrated that prediction of knee joint loading maxima may be possible without laboratory-measured motion capture data. This is a promising step in facilitating knee joint loading predictions in simple environments, such as a physician's appointment. In future, the rapid measurement and analysis setup could be utilized to guide patients in rehabilitation to slow development of joint disorders, such as osteoarthritis.

Sensitivity of simulated knee joint mechanics to selected human and bovine fibril-reinforced poroelastic material properties

Fibril-reinforced poroviscoelastic material models are considered state-of-the-art in modeling articular cartilage biomechanics. Yet, cartilage material parameters are often based on bovine tissue properties in computational knee joint models, although bovine properties are distinctly different from those of humans. Thus, we aimed to investigate how cartilage mechanical responses are affected in the knee joint model during walking when fibril-reinforced poroviscoelastic properties of cartilage are based on human data instead of bovine. We constructed a finite element knee joint model in which tibial and femoral cartilages were modeled as fibril-reinforced poroviscoelastic material using either human or bovine data. Joint loading was based on subject-specific gait data. The resulting mechanical responses of knee cartilage were compared between the knee joint models with human or bovine fibril-reinforced poroviscoelastic cartilage properties. Furthermore, we conducted a sensitivity analysis to determine which fibril-reinforced poroviscoelastic material parameters have the greatest impact on cartilage mechanical responses in the knee joint during walking. In general, bovine cartilage properties yielded greater maximum principal stresses and fluid pressures (both up to 30%) when compared to the human cartilage properties during the loading response in both femoral and tibial cartilage sites. Cartilage mechanical responses were very sensitive to the collagen fibril-related material parameter variations during walking while they were unresponsive to proteoglycan matrix or fluid flow-related material parameter variations. Taken together, human cartilage material properties should be accounted for when the goal is to compare absolute mechanical responses of knee joint cartilage as bovine material parameters lead to substantially different cartilage mechanical responses.

Validation and evaluation of subject-specific finite element models of the pediatric knee

Finite element (FE) models have been widely used to investigate knee joint biomechanics. Most of these models have been developed to study adult knees, neglecting pediatric populations. In this study, an atlas-based approach was employed to develop subject-specific FE models of the knee for eight typically developing pediatric individuals. Initially, validation simulations were performed at four passive tibiofemoral joint (TFJ) flexion angles, and the resulting TFJ and patellofemoral joint (PFJ) kinematics were compared to corresponding patient-matched measurements derived from magnetic resonance imaging (MRI). A neuromusculoskeletal-(NMSK)-FE pipeline was then used to simulate knee biomechanics during stance phase of walking gait for each participant to evaluate model simulation of a common motor task. Validation simulations demonstrated minimal error and strong correlations between FE-predicted and MRI-measured TFJ and PFJ kinematics (ensemble average of root mean square errors < 5 mm for translations and < 4.1° for rotations). The FE-predicted kinematics were strongly correlated with published reports (ensemble average of Pearson's correlation coefficients (ρ) > 0.9 for translations and ρ > 0.8 for rotations), except for TFJ mediolateral translation and abduction/adduction rotation. For walking gait, NMSK-FE model-predicted knee kinematics, contact areas, and contact pressures were consistent with experimental reports from literature. The strong agreement between model predictions and experimental reports underscores the capability of sequentially linked NMSK-FE models to accurately predict pediatric knee kinematics, as well as complex contact pressure distributions across the TFJ articulations. These models hold promise as effective tools for parametric analyses, population-based clinical studies, and enhancing our understanding of various pediatric knee injury mechanisms. They also support intervention design and prediction of surgical outcomes in pediatric populations.

Comparison of constitutive models for meniscus and their effect on the knee joint biomechanics during gait

Mechanical behavior of meniscus can be modeled using constitutive material models of varying complexity, such as isotropic elastic or fibril reinforced poroelastic (FRPE). However, the FRPE material is complex to implement, computationally demanding in 3D geometries, and simulation is time-consuming. Hence, we aimed to quantify the most suitable and efficient constitutive model of meniscus for simulation of cartilage responses in the knee joint during walking. We showed that simpler constitutive material models can reproduce similar cartilage responses to a knee model with the FRPE meniscus, but only knee models that consider orthotropic elastic meniscus can also reproduce meniscus responses adequately.

Towards a Transferable Modeling Method of the Knee to Distinguish Between Future Healthy Joints from Osteoarthritic Joints: Data from the Osteoarthritis Initiative

Computational models can be used to predict the onset and progression of knee osteoarthritis. Ensuring the transferability of these approaches among computational frameworks is urgent for their reliability. In this work, we assessed the transferability of a template-based modeling strategy, based on the finite element (FE) method, by implementing it on two different FE softwares and comparing their results and conclusions. For that, we simulated the knee joint cartilage biomechanics of 154 knees using healthy baseline conditions and predicted the degeneration that occurred after 8 years of follow-up. For comparisons, we grouped the knees using their Kellgren-Lawrence grade at the 8-year follow-up time and the simulated volume of cartilage tissue that exceeded age-dependent thresholds of maximum principal stress. We considered the medial compartment of the knee in the FE models and used ABAQUS and FEBio FE softwares for simulations. The two FE softwares detected different volumes of overstressed tissue in corresponding knee samples (p < 0.01). However, both programs correctly distinguished between the joints that remained healthy and those that developed severe osteoarthritis after the follow-up (AUC = 0.73). These results indicate that different software implementations of a template-based modeling method similarly classify future knee osteoarthritis grades, motivating further evaluations using simpler cartilage constitutive models and additional studies on the reproducibility of these modeling strategies.

Visible and Near-Infrared Spectroscopy Enables Differentiation of Normal and Early Osteoarthritic Human Knee Joint Articular Cartilage

Osteoarthritis degenerates cartilage and impairs joint function. Early intervention opportunities are missed as current diagnostic methods are insensitive to early tissue degeneration. We investigated the capability of visible light-near-infrared spectroscopy (Vis-NIRS) to differentiate normal human cartilage from early osteoarthritic one. Vis-NIRS spectra, biomechanical properties and the state of osteoarthritis (OARSI grade) were quantified from osteochondral samples harvested from different anatomical sites of human cadaver knees. Two support vector machines (SVM) classifiers were developed based on the Vis-NIRS spectra and OARSI scores. The first classifier was designed to distinguish normal (OARSI: 0-1) from general osteoarthritic cartilage (OARSI: 2-5) to check the general suitability of the approach yielding an average accuracy of 75% (AUC = 0.77). Then, the second classifier was designed to distinguish normal from early osteoarthritic cartilage (OARSI: 2-3) yielding an average accuracy of 71% (AUC = 0.73). Important wavelength regions for differentiating normal from early osteoarthritic cartilage were related to collagen organization (wavelength region: 400-600 nm), collagen content (1000-1300 nm) and proteoglycan content (1600-1850 nm). The findings suggest that Vis-NIRS allows objective differentiation of normal and early osteoarthritic tissue, e.g., during arthroscopic repair surgeries.

A novel knee joint model in FEBio with inhomogeneous fibril-reinforced biphasic cartilage simulating tissue mechanical responses during gait: data from the osteoarthritis initiative

We developed a novel knee joint model in FEBio to simulate walking. Knee cartilage was modeled using a fibril-reinforced biphasic (FRB) formulation with depth-wise collagen architecture and split-lines to account for cartilage structure. Under axial compression, the knee model with FRB cartilage yielded contact pressures, similar to reported experimental data. Furthermore, gait analysis with FRB cartilage simulated spatial and temporal trends in cartilage fluid pressures, stresses, and strains, comparable to those of a fibril-reinforced poroviscoelastic (FRPVE) material in Abaqus. This knee joint model in FEBio could be used for further studies of knee disorders using physiologically relevant loading.

Signatures of disease progression in knee osteoarthritis: insights from an integrated multi-scale modeling approach, a proof of concept

Introduction: Knee osteoarthritis (KOA) is characterized by articular cartilage degeneration. It has been widely accepted that the mechanical joint environment plays a significant role in the onset and progression of this disease. In silico models have been used to study the interplay between mechanical loading and cartilage degeneration, hereby relying mainly on two key mechanoregulatory factors indicative of collagen degradation and proteoglycans depletion. These factors are the strain in collagen fibril direction (SFD) and maximum shear strain (MSS) respectively. Methods: In this study, a multi-scale in silico modeling approach was used based on a synergy between musculoskeletal and finite element modeling to evaluate the SFD and MSS. These strains were evaluated during gait based on subject-specific gait analysis data collected at baseline (before a 2-year follow-up) for a healthy and progressive early-stage KOA subject with similar demographics. Results: The results show that both SFD and MSS factors allowed distinguishing between a healthy subject and a KOA subject, showing progression at 2 years follow-up, at the instance of peak contact force as well as during the stance phase of the gait cycle. At the peak of the stance phase, the SFD were found to be more elevated in the KOA patient with the median being 0.82% higher in the lateral and 0.4% higher in the medial compartment of the tibial cartilage compared to the healthy subject. Similarly, for the MSS, the median strains were found to be 3.6% higher in the lateral and 0.7% higher in the medial tibial compartment of the KOA patient compared to the healthy subject. Based on these intersubject SFD and MSS differences, we were additionally able to identify that the tibial compartment of the KOA subject at risk of progression. Conclusion/discussion: We confirmed the mechanoregulatory factors as potential biomarkers to discriminate patients at risk of disease progression. Future studies should evaluate the sensitivity of the mechanoregulatory factors calculated based on this multi-scale modeling workflow in larger patient and control cohorts.

Load-induced fluid pressurisation in hydrogel systems before and after reinforcement by melt-electrowritten fibrous meshes

Fluid pressure develops transiently within mechanically-loaded, cell-embedding hydrogels, but its magnitude depends on the intrinsic material properties of the hydrogel and cannot be easily altered. The recently developed melt-electrowriting (MEW) technique enables three-dimensional printing of structured fibrous mesh with small fibre diameter (20 μm). The MEW mesh with 20 μm fibre diameter can synergistically increase the instantaneous mechanical stiffness of soft hydrogels. However, the reinforcing mechanism of the MEW meshes is not well understood, and may involve load-induced fluid pressurisation. Here, we examined the reinforcing effect of MEW meshes in three hydrogels: gelatin methacryloyl (GelMA), agarose and alginate, and the role of load-induced fluid pressurisation in the MEW reinforcement. We tested the hydrogels with and without MEW mesh (i.e., hydrogel alone, and MEW-hydrogel composite) using micro-indentation and unconfined compression, and analysed the mechanical data using biphasic Hertz and mixture models. We found that the MEW mesh altered the tension-to-compression modulus ratio differently for hydrogels that are cross-linked differently, which led to a variable change to their load-induced fluid pressurisation. MEW meshes only enhanced the fluid pressurisation for GelMA, but not for agarose or alginate. We speculate that only covalently cross-linked hydrogels (GelMA) can effectively tense the MEW meshes, thereby enhancing the fluid pressure developed during compressive loading. In conclusion, load-induced fluid pressurisation in selected hydrogels was enhanced by MEW fibrous mesh, and may be controlled by MEW mesh of different designs in the future, thereby making fluid pressure a tunable cell growth stimulus for tissue engineering involving mechanical stimulation.

Morphology of proximal and distal human semitendinosus compartments and the effects of distal tendon harvesting for anterior cruciate ligament reconstruction

The human semitendinosus muscle is characterized by a tendinous inscription separating proximal and distal neuromuscular compartments. As each compartment is innervated by separate nerve branches, potential exists for independent operation and control of compartments. However, the morphology and function of each compartment have not been thoroughly examined in an adult human population. Further, the distal semitendinosus tendon is typically harvested for use in anterior cruciate ligament reconstruction surgery, which induces long-term morphological changes to the semitendinosus muscle-tendon unit. It remains unknown if muscle morphological alterations following anterior cruciate ligament reconstruction are uniform between proximal and distal semitendinosus compartments. Here, we performed magnetic resonance imaging on 10 individuals who had undergone anterior cruciate ligament reconstruction involving an ipsilateral distal semitendinosus tendon graft 14 ± 4 months prior, extracting morphological parameters of the whole semitendinosus muscle and each individual compartment from both the (non-injured) contralateral and surgical legs. In the contralateral leg, volume and length of the proximal compartment were smaller than the distal compartment. No between-compartment differences in volume or length were found for anterior cruciate ligament reconstructed legs, likely due to greater shortening of the distal compared to the proximal compartment after anterior cruciate ligament reconstruction. The maximal anatomical cross-sectional area of both compartments was substantially smaller on the anterior cruciate ligament reconstructed leg but did not differ between compartments on either leg. The absolute and relative between-leg differences in proximal compartment morphology on the anterior cruciate ligament reconstructed leg were strongly correlated with the corresponding between-leg differences in distal compartment morphological parameters. Specifically, greater between-leg morphological differences in one compartment were highly correlated with large between-leg differences in the other compartment, and vice versa for smaller differences. These relationships indicate that despite the heterogeneity in compartment length and volume, compartment atrophy is not independent or random. Further, the tendinous inscription endpoints were generally positioned at the same proximodistal level as the compartment maximal anatomical cross-sectional areas, providing a wide area over which the tendinous inscription could mechanically interact with compartments. Overall, results suggest the two human semitendinosus compartments are not mechanically independent.

Micro-computed Tomography-Based Collagen Orientation and Anisotropy Analysis of Rabbit Articular Cartilage

The collagen network is the highly organized backbone of articular cartilage providing tissue tensile stiffness and restricting proteoglycan bleaching out of the tissue. Osteoarthritis (OA) diminishes proper collagen network adaptation. Our aim was to provide quantitative three-dimensional (3D) information of the cartilage collagen network adaptation in early osteoarthritis using high resolution micro-computed tomography (µCT)-imaging. Osteochondral samples from the femoral condyles were collected from healthy (N = 8, both legs) and experimental OA rabbit model with anterior cruciate ligament transection (N = 14, single leg). Samples were processed for cartilage µCT-imaging and histological evaluation with polarized light microscopy (PLM). Structure tensor analysis was used to analyse the collagen fibre orientation and anisotropy of the µCT-images, and PLM was used as a validation for structural changes. Depth-wise comparison of collagen fibre orientation acquired with µCT-imaging and PLM correlated well, but the values obtained with PLM were systematically greater than those measured with µCT-imaging. Structure tensor analysis allowed for 3D quantification of collagen network anisotropy. Finally, µCT-imaging revealed only minor differences between the control and experimental groups.

Adaptation of Fibril-Reinforced Poroviscoelastic Properties in Rabbit Collateral Ligaments 8 Weeks After Anterior Cruciate Ligament Transection

Ligaments of the knee provide stability and prevent excessive motions of the joint. Rupture of the anterior cruciate ligament (ACL), a common sports injury, results in an altered loading environment for other tissues in the joint, likely leading to their mechanical adaptation. In the collateral ligaments, the patterns and mechanisms of biomechanical adaptation following ACL transection (ACLT) remain unknown. We aimed to characterize the adaptation of elastic and viscoelastic properties of the lateral and medial collateral ligaments eight weeks after ACLT. Unilateral ACLT was performed in six rabbits, and collateral ligaments were harvested from transected and contralateral knee joints after eight weeks, and from an intact control group (eight knees from four animals). The cross-sectional areas were measured with micro-computed tomography. Stepwise tensile stress-relaxation testing was conducted up to 6% final strain, and the elastic and viscoelastic properties were characterized with a fibril-reinforced poroviscoelastic material model. We found that the cross-sectional area of the collateral ligaments in the ACL transected knees increased, the nonlinear elastic collagen network modulus of the LCL decreased, and the amount of fast relaxation in the MCL decreased. Our results indicate that rupture of the ACL leads to an early adaptation of the elastic and viscoelastic properties of the collagen fibrillar network in the collateral ligaments. These adaptations may be important to consider when evaluating whole knee joint mechanics after ACL rupture, and the results aid in understanding the consequences of ACL rupture on other tissues.

Correction: A musculoskeletal finite element model of rat knee joint for evaluating cartilage biomechanics during gait

[This corrects the article DOI: 10.1371/journal.pcbi.1009398.].

Dual-contrast micro-CT enables cartilage lesion detection and tissue condition evaluation ex vivo

BACKGROUND

Post-traumatic osteoarthritis is a frequent joint disease in the horse. Currently, equine medicine lacks effective methods to diagnose the severity of chondral defects after an injury.

OBJECTIVES

To investigate the capability of dual-contrast-enhanced computed tomography (dual-CECT) for detection of chondral lesions and evaluation of the severity of articular cartilage degeneration in the equine carpus ex vivo.

STUDY DESIGN

Pre-clinical experimental study.

METHODS

In nine Shetland ponies, blunt and sharp grooves were randomly created (in vivo) in the cartilage of radiocarpal and middle carpal joints. The contralateral joint served as control. The ponies were subjected to an 8-week exercise protocol and euthanised 39 weeks after surgery. CECT scanning (ex vivo) of the joints was performed using a micro-CT scanner 1 hour after an intra-articular injection of a dual-contrast agent. The dual-contrast agent consisted of ioxaglate (negatively charged, q = -1) and bismuth nanoparticles (BiNPs, q = 0, diameter ≈ 0.2 µm). CECT results were compared to histological cartilage proteoglycan content maps acquired using digital densitometry.

RESULTS

BiNPs enabled prolonged visual detection of both groove types as they are too large to diffuse into the cartilage. Furthermore, proportional ioxaglate diffusion inside the tissue allowed differentiation between the lesion and ungrooved articular cartilage (3 mm from the lesion and contralateral joint). The mean ioxaglate partition in the lesion was 19 percentage points higher (P < 0.001) when compared with the contralateral joint. The digital densitometry and the dual-contrast CECT findings showed good subjective visual agreement.

MAIN LIMITATIONS

Ex vivo study protocol and a low number of investigated joints.

CONCLUSIONS

The dual-CECT methodology, used in this study for the first time to image whole equine joints, is capable of effective lesion detection and simultaneous evaluation of the condition of the articular cartilage.

Dependence of crack shape in loaded articular cartilage on the collagenous structure

Cartilage cracks disrupt tissue mechanics, alter cell mechanobiology, and often trigger tissue degeneration. Yet, some tissue cracks heal spontaneously. A primary factor determining the fate of tissue cracks is the compression-induced mechanics, specifically whether a crack opens or closes when loaded. Crack deformation is thought to be affected by tissue structure, which can be probed by quantitative polarized light microscopy (PLM). It is unclear how the PLM measures are related to deformed crack morphology. Here, we investigated the relationship between PLM-derived cartilage structure and mechanical behavior of tissue cracks by testing if PLM-derived structural measures correlated with crack morphology in mechanically indented cartilages.

METHODS

Knee joint cartilages harvested from mature and immature animals were used for their distinct collagenous fibrous structure and composition. The cartilages were cut through thickness, indented over the cracked region, and processed histologically. Sample-specific birefringence was quantified as two-dimensional (2D) maps of azimuth and retardance, two measures related to local orientation and degree of alignment of the collagen fibers, respectively. The shape of mechanically indented tissue cracks, measured as depth-dependent crack opening, were compared with azimuth, retardance, or "PLM index," a new parameter derived by combining azimuth and retardance.

RESULTS

Of the three parameters, only the PLM index consistently correlated with the crack shape in immature and mature tissues.

CONCLUSION

In conclusion, we identified the relative roles of azimuth and retardance on the deformation of tissue cracks, with azimuth playing the dominant role. The applicability of the PLM index should be tested in future studies using naturally-occurring tissue cracks.

Injury-related cell death and proteoglycan loss in articular cartilage: Numerical model combining necrosis, reactive oxygen species, and inflammatory cytokines

Osteoarthritis (OA) is a common musculoskeletal disease that leads to deterioration of articular cartilage, joint pain, and decreased quality of life. When OA develops after a joint injury, it is designated as post-traumatic OA (PTOA). The etiology of PTOA remains poorly understood, but it is known that proteoglycan (PG) loss, cell dysfunction, and cell death in cartilage are among the first signs of the disease. These processes, influenced by biomechanical and inflammatory stimuli, disturb the normal cell-regulated balance between tissue synthesis and degeneration. Previous computational mechanobiological models have not explicitly incorporated the cell-mediated degradation mechanisms triggered by an injury that eventually can lead to tissue-level compositional changes. Here, we developed a 2-D mechanobiological finite element model to predict necrosis, apoptosis following excessive production of reactive oxygen species (ROS), and inflammatory cytokine (interleukin-1)-driven apoptosis in cartilage explant. The resulting PG loss over 30 days was simulated. Biomechanically triggered PG degeneration, associated with cell necrosis, excessive ROS production, and cell apoptosis, was predicted to be localized near a lesion, while interleukin-1 diffusion-driven PG degeneration was manifested more globally. Interestingly, the model also showed proteolytic activity and PG biosynthesis closer to the levels of healthy tissue when pro-inflammatory cytokines were rapidly inhibited or cleared from the culture medium, leading to partial recovery of PG content. The numerical predictions of cell death and PG loss were supported by previous experimental findings. Furthermore, the simulated ROS and inflammation mechanisms had longer-lasting effects (over 3 days) on the PG content than localized necrosis. The mechanobiological model presented here may serve as a numerical tool for assessing early cartilage degeneration mechanisms and the efficacy of interventions to mitigate PTOA progression.

Contribution of collagen degradation and proteoglycan depletion to cartilage degeneration in primary and secondary osteoarthritis: an in silico study

OBJECTIVES

Current experimental approaches cannot elucidate the effect of maladaptive changes on the main cartilage constituents during the degeneration process in osteoarthritis (OA). In silico approaches, however, allow creating 'virtual knock-out' cases to elucidate these effects in a constituent-specific manner. We used such an approach to study the main mechanisms of cartilage degeneration in different mechanical loadings associated with the following OA etiologies: (1) physiological loading of degenerated cartilage, (2) injurious loading of healthy intact cartilage and (3) physiological loading of cartilage with a focal defect.

METHODS

We used the recently developed Cartilage Adaptive REorientation Degeneration (CARED) framework to simulate cartilage degeneration associated with primary and secondary OA (OA cases (1)-(3)). CARED incorporates numerical description of tissue-level cartilage degeneration mechanisms in OA, namely, collagen degradation, collagen reorientation, fixed charged density loss and tissue hydration increase following mechanical loading. We created 'virtual knock-out' scenarios by deactivating these degenerative processes one at a time in each of the three OA cases.

RESULTS

In the injurious loading of intact and physiological loading of degenerated cartilage, collagen degradation drives degenerative changes through fixed charge density loss and tissue hydration rise. In contrast, the two later mechanisms were more prominent in the focal defect cartilage model.

CONCLUSION

The virtual knock-out models reveal that injurious loading to intact cartilage and physiological loading to degenerated cartilage induce initial degenerative changes in the collagen network, whereas, in the presence of a focal cartilage defect, mechanical loading initially causes proteoglycans (PG) depletion, before changes in the collagen fibril network occur.

Multiscale In Silico Modeling of Cartilage Injuries

Injurious loading of the joint can be accompanied by articular cartilage damage and trigger inflammation. However, it is not well-known which mechanism controls further cartilage degradation, ultimately leading to post-traumatic osteoarthritis. For personalized prognostics, there should also be a method that can predict tissue alterations following joint and cartilage injury. This chapter gives an overview of experimental and computational methods to characterize and predict cartilage degradation following joint injury. Two mechanisms for cartilage degradation are proposed. In (1) biomechanically driven cartilage degradation, it is assumed that excessive levels of strain or stress of the fibrillar or non-fibrillar matrix lead to proteoglycan loss or collagen damage and degradation. In (2) biochemically driven cartilage degradation, it is assumed that diffusion of inflammatory cytokines leads to degradation of the extracellular matrix. When implementing these two mechanisms in a computational in silico modeling workflow, supplemented by in vitro and in vivo experiments, it is shown that biomechanically driven cartilage degradation is concentrated on the damage environment, while inflammation via synovial fluid affects all free cartilage surfaces. It is also proposed how the presented in silico modeling methodology may be used in the future for personalized prognostics and treatment planning of patients with a joint injury.

Associations of human femoral condyle cartilage structure and composition with viscoelastic and constituent-specific material properties at different stages of osteoarthritis

The relationships between structure and function in human knee femoral cartilage are not well-known at different stages of osteoarthritis. Thus, our aim was to characterize the depth-dependent composition and structure (proteoglycan content, collagen network organization and collagen content) of normal and osteoarthritic human femoral condyle cartilage (n = 47) and relate them to their viscoelastic and constituent-specific mechanical properties that are obtained through dynamic sinusoidal testing and fibril-reinforced poroelastic material modeling of stress-relaxation testing, respectively. We characterized the proteoglycan content using digital densitometry, collagen network organization (orientation angle and anisotropy) using polarized light microscopy and collagen content using Fourier transform infrared spectroscopy. In the superficial cartilage (0-10 % of thickness), the collagen network disorganization and proteoglycan loss were associated with the smaller initial fibril network modulus - a parameter representing the pretension of the collagen network. Furthermore, the proteoglycan loss was associated with the greater strain-dependent fibril network modulus - a measure of nonlinear mechanical behavior. The proteoglycan loss was also associated with greater cartilage viscosity at a low loading frequency (0.005 Hz), while the collagen network disorganization was associated with greater cartilage viscosity at a high loading frequency (1 Hz). Our results suggest that proteoglycan loss and collagen network disorganization reduce the pretension of the collagen network while proteoglycan degradation also increases the nonlinear mechanical behavior of the collagen network. Further, the results also highlight that proteoglycan loss and collagen disorganization increase the viscosity of femoral cartilage, but their contribution to increased viscosity occurs in completely different loading frequencies.

Effect of cells on spatial quantification of proteoglycans in articular cartilage of small animals

OBJECTIVE

Histochemical characterization of proteoglycan content in articular cartilage is important for the understanding of osteoarthritis pathogenesis. However, cartilage cells may interfere with the measurement of matrix proteoglycan content in small animal models (e.g. mice and rats) due to the high cell volume fraction (38%) in mice compared to human tissue (~1%). We investigated whether excluding the cells from image analysis affects the histochemically measured proteoglycan content of rat knee joint cartilage and assessed the effectiveness of a deep learning algorithm-based tool named U-Net in cell segmentation.

DESIGN

Histological sections were stained with Safranin-O, after which optical densities were measured using digital densitometry to estimate proteoglycan content. U-Net was trained with 600 annotated Safranin-O cartilage images for exclusion of cells from the cartilage extracellular matrix. Optical densities of the ECM with and without cells were compared as a function of normalized tissue depth.

RESULTS

U-Net cell segmentation was accurate, with the measured cell area fraction following largely that of ground-truth images (average difference: 4.3%). Cell area fraction varied as a function of tissue depth and took up 8-21% of the tissue area. The exclusion of cells from the analysis led to an increase in the analyzed depth-dependent optical density of cartilage by approximately 0.6-1.8% (p < 0.01).

CONCLUSIONS

Although the effect of cells on the analyzed proteoglycan content is small, it should be considered for improved sensitivity, especially at the onset of the disease during which cells may proliferate in small animals.

The intrinsic quality of proteoglycans, but not collagen fibres, degrades in osteoarthritic cartilage

The function of articular cartilage as a load-bearing connective tissue is derived primarily from a balanced interaction between the swelling proteoglycan (PG) matrix and tension-resistant collagen fibrous network. Such balance is compromised during joint disease such as osteoarthritis (OA) due to degradation to PGs and/or collagens. While the PG degradation is generally thought to be related to a loss of protein abundance, the collagenous degradation is more complex as it can be caused independently by a decrease of collagen content, disorganisation of fibrous structure and softening of individual collagen fibrils. A comprehensive understanding of the initial trajectories of degradation of PGs and collagen network can improve our chance of finding potential therapeutic solutions for OA. Here, we developed geometrically, structurally, and compositionally realistic and sample-specific Finite Element (FE) models under the framework of multiphasic mixture theory, from which the elastic moduli of collagen fibres and the PG load-bearing quality in healthy and diseased cartilages were estimated by numerical optimisation of the multi-step indentation stress relaxation force-time curves. We found the intrinsic quality of collagen fibres, measured by their elastic moduli, to stay constant for healthy and diseased cartilages. Combining with previous findings which show unaltered collagen content during early stages of OA, our results suggest the disorganisation of collagen fibrous network as the first form of collagenous degradation in osteoarthritic cartilage. We also found that PG degradation involves not only a loss of protein abundance, but also the quality of the remaining PGs in generating sufficient osmotic pressure for load bearing. This study sheds light on the mechanism of OA pathogenesis and highlights the restoration of collageneous organisation in cartilage as key medical intervention for OA. STATEMENT OF SIGNIFICANCE: Collagen network in articular cartilage consists of individual fibres that are organised into depth-dependent structure specialised for joint load-bearing and lubrication. During osteoarthritis, the collagen network undergoes mechanical degradation, but it is unclear if a loss of content, disorganisation of fibrous structure, or softening of individual fibres causes this degeneration. Using mechanical indentation, Finite Element modelling, and numerical optimisation methods, we determined that individual fibres did not soften in early disease stage. Together with previous findings showing unaltered collagen content, our results pinpoint the disorganisation of collagen structure as the main culprit for early collagenous degradation in osteoarthritic cartilage. Thus, early restoration in cartilage of collagen organisation, instead of individual fibre quality, may be key to slow osteoarthritis development.

Deformation behaviors and mechanical impairments of tissue cracks in immature and mature cartilages

Degeneration of articular cartilage is often triggered by a small tissue crack. As cartilage structure and composition change with age, the mechanics of cracked cartilage may depend on the tissue age, but this relationship is poorly understood. Here, we investigated cartilage mechanics and crack deformation in immature and mature cartilage exposed to a full-thickness tissue crack using indentation testing and histology, respectively. When a cut was introduced, tissue cracks opened wider in the mature cartilage compared to the immature cartilage. However, the opposite occurred upon mechanical indentation over the cracked region. Functionally, the immature-cracked cartilages stress-relaxed faster, experienced increased tissue strain, and had reduced instantaneous stiffness, compared to the mature-cracked cartilages. Taken together, mature cartilage appears to withstand surface cracks and maintains its mechanical properties better than immature cartilage and these superior properties can be explained by the structure of their collagen fibrous network.

Subject-specific biomechanical analysis to estimate locations susceptible to osteoarthritis-Finite element modeling and MRI follow-up of ACL reconstructed patients

The aims of this case-control study were to: (1) Identify cartilage locations and volumes at risk of osteoarthritis (OA) using subject-specific finite element (FE) models; (2) Quantify the relationships between the simulated biomechanical parameters and T₂ and T_1ρ relaxation times of magnetic resonance imaging (MRI). We created subject-specific FE models for seven patients with anterior cruciate ligament (ACL) reconstruction and six controls based on a previous proof-of-concept study. We identified locations and cartilage volumes susceptible to OA, based on maximum principal stresses and absolute maximum shear strains in cartilage exceeding thresholds of 7 MPa and 32%, respectively. The locations and volumes susceptible to OA were compared qualitatively and quantitatively against 2-year longitudinal changes in T₂ and T_1ρ relaxation times. The degeneration volumes predicted by the FE models, based on excessive maximum principal stresses, were significantly correlated (r = 0.711, p < 0.001) with the degeneration volumes determined from T₂ relaxation times. There was also a significant correlation between the predicted stress values and changes in T₂ relaxation time (r = 0.649, p < 0.001). Absolute maximum shear strains and changes in T_1ρ relaxation time were not significantly correlated. Five out of seven patients with ACL reconstruction showed excessive maximum principal stresses in either one or both tibial cartilage compartments, in agreement with follow-up information from MRI. Expectedly, for controls, the FE models and follow-up information showed no degenerative signs. Our results suggest that the presented modelling methodology could be applied to prospectively identify ACL reconstructed patients at risk of biomechanically driven OA, particularly by the analysis of maximum principal stresses of cartilage.

Assessment of Ligament Viscoelastic Properties Using Raman Spectroscopy

Injuries to the ligaments of the knee commonly impact vulnerable and physically active individuals. These injuries can lead to the development of degenerative diseases such as post-traumatic osteoarthritis (PTOA). Non-invasive optical modalities, such as infrared and Raman spectroscopy, provide means for quantitative evaluation of knee joint tissues and have been proposed as potential quantitative diagnostic tools for arthroscopy. In this study, we evaluate Raman spectroscopy as a viable tool for estimating functional properties of collateral ligaments. Artificial trauma was induced by anterior cruciate ligament transection (ACLT) in the left or right knee joint of skeletally mature New Zealand rabbits. The corresponding contralateral (CL) samples were extracted from healthy unoperated joints along with a separate group of control (CNTRL) animals. The rabbits were sacrificed at 8 weeks after ACLT. The ligaments were then harvested and measured using Raman spectroscopy. A uniaxial tensile stress-relaxation testing protocol was adopted for determining several biomechanical properties of the samples. Partial least squares (PLS) regression models were then employed to correlate the spectral data with the biomechanical properties. Results show that the capacity of Raman spectroscopy for estimating the biomechanical properties of the ligament samples varies depending on the target property, with prediction error ranging from 15.78% for tissue cross-sectional area to 30.39% for stiffness. The hysteresis under cyclic loading at 2 Hz (RMSE = 6.22%, Normalized RMSE = 22.24%) can be accurately estimated from the Raman data which describes the viscous damping properties of the tissue. We conclude that Raman spectroscopy has the potential for non-destructively estimating ligament biomechanical properties in health and disease, thus enhancing the diagnostic value of optical arthroscopic evaluations of ligament integrity.

Shear strain and inflammation-induced fixed charge density loss in the knee joint cartilage following ACL injury and reconstruction: A computational study

Excessive tissue deformation near cartilage lesions and acute inflammation within the knee joint after anterior cruciate ligament (ACL) rupture and reconstruction surgery accelerate the loss of fixed charge density (FCD) and subsequent cartilage tissue degeneration. Here, we show how biomechanical and biochemical degradation pathways can predict FCD loss using a patient-specific finite element model of an ACL reconstructed knee joint exhibiting a chondral lesion. Biomechanical degradation was based on the excessive maximum shear strains that may result in cell apoptosis, while biochemical degradation was driven by the diffusion of pro-inflammatory cytokines. We found that the biomechanical model was able to predict substantial localized FCD loss near the lesion and on the medial areas of the lateral tibial cartilage. In turn, the biochemical model predicted FCD loss all around the lesion and at intact areas; the highest FCD loss was at the cartilage-synovial fluid-interface and decreased toward the deeper zones. Interestingly, simulating a downturn of an acute inflammatory response by reducing the cytokine concentration exponentially over time in synovial fluid led to a partial recovery of FCD content in the cartilage. Our novel numerical approach suggests that in vivo FCD loss can be estimated in injured cartilage following ACL injury and reconstruction. Our novel modeling platform can benefit the prediction of PTOA progression and the development of treatment interventions such as disease-modifying drug testing and rehabilitation strategies.

Towards Tailored Rehabilitation by Implementation of a Novel Musculoskeletal Finite Element Analysis Pipeline

Tissue-level mechanics (e.g., stress and strain) are important factors governing tissue remodeling and development of knee osteoarthritis (KOA), and hence, the success of physical rehabilitation. To date, no clinically feasible analysis toolbox has been introduced and used to inform clinical decision making with subject-specific in-depth joint mechanics of different activities. Herein, we utilized a rapid state-of-the-art electromyography-assisted musculoskeletal finite element analysis toolbox with fibril-reinforced poro(visco)elastic cartilages and menisci to investigate knee mechanics in different activities. Tissue mechanical responses, believed to govern collagen damage, cell death, and fixed charge density loss of proteoglycans, were characterized within 15 patients with KOA while various daily activities and rehabilitation exercises were performed. Results showed more inter-participant variation in joint mechanics during rehabilitation exercises compared to daily activities. Accordingly, the devised workflow may be used for designing subject-specific rehabilitation protocols. Further, results showed the potential to tailor rehabilitation exercises, or assess capacity for daily activity modifications, to optimally load knee tissue, especially when mechanically-induced cartilage degeneration and adaptation are of interest.

An EMG-assisted muscle-force driven finite element analysis pipeline to investigate joint- and tissue-level mechanical responses in functional activities: towards a rapid assessment toolbox

Joint tissue mechanics (e.g., stress and strain) are believed to have a major involvement in the onset and progression of musculoskeletal disorders, e.g., knee osteoarthritis (KOA). Accordingly, considerable efforts have been made to develop musculoskeletal finite element (MS-FE) models to estimate highly detailed tissue mechanics that predict cartilage degeneration. However, creating such models is time-consuming and requires advanced expertise. This limits these complex, yet promising MS-FE models to research applications with few participants and makes the models impractical for clinical assessments. Also, these previously developed MS-FE models have not been used to assess activities other than gait. This study introduces and verifies a semi-automated rapid state-of-the-art MS-FE modeling and simulation toolbox incorporating an electromyography- (EMG) assisted MS model and a muscle-force driven FE model of the knee with fibril-reinforced poro(visco)elastic cartilages and menisci. To showcase the usability of the pipeline, we estimated joint- and tissue-level knee mechanics in 15 KOA individuals performing different daily activities. The pipeline was verified by comparing the estimated muscle activations and joint mechanics to existing experimental data. To determine the importance of EMG-assisted MS approach, results were compared to those from the same FE models but driven by static-optimization-based MS models. The EMG-assisted MS-FE pipeline bore a closer resemblance to experiments compared to the static-optimization-based MS-FE pipeline. Importantly, the developed pipeline showed great potential as a rapid MS-FE analysis toolbox to investigate multiscale knee mechanics during different activities of individuals with KOA. The template FE model of the study is freely available here.

Near infrared spectroscopic evaluation of biochemical and crimp properties of knee joint ligaments and patellar tendon

Knee ligaments and tendons play an important role in stabilizing and controlling the motions of the knee. Injuries to the ligaments can lead to abnormal mechanical loading of the other supporting tissues (e.g., cartilage and meniscus) and even osteoarthritis. While the condition of knee ligaments can be examined during arthroscopic repair procedures, the arthroscopic evaluation suffers from subjectivity and poor repeatability. Near infrared spectroscopy (NIRS) is capable of non-destructively quantifying the composition and structure of collagen-rich connective tissues, such as articular cartilage and meniscus. Despite the similarities, NIRS-based evaluation of ligament composition has not been previously attempted. In this study, ligaments and patellar tendon of ten bovine stifle joints were measured with NIRS, followed by chemical and histological reference analysis. The relationship between the reference properties of the tissue and NIR spectra was investigated using partial least squares regression. NIRS was found to be sensitive towards the water (R2_CV = .65) and collagen (R2_CV = .57) contents, while elastin, proteoglycans, and the internal crimp structure remained undetectable. As collagen largely determines the mechanical response of ligaments, we conclude that NIRS demonstrates potential for quantitative evaluation of knee ligaments.

Site- and Zone-Dependent Changes in Proteoglycan Content and Biomechanical Properties of Bluntly and Sharply Grooved Equine Articular Cartilage

In this study, we mapped and quantified changes of proteoglycan (PG) content and biomechanical properties in articular cartilage in which either blunt or sharp grooves had been made, both close to the groove and more remote of it, and at the opposing joint surface (kissing site) in equine carpal joints. In nine adult Shetland ponies, standardized blunt and sharp grooves were surgically made in the radiocarpal and middle carpal joints of a randomly chosen front limb. The contralateral control limb was sham-operated. At 39 weeks after surgery, ponies were euthanized. In 10 regions of interest (ROIs) (six remote from the grooves and four directly around the grooves), PG content as a function of tissue-depth and distance-to-groove was estimated using digital densitometry. Biomechanical properties of the cartilage were evaluated in the six ROIs remote from the grooves. Compared to control joints, whole tissue depth PG loss was found in sites adjacent to sharp and, to a larger extent, blunt grooves. Also, superficial PG loss of the surgically untouched kissing cartilage layers was observed. Significant PG loss was observed up to 300 µm (sharp) and at 500 µm (blunt) from the groove into the surrounding tissue. Equilibrium modulus was lower in grooved cartilage than in controls. Grooves, in particular blunt grooves, gave rise to severe PG loss close to the grooved sites and to mild degeneration more remote from the grooves in both sharply and bluntly grooved cartilage and at the kissing sites, resulting in loss of mechanical strength over the 9-month period.

In vivo assessment of the passive stretching response of the bi-compartmental human semitendinosus muscle using shear wave elastography

The semitendinosus muscle contains distinct proximal and distal compartments arranged anatomically in-series but separated by a tendinous inscription, with each compartment innervated by separate nerve branches. Although extensively investigated in other mammals, compartment-specific mechanical properties within the human semitendinosus have scarcely been assessed in vivo. Experimental data obtained during muscle-tendon unit stretching (e.g., slack angle) can also be used to validate and/or improve musculoskeletal model estimates of semitendinosus muscle force. The purpose of this study was to investigate the passive stretching response of proximal and distal humans semitendinosus compartments to distal joint extension. Using two-dimensional shear wave elastography, we bilaterally obtained shear moduli of both semitendinosus compartments from 14 prone-positioned individuals at ten knee flexion angles (from 90° to 0° [full extension] at 10° intervals). Passive muscle mechanical characteristics (slack angle, slack shear modulus, and the slope of the increase in shear modulus) were determined for each semitendinosus compartment by fitting a piecewise exponential model to the shear modulus-joint angle curves. We found no differences between compartments or legs for slack angle, slack shear modulus, or the slope of the increase in shear modulus. We also found the experimentally determined slack angle occurred at ~15-80° higher knee flexion angles compared to estimates from two commonly used musculoskeletal models, depending on participant and model used. Overall, these findings demonstrate that passive shear modulus-joint angle curves do not differ between proximal and distal human semitendinosus compartments, and provide experimental data to improve semitendinosus force estimates derived from musculoskeletal models.

Early changes in osteochondral tissues in a rabbit model of post-traumatic osteoarthritis

Concurrent osteoarthritic (OA) manifestations in bone and cartilage are poorly known. To shed light on this issue, this study aims to investigate changes in subchondral bone and articular cartilage at two time points after anterior cruciate ligament transection (ACLT) in a rabbit model. 2 (N = 16) and 8 (N = 10) weeks after ACLT, the subchondral bone structure, cartilage thickness, Osteoarthritis Research Society International (OARSI) score, fixed charged density (FCD), and collagen orientation angle were analyzed. OA related changes were evaluated by comparing the ACLT to the contralateral (C-L) and control knees. Already 2 weeks after ACLT, higher trabecular number in the medial femoral condyle and femoral groove, greater OARSI score in the femoral condyles, and thinner trabeculae in the lateral tibial plateau and femoral groove were observed in ACLT compared to C-L knees. Only minor changes of cartilage collagen orientation in the femoral condyles and femoral groove and smaller FCD in the femoral condyles, medial tibial plateau, femoral groove and patella were observed. 8 weeks post-ACLT, the surgical knees had thinner subchondral plate and trabeculae, and smaller trabecular bone volume fraction in most of the knee locations. OARSI score was greater in the femoral condyle and lateral tibial plateau cartilage. FCD loss was progressive only in the femoral condyle, femoral groove, and patellar cartilage, and minor changes of cartilage collagen orientation angle were present in the femoral condyles, femoral groove, and lateral tibial plateau. We conclude that ACLT induces progressive subchondral bone loss, during which proteoglycan loss occurs followed by their partly recovery, as indicated by FCD results.

Structural, compositional, and functional effects of blunt and sharp cartilage damage on the joint: A 9-month equine groove model study

This study aimed to quantify the long-term progression of blunt and sharp cartilage defects and their effect on joint homeostasis and function of the equine carpus. In nine adult Shetland ponies, the cartilage in the radiocarpal and middle carpal joint of one front limb was grooved (blunt or sharp randomized). The ponies were subjected to an 8-week exercise protocol and euthanized at 39 weeks. Structural and compositional alterations in joint tissues were evaluated in vivo using serial radiographs, synovial biopsies, and synovial fluid samples. Joint function was monitored by quantitative gait analysis. Macroscopic, microscopic, and biomechanical evaluation of the cartilage and assessment of subchondral bone parameters were performed ex vivo. Grooved cartilage showed higher OARSI microscopy scores than the contra-lateral sham-operated controls (p < 0.0001). Blunt-grooved cartilage scored higher than sharp-grooved cartilage (p = 0.007) and fixed charge density around these grooves was lower (p = 0.006). Equilibrium and instantaneous moduli trended lower in grooved cartilage than their controls (significant for radiocarpal joints). Changes in other tissues included a threefold to sevenfold change in interleukin-6 expression in synovium from grooved joints at week 23 (p = 0.042) and an increased CPII/C2C ratio in synovial fluid extracted from blunt-grooved joints at week 35 (p = 0.010). Gait analysis outcome revealed mild, gradually increasing lameness. In conclusion, blunt and, to a lesser extent, sharp grooves in combination with a period of moderate exercise, lead to mild degeneration in equine carpal cartilage over a 9-month period, but the effect on overall joint health remains limited.

High-resolution infrared microspectroscopic characterization of cartilage cell microenvironment

The lateral resolution of infrared spectroscopy has been inadequate for accurate biochemical characterization of the cell microenvironment, a region regulating biochemical and biomechanical signals to cells. In this study, we demonstrate the capacity of a high-resolution Fourier transform infrared microspectroscopy (HR-FTIR-MS) to characterize the collagen content of this region. Specifically, we focus on the collagen content in the cartilage cell (chondrocyte) microenvironment of healthy and osteoarthritic (OA) cartilage. Human tibial cartilage samples (N = 28) were harvested from 7 cadaveric donors and graded for OA severity (healthy, early OA, advanced OA). HR-FTIR-MS was used to analyze the collagen content of the chondrocyte microenvironment of five distinct zones across the tissue depth. HR-FTIR-MS successfully showed collagen content distribution across chondrocytes and their environment. In zones 2 and 3 (10 - 50% of the tissue thickness), we observed that collagen content was smaller (P < 0.05) in early OA compared to the healthy tissue in the vicinity of cells (pericellular region). The collagen content loss was extended to the extracellular matrix in advanced OA tissue. No significant differences in the collagen content of the chondrocyte microenvironment were observed between the groups in the most superficial (0-10%) and deep zones (50-100%). HR-FTIR-MS revealed collagen loss in the early OA cartilage pericellular region before detectable changes in the extracellular matrix in advanced OA. HR-FTIR-MS-based compositional assessment enables a better understanding of OA-related changes in tissues. This technique can be used to identify new disease mechanisms enabling better intervention strategies. STATEMENT OF SIGNIFICANCE: Osteoarthritis (OA) is the most common degenerative joint disease causing pain and disability. While significant progress has been made in OA research, OA pathogenesis is still poorly understood and current OA treatments are mainly palliative. This study demonstrates that high-resolution FTIR microspectroscopy (HR-FTIR-MS) can characterize OA-induced compositional changes in the cell microenvironment (pericellular matrix) during the early disease stages before tissue changes in the extracellular matrix become apparent. This technique may further enable the identification of new OA mechanisms and improve our current understanding of OA pathogenesis, thus, enabling the development of better treatment methods.

Functional and structural properties of human patellar articular cartilage in osteoarthritis

Changes in the fibril-reinforced poroelastic (FRPE) mechanical material parameters of human patellar cartilage at different stages of osteoarthritis (OA) are not known. Further, the patellofemoral joint loading is thought to include more sliding and shear compared to other knee joint locations, thus, the relations between structural and functional changes may differ in OA. Thus, our aim was to determine the patellar cartilage FRPE properties followed by associating them with the structure and composition. Osteochondral plugs (n = 14) were harvested from the patellae of six cadavers. Then, the FRPE material properties were determined, and those properties were associated with proteoglycan content, collagen fibril orientation angle, optical retardation (fibril parallelism), and the state of OA of the samples. The initial fibril network modulus and permeability strain-dependency factor were 72% and 63% smaller in advanced OA samples when compared to early OA samples. Further, we observed a negative association between the initial fibril network modulus and optical retardation (r = -0.537, p < 0.05). We also observed positive associations between 1) the initial permeability and optical retardation (r = 0.547, p < 0.05), and 2) the initial fibril network modulus and optical density (r = 0.670, p < 0.01).These results suggest that the reduced pretension of the collagen fibrils, as shown by the reduced initial fibril network modulus, is linked with the loss of proteoglycans and cartilage swelling in human patellofemoral OA. The characterization of these changes is important to improve the representativeness of knee joint models in tissue and cell scale.

Guide to mechanical characterization of articular cartilage and hydrogel constructs based on a systematic in silico parameter sensitivity analysis

Osteoarthritis is a whole joint disease with cartilage degeneration being an important manifestation. Tissue engineering treatment is a solution for repairing cartilage defects by implantation of chondrocyte-laden hydrogel constructs within the defect. In silico models have recently been introduced to simulate and optimize the design of these constructs. These models require accurate knowledge on the mechanical properties of the hydrogel constructs and cartilage explants, which are challenging to obtain due to their anisotropic structure and time-dependent behaviour. We performed a systematic in silico parameter sensitivity analysis to find the most efficient unconfined compression testing protocols for mechanical characterization of hydrogel constructs and cartilage explants, with a minimum number of tests but maximum identifiability of the material parameters. The construct and explant were thereby modelled as porohyperelastic and fibril-reinforced poroelastic materials, respectively. Three commonly used loading regimes were simulated in Abaqus (ramp, relaxation and dynamic loading) with varying compressive strain magnitudes and rates. From these virtual experiments, the resulting material parameters were obtained for each combination using a numerical inverse identification scheme. For hydrogels, maximum sensitivity to the different material parameters was found when using a single step ramp loading (20% compression with 10%/s rate) followed by 15 min relaxation. For cartilage explants, a two-stepped ramp loading (10% compression with 10%/s rate and 10% compression with 1%/s rate), each step followed by 15 min relaxation, yielded the maximum sensitivity to the different material parameters. With these protocols, the material parameters could be retrieved with the lowest amount of uncertainty (hydrogel: < 2% and cartilage: < 6%). These specific results and the overall methodology can be used to optimize mechanical testing protocols to yield reliable material parameters for in silico models of cartilage and hydrogel constructs.