Causes of mechanically induced collagen damage in articular cartilage

Sensitivity of knee cartilage biomechanics in finite element analysis to selected Musculoskeletal models

Knee joint kinematics and kinetics analyzed by musculoskeletal (MS) modeling are often utilized in finite element (FE) models, estimating tissue-level mechanical responses. We compared knee cartilage stresses, strains, and centers of pressure of FE models driven by two widely used MS models, implemented in AnyBody and OpenSim. Minor discrepancies in the results were observed between the models. AnyBody-driven FE models showed slightly higher stresses in the medial tibial cartilage, while OpenSim-driven FE models estimated more anterior and lateral center of pressure. Recognizing these differences in the MS-FE models is important to ensure reliable analysis of cartilage mechanics and failure and simulation of rehabilitation.

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.

Charting Aging Trajectories of Knee Cartilage Thickness for Early Osteoarthritis Risk Prediction: An MRI Study from the Osteoarthritis Initiative Cohort

Knee osteoarthritis (OA), a prevalent joint disease in the U.S., poses challenges in terms of predicting of its early progression. Although high-resolution knee magnetic resonance imaging (MRI) facilitates more precise OA diagnosis, the heterogeneous and multifactorial aspects of OA pathology remain significant obstacles for prognosis. MRI-based scoring systems, while standardizing OA assessment, are both time-consuming and labor-intensive. Current AI technologies facilitate knee OA risk scoring and progression prediction, but these often focus on the symptomatic phase of OA, bypassing initial-stage OA prediction. Moreover, their reliance on complex algorithms can hinder clinical interpretation. To this end, we make this effort to construct a computationally efficient, easily-interpretable, and state-of-the-art approach aiding in the radiographic OA (rOA) auto-classification and prediction of the incidence and progression, by contrasting an individual's cartilage thickness with a similar demographic in the rOA-free cohort. To better visualize, we have developed the toolset for both prediction and local visualization. A movie demonstrating different subtypes of dynamic changes in local centile scores during rOA progression is available at https://tli3.github.io/KneeOA/. Specifically, we constructed age-BMI-dependent reference charts for knee OA cartilage thickness, based on MRI scans from 957 radiographic OA (rOA)-free individuals from the Osteoarthritis Initiative cohort. Then we extracted local and global centiles by contrasting an individual's cartilage thickness to the rOA-free cohort with a similar age and BMI. Using traditional boosting approaches with our centile-based features, we obtain rOA classification of KLG ≤ 1 versus KLG = 2 (AUC = 0.95, F1 = 0.89), KLG ≤ 1 versus KLG ≥ 2 (AUC = 0.90, F1 = 0.82) and prediction of KLG2 progression (AUC = 0.98, F1 = 0.94), rOA incidence (KLG increasing from < 2 to ≥ 2; AUC = 0.81, F1 = 0.69) and rOA initial transition (KLG from 0 to 1; AUC = 0.64, F1 = 0.65) within a future 48-month period. Such performance in classifying KLG ≥ 2 matches that of deep learning methods in recent literature. Furthermore, its clinical interpretation suggests that cartilage changes, such as thickening in lateral femoral and anterior femoral regions and thinning in lateral tibial regions, may serve as indicators for prediction of rOA incidence and early progression. Meanwhile, cartilage thickening in the posterior medial and posterior lateral femoral regions, coupled with a reduction in the central medial femoral region, may signify initial phases of rOA transition.

Computational assessment on the impact of collagen fiber orientation in cartilages on healthy and arthritic knee kinetics/kinematics

BACKGROUND

The inhomogeneous distribution of collagen fiber in cartilage can substantially influence the knee kinematics. This becomes vital for understanding the mechanical response of soft tissues, and cartilage deterioration including osteoarthritis (OA). Though the conventional computational models consider geometrical heterogeneity along with fiber reinforcements in the cartilage model as material heterogeneity, the influence of fiber orientation on knee kinetics and kinematics is not fully explored. This work examines how the collagen fiber orientation in the cartilage affects the healthy (intact knee) and arthritic knee response over multiple gait activities like running and walking.

METHODS

A 3D finite element knee joint model is used to compute the articular cartilage response during the gait cycle. A fiber-reinforced porous hyper elastic (FRPHE) material is used to model the soft tissue. A split-line pattern is used to implement the fiber orientation in femoral and tibial cartilage. Four distinct intact cartilage models and three OA models are simulated to assess the impact of the orientation of collagen fibers in a depth wise direction. The cartilage models with fibers oriented in parallel, perpendicular, and inclined to the articular surface are investigated for multiple knee kinematics and kinetics.

FINDINGS

The comparison of models with fiber orientation parallel to articulating surface for walking and running gait has the highest elastic stress and fluid pressure compared with inclined and perpendicular fiber-oriented models. Also, the maximum contact pressure is observed to be higher in the case of intact models during the walking cycle than for OA models. In contrast, the maximum contact pressure is higher during running in OA models than in intact models. Additionally, parallel-oriented models produce higher maximum stresses and fluid pressure for walking and running gait than proximal-distal-oriented models. Interestingly, during the walking cycle, the maximum contact pressure with intact models is approximately three times higher than on OA models. In contrast, the OA models exhibit higher contact pressure during the running cycle.

INTERPRETATION

Overall, the study indicates that collagen orientation is crucial for tissue responsiveness. This investigation provides insights into the development of tailored implants.

Earlier and More Severe Cartilage Degeneration Occurs After Meniscal Tears in Juvenile Rabbits Compared with Adults

OBJECTIVE

To compare the severity of cartilage degeneration after meniscal tears between juvenile and adult rabbits.

DESIGN

This study included 20 juvenile rabbits (2 weeks after birth) and 20 adult rabbits (6 months after birth). Meniscal tears were prepared in the anterior horn of medial menisci of right knees. Rabbits were sacrificed at 1, 3, 6, and 12 weeks postoperatively. Cartilage degenerations in the medial femoral condyle and medial tibial plateau were evaluated macroscopically and histologically. The semiquantitative assessment of cartilage degeneration was graded by macroscopic Outerbridge scoring system and histological Osteoarthritis Research Society International (OARSI) scoring system.

RESULTS

In juvenile rabbits, the morphologically intact cartilage and normal extracellular matrix architecture were observed at the first week postoperatively. Mild uneven cartilage surface and toluidine blue depletion in the medial femoral condyle were observed on histological assessment at 3 weeks postoperatively. The worsened cartilage deterioration demonstrating chondral fibrillation, prominent cell death, and glycosaminoglycan (GAG) release was observed at 6 and 12 weeks postoperatively. In adult rabbits, only mild cartilage degeneration was observed in the medial femoral condyle at 12 weeks postoperatively. The outcomes of Outerbridge and OARSI scores were consistent with the aforementioned findings in juvenile and adult rabbits.

CONCLUSIONS

Our study validated that earlier and more severe cartilage degenerations were observed in juvenile rabbits after meniscal tears compared with adult rabbits. Moreover, the post-tear cartilage degeneration demonstrated regional specificity corresponded to the tear position. However, caution is warranted when extrapolating results of animal models to humans.

Comparative Morphological and Morphometric Study between Medial and Lateral Menisci in Aged Male and Female Human Cadavers

Background

The meniscal cartilages are fibrous discs that are important for knee structures and have the ability to bear weight and stabilize joints. However, morphological and standard data for the meniscus are limited. Therefore, this work will compare anatomical and histological parameters of meniscal cartilages. The results will be important for the different measurements that are necessary for knee joint surgery.

Materials and Methods

A total of 24 aged cadavers (12 males and 12 females) were included. Knee joints were dissected and the menisci were excised and labeled as medial or lateral, right or left, male or female. Then, the menisci were kept in 10% formalin solution. Morphological variations of the meniscal shapes were macroscopically categorized. Different measurements, including the distance between anterior and posterior horns, outer and inner circumferences, width (breadth), and thickness, were done using a digital Vernier caliper and recorded manually.

Results

48 medial menisci (MMi) cartilages were studied, they were 54.6% crescent-shaped, 34.6% V-shaped, and 10.8% U-shaped. 48 lateral menisci (LMi) cartilages were studied, 41.6% were crescent-shaped, 56.4% were C-shaped, and 2% were disc-shaped articular cartilage. Findings included differences in their lengths and thickness.

Conclusion

The findings of this study were significant in providing new information on various morphological and morphometric parameters of the MMi and LMi in aged males and females, which are necessary to require more precise and comprehensive fundamental data that will be helpful for many specialists for better diagnostic and therapeutic approaches; aiming to restore normal joint conditions in senile people complaining of different meniscal pathologies.

Topographic Orientation of Scaffolds for Tissue Regeneration: Recent Advances in Biomaterial Design and Applications

Tissue engineering to develop alternatives for the maintenance, restoration, or enhancement of injured tissues and organs is gaining more and more attention. In tissue engineering, the scaffold used is one of the most critical elements. Its characteristics are expected to mimic the native extracellular matrix and its unique topographical structures. Recently, the topographies of scaffolds have received increasing attention, not least because different topographies, such as aligned and random, have different repair effects on various tissues. In this review, we have focused on various technologies (electrospinning, directional freeze-drying, magnetic freeze-casting, etching, and 3-D printing) to fabricate scaffolds with different topographic orientations, as well as discussed the physicochemical (mechanical properties, porosity, hydrophilicity, and degradation) and biological properties (morphology, distribution, adhesion, proliferation, and migration) of different topographies. Subsequently, we have compiled the effect of scaffold orientation on the regeneration of vessels, skin, neural tissue, bone, articular cartilage, ligaments, tendons, cardiac tissue, corneas, skeletal muscle, and smooth muscle. The compiled information in this review will facilitate the future development of optimal topographical scaffolds for the regeneration of certain tissues. In the majority of tissues, aligned scaffolds are more suitable than random scaffolds for tissue repair and regeneration. The underlying mechanism explaining the various effects of aligned and random orientation might be the differences in “contact guidance”, which stimulate certain biological responses in cells.

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.

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.

Rapid X-Ray-Based 3-D Finite Element Modeling of Medial Knee Joint Cartilage Biomechanics During Walking

Finite element (FE) modeling is becoming an increasingly popular method for analyzing knee joint mechanics and biomechanical mechanisms leading to osteoarthritis (OA). The most common and widely available imaging method for knee OA diagnostics is planar X-ray imaging, while more sophisticated imaging methods, e.g., magnetic resonance imaging (MRI) and computed tomography (CT), are seldom used. Hence, the capability to produce accurate biomechanical knee joint models directly from X-ray imaging would bring FE modeling closer to clinical use. Here, we extend our atlas-based framework by generating FE knee models from X-ray images (N = 28). Based on measured anatomical landmarks from X-ray and MRI, knee joint templates were selected from the atlas library. The cartilage stresses and strains of the X-ray-based model were then compared with the MRI-based model during the stance phase of the gait. The biomechanical responses were statistically not different between MRI- vs. X-ray-based models when the template obtained from X-ray imaging was the same as the MRI template. However, if this was not the case, the peak values of biomechanical responses were statistically different between X-ray and MRI models. The developed X-ray-based framework may pave the way for a clinically feasible approach for knee joint FE modeling.

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.

A Comprehensive Framework to Evaluate the Effects of Anterior Cruciate Ligament Injury and Reconstruction on Graft and Cartilage Status through the Analysis of MRI T2 Relaxation Time and Knee Laxity: A Pilot Study

BACKGROUND

Anterior cruciate ligament (ACL) tear represents a common orthopedic traumatic issue that often leads to an early development of osteoarthritis. To improve the diagnostic and prognostic techniques involved in the assessment of the joint after the trauma and during the healing process, the present work proposes a multi-parametric approach that aims to investigate the relationship between joint function and soft tissue status before and after ACL reconstruction.

METHODS

Thirteen consecutive patients who underwent ACL reconstruction were preliminarily enrolled in this study. Joint laxity assessment as well as magnetic resonance imaging with T2 mapping were performed in the pre-operative stage, at four and 18 months after surgery to acquire objective information to correlate knee function and soft tissue condition.

RESULTS

Correlations were found between graft and cartilage T2 signal, suggesting an interplay between these tissues within the knee joint. Moreover, graft maturation resulted in being connected to joint laxity, as underlined by the correlation between the graft T2 signal and the temporal evolution of knee function.

CONCLUSIONS

This preliminary study represents a step forward in assessing the effects of ACL graft maturation on knee biomechanics, and vice versa. The presented integrated framework underlines the possibility to quantitatively assess the impact of ACL reconstruction on trauma recovery and cartilage homeostasis. Moreover, the reported findings-despite the preliminary nature of the clinical impacts-evidence the possibility of monitoring the surgery outcomes using a multi-parametric prognostic investigation tool.

Chondral Defects Cause Kissing Lesions in a Porcine Model

OBJECTIVE

To assess the development of kissing lesions 12 months after the generation of full-thickness chondral defects.

DESIGN

Eight minipigs were randomized into 2 groups: the Φ8.5 mm full-thickness chondral defect group (8.5FT group) and the Φ6.5 mm full-thickness chondral defect group (6.5FT group). The Φ8.5 mm or Φ6.5 mm full-thickness chondral defects were prepared in the medial femoral condyle. Knee magnetic resonance imaging (MRI) was performed before sacrifice. India ink staining was performed to macroscopically assess kissing lesions. Histologic staining (hematoxylin-eosin [HE], safranin O/fast green, toluidine blue staining) and immunohistochemistry (collagen I, collagen II, collagen X, MMP-3) were performed. Microcomputed tomography analysis was completed to assess subchondral bone alterations.

RESULTS

Obvious kissing lesions were observed on the tibial plateau. Knee MRI demonstrated high cartilage signal intensity in the medial femoral condyle and opposite tibial plateau. HE staining demonstrated cartilage fibrillation and prominent cell death. The depletion of safranin O, toluidine blue staining, and collagen II was observed in the kissing lesion areas. The kissing lesion areas demonstrated increased collagen I, Collagen X, and MMP-3 expression. The 8.5FT group showed a significantly lower mean trabecular number (2.80 1/mm) than the control group (3.26 1/mm). The 6.5FT group showed a significantly increased mean trabecular thickness (0.54 mm) and a decreased mean trabecular number (2.71 1/mm) compared to the control group (0.32 mm; 3.26 1/mm).

CONCLUSIONS

Obvious kissing lesions were observed on the tibial plateau. Knee MRI demonstrated high cartilage signal The presented findings support the development of kissing lesions caused by full-thickness chondral defects.

Advances in the Mechanisms Affecting Meniscal Avascular Zone Repair and Therapies

Injuries to menisci are the most common disease among knee joint-related morbidities and cover a widespread population ranging from children and the general population to the old and athletes. Repair of the injuries in the meniscal avascular zone remains a significant challenge due to the limited intrinsic healing capacity compared to the peripheral vascularized zone. The current surgical strategies for avascular zone injuries remain insufficient to prevent the development of cartilage degeneration and the ultimate emergence of osteoarthritis (OA). Due to the drawbacks of current surgical methods, the research interest has been transferred toward facilitating meniscal avascular zone repair, where it is expected to maintain meniscal tissue integrity, prevent secondary cartilage degeneration and improve knee joint function, which is consistent with the current prevailing management idea to maintain the integrity of meniscal tissue whenever possible. Biological augmentations have emerged as an alternative to current surgical methods for meniscal avascular zone repair. However, understanding the specific biological mechanisms that affect meniscal avascular zone repair is critical for the development of novel and comprehensive biological augmentations. For this reason, this review firstly summarized the current surgical techniques, including meniscectomies and meniscal substitution. We then discuss the state-of-the-art biological mechanisms, including vascularization, inflammation, extracellular matrix degradation and cellular component that were associated with meniscal avascular zone healing and the advances in therapeutic strategies. Finally, perspectives for the future biological augmentations for meniscal avascular zone injuries will be given.

An in silico Framework of Cartilage Degeneration That Integrates Fibril Reorientation and Degradation Along With Altered Hydration and Fixed Charge Density Loss

Injurious mechanical loading of articular cartilage and associated lesions compromise the mechanical and structural integrity of joints and contribute to the onset and progression of cartilage degeneration leading to osteoarthritis (OA). Despite extensive in vitro and in vivo research, it remains unclear how the changes in cartilage composition and structure that occur during cartilage degeneration after injury, interact. Recently, in silico techniques provide a unique integrated platform to investigate the causal mechanisms by which the local mechanical environment of injured cartilage drives cartilage degeneration. Here, we introduce a novel integrated Cartilage Adaptive REorientation Degeneration (CARED) algorithm to predict the interaction between degenerative variations in main cartilage constituents, namely collagen fibril disorganization and degradation, proteoglycan (PG) loss, and change in water content. The algorithm iteratively interacts with a finite element (FE) model of a cartilage explant, with and without variable depth to full-thickness defects. In these FE models, intact and injured explants were subjected to normal (2 MPa unconfined compression in 0.1 s) and injurious mechanical loading (4 MPa unconfined compression in 0.1 s). Depending on the mechanical response of the FE model, the collagen fibril orientation and density, PG and water content were iteratively updated. In the CARED model, fixed charge density (FCD) loss and increased water content were related to decrease in PG content. Our model predictions were consistent with earlier experimental studies. In the intact explant model, minimal degenerative changes were observed under normal loading, while the injurious loading caused a reorientation of collagen fibrils toward the direction perpendicular to the surface, intense collagen degradation at the surface, and intense PG loss in the superficial and middle zones. In the injured explant models, normal loading induced intense collagen degradation, collagen reorientation, and PG depletion both on the surface and around the lesion. Our results confirm that the cartilage lesion depth is a crucial parameter affecting tissue degeneration, even under physiological loading conditions. The results suggest that potential fibril reorientation might prevent or slow down fibril degradation under conditions in which the tissue mechanical homeostasis is perturbed like the presence of defects or injurious loading.

Proof-of-concept for the detection of early osteoarthritis pathology by clinically applicable endomicroscopy and quantitative AI-supported optical biopsy

OBJECTIVE

Clinical trials for osteoarthritis (OA), the leading cause of global disability, are unable to pinpoint the early, potentially reversible disease with clinical technology. Hence, disease-modifying drug candidates cannot be tested early in the disease. To overcome this obstacle, we asked whether early OA-pathology detection is possible with current clinical technology.

METHODS

We determined the relationship between two sensitive early OA markers, atomic force microscopy (AFM)-measured human articular cartilage (AC) surface stiffness, and location-matched superficial zone chondrocyte spatial organizations (SCSOs), asking whether a significant loss of surface stiffness can be detected in early OA SCSO stages. We then tested whether current clinical technology can visualize and accurately diagnose the SCSOs using an approved probe-based confocal laser-endomicroscope and a random forest (RF) model.

RESULTS

We demonstrated a correlation between AC surface stiffness and the SCSO (r_rm = -0.91; 95%CI: -0.97, -0.73), and an extensive loss of surface stiffness specifically in those ACs with early OA-typical SCSO (95%CIs: string SCSO: 269-173 kPa, double string SCSO: 77-46 kPa). This established the SCSO as a visualizable, functionally relevant surrogate marker of early OA AC surface pathology. Moreover, SCSO-based stiffness discrimination worked well in each patient's AC. We then demonstrated feasibility of visualizing the SCSO by clinical laser-endomicroscopy and, importantly, accurate SCSO diagnosis using RF.

CONCLUSION

We present the proof-of-concept of early OA-pathology detection with available clinical technology, introducing a future-oriented, AI-supported, non-destructive quantitative optical biopsy for early disease detection. Operationalizing SCSO recognition, this approach allows testing for correlations between local tissue architectures with other experimental and clinical read-outs, but needs clinical validation and a larger sample size for defining diagnostic thresholds.

Characterization and finite element validation of transchondral strain in the human hip during static and dynamic loading

Distribution of strain through the thickness of articular cartilage, or transchondral strain, is highly dependent on the geometry of the joint involved. Excessive transchondral strain can damage the solid matrix and ultimately lead to osteoarthritis. Currently, high-resolution transchondral strain distribution is unknown in the human hip. Thus, knowledge of transchondral strain patterns is of fundamental importance to interpreting the patterns of injury that occur in prearthritic hip joints. This study had three main objectives. We sought to 1) quantify high-resolution transchondral strain in the native human hip, 2) determine differences in transchondral strain between static and dynamic loading conditions to better understand recovery and repressurization of cartilage in the hip, and 3) create finite element (FE) models of the experimental testing to validate a modeling framework for future analysis. The transchondral strain patterns found in this study provide insight on the localization of strain within cartilage of the hip. Most notably, the chondrolabral junction experienced high tensile and shear strain across all samples, which explains clinical data reporting it as the most common region of damage in cartilage of the hip. Further, the representative FE framework was able to match the experimental static results and predict the dynamic results with very good agreement. This agreement provides confidence for both experimental and computational measurement methods and demonstrates that the specific anisotropic biphasic FE framework used in this study can both describe and predict the experimental results.

Rheometer enabled study of cartilage frequency-dependent properties

Despite the well-established dependence of cartilage mechanical properties on the frequency of the applied load, most research in the field is carried out in either load-free or constant load conditions because of the complexity of the equipment required for the determination of time-dependent properties. These simpler analyses provide a limited representation of cartilage properties thus greatly reducing the impact of the information gathered hindering the understanding of the mechanisms involved in this tissue replacement, development and pathology. More complex techniques could represent better investigative methods, but their uptake in cartilage research is limited by the highly specialised training required and cost of the equipment. There is, therefore, a clear need for alternative experimental approaches to cartilage testing to be deployed in research and clinical settings using more user-friendly and financial accessible devices. Frequency dependent material properties can be determined through rheometry that is an easy to use requiring a relatively inexpensive device; we present how a commercial rheometer can be adapted to determine the viscoelastic properties of articular cartilage. Frequency-sweep tests were run at various applied normal loads on immature, mature and trypsinased (as model of osteoarthritis) cartilage samples to determine the dynamic shear moduli (G*, G′ G″) of the tissues. Moduli increased with increasing frequency and applied load; mature cartilage had generally the highest moduli and GAG depleted samples the lowest. Hydraulic permeability (K_H) was estimated from the rheological data and decreased with applied load; GAG depleted cartilage exhibited higher hydraulic permeability than either immature or mature tissues. The rheometer-based methodology developed was validated by the close comparison of the rheometer-obtained cartilage characteristics (G*, G′, G″, K_H) with results obtained with more complex testing techniques available in literature. Rheometry is relatively simpler and does not require highly capital intensive machinery and staff training is more accessible; thus the use of a rheometer would represent a cost-effective approach for the determination of frequency-dependent properties of cartilage for more comprehensive and impactful results for both healthcare professional and R&D.

12 Degrees of Freedom Muscle Force Driven Fibril-Reinforced Poroviscoelastic Finite Element Model of the Knee Joint

Accurate knowledge of the joint kinematics, kinetics, and soft tissue mechanical responses is essential in the evaluation of musculoskeletal (MS) disorders. Since in vivo measurement of these quantities requires invasive methods, musculoskeletal finite element (MSFE) models are widely used for simulations. There are, however, limitations in the current approaches. Sequentially linked MSFE models benefit from complex MS and FE models; however, MS model's outputs are independent of the FE model calculations. On the other hand, due to the computational burden, embedded (concurrent) MSFE models are limited to simple material models and cannot estimate detailed responses of the soft tissue. Thus, first we developed a MSFE model of the knee with a subject-specific MS model utilizing an embedded 12 degrees of freedom (DoFs) knee joint with elastic cartilages in which included both secondary kinematic and soft tissue deformations in the muscle force estimation (inverse dynamics). Then, a muscle-force-driven FE model with fibril-reinforced poroviscoelastic cartilages and fibril-reinforced poroelastic menisci was used in series to calculate detailed tissue mechanical responses (forward dynamics). Second, to demonstrate that our workflow improves the simulation results, outputs were compared to results from the same FE models which were driven by conventional MS models with a 1 DoF knee, with and without electromyography (EMG) assistance. The FE model driven by both the embedded and the EMG-assisted MS models estimated similar results and consistent with experiments from literature, compared to the results estimated by the FE model driven by the MS model with 1 DoF knee without EMG assistance.

Identification of locations susceptible to osteoarthritis in patients with anterior cruciate ligament reconstruction: Combining knee joint computational modelling with follow-up T1ρ and T2 imaging

BACKGROUND

Finite element modelling can be used to evaluate altered loading conditions and failure locations in knee joint tissues. One limitation of this modelling approach has been experimental comparison. The aims of this proof-of-concept study were: 1) identify areas susceptible to osteoarthritis progression in anterior cruciate ligament reconstructed patients using finite element modelling; 2) compare the identified areas against changes in T₂ and T_1ρ values between 1-year and 3-year follow-up timepoints.

METHODS

Two patient-specific finite element models of knee joints with anterior cruciate ligament reconstruction were created. The knee geometry was based on clinical magnetic resonance imaging and joint loading was obtained via motion capture. We evaluated biomechanical parameters linked with cartilage degeneration and compared the identified risk areas against T₂ and T_1ρ maps.

FINDINGS

The risk areas identified by the finite element models matched the follow-up magnetic resonance imaging findings. For Patient 1, excessive values of maximum principal stresses and shear strains were observed in the posterior side of the lateral tibial and femoral cartilage. For Patient 2, high values of maximum principal stresses and shear strains of cartilage were observed in the posterior side of the medial joint compartment. For both patients, increased T₂ and T_1ρ values between the follow-up times were observed in the same areas.

INTERPRETATION

Finite element models with patient-specific geometries and motions and relatively simple material models of tissues were able to identify areas susceptible to post-traumatic knee osteoarthritis. We suggest that the methodology presented here may be applied in large cohort studies.

Pathogenesis and clinical management of obesity-related knee osteoarthritis: Impact of mechanical loading

Obesity-related osteoarthritis (OA) is a complex, multifactorial condition that can cause significant impact on patients' quality of life. Whilst chronic inflammation, adipocytokines and metabolic factors are considered to be important pathogenic factors in obesity related OA, there has been limited investigation into the biomechanical impact of obesity on OA development. This review aims to demonstrate that mechanical factors are the major pathological cause of obesity-related OA. The effect of obesity on pathological changes to the osteochondral unit and surrounding connective tissues in OA is summarized, as well as the impact of obesity-related excessive and abnormal joint loading, concomitant joint malalignment and muscle weakness. An integrated therapeutic strategy based on this multi-factorial presentation is presented, to assist in the management of obesity related OA.

THE TRANSLATIONAL POTENTIAL OF THIS ARTICLE

Despite the high prevalence of obesity-related OA, there is no specific guideline available for obesity-related OA management. In this review, we demonstrated the pathological changes of obesity-related OA and summarized the impact of biomechanical factors by proposing a hypothetical model of obesity-related OA change. Therapeutic strategies based on adjusting abnormal mechanical effects are presented to assist in the management of obesity-related OA.

Clinical Contrast-Enhanced Computed Tomography With Semi-Automatic Segmentation Provides Feasible Input for Computational Models of the Knee Joint

Computational models can provide information on joint function and risk of tissue failure related to progression of osteoarthritis (OA). Currently, the joint geometries utilized in modeling are primarily obtained via manual segmentation, which is time-consuming and hence impractical for direct clinical application. The aim of this study was to evaluate the applicability of a previously developed semi-automatic method for segmenting tibial and femoral cartilage to serve as input geometry for finite element (FE) models. Knee joints from seven volunteers were first imaged using a clinical computed tomography (CT) with contrast enhancement and then segmented with semi-automatic and manual methods. In both segmentations, knee joint models with fibril-reinforced poroviscoelastic (FRPVE) properties were generated and the mechanical responses of articular cartilage were computed during physiologically relevant loading. The mean differences in the absolute values of maximum principal stress, maximum principal strain, and fibril strain between the models generated from semi-automatic and manual segmentations were <1 MPa, <0.72% and <0.40%, respectively. Furthermore, contact areas, contact forces, average pore pressures, and average maximum principal strains were not statistically different between the models (p >0.05). This semi-automatic method speeded up the segmentation process by over 90% and there were only negligible differences in the results provided by the models utilizing either manual or semi-automatic segmentations. Thus, the presented CT imaging-based segmentation method represents a novel tool for application in FE modeling in the clinic when a physician needs to evaluate knee joint function.

Use of Computational Modeling to Study Joint Degeneration: A Review

Osteoarthritis (OA), a degenerative joint disease, is the most common chronic condition of the joints, which cannot be prevented effectively. Computational modeling of joint degradation allows to estimate the patient-specific progression of OA, which can aid clinicians to estimate the most suitable time window for surgical intervention in osteoarthritic patients. This paper gives an overview of the different approaches used to model different aspects of joint degeneration, thereby focusing mostly on the knee joint. The paper starts by discussing how OA affects the different components of the joint and how these are accounted for in the models. Subsequently, it discusses the different modeling approaches that can be used to answer questions related to OA etiology, progression and treatment. These models are ordered based on their underlying assumptions and technologies: musculoskeletal models, Finite Element models, (gene) regulatory models, multiscale models and data-driven models (artificial intelligence/machine learning). Finally, it is concluded that in the future, efforts should be made to integrate the different modeling techniques into a more robust computational framework that should not only be efficient to predict OA progression but also easily allow a patient’s individualized risk assessment as screening tool for use in clinical practice.

EMG-Assisted Muscle Force Driven Finite Element Model of the Knee Joint with Fibril-Reinforced Poroelastic Cartilages and Menisci

Abnormal mechanical loading is essential in the onset and progression of knee osteoarthritis. Combined musculoskeletal (MS) and finite element (FE) modeling is a typical method to estimate load distribution and tissue responses in the knee joint. However, earlier combined models mostly utilize static-optimization based MS models and muscle force driven FE models typically use elastic materials for soft tissues or analyze specific time points of gait. Therefore, here we develop an electromyography-assisted muscle force driven FE model with fibril-reinforced poro(visco)elastic cartilages and menisci to analyze knee joint loading during the stance phase of gait. Moreover, since ligament pre-strains are one of the important uncertainties in joint modeling, we conducted a sensitivity analysis on the pre-strains of anterior and posterior cruciate ligaments (ACL and PCL) as well as medial and lateral collateral ligaments (MCL and LCL). The model produced kinematics and kinetics consistent with previous experimental data. Joint contact forces and contact areas were highly sensitive to ACL and PCL pre-strains, while those changed less cartilage stresses, fibril strains, and fluid pressures. The presented workflow could be used in a wide range of applications related to the aetiology of cartilage degeneration, optimization of rehabilitation exercises, and simulation of knee surgeries.

EMG-Assisted Muscle Force Driven Finite Element Model of the Knee Joint with Fibril-Reinforced Poroelastic Cartilages and Menisci

Abnormal mechanical loading is essential in the onset and progression of knee osteoarthritis. Combined musculoskeletal (MS) and finite element (FE) modeling is a typical method to estimate load distribution and tissue responses in the knee joint. However, earlier combined models mostly utilize static-optimization based MS models and muscle force driven FE models typically use elastic materials for soft tissues or analyze specific time points of gait. Therefore, here we develop an electromyography-assisted muscle force driven FE model with fibril-reinforced poro(visco)elastic cartilages and menisci to analyze knee joint loading during the stance phase of gait. Moreover, since ligament pre-strains are one of the important uncertainties in joint modeling, we conducted a sensitivity analysis on the pre-strains of anterior and posterior cruciate ligaments (ACL and PCL) as well as medial and lateral collateral ligaments (MCL and LCL). The model produced kinematics and kinetics consistent with previous experimental data. Joint contact forces and contact areas were highly sensitive to ACL and PCL pre-strains, while those changed less cartilage stresses, fibril strains, and fluid pressures. The presented workflow could be used in a wide range of applications related to the aetiology of cartilage degeneration, optimization of rehabilitation exercises, and simulation of knee surgeries.

The Effect of Articular Cartilage Focal Defect Size and Location in Whole Knee Biomechanics Models

Articular cartilage focal defects are common soft tissue injuries potentially linked to osteoarthritis (OA) development. Although several defect characteristics likely contribute to osteoarthritis, their relationship to local tissue deformation remains unclear. Using finite element models with various femoral cartilage geometries, we explore how defects change cartilage deformation and joint kinematics assuming loading representative of the maximum joint compression during the stance phase of gait. We show how defects, in combination with location-dependent cartilage mechanics, alter deformation in affected and opposing cartilages, as well as joint kinematics. Small and average sized defects increased maximum compressive strains by approximately 50% and 100%, respectively, compared to healthy cartilage. Shifts in the spatial locations of maximum compressive strains of defect containing models were also observed, resulting in loading of cartilage regions with reduced initial stiffnesses supporting the new, elevated loading environments. Simulated osteoarthritis (modeled as a global reduction in mean cartilage stiffness) did not significantly alter joint kinematics, but exacerbated tissue deformation. Femoral defects were also found to affect healthy tibial cartilage deformations. Lateral femoral defects increased tibial cartilage maximum compressive strains by 25%, while small and average sized medial defects exhibited decreases of 6% and 15%, respectively, compared to healthy cartilage. Femoral defects also affected the spatial distributions of deformation across the articular surfaces. These deviations are especially meaningful in the context of cartilage with location-dependent mechanics, leading to increases in peak contact stresses supported by the cartilage of between 11% and 34% over healthy cartilage.

Detrimental effects of long sedentary bouts on the biomechanical response of cartilage to sliding

Purpose/Aim: Epidemiological evidence suggests, contrary to popular mythos, that increased exercise/joint activity does not place articular cartilage at increased risk of disease, but instead promotes joint health. One explanation for this might be activity-induced cartilage rehydration; where joint articulation drives restoration of tissue hydration, thickness, and dependent tribomechanical outcomes (e.g., load support, stiffness, and lubricity) lost to joint loading. However, there have been no studies investigating how patterning of intermittent articulation influences the hydration and biomechanical functions of cartilage.Materials and Methods: Here we leveraged the convergent stationary contact area (cSCA) testing configuration and its unique ability to drive tribological rehydration, to elucidate how intermittency of activity affects the biomechanical functions of bovine stifle cartilage under well-controlled sliding conditions that have been designed to model a typical "day" of human joint activity.Results: For a fixed volume of "daily" activity (30 min) and sedentary time (60 min), breaking up intermittent activity into longer and less-frequent bouts (corresponding to longer continuous sedentary periods) resulted in the exposure of articular cartilage to markedly greater strains, losses of interstitial pressure, and friction coefficients.Conclusions: These results demonstrated that the regularity of ex vivo activity regimens, specifically the duration of sedentary bouts, had a dominant effect on the biomechanical functions of articular cartilage. In more practical terms, the results suggest that brief but regular movement patterns (e.g., every hour) may be biomechanically preferred to long and infrequent movement patterns (e.g., a long walk after a sedentary day) when controlling for daily activity volume (e.g., 30 min).

Chondrocyte viability is lost during high-rate impact loading by transfer of amplified strain, but not stress, to pericellular and cellular regions

OBJECTIVE

Deleterious impact loading to cartilage initiates post-traumatic osteoarthritis (OA). While cytokine and enzyme levels regulate disease progression, specific mechanical cues that elucidate cellular OA origins merit further investigation. We defined the dominant pericellular and cellular strain/stress transfer mechanisms following bulk-tissue injury associated with cell death.

METHOD

Using an in vitro model, we investigated rate-dependent loading and spatial localization of cell viability in acute indentation and time-course studies. Atomic force microscopy (AFM) and magnetic resonance imaging (MRI) confirmed depth-wise changes in cartilage micro-/macro-mechanics and structure post-indentation. To understand the transfer of loading to cartilage domains, we computationally modeled full-field strain and stress measures in interstitial matrix, pericellular and cellular regions.

RESULTS

Chondrocyte viability decreased following rapid impact (80%/s) vs slow loading (0.1%/s) or unloaded controls. Viability was lost immediately during impact within regions near the indenter-tissue contact but did not change over 7 days of tissue culture. AFM studies revealed a loss of stiffness following 80%/s loading, and MRI studies confirmed an increased tensile and shear strain, but not relaxometry. Image-based patterns of chondrocyte viability closely matched computational estimates of amplified maximum principal and shear strain in interstitial matrix, pericellular and cellular regions.

CONCLUSION

Rapid indentation worsens chondrocyte death and degrades cartilage matrix stiffness in indentation regions. Cell death at high strain rates may be driven by elevated tensile strains, but not matrix stress. Strain amplification beyond critical thresholds in the pericellular matrix and cells may define a point of origin for early damage in post-traumatic OA.

Femorotibial kinematics in dogs treated with tibial plateau leveling osteotomy for cranial cruciate ligament insufficiency: An in vivo fluoroscopic analysis during walking

OBJECTIVE

To determine the ability of tibial plateau leveling osteotomy (TPLO) to address abnormal femorotibial kinematics caused by cranial cruciate ligament (CCL) rupture during walking in dogs.

STUDY DESIGN

Prospective, clinical.

ANIMALS

Sixteen dogs (20-40 kg) with unilateral complete CCL rupture.

METHODS

Lateral view fluoroscopy was performed during treadmill walking preoperatively and 6 months after TPLO. Digital three-dimensional (3D) models of the femora and tibiae were created from computed tomographic (CT) images. Gait cycles were analyzed by using a 3D-to-2D image registration process. Craniocaudal translation, internal/external rotation, and flexion/extension of the femorotibial joint were compared between preoperative and 6-month postoperative time points for the affected stifle and 6-month postoperative unaffected contralateral (control) stifles.

RESULTS

In the overall population, CCL rupture resulted in 10 ± 2.2 mm (mean ± SD) cranial tibial translation at midstance phase, which was converted to 2.1 ± 4.3 mm caudal tibial translation after TPLO. However, five of 16 TPLO-treated stifles had 4.1 ± 0.3 mm of cranial tibial subluxation during mid-to-late stance phase, whereas 10 of 16 TPLO-treated stifles had 4.3 ± 0.4 mm of caudal tibial subluxation throughout the gait cycle. Overall, postoperative axial rotational and flexion/extension patterns were not different from control, but stifles with caudal tibial subluxation had more external tibial rotation during mid-to-late stance phase compared with stifles with cranial tibial subluxation.

CONCLUSION

TPLO mitigated abnormal femorotibial kinematics but did not restore kinematics to control values in 15 of 16 dogs during walking.

CLINICAL SIGNIFICANCE

Tibial plateau leveling osteotomy reduces cranial tibial subluxation during walking, but persistent instability is common.

Conditional deletion of E11/podoplanin in bone protects against load-induced osteoarthritis

Background

Subchondral bone (SCB) thickening is one of the earliest detectable changes in osteoarthritic joints and is considered a potential trigger for subsequent articular cartilage degeneration. In this manuscript, we examine whether disruption to the SCB osteocyte network contributes to the initiation and pathogenesis of osteoarthritis.

Methods

We examined expression patterns of the glycoprotein E11/podoplanin by immunohistochemical labelling in murine, human and canine osteoarthritis models. We also examined the effects of twice-weekly administration of Bortezomib, a proteasome inhibitor which stabilises osteocyte E11 levels, to C57/BL6 wild-type male mice (1 mg/kg/day) for 8 weeks after surgical destabilisation of the medial meniscus. By inducing osteoarthritis-like changes in the right knee joint of 12-week-old male E11 hypomorphic mice (and corresponding controls) using a post-traumatic joint loading model, we also investigated whether a bone-specific E11 deletion in mice increases joint vulnerability to osteoarthritis. Articular cartilage degradation and osteophyte formation were assessed by histology and in line with the OARSI grading system.

Results

Our studies reveal increased E11 expression in osteocytes of human and canine osteoarthritic SCB. We found that Bortezomib administration had no effect on surgically-induced osteoarthritis, potentially due to a lack of the expected stabilisation of E11 in the SCB. We also found, in concordance with our previous work, wild-type mice exhibited significant load-induced articular cartilage lesions on the lateral femoral condyle (p < 0.01) and osteophyte formation. In contrast, E11 hypomorphic mice did not develop osteophytes or any corresponding articular lesions.

Conclusions

Overall, these data suggest that an intact osteocyte network in the SCB contributes to the development of mechanically-driven osteoarthritis. Further, the data presented here indicate that the molecular pathways that preserve the osteocyte network, such as those driven by E11, may be targeted to limit osteoarthritis pathogenesis.

Electronic supplementary material

The online version of this article (10.1186/s12891-019-2731-9) contains supplementary material, which is available to authorized users.

Computational evaluation of altered biomechanics related to articular cartilage lesions observed in vivo

Chondral lesions provide a potential risk factor for development of osteoarthritis. Despite the variety of in vitro studies on lesion degeneration, in vivo studies that evaluate relation between lesion characteristics and the risk for the possible progression of OA are lacking. Here, we aimed to characterize different lesions and quantify biomechanical responses experienced by surrounding cartilage tissue. We generated computational knee joint models with nine chondral injuries based on clinical in vivo arthrographic computed tomography images. Finite element models with fibril-reinforced poro(visco)elastic cartilage and menisci were constructed to simulate physiological loading. Systematically, the lesions experienced increased peak values of maximum principal strain, maximum shear strain, and minimum principal strain in the surrounding chondral tissue (p < 0.01) compared with intact tissue. Depth, volume, and area of the lesion correlated with the maximum shear strain (p < 0.05, Spearman rank correlation coefficient ρ = 0.733-0.917). Depth and volume of the lesion correlated also with the maximum principal strain (p < 0.05, ρ = 0.767, and ρ = 0.717, respectively). However, the lesion area had non-significant correlation with this strain parameter (p = 0.06, ρ = 0.65). Potentially, the introduced approach could be developed for clinical evaluation of biomechanical risks of a chondral lesion and planning an intervention. Statement of Clinical Relevance: In this study, we computationally characterized different in vivo chondral lesions and evaluated their risk of cartilage degeneration. This information is vital in decision-making for intervention in order to prevent post-traumatic osteoarthritis. © 2019 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res.

Microplate assay for denatured collagen using collagen hybridizing peptides

The purpose of this study was to develop a microplate assay for quantifying denatured collagen by measuring the fluorescence of carboxyfluorescein bound collagen hybridizing peptides (F-CHP). We have shown that F-CHP binds selectively with denatured collagen, and that mechanical overload of tendon fascicles causes collagen denaturation. Proteinase K was used to homogenize tissue samples after F-CHP staining, allowing fluorescence measurement using a microplate reader. We compared our new assay to our previous image analysis method and the trypsin-hydroxyproline assay, which is the only other available method to directly quantify denatured collagen. Relative quantification of denatured collagen was performed in rat tail tendon fascicles subjected to incremental tensile overload, and normal and ostoeoarthritic guinea pig cartilage. In addition, the absolute amount of denatured collagen was determined in rat tail tendon by correlating F-CHP fluorescence with percent denatured collagen as determined by the trypsin-hydroxyproline assay. Rat tail tendon fascicles stretched to low strains (<7.5%) exhibited minimal denatured collagen, but values rapidly increased at medium strains (7.5-10.5%) and plateaued at high strains (≥12%). Osteoarthritic cartilage had higher F-CHP fluorescence than healthy cartilage. Both of these outcomes are consistent with previous studies. With the calibration curve, the microplate assay was able to absolutely quantify denatured collagen in mechanically damaged rat tail tendon fascicles as reliably as the trypsin-hydroxyproline assay. Further, we achieved these results more efficiently than current methods in a rapid, high-throughput manner, with multiple types of collagenous tissue while maintaining accuracy. © 2018 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 37:431-438, 2019.

Maximum shear strain-based algorithm can predict proteoglycan loss in damaged articular cartilage

Post-traumatic osteoarthritis (PTOA) is a common disease, where the mechanical integrity of articular cartilage is compromised. PTOA can be a result of chondral defects formed due to injurious loading. One of the first changes around defects is proteoglycan depletion. Since there are no methods to restore injured cartilage fully back to its healthy state, preventing the onset and progression of the disease is advisable. However, this is problematic if the disease progression cannot be predicted. Thus, we developed an algorithm to predict proteoglycan loss of injured cartilage by decreasing the fixed charge density (FCD) concentration. We tested several mechanisms based on the local strains or stresses in the tissue for the FCD loss. By choosing the degeneration threshold suggested for inducing chondrocyte apoptosis and cartilage matrix damage, the algorithm driven by the maximum shear strain showed the most substantial FCD losses around the lesion. This is consistent with experimental findings in the literature. We also observed that by using coordinate system-independent strain measures and selecting the degeneration threshold in an ad hoc manner, all the resulting FCD distributions would appear qualitatively similar, i.e., the greatest FCD losses are found at the tissue adjacent to the lesion. The proposed strain-based FCD degeneration algorithm shows a great potential for predicting the progression of PTOA via biomechanical stimuli. This could allow identification of high-risk defects with an increased risk of PTOA progression.

Articular cartilage gene expression patterns in the tissue surrounding the impact site following applications of shear and axial loads

BACKGROUND

Osteoarthritis is a degradative joint disease found in humans and commercial swine which can develop from a number of factors, including prior joint trauma. An impact injury model was developed to deliver in vitro loads to disease-free porcine patellae in a model of OA.

METHODS

Axial impactions (2000 N normal) and shear impactions (500 N normal with induced shear forces) were delivered to 48 randomly assigned patellae. The patellae were then cultured for 0, 3, 7, or 14 days following the impact. Specimens in the tissue surrounding the loading site were harvested and expression of 18 OA related genes was studied via quantitative PCR. The selected genes were previously identified from published work and fell into four categories: cartilage matrix, degradative enzymes, inflammatory response, and apoptosis.

RESULTS

Type II collagen (Col2a1) showed significantly lower expression in shear vs. axial adjacent tissue at day 0 and 7 (fold changes of 0.40 & 0.19, respectively). In addition, higher expression of degradative enzymes and Fas, an apoptosis gene, was observed in the shear specimens.

CONCLUSIONS

The results suggest that a more physiologically valid shear load may induce more damage to surrounding articular cartilage than a normal load alone.

Osteochondral resurfacing implantation angle is more important than implant material stiffness

Osteochondral resurfacing implants are a promising treatment for focal cartilage defects. Several implant-factors may affect the clinical outcome of this treatment, such as the implant material stiffness and the accuracy of implant placement, known to be challenging. In general, softer implants are expected to be more accommodating for implant misalignment than stiffer implants, and motion is expected to increase effects from implant misalignment and stiffness. 3D finite element models of cartilage/cartilage contact were employed in which implantation angle (0°, 5°, 10°) and implant material stiffness (E = 5 MPa, 100 MPa, 2 GPa) were varied. A creep loading (0.6 MPa) was simulated, followed by a sliding motion. Creep loading resulted in low maximum collagen strains of 2.5% in the intact case compared to 11.7% with an empty defect. Implants mostly positively affected collagen strains, deviatoric strains, and hydrostatic pressures in the adjacent cartilage, but these effects were superior for correct alignment (0°). The main effect of implant misalignment was bulging of opposing cartilage tissue into the gap caused by the misalignment. This increased collagen strains and hydrostatic pressures. Deviatoric strains were increased adjacent to the gap. Subsequent sliding initially increased strains for a stiff, misaligned implant, but generally sliding decreased strains. In conclusion, implants can decrease the detrimental effect of defects, but correct implant alignment is crucial, more than implant material stiffness. Implant misalignment causes a gap, causing potentially damaging cartilage deformation during prolonged loading, for example, standing, even for soft implants. Mild motion may positively affect the cartilage. © 2018 The Authors. Journal of Orthopaedic Research® published by Wiley Periodicals, Inc. on behalf of Orthopaedic Research Society. J Orthop Res 36:2911-2922, 2018.

A novel mechanobiological model can predict how physiologically relevant dynamic loading causes proteoglycan loss in mechanically injured articular cartilage

Cartilage provides low-friction properties and plays an essential role in diarthrodial joints. A hydrated ground substance composed mainly of proteoglycans (PGs) and a fibrillar collagen network are the main constituents of cartilage. Unfortunately, traumatic joint loading can destroy this complex structure and produce lesions in tissue, leading later to changes in tissue composition and, ultimately, to post-traumatic osteoarthritis (PTOA). Consequently, the fixed charge density (FCD) of PGs may decrease near the lesion. However, the underlying mechanisms leading to these tissue changes are unknown. Here, knee cartilage disks from bovine calves were injuriously compressed, followed by a physiologically relevant dynamic compression for twelve days. FCD content at different follow-up time points was assessed using digital densitometry. A novel cartilage degeneration model was developed by implementing deviatoric and maximum shear strain, as well as fluid velocity controlled algorithms to simulate the FCD loss as a function of time. Predicted loss of FCD was quite uniform around the cartilage lesions when the degeneration algorithm was driven by the fluid velocity, while the deviatoric and shear strain driven mechanisms exhibited slightly discontinuous FCD loss around cracks. Our degeneration algorithm predictions fitted well with the FCD content measured from the experiments. The developed model could subsequently be applied for prediction of FCD depletion around different cartilage lesions and for suggesting optimal rehabilitation protocols.

Collagen Damage Location in Articular Cartilage Differs if Damage is Caused by Excessive Loading Magnitude or Rate

Collagen damage in articular cartilage is considered nearly irreversible and may be an early indication of cartilage degeneration. Surface fibrillation and internal collagen damage may both develop after overloading. This study hypothesizes that damage develops at these different locations, because the distribution of excessive strains varies with loading rate as a consequence of time-dependent cartilage properties. The objective is to explore whether collagen damage could preferentially occur superficially or internally, depending on the magnitude and rate of overloading. Bovine osteochondral plugs were compressed with a 2 mm diameter indenter to 15, 25, 35 and 45 N, and at 5, 60 and 120 mm/min. Surface fibrillation and internal collagen damage were graded by four observers, based on histology and staining of collagen damage. Results show that loading magnitude affects the degree of collagen damage, while loading rate dominates the location of network damage: low rates predominantly damage superficial collagen, while at high rates, internal collagen damage occurs. The proposed explanation for the rate-dependent location is that internal fluid flows govern the time-dependent internal tissue deformation and therewith the location of overstained and damaged areas. This supports the hypothesis that collagen damage development is influenced by the time-dependent material behaviour of cartilage.

The critical size of focal articular cartilage defects is associated with strains in the collagen fibers

The size of full-thickness focal cartilage defect is accepted to be predictive of its fate, but at which size threshold treatment is required is unclear. Clarification of the mechanism behind this threshold effect will help determining when treatment is required. The objective was to investigate the effect of defect size on strains in the collagen fibers and the non-fibrillar matrix of surrounding cartilage. These strains may indicate matrix disruption. Tissue deformation into the defect was expected, stretching adjacent superficial collagen fibers, while an osteochondral implant was expected to prevent these deformations. Finite element simulations of cartilage/cartilage contact for intact, 0.5 to 8mm wide defects and 8mm implant cases were performed. Impact, a load increase to 2MPa in 1ms, and creep loading, a constant load of 0.5MPa for 900s, scenarios were simulated. A composition-based material model for articular cartilage was employed. Impact loading caused low strain levels for all models. Creep loading increased deviatoric strains and collagen strains in the surrounding cartilage. Deviatoric strains increased gradually with defect size, but the surface area at which collagen fiber strains exceeded failure thresholds, abruptly increased for small increases of defect size. This was caused by a narrow distribution of collagen fiber strains resulting from the non-linear stiffness of the fibers. We postulate this might be the mechanism behind the existence of a critical defect size. Filling of the defect with an implant reduced deviatoric and collagen fiber strains towards values for intact cartilage.

Development of a Cartilage Shear-Damage Model to Investigate the Impact of Surface Injury on Chondrocytes and Extracellular Matrix Wear

Background Many i n vitro damage models investigate progression of cartilage degradation after a supraphysiologic, compressive impact at the surface and do not model shear-induced damage processes. Models also neglect the response to uninterrupted tribological stress after damage. It was hypothesized that shear-induced removal of the superficial zone would accelerate matrix degradation when damage was followed by continued load and articulation. Methods Bovine cartilage underwent a 5-day test. Shear-damaged samples experienced 2 days of damage induction with articulation against polyethylene and then continued articulation against cartilage (CoC), articulation against metal (MoC), or rest as free-swelling control (FSC). Surface-intact samples were randomized to CoC, MoC, or FSC for the entire 5-day test. Samples were evaluated for chondrocyte viability, GAG (glycosaminoglycan) release (matrix wear surrogate), and histological integrity. Results Shear induction wore away the superficial zone. Damaged samples began continued articulation with collagen matrix disruption and increased cell death compared to intact samples. In spite of the damaged surface, these samples did not exhibit higher GAG release than intact samples articulating against the same counterface ( P = 0.782), contrary to our hypothesis. Differences in GAG release were found to be due to tribological testing against metal ( P = 0.003). Conclusion Shear-induced damage lowers chondrocyte viability and affects extracellular matrix integrity. Continued motion of either cartilage or metal against damaged surfaces did not increase wear compared with intact samples. We conjecture that favorable reorganization of the surface collagen fibers during articulation protected the underlying matrix. This finding suggests a potential window for clinical interventions to slow matrix degradation after traumatic incidents.

Experimental Study on the Mechanical Properties of Porcine Cartilage with Microdefect under Rolling Load

Objectives

To investigate the mechanical responses of microdefect articular cartilage under rolling load and find out the failure rule.

Methods

Rolling load was applied to the porcine articular cartilage samples with rectangular notches of different depths. The displacement and strain near the notches were obtained by the noncontact digital image correlation technique.

Results

The strain value and peak frequency around the notch increased; the maximum equivalent strain value could be observed at both bottom corners of the notch; the equivalent strain value first increased and then decreased at the points in the superficial and middle layers with the increase of rolling velocity; the points in the deep layer were less affected by rolling velocity; the equivalent strain value of the points in the superficial layer declined after rising with the increase of defect depth, while a decreased trend could be found for the points in the middle and deep layers.

Conclusions

The shear strain, which rose with the increase in defect depth, was the main factor in cartilage destruction. The cartilage tended to be destructed firstly at the bottom corner of the defect. Rolling velocity showed significant effects on superficial and middle layers. Cartilage had the ability to resist destruction.

The initial repair response of articular cartilage after mechanically induced damage

The regenerative potential of articular cartilage (AC) defects is limited and depends on defect size, biomechanical conditions, and age. Early events after overloading might be predictive for cartilage degeneration in the long term. Therefore, the present aim is to investigate the temporal response of cartilage to overloading at cell, matrix, and tissue level during the first period after mechanical overloading. In the present study, the effect of high loading (∼8 MPa) at a high rate (∼14 MPa/s) at day 0 during a 9 day period on collagen damage, gene expression, cell death, and biochemical composition in AC was investigated. A model system was developed which enabled culturing osteochondral explants after loading. Proteoglycan content was repeatedly monitored over time using μCT, whereas other evaluations required destructive measurements. Changes in matrix related gene expressions indicated a degenerative response during the first 6 h after loading. After 24 h, this was restored and data suggested an initial repair response. Cell death and microscopic damage increased after 24 h following loading. These degradative changes were not restored within the 9 day culture period, and were accompanied by a slight loss of proteoglycans at the articular surface that extended into the middle zones. The combined findings indicate that high magnitude loading of articular cartilage at a high rate induces an initial damage that later initiates a healing response that can probably not be retained due to loss of cell viability. Consequently, the matrix cannot be restored in the short term. © 2016 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 35:1265-1273, 2017.

On mechanical mechanism of damage evolution in articular cartilage

Superficial lesions of cartilage are the direct indication of osteoarthritis. To investigate the mechanical mechanism of cartilage with micro-defect under external loading, a new plain strain numerical model with micro-defect was proposed and damage evolution progression in cartilage over time has been simulated, the parameter were studied including load style, velocity of load and degree of damage. The new model consists of the hierarchical structure of cartilage and depth-dependent arched fibers. The numerical results have shown that not only damage of the cartilage altered the distribution of the stress but also matrix and fiber had distinct roles in affecting cartilage damage, and damage in either matrix or fiber could promote each other. It has been found that the superficial cracks in cartilage spread preferentially along the tangent direction of the fibers. It is the arched distribution form of fibers that affects the crack spread of cartilage, which has been verified by experiment. During the process of damage evolution, its extension direction and velocity varied constantly with the damage degree. The rolling load could cause larger stress and strain than sliding load. Strain values of the matrix initially increased and then decreased gradually with the increase of velocity, and velocity had a greater effect on matrix than fibers. Damage increased steadily before reaching 50%, sharply within 50 to 85%, and smoothly and slowly after 85%. The finding of the paper may help to understand the mechanical mechanism why the cracks in cartilage spread preferentially along the tangent direction of the fibers.

Genipin crosslinking decreases the mechanical wear and biochemical degradation of impacted cartilage in vitro

High energy trauma to cartilage causes surface fissures and microstructural damage, but the degree to which this damage renders the tissue more susceptible to wear and contributes to the progression of post-traumatic osteoarthritis (PTOA) is unknown. Additionally, no treatments are currently available to strengthen cartilage after joint trauma and to protect the tissue from subsequent degradation and wear. The purposes of this study were to investigate the role of mechanical damage in the degradation and wear of cartilage, to evaluate the effects of impact and subsequent genipin crosslinking on the changes in the viscoelastic parameters of articular cartilage, and to test the hypothesis that genipin crosslinking is an effective treatment to enhance the resistance to biochemical degradation and mechanical wear. Results demonstrate that cartilage stiffness decreases after impact loading, likely due to the formation of fissures and microarchitectural damage, and is partially or fully restored by crosslinking. The wear resistance of impacted articular cartilage was diminished compared to undamaged cartilage, suggesting that mechanical damage that is directly induced by the impact may contribute to the progression of PTOA. However, the decrease in wear resistance was completely reversed by the crosslinking treatments. Additionally, the crosslinking treatments improved the resistance to collagenase digestion at the impact-damaged articular surface. These results highlight the potential therapeutic value of collagen crosslinking via genipin in the prevention of cartilage degeneration after traumatic injury. © 2016 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 35:558-565, 2017.

Quantitative Evaluation of the Mechanical Risks Caused by Focal Cartilage Defects in the Knee

Focal cartilage lesions can proceed to severe osteoarthritis or remain unaltered even for years. A method to identify high risk defects would be of utmost importance to guide clinical decision making and to identify the patients that are at the highest risk for the onset and progression of osteoarthritis. Based on cone beam computed tomography arthrography, we present a novel computational model for evaluating changes in local mechanical responses around cartilage defects. Our model, based on data obtained from a human knee in vivo, demonstrated that the most substantial alterations around the defect, as compared to the intact tissue, were observed in minimum principal (compressive) strains and shear strains. Both strain values experienced up to 3-fold increase, exceeding levels previously associated with chondrocyte apoptosis and failure of collagen crosslinks. Furthermore, defects at the central regions of medial tibial cartilage with direct cartilage-cartilage contact were the most vulnerable to loading. Also locations under the meniscus experienced substantially increased minimum principal strains. We suggest that during knee joint loading particularly minimum principal and shear strains are increased above tissue failure limits around cartilage defects which might lead to osteoarthritis. However, this increase in strains is highly location-specific on the joint surface.

The effect of loading rate on the development of early damage in articular cartilage

Experimental reports suggest that cartilage damage depends on strain magnitude. Additionally, because of its poro-viscoelastic nature, strain magnitude in cartilage can depend on strain rate. The present study explores whether cartilage damage may develop dependent on strain rate, even when the presented damage numerical model is strain-dependent but not strain-rate-dependent. So far no experiments have been distinguished whether rate-dependent cartilage damage occurs in the collagen or in the non-fibrillar network. Thus, this research presents a finite element analysis model where, among others, collagen and non-fibrillar matrix are incorporated as well as a strain-dependent damage mechanism for these components. Collagen and non-fibrillar matrix stiffness decrease when a given strain is reached until complete failure upon reaching a maximum strain. With such model, indentation experiments at increasing strain rates were simulated on cartilage plugs and damage development was monitored over time. Collagen damage increased with increasing strain rate from 21 to 42 %. In contrast, damage in the non-fibrillar matrix decreased with increasing strain rates from 72 to 34 %. Damage started to develop at a depth of approximately 20 % of the sample height, and this was more pronounced for the slow and modest loading rates. However, the most severe damage at the end of the compression step occurred at the surface for the plugs subjected to 120 mm/min strain rate. In conclusion, the present study confirms that the location and magnitude of damage in cartilage may be strongly dependent on strain rate, even when damage occurs solely through a strain-dependent damage mechanism.

The acetabular labrum

The acetabular labrum is a soft-tissue structure which lines the acetabular rim of the hip joint. Its role in hip joint biomechanics and joint health has been of particular interest over the past decade. In normal hip joint biomechanics, the labrum is crucial in retaining a layer of pressurised intra-articular fluid for joint lubrication and load support/distribution. Its seal around the femoral head is further regarded as a contributing to hip stability through its suction effect. The labrum itself is also important in increasing contact area thereby reducing contact stress. Given the labrum’s role in normal hip joint biomechanics, surgical techniques for managing labral damage are continuously evolving as our understanding of its anatomy and function continue to progress. The current paper aims to review the anatomy and biomechanical function of the labrum and how they are affected by differing surgical techniques.

Take home message: The acetabular labrum plays a critical role in hip function and maintaining and restoring its function during surgical intervention remain an essential goal.

Cite this article: Bone Joint J 2016;98-B:730–5.