In vivo cartilage contact deformation in the healthy human tibiofemoral joint

Ultrasound Measurement of Femoral Articular Cartilage Thickness Before and After Marathon Running

OBJECTIVE

The purpose of this study was to use ultrasonography to measure femoral articular cartilage thickness changes during marathon running, which could support MRI studies showing that deformation of knee cartilage during long-distance running is no greater than that for other weight-bearing activities.

MATERIALS AND METHODS

Participants included 38 marathon runners with no knee pain or history of knee injury, aged 18-39. Ultrasound images of the femoral articular cartilage were taken two hours before and immediately after the race. Femoral articular cartilage thickness was measured at both the medial and lateral femoral condyles.

RESULTS

The maximum change in femoral articular cartilage thickness, measured at the left outer lateral femoral condyle, was 6.94% (P=.006). All other femoral articular cartilage thickness changes were not significant.

CONCLUSION

A change in femoral articular cartilage thickness of 6.94% supports our hypothesis that long-distance running does not induce deformational changes greater than that of regular daily activities. This study using ultrasonography supports MRI evidence that knee cartilage tolerates marathon running well.

Revisiting the flexion-valgus type unicondylar posterolateral tibial plateau depression fracture pattern: classification and treatment

OBJECTIVE

This study aimed to reclassify posterolateral tibial plateau fractures caused by a flexion-valgus force and describe this fracture pattern to provide a relatively programmed surgical treatment based on morphological characteristics that may improve reduction and stabilization.

METHODS

We retrospectively reviewed the fracture pattern and injury mechanism of patients with posterolateral tibial plateau fractures who underwent surgery at the First Affiliated Hospital of Nanjing Medical University between January 2014 and April 2020. The cohort was divided into three types. Type I was a depression fracture of the posterolateral platform with an intact posterolateral cortex. Type II was a depression fracture of the posterolateral platform with a disrupted posterolateral cortex. Type III was a depression fracture of the posterolateral platform in combination with anterior cruciate ligament (ACL) rupture or tibial insertion avulsion fracture of the ACL. The lateral window of the modified Frosch approach with an L-type locking plate was used for patients with type I and type III fractures. For patients with type II fractures, both lateral and posterolateral windows of the modified Frosch approach were used for surgery, and a T-plate on the posterior side with an L-plate on the lateral side were used for fixation. The Rasmussen radiology scoring was used to evaluate the quality of surgical reduction and the Rasmussen functional scoring evaluation standard was used to evaluate knee joint function.

RESULTS

A total of 69 tibial plateau fractures (36 male, 33 female) involving the posterolateral platform were discovered and included in this study. All patients suffered flexion-valgus force at the moment of the accident. There were 32 cases of Type I fracture, 28 cases of Type II fracture, and 9 cases of Type III fracture. The patients were followed up for 12-30 (mean 20.8 ± 9.4) months. The postoperative Rasmussen radiological scores for the three types of fractures were 15-17 (mean 16.31 ± 0.78), 14-17 (mean 15.93 ± 0.94), and 14-17 (mean 16.22 ± 0.97), respectively. The postoperative Rasmussen functional scores for the three types of fractions were 27-30 (mean 27.97 ± 0.90), 27-29 (mean 27.36 ± 0.56), and 27-29 (mean 27.56 ± 0.73), respectively.

CONCLUSION

Flexion-valgus posterolateral tibial plateau fractures were divided into three types based on the integrity of the posterolateral wall and ACL injuries. We hope the classification can play a certain reference role in recognizing and treating flexion-valgus-type posterolateral tibial plateau fractures.

Symmetry in knee arthrokinematics in healthy collegiate athletes during fast running and drop jump revealed through dynamic biplane radiography

OBJECTIVE

Changes in cartilage contact area and/or contact location after knee injury can initiate and exacerbate cartilage degeneration. Typically, the contralateral knee is used as a surrogate for native cartilage contact patterns on the injured knee. However, symmetry in cartilage contact patterns between healthy knees during high-impact activities is unknown.

METHOD

Tibiofemoral kinematics were measured on 19 collegiate athletes during fast running and drop jump using dynamic biplane radiography and a validated registration process that matched computed tomography (CT)-based bone models to the biplane radiographs. Cartilage contact area and location were measured with participant-specific magnetic resonance imaging (MRI)-based cartilage models superimposed on the CT-based bone models. Symmetry in cartilage contact area and location was assessed by the absolute side-to-side differences (SSD) within participants.

RESULTS

The SSD in contact area during running (7.7 ± 6.1% and 8.0 ± 4.6% in the medial and lateral compartments, respectively) was greater than during drop jump (4.2 ± 3.7% and 5.7 ± 2.6%, respectively) (95% CI of the difference: medial [2.4%, 6.6%], lateral [1.5%, 4.9%]). The average SSD in contact location was 3.5 mm or less in the anterior-posterior (AP) direction and 2.1 mm or less in the medial-lateral (ML) direction on the femur and tibia for both activities. The SSD in AP contact location on the femur was greater during running than during drop jump (95% CI of the difference: medial [1.6 mm, 3.6 mm], lateral [0.6 mm, 1.9 mm]).

CONCLUSION

This study provides context for interpreting results from previous studies on tibiofemoral arthrokinematics. Previously reported differences between ligament-repaired and contralateral knee arthrokinematics fall within the range of typical SSDs observed in healthy athletes. Previously reported arthrokinematics differences that exceed SSDs found in these healthy athletes occur only in the presence of anterior cruciate ligament (ACL) deficiency or meniscectomy.

Effect of surface matching mismatch of focal knee articular prosthetic on tibiofemoral contact stress using finite element analysis

Aims

Focal knee arthroplasty is an attractive alternative to knee arthroplasty for young patients because it allows preservation of a large amount of bone for potential revisions. However, the mechanical behaviour of cartilage has not yet been investigated because it is challenging to evaluate in vivo contact areas, pressure, and deformations from metal implants. Therefore, this study aimed to determine the contact pressure in the tibiofemoral joint with a focal knee arthroplasty using a finite element model.

Methods

The mechanical behaviour of the cartilage surrounding a metal implant was evaluated using finite element analysis. We modelled focal knee arthroplasty with placement flush, 0.5 mm deep, or protruding 0.5 mm with regard to the level of the surrounding cartilage. We compared contact stress and pressure for bone, implant, and cartilage under static loading conditions.

Results

Contact stress on medial and lateral femoral and tibial cartilages increased and decreased, respectively, the most and the least in the protruding model compared to the intact model. The deep model exhibited the closest tibiofemoral contact stress to the intact model. In addition, the deep model demonstrated load sharing between the bone and the implant, while the protruding and flush model showed stress shielding. The data revealed that resurfacing with a focal knee arthroplasty does not cause increased contact pressure with deep implantation. However, protruding implantation leads to increased contact pressure, decreased bone stress, and biomechanical disadvantage in an in vivo application.

Conclusion

These results show that it is preferable to leave an edge slightly deep rather than flush and protruding.

Auto-segmentation of the tibia and femur from knee MR images via deep learning and its application to cartilage strain and recovery

The ability to efficiently and reproducibly generate subject-specific 3D models of bone and soft tissue is important to many areas of musculoskeletal research. However, methodologies requiring such models have largely been limited by lengthy manual segmentation times. Recently, machine learning, and more specifically, convolutional neural networks, have shown potential to alleviate this bottleneck in research throughput. Thus, the purpose of this work was to develop a modified version of the convolutional neural network architecture U-Net to automate segmentation of the tibia and femur from double echo steady state knee magnetic resonance (MR) images. Our model was trained on a dataset of over 4,000 MR images from 34 subjects, segmented by three experienced researchers, and reviewed by a musculoskeletal radiologist. For our validation and testing sets, we achieved dice coefficients of 0.985 and 0.984, respectively. As further testing, we applied our trained model to a prior study of tibial cartilage strain and recovery. In this analysis, across all subjects, there were no statistically significant differences in cartilage strain between the machine learning and ground truth bone models, with a mean difference of 0.2 ± 0.7 % (mean ± 95 % confidence interval). This difference is within the measurement resolution of previous cartilage strain studies from our lab using manual segmentation. In summary, we successfully trained, validated, and tested a machine learning model capable of segmenting MR images of the knee, achieving results that are comparable to trained human segmenters.

Empirical joint contact mechanics: A comprehensive review

Empirical joint contact mechanics measurement (EJCM; e.g. contact area or force, surface velocities) enables critical investigations of the relationship between changing joint mechanics and the impact on surface-to-surface interactions. In orthopedic biomechanics, understanding the changes to cartilage contact mechanics following joint pathology or aging is critical due to its suggested role in the increased risk of osteoarthritis (OA), which might be due to changed kinematics and kinetics that alter the contact patterns within a joint. This article reviews and discusses EJCM approaches that have been applied to articulating joints such that readers across different disciplines will be informed of the various measurement and analysis techniques used in this field. The approaches reviewed include classical measurement approaches (radiographic and sectioning, dye staining, casting, surface proximity, and pressure measurement), stereophotogrammetry/motion analysis, computed tomography (CT), magnetic resonance imaging (MRI), and high-speed videoradiography. Perspectives on approaches to advance this field of EJCM are provided, including the value of considering relative velocity in joints, tractional stress, quantification of joint contact area shape, consideration of normalization techniques, net response (superposition) of multiple input variables, and establishing linkages to regional cartilage health status. EJCM measures continue to provide insights to advance our understanding of cartilage health and degeneration and provide avenues to assess the efficacy and guide future directions of developing interventions (e.g. surgical, biological, rehabilitative) to optimize joint's health and function long term.

Acetabular labrum and cartilage contact mechanics during pivoting and walking tasks in individuals with cam femoroacetabular impingement syndrome

Femoroacetabular impingement syndrome (FAIS) is a motion-related pathology of the hip characterized by pain, morphological abnormalities of the proximal femur, and an elevated risk of joint deterioration and hip osteoarthritis. Activities that require deep flexion are understood to induce impingement in cam FAIS patients, however, less demanding activities such as walking and pivoting may induce pain as well as alterations in kinematics and joint stability. Still, the paucity of quantitative descriptions of cam FAIS has hindered understanding underlying hip joint mechanics during such activities. Previous in silico studies have employed generalized model geometry or kinematics to simulate impingement between the femur and acetabulum, which may not accurately capture the interplay between morphology and motion. In this study, we utilized models with participant-specific bone and articular soft tissue anatomy and kinematics measured by dual-fluoroscopy to compare hip contact mechanics of cam FAIS patients to controls during four activities of daily living (internal/external pivoting and level/incline walking). Averaged across the gait cycle during incline walking, patients displayed increased strain in the anterior joint (labrum strain: p-value = 0.038, patients: 11.7 ± 6.7 %, controls: 5.0 ± 3.6 %; cartilage strain: p-value = 0.029, patients: 9.1 ± 3.3 %, controls: 4.2 ± 2.3). Patients also exhibited increased average anterior cartilage strains during external pivoting (p-value = 0.039; patients: 13.0 ± 9.2 %, controls: 3.9 ± 3.2 %]). No significant differences between patient and control contact area and strain were found for level walking and internal pivoting. Our study provides new insights into the biomechanics of cam FAIS, including spatiotemporal hip joint contact mechanics during activities of daily living.

Cartilage biomechanics: From the basic facts to the challenges of tissue engineering

Articular cartilage (AC) is the thin tissue that covers the long bone ends in the joints and that ensures the transmission of forces between adjacent bones while allowing nearly frictionless movements between them. AC repair is a technologic and scientific challenge that has been addressed with numerous approaches. A major deadlock is the capacity to take in account its complex mechanical properties in repair strategies. In this review, we first describe the major mechanical behaviors of AC for the non-specialists. Then, we show how researchers have progressively identified specific mechanical parameters using mathematical models. There are still gaps in our understanding of some of the observations concerning AC biomechanical properties, particularly the differences in extracellular matrix stiffness measured at the microscale and at the millimetric scale. Nevertheless, for bioengineering applications, AC repair strategies must take into account what are commonly considered the main mechanical features of cartilage: its ability to withstand high stresses through three main behaviors (elasticity, poroelasticity and swelling). Finally, we emphasize that future studies need to investigate AC mechanical properties at different scales, particularly the gradient of mechanical properties around cells and across the cartilage depth, and the differences in mechanical properties at different scales. This multi-scale approach could greatly enhance the success of AC restorative approaches.

Small Chondral Defects Affect Tibiofemoral Contact Area and Stress: Should a Lower Threshold Be Used for Intervention?

BACKGROUND

Chondral defects in the knee have biomechanical differences because of defect size and location. Prior literature only compares the maximum stress experienced with large defects.

HYPOTHESIS

It was hypothesized that pressure surrounding the chondral defect would increase with size and vary in location, such that a size cutoff exists that suggests surgical intervention.

STUDY DESIGN

Controlled laboratory study.

METHODS

Isolated chondral defects from 0.09 to 1.0 cm2 were created on the medial and lateral femoral condyles of 6 human cadaveric knees. The knees were fixed to a uniaxial load frame and loaded from 0 to 600 N at full extension. Another defect was created at the point of tibiofemoral contact at 30° of flexion. Tibiofemoral contact pressures were measured. Peak contact pressure was the highest value in the area delimited within a 3-mm rim around the defect. The location of the peak contact pressure was determined.

RESULTS

At full extension, the mean maximum pressures on the medial femoral condyle ranged from 4.30 to 6.91 MPa at 0.09 and 1.0 cm2, respectively (P < .01). The location of the peak pressure was found posteromedial in defects between 0.09 and 0.25 cm2, shifting anterolaterally at sizes 0.49 and 1.0 cm2 (P < .01). The maximum pressures on the lateral femoral condyle ranged from 3.63 to 5.81 MPa at 0.09 and 1.0 cm2, respectively (P = .02). The location of the peak contact pressure point was anterolateral in defects between 0.09 and 0.25 cm2, shifting posterolaterally at 0.49 and 1.0 cm2 (P < .01). No differences in contact pressure between full extension and 30° of flexion were found for either the lateral or medial condyles.

CONCLUSION

Full-thickness chondral defects bilaterally had a significant increase in contact pressure between defect sizes of 0.49 and 1.0 cm2. The location of the maximum contact pressures surrounding the lesion also varied with larger defects. Contact area redistribution and cartilage stress change may affect adjacent cartilage integrity.

CLINICAL RELEVANCE

Size cutoffs may exist earlier in the natural history of chondral defects than previously realized, suggesting a lower threshold for intervention.

Articular contact motion at the knee during daily activities

We combined mobile biplane X-ray imaging and magnetic resonance imaging to measure the regions of articular cartilage contact and cartilage thickness at the tibiofemoral and patellofemoral joints during six functional activities: standing, level walking, downhill walking, stair ascent, stair descent, and open-chain (non-weight-bearing) knee flexion. The contact centers traced similar paths on the medial and lateral femoral condyles, femoral trochlea, and patellar facet in all activities while their locations on the tibial plateau were more varied. The translations of the contact centers on the femur and patella were tightly coupled to the tibiofemoral flexion angle in all activities (r²  > 0.95) whereas those on the tibia were only moderately related to the flexion angle (r²  > 0.62). The regions of contacting cartilage were significantly thicker than the regions of non-contacting cartilage on the patella, femoral trochlea, and the medial and lateral tibial plateaus in all activities (p < 0.001). There were no significant differences in thickness between contacting and non-contacting cartilage on the medial and lateral femoral condyles in all activities, except open-chain knee flexion. Our results provide partial support for the proposition that cartilage thickness is adapted to joint load and do not exclude the possibility that other factors, such as joint congruence, also play a role in regulating the structure and organization of healthy cartilage. The data obtained in this study may serve as a guide when evaluating articular contact motion in osteoarthritic and reconstructed knees.

Design and validation of a semi-automatic bone segmentation algorithm from MRI to improve research efficiency

Segmentation of medical images into different tissue types is essential for many advancements in orthopaedic research; however, manual segmentation techniques can be time- and cost-prohibitive. The purpose of this work was to develop a semi-automatic segmentation algorithm that leverages gradients in spatial intensity to isolate the patella bone from magnetic resonance (MR) images of the knee that does not require a training set. The developed algorithm was validated in a sample of four human participants (in vivo) and three porcine stifle joints (ex vivo) using both magnetic resonance imaging (MRI) and computed tomography (CT). We assessed the repeatability (expressed as mean ± standard deviation) of the semi-automatic segmentation technique on: (1) the same MRI scan twice (Dice similarity coefficient = 0.988 ± 0.002; surface distance = − 0.01 ± 0.001 mm), (2) the scan/re-scan repeatability of the segmentation technique (surface distance = − 0.02 ± 0.03 mm), (3) how the semi-automatic segmentation technique compared to manual MRI segmentation (surface distance = − 0.02 ± 0.08 mm), and (4) how the semi-automatic segmentation technique compared when applied to both MRI and CT images of the same specimens (surface distance = − 0.02 ± 0.06 mm). Mean surface distances perpendicular to the cartilage surface were computed between pairs of patellar bone models. Critically, the semi-automatic segmentation algorithm developed in this work reduced segmentation time by approximately 75%. This method is promising for improving research throughput and potentially for use in generating training data for deep learning algorithms.

A Method to Compare Heterogeneous Types of Bone and Cartilage Meshes

Accurately capturing the bone and cartilage morphology and generating a mesh remains a critical step in the workflow of computational knee joint modeling. Currently, there is no standardized method to compare meshes of different element types and nodal densities, making comparisons across research teams a significant challenge. The aim of this paper is to describe a method to quantify differences in knee joint bone and cartilages meshes, independent of bone and cartilage mesh topology. Bone mesh-to-mesh distances, subchondral bone boundaries, and cartilage thicknesses from meshes of any type of mesh are obtained using a series of steps involving registration, resampling, and radial basis function fitting after which the comparisons are performed. Subchondral bone boundaries and cartilage thicknesses are calculated and visualized in a common frame of reference for comparison. The established method is applied to models developed by five modeling teams. Our approach to obtain bone mesh-to-mesh distances decreased the divergence seen in selecting a reference mesh (i.e., comparing mesh A-to-B versus mesh B-to-A). In general, the bone morphology was similar across teams. The cartilage thicknesses for all models were calculated and the mean absolute cartilage thickness difference was presented, the articulating areas had the best agreement across teams. The teams showed disagreement on the subchondral bone boundaries. The method presented in this paper allows for objective comparisons of bone and cartilage geometry that is agnostic to mesh type and nodal density.

Tibiofemoral contact and alignment in patients with anterior cruciate ligament rupture treated nonoperatively versus reconstruction : an upright, open MRI study

AIMS

Anterior cruciate ligament (ACL) rupture commonly leads to post-traumatic osteoarthritis, regardless of surgical reconstruction. This study uses standing MRI to investigate changes in contact area, contact centroid location, and tibiofemoral alignment between ACL-injured knees and healthy controls, to examine the effect of ACL reconstruction on these parameters.

METHODS

An upright, open MRI was used to directly measure tibiofemoral contact area, centroid location, and alignment in 18 individuals with unilateral ACL rupture within the last five years. Eight participants had been treated nonoperatively and ten had ACL reconstruction performed within one year of injury. All participants were high-functioning and had returned to sport or recreational activities. Healthy contralateral knees served as controls. Participants were imaged in a standing posture with knees fully extended.

RESULTS

Participants' mean age was 28.4 years (SD 7.3), the mean time since injury was 2.7 years (SD 1.6), and the mean International Knee Documentation Subjective Knee Form score was 84.4 (SD 13.5). ACL injury was associated with a 10% increase (p = 0.001) in contact area, controlling for compartment, sex, posture, age, body mass, and time since injury. ACL injury was associated with a 5.2% more posteriorly translated medial centroid (p = 0.001), equivalent to a 2.6 mm posterior translation on a representative tibia with mean posteroanterior width of 49.4 mm. Relative to the femur, the tibiae of ACL ruptured knees were 2.3 mm more anteriorly translated (p = 0.003) and 2.6° less externally rotated (p = 0.010) than healthy controls. ACL reconstruction was not associated with an improvement in any measure.

CONCLUSION

ACL rupture was associated with an increased contact area, posteriorly translated medial centroid, anterior tibial translation, and reduced tibial external rotation in full extension. These changes were present 2.7 years post-injury regardless of ACL reconstruction status. Cite this article: Bone Joint J 2021;103-B(9):1505-1513.

Multifaceted Design and Emerging Applications of Tissue Adhesives

Tissue adhesives can form appreciable adhesion with tissues and have found clinical use in a variety of medical settings such as wound closure, surgical sealants, regenerative medicine, and device attachment. The advantages of tissue adhesives include ease of implementation, rapid application, mitigation of tissue damage, and compatibility with minimally invasive procedures. The field of tissue adhesives is rapidly evolving, leading to tissue adhesives with superior mechanical properties and advanced functionality. Such adhesives enable new applications ranging from mobile health to cancer treatment. To provide guidelines for the rational design of tissue adhesives, here, existing strategies for tissue adhesives are synthesized into a multifaceted design, which comprises three design elements: the tissue, the adhesive surface, and the adhesive matrix. The mechanical, chemical, and biological considerations associated with each design element are reviewed. Throughout the report, the limitations of existing tissue adhesives and immediate opportunities for improvement are discussed. The recent progress of tissue adhesives in topical and implantable applications is highlighted, and then future directions toward next-generation tissue adhesives are outlined. The development of tissue adhesives will fuse disciplines and make broad impacts in engineering and medicine.

Reliability of tibiofemoral contact area and centroid location in upright, open MRI

BACKGROUND

Imaging cannot be performed during natural weightbearing in biomechanical studies using conventional closed-bore MRI, which has necessitated simulating weightbearing load on the joint. Upright, open MRI (UO-MRI) allows for joint imaging during natural weightbearing and may have the potential to better characterize the biomechanical effect of tibiofemoral pathology involving soft tissues. However open MRI scanners have lower field strengths than closed-bore scanners, which limits the image quality that can be obtained. Thus, there is a need to establish the reliability of measurements in upright weightbearing postures obtained using UO-MRI.

METHODS

Knees of five participants with prior anterior cruciate ligament (ACL) rupture were scanned standing in a 0.5 T upright open MRI scanner using a 3D DESS sequence. Manual segmentation of cartilage regions in contact was performed and centroids of these contact areas were automatically determined for the medial and lateral tibiofemoral compartments. Inter-rater, test-retest, and intra-rater reliability were determined and quantified using intra-class correlation (ICC_3,1), standard error of measurement (SEM), and smallest detectable change with 95% confidence (SDC₉₅). Accuracy was assessed by using a high-resolution 7 T MRI as a reference.

RESULTS

Contact area and centroid location reliability (inter-rater, test-retest, and intra-rater) for sagittal scans in the medial compartment had ICC_3,1 values from 0.95-0.99 and 0.98-0.99 respectively. In the lateral compartment, contact area and centroid location reliability ICC_3,1 values ranged from 0.83-0.91 and 0.95-1.00 respectively. The smallest detectable change in contact area was 1.28% in the medial compartment and 0.95% in the lateral compartment. Contact area and centroid location reliability for coronal scans in the medial compartment had ICC_3,1 values from 0.90-0.98 and 0.98-1.00 respectively, and in the lateral compartment ICC_3,1 ranged from 0.76-0.94 and 0.93-1.00 respectively. The smallest detectable change in contact area was 0.65% in the medial compartment and 1.41% in the lateral compartment. Contact area was accurate to within a mean absolute error of 11.0 mm².

CONCLUSIONS

Knee contact area and contact centroid location can be assessed in upright weightbearing MRI with good to excellent reliability. The lower field strength used in upright, weightbearing MRI does not compromise the reliability of tibiofemoral contact area and centroid location measures.

An MRI-compatible varus-valgus loading device for whole-knee joint functionality assessment based on compartmental compression: a proof-of-concept study

OBJECTIVE

Beyond static assessment, functional techniques are increasingly applied in magnetic resonance imaging (MRI) studies. Stress MRI techniques bring together MRI and mechanical loading to study knee joint and tissue functionality, yet prototypical axial compressive loading devices are bulky and complex to operate. This study aimed to design and validate an MRI-compatible pressure-controlled varus-valgus loading device that applies loading along the joint line.

METHODS

Following the device's thorough validation, we demonstrated proof of concept by subjecting a structurally intact human cadaveric knee joint to serial imaging in unloaded and loaded configurations, i.e. to varus and valgus loading at 7.5 kPa (= 73.5 N), 15 kPa (= 147.1 N), and 22.5 kPa (= 220.6 N). Following clinical standard (PDw fs) and high-resolution 3D water-selective cartilage (WATSc) sequences, we performed manual segmentations and computations of morphometric cartilage measures. We used CT and radiography (to quantify joint space widths) and histology and biomechanics (to assess tissue quality) as references.

RESULTS

We found (sub)regional decreases in cartilage volume, thickness, and mean joint space widths reflective of areal pressurization of the medial and lateral femorotibial compartments.

DISCUSSION

Once substantiated by larger sample sizes, varus-valgus loading may provide a powerful alternative stress MRI technique.

Cartilage and collagen mechanics under large-strain shear within in vivo and at supraphysiogical temperatures

Human joints, particularly those of extremities, experience a significant range of temperatures in vivo. Joint temperature influences the mechanics of both joint and cartilage, and the mechanics of cartilage can affect the temperature of both joint and cartilage. Thermal treatments and tissue repairs, such as thermal chondroplasty, and ex vivo tissue engineering may also expose cartilage to supraphysiological temperatures. Furthermore, although cartilage undergoes principally compressive loads in vivo, shear strain plays a significant role at larger compressive strains. Thus, we aimed to determine whether and how the bulk mechanical responses of cartilage undergoing large-strain shear change (1) within the range of temperatures relevant in vivo, and (2) both during and after supraphysiological thermal treatments. We completed large-strain shear tests (10 and 15%) at four thermal conditions: 24^∘C and 40^∘C to span the in vivo range, and 70^∘C and 24^∘C repeated after 70^∘C to explore mechanics during and after potential treatments. We calculated the bulk mechanical responses (strain-energy dissipation densities, peak-to-peak shear stresses, and peak-effective shear moduli) as of function of temperature and used statistical methods to probe significant differences. To probe the mechanisms underlying differences we assessed specimens, principally the type II collagen, with imaging (second harmonic generation and transmission electron microscopies, and histology) and assessed the temperature-dependent mechanics of type II collagen molecules within cartilage using steered molecular dynamics simulations. Our results suggest that the bulk mechanical responses of cartilage depend significantly on temperature both within the in vivo range and at supraphysiological temperatures, showing significant reductions in all mechanical measures with increasing temperature. Using imaging and simulations we determined that one underlying mechanism explaining our results may be changes in the molecular deformation profiles of collagen molecules versus temperature, likely compounded at larger length scales. These new insights into the mechanics of cartilage and collagen may suggest new treatment targets for damaged or osteoarthritic cartilage.

Through-thickness patterns of shear strain evolve in early osteoarthritis

OBJECTIVE

Given the structural changes associated with the progression of Osteoarthritis (OA), we hypothesized that patterns of through-thickness, large-strain shear evolve with early-stage OA. We therefore aimed to determine whether and how patterns of shear strains change during early-stage OA to 1) gain insight into the progression of OA by quantifying changes in local deformations; 2) gauge the potential of patterns in shear strain to serve as image-based biomarkers of early-stage OA; and 3) provide high-resolution, through-thickness data for proposing, fitting, and validating constitutive models for cartilage.

DESIGN

We completed displacement-driven, large-strain shear tests (5, 10, 15%) on 44 specimens of variably advanced osteoarthritic human articular cartilage as determined by both Osteoarthritis Research Society International (OARSI) grade and PLM-CO score. We recorded the through-thickness deformations with a stereo-camera system and processed these data using digital image correlation (DIC) to determine full-thickness patterns of strains and relative zonal recruitments, i.e., the average shear strain in a through-thickness zone weighted by its relative thickness and normalized by the applied strain.

RESULTS

We observed three general shapes for the curves of averaged through-thickness, Green-Lagrange shear strains during progression of OA. We also observed that during the progression of OA only the deep zone is recruited differently under shear in a statistically significant way.

CONCLUSIONS

We propose that changes in through-thickness patterns of shear strain could provide sensitive biomarkers for early clinical detection of OA. The relative zonal recruitment of the deep zone decreases with progressing OA (OARSI grade) and microstructural remodeling (PLM-CO score), which do not consistently affect recruitment of the superficial and middle zones.

Tibiofemoral Cartilage Contact Differences Between Level Walking and Downhill Running

Background:

Some studies have suggested that altered tibiofemoral cartilage contact behavior (arthrokinematics) may contribute to long-term cartilage degeneration, potentially leading to tibiofemoral osteoarthritis. However, few studies have assessed normal tibiofemoral arthrokinematics during dynamic activities.

Purpose:

To characterize tibiofemoral arthrokinematics during the impact phase of level walking and downhill running.

Study Design:

Descriptive laboratory study.

Methods:

Arthrokinematic data were collected on uninjured knees of 44 participants (mean age, 20.7 ± 6.6 years). Using a dynamic stereoradiographic imaging system with superimposed 3-dimensional bone models from computed tomography and magnetic resonance imaging of participant-specific tibiofemoral joints, arthrokinematics were assessed during the first 15% of the gait cycle during level walking and the first 10% of the gait cycle during downhill running.

Results:

During level walking and downhill running, the medial compartment had a greater cartilage contact area versus the lateral compartment. Both compartments had a significantly less cartilage contact area during running versus walking (medial compartment gait cycle affected: 8%-10%; lateral compartment gait cycle affected: 5%-10%). Further, medial and lateral compartment tibiofemoral contact paths were significantly more posterior and longer during downhill running.

Conclusion:

There was a decreased tibiofemoral cartilage contact area during downhill running compared with level walking, suggesting that underlying bone morphology may play a key role in determining the size of cartilage contact regions.

Clinical Relevance:

This study provides the first data characterizing tibiofemoral cartilage contact patterns during level walking and downhill running. These results provide evidence in support of performing biomechanical assessments during both level walking and downhill running to obtain a comprehensive picture of tibiofemoral cartilage behavior after clinical interventions.

Low field magnetic resonance imaging of the equine distal interphalangeal joint: Comparison between weight-bearing and non-weight-bearing conditions

This descriptive study aimed to compare the magnetic resonance appearance of the distal interphalangeal joint articular cartilage between standing weight-bearing and non-weight-bearing conditions. Ten forefeet of live horses were scanned in a standing low-field magnetic resonance system (0.27 T). After euthanasia for reasons unrelated to the study, the non-weight-bearing isolated feet were scanned in a vertical positioning reproducing limb orientation in live horses. The same acquisition settings as during the weight-bearing examination were used. Thickness and cross-sectional area of the distal interphalangeal articular cartilage and joint space were measured on tridimensional T1-weighted gradient echo high resolution frontal and sagittal images at predetermined landmarks in both conditions and were compared using a linear mixed-effects model. Frontal images were randomized and submitted to 9 blinded readers with 3 different experience levels for identification of weight-bearing versus non-weight-bearing acquisitions based on cartilage appearance. Weight-bearing limbs had significantly thinner distal interphalangeal cartilage (p = 0.0001) than non-weight-bearing limbs. This change was greater in the distal phalanx cartilage than that of the middle phalanx. Blinded readers correctly identified 83% (range 65 to 95%) of the images as weight-bearing or non-weight-bearing acquisitions, with significantly different results observed among the different readers (p < 0.001) and groups (p < 0.001). These results indicate that distal interphalangeal articular cartilage and particularly cartilage of the distal phalanx thins when weight-bearing compared to the non-weight-bearing standing postmortem conditions and suggest that cartilage abnormalities may be more difficult to identify on weight-bearing standing magnetic resonance imaging.

Cartilage defect location and stiffness predispose the tibiofemoral joint to aberrant loading conditions during stance phase of gait

OBJECTIVES

The current study quantified the influence of cartilage defect location on the tibiofemoral load distribution during gait. Furthermore, changes in local mechanical stiffness representative for matrix damage or bone ingrowth were investigated. This may provide insights in the mechanical factors contributing to cartilage degeneration in the presence of an articular cartilage defect.

METHODS

The load distribution following cartilage defects was calculated using a musculoskeletal model that included tibiofemoral and patellofemoral joints with 6 degrees-of-freedom. Circular cartilage defects of 100 mm2 were created at different locations in the tibiofemoral contact geometry. By assigning different mechanical properties to these defect locations, softening and hardening of the tissue were evaluated.

RESULTS

Results indicate that cartilage defects located at the load-bearing area only affect the load distribution of the involved compartment. Cartilage defects in the central part of the tibia plateau and anterior-central part of the medial femoral condyle present the largest influence on load distribution. Softening at the defect location results in overloading, i.e., increased contact pressure and compressive strains, of the surrounding tissue. In contrast, inside the defect, the contact pressure decreases and the compressive strain increases. Hardening at the defect location presents the opposite results in load distribution compared to softening. Sensitivity analysis reveals that the surrounding contact pressure, contact force and compressive strain alter significantly when the elastic modulus is below 7 MPa or above 18 MPa.

CONCLUSION

Alterations in local mechanical behavior within the high load bearing area resulted in aberrant loading conditions, thereby potentially affecting the homeostatic balance not only at the defect but also at the tissue surrounding and opposing the defect. Especially, cartilage softening predisposes the tissue to loads that may contribute to accelerated risk of cartilage degeneration and the initiation or progression towards osteoarthritis of the whole compartment.

On incorporating osmotic prestretch/prestress in image-driven finite element simulations of cartilage

Medical imaging performed in vivo captures geometries under Donnan osmotic loading, even when the articulating joint is otherwise mechanically unloaded. Hence patient-specific finite element (FE) models constructed from such medical images of cartilage represent osmotically induced prestretched/prestressed states. When applying classical modeling approaches to patient-specific simulations of cartilage a theoretical inconsistency arises: the in-vivo imaged geometry (used to construct the model) is not an unloaded, stress-free reference configuration. Furthermore when fitting nonlinear constitutive models that include osmotic swelling (to obtain material parameters), if one assumes that experimental data-generated from osmotically loaded cartilage-begin from a stress-free reference configuration the fitted stress-stretch relationship (parameters) obtained will actually describe a different behavior. In this study we: (1) establish a practical computational method to include osmotically induced prestretch/prestress in image-driven simulations of cartilage; and (2) investigate the influence of considering the prestretched/prestressed state both when fitting fiber-reinforced, biphasic constitutive models of cartilage that include osmotic swelling and when simulating cartilage responses. Our results highlight the importance of determining the prestretched/prestressed state within cartilage induced by osmotic loading in the imaged configuration prior to solving boundary value problems of interest. With our new constitutive model and modeling methods, we aim to improve the fidelity of FE-based, patient-specific biomechanical simulations of joints and cartilage. Improved simulations can provide medical researchers with new information often unavailable in a clinical setting, information that may contribute to better insight into the pathophysiology of cartilage diseases.

Analysis of in-vivo articular cartilage contact surface of the knee during a step-up motion

BACKGROUND

Numerous studies have reported on the tibiofemoral articular cartilage contact kinematics, however, no data has been reported on the articular cartilage geometry at the contact area. This study investigated the in-vivo tibiofemoral articular cartilage contact biomechanics during a dynamic step-up motion.

METHODS

Ten healthy subjects were imaged using a validated magnetic resonance and dual fluoroscopic imaging technique during a step-up motion. Three-dimensional bone and cartilage models were constructed from the magnetic resonance images. The cartilage contact along the motion path was analyzed, including cartilage contact location and the cartilage surface geometry at the contact area.

FINDINGS

The cartilage contact excursions were similar in anteroposterior and mediolateral directions in the medial and lateral compartments of the tibia plateau (P>0.05). Both medial and lateral compartments were under convex (femur) to convex (tibia) contact in the sagittal plane, and under convex (femur) to concave (tibia) contact in the coronal plane. The medial tibial articular contact radius was larger than the lateral side in the sagittal plane along the motion path (P<0.001).

INTERPRETATIONS

These data revealed that both the medial and lateral compartments of the knee experienced convex (femur) to convex (tibia) contact in sagittal plane (or anteroposterior direction) during the dynamic step-up motion. These data could provide new insight into the in-vivo cartilage contact biomechanics research, and may provide guidelines for development of anatomical total knee arthroplasties that are aimed to reproduce normal knee joint kinematics.

Comparison of mechanical compressive properties of commercial and autologous fibrin glues for tissue engineering applications

BACKGROUND

Fibrin glues are widely used in orthopedic surgery as adhesives and hemostatic agents. We evaluated the compressive properties of selected fibrin glues in order to identify which are appropriate for tissue regeneration applications subject to compression.

METHODS

Uniaxial unconfined compression tests were performed on fibrin gels prepared from commercial and autologous products: (1) Evicel (Ethicon), (2) Tisseel (Baxter), (3) Angel (Arthrex), and (4) ProPlaz (Biorich). Cyclic loads were applied from 0 to 30% strain for 100cycles at 0.5Hz. Following cyclic testing, specimens were subjected to ramp displacement of 1% strain per second to 80% strain.

FINDINGS

Throughout cyclic loading, Evicel and Tisseel deformed (shortened) less than Angel at all but one time point, and deformed less than ProPlaz at cycles 10 and 20. The dynamic moduli, peak stress, and strain energy were significantly greater in Tisseel than all other groups. Evicel displayed significantly greater dynamic moduli, peak stress, and strain energy than Angel and ProPlaz. Following cyclic testing, Tisseel and Evicel were significantly less deformed than Angel. No specimens exhibited gross failure during ramp loading to 80% strain. Ramp loading trends mirrored those of cyclic loading.

INTERPRETATION

The tested commercial glues were significantly more resistant to compression than the autologous products. The compressive properties of Tisseel were approximately twice those of Evicel. All preparations displayed moduli multiple orders of magnitude less than that of native articular cartilage. We conclude that in knee surgeries requiring fibrin glue to undergo compression of daily activity, commercial products are preferable to autologous preparations from platelet-poor plasma, though both will deform significantly.

Coculture of bovine cartilage with synovium and fibrous joint capsule increases aggrecanase and matrix metalloproteinase activity

Background

A hallmark of osteoarthritis is increased proteolytic cleavage of aggrecan. Cross talk between cartilage and the synovium + joint capsule (SJC) can drive cartilage degradation by activating proteases in both tissues. We investigated aggrecan proteolysis patterns in cartilage explants using a physiologically relevant explant model of joint injury combining cartilage mechanical compression and coincubation with SJC.

Methods

Bovine cartilage explants were untreated; coincubated with SJC; or subjected to mechanical injury and coincubated with SJC, mechanical injury alone, or mechanical injury and incubated with tumor necrosis factor-α (TNF-α). To compare the patterns of aggrecan proteolysis between 6 h and 16 days, release of sulfated glycosaminoglycans and specific proteolytic aggrecan fragments into medium or remaining in cartilage explants was measured by dimethylmethylene blue and Western blot analysis.

Results

Aggrecanase activity toward aggrecan was observed in all conditions, but it was directed toward the TEGE↓ARGS interglobular domain (IGD) site only when cartilage was coincubated with SJC or TNF-α. Matrix metalloproteinase (MMP) activity at the aggrecan IGD site (IPES↓FFGV) was not detected when cartilage was exposed to TNF-α (up to 6 days), but it was in all other conditions. Compared with when bovine cartilage was left untreated or subjected to mechanical injury alone, additional aggrecan fragment types were released into medium and proteolysis of aggrecan started at an earlier time when SJC was present.

Conclusions

Indicative of different proteolytic pathways for aggrecan degradation, the SJC increases both aggrecanase and MMP activity toward aggrecan, whereas TNF-α inhibits MMP activity against the IGD of aggrecan.

Electronic supplementary material

The online version of this article (doi:10.1186/s13075-017-1318-9) contains supplementary material, which is available to authorized users.

Tibio-femoral joint contact in healthy and osteoarthritic knees during quasi-static squat: A bi-planar X-ray analysis

The aim of this study was to quantify the tibio-femoral contact point (CP) locations in healthy and osteoarthritic (OA) subjects during a weight-bearing squat using stand-alone biplanar X-ray images. Ten healthy and 9 severe OA subjects performed quasi-static squats. Bi-planar X-ray images were recorded at 0°, 15°, 30°, 45°, and 70° of knee flexion. A reconstruction/registration process was used to create 3D models of tibia, fibula, and femur from bi-planar X-rays and to measure their positions at each posture. A weighted centroid of proximity algorithm was used to calculate the tibio-femoral CP locations. The accuracy of the reconstruction/registration process in measuring the quasi-static kinematics and the contact parameters was evaluated in a validation study. The quasi-static kinematics data revealed that in OA knees, adduction angles were greater (p<0.01), and the femur was located more medially relative to the tibia (p<0.01). Similarly, the average CP locations on the medial and lateral tibial plateaus of the OA patients were shifted (6.5±0.7mm; p<0.01) and (9.6±3.1mm; p<0.01) medially compared to the healthy group. From 0° to 70° flexion, CPs moved 8.1±5.3mm and 8.9±5.3mm posteriorly on the medial and lateral plateaus of healthy knees; while in OA joints CPs moved 10.1±8.4mm and 3.6±2.8mm posteriorly. The average minimum tibio-femoral bone-to-bone distances of the OA joints were lower in both compartments (p<0.01). The CPs in the OA joints were located more medially and displayed a higher ratio of medial to lateral posterior translations compared to healthy joints.

Osteochondral Allograft Donor-Host Matching by the Femoral Condyle Radius of Curvature

BACKGROUND

Conventional osteochondral allograft (OCA) matching, requiring orthotopic, size-matched condyles, and narrow surgical time windows often prohibit timely transplantation.

HYPOTHESIS

The femoral condyle radius of curvature (RoC) is an appropriate, isolated criterion for donor-host matching in fresh OCAs, potentially enhancing matching efficiency when compared with conventional matching techniques.

STUDY DESIGN

Descriptive laboratory study.

METHODS

In part 1 of this study, 3-dimensional digital reconstructions of 14 randomly selected, cadaveric distal femoral hemicondyles were performed. Each condyle was divided into anterior, middle, and posterior zones. A virtual best-fit grid was applied to each, and each zone's sagittal- and coronal-plane RoCs were determined. Seven nonorthotopic OCA transplantations were performed based on RoC matching with 1-mm tolerance, and the preoperative and postoperative surface geometry were quantified to assess the accuracy of articular surface restoration. Of note, each donor-host pair did not match by the conventional method. In part 2 of this study, 12 cadaveric distal femora were categorized by size and digitized in the aforementioned manner. Simulated circular defects measuring 20, 25, and 30 mm in diameter were introduced into each zone. OCA matches were determined based on donor and host RoCs, and the total number of potential matches (of 71 total comparisons) was recorded as a percentage for each simulated defect. Finally, the results of RoC matching were compared with the conventional method for simulated defects in all zones of both the medial and lateral femoral condyles.

RESULTS

Part 1: The mean surface deviation after OCA transplantation was -0.09 mm, with a mean maximum protrusion at any point of 0.59 mm. Part 2: Using the RoC, 20-mm defects had a 100% chance of being matched. Defects of 25 and 30 mm had a 91% and 64% chance of being matched, respectively. Compared with the conventional method, the RoC method yielded a 3.2-fold greater match rate for lesions of the medial and lateral femoral condyles ( P = .02).

CONCLUSION

This investigation shows that femoral condyle RoCs in the sagittal and coronal planes may be useful, alternative matching criteria, expanding on current standards.

CLINICAL RELEVANCE

These matching criteria may increase the number of available matches, reduce wait times for patients, and reduce the number of wasted grafts.

T2* mapping of subtalar cartilage: Precision and association between anatomical variants and cartilage composition

Hindfoot arthritis is an important contributor to foot pain and physical disability. While the subtalar joint (STJ) is most frequently affected, anatomical variants such as facet configuration were suggested to further STJ cartilage deterioration. T2* mapping enables detection of ultra-structural cartilage change, particularly in thin cartilage layers, but its feasibility in the STJ has not yet been evaluated. The purpose of this study was to evaluate segmentation consistency and inter-scan short-term precision error of T2* mapping of talocalcaneal cartilage and to investigate the relationship between facet configuration and STJ T2* values. Using 3Tesla morphological magnetic resonance imaging, STJ configuration was categorized according to the degree of fusion between anterior, medial, or posterior facets. Subsequently, two repeats of multi-echo gradient recalled echo sequences were performed to obtain T2* maps with repositioning. Segmentation consistency of T2* values attained an ICC of 0.90 (95%CI 0.69-0.99). Precision errors comprised a coefficient of variation (CV) ranging 0.01-0.05, corresponding to a root mean square CV of 0.03-0.04. A 2-joint configuration type (i.e., fused anterior-medial facets) was significantly associated with a decrease in posterior facet T2* values (β = -0.6, p = 0.046). STJ T2* mapping is a reliable method requiring at least a 4% difference within people to enable detection of significant change. Anatomical variants in STJ configuration were associated with T2* values with the more stable 3-joint types exhibiting more favorable cartilage outcomes. Longer-term larger-scaled studies focusing on arthritis pathology are needed to further support the use of T2* mapping in hindfoot disease monitoring. © 2016 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 34:1969-1976, 2016.

Shear deformations of human articular cartilage: Certain mechanical anisotropies apparent at large but not small shear strains

Articular cartilage has pronounced through-the-thickness heterogeneity in both ultrastructure and mechanical function. The tissue undergoes a combination of large deformations in vivo, where shear is critical in both failure and chondrocyte death. Yet the microstructure mechanical response of cartilage to multi-axial large shear deformations is unknown. We harvested a total of 42 cartilage specimens from seven matched locations across the lateral femoral condyles and patellofemoral grooves of six human male donors (30.2±8.8yrs old, M±SD). With each specimen we applied a range of quasi-static, multi-axial large (simple) shear displacements both parallel and perpendicular to the local split-line direction (SLD). Shear stresses in cartilage specimens from the patellofemoral grooves were higher, and more energy was dissipated, at all applied strains under loading parallel to the local SLD versus perpendicular, while specimens from the lateral condyles were mechanically anisotropic only under larger strains of 20% and 25%. Cartilage also showed significant intra-donor variability at larger shear strains but no significant inter-donor variability. Overall, shear strain-energy dissipation was almost constant at 5% applied shear strain and increased nonlinearly with increasing shear magnitude. Our results suggest that full understanding of cartilage mechanics requires large-strain analyses to account for nonlinear, anisotropic and location-dependent effects not fully realized at small strains.

Effect of normal gait on in vivo tibiofemoral cartilage strains

Altered cartilage loading is believed to be associated with osteoarthritis development. However, there are limited data regarding the influence of normal gait, an essential daily loading activity, on cartilage strains. In this study, 8 healthy subjects with no history of knee surgery or injury underwent magnetic resonance imaging of a single knee prior to and following a 20-min walking activity at approximately 1.1m/s. Bone and cartilage surfaces were segmented from these images and compiled into 3-dimensional models of the tibia, femur, and associated cartilage. Thickness changes were measured across a grid of evenly spaced points spanning the models of the articular surfaces. Averaged compartmental strains and local strains were then calculated. Overall compartmental strains after the walking activity were found to be significantly different from zero in all four tibiofemoral compartments, with tibial cartilage strain being significantly larger than femoral cartilage strain. These results provide baseline data regarding the normal tibiofemoral cartilage strain response to gait. Additionally, the technique employed in this study has potential to be used as a "stress test" to understand how factors including age, weight, and injury influence tibiofemoral cartilage strain response, essential information in the development of potential treatment strategies for the prevention of osteoarthritis.

Importance of Patella, Quadriceps Forces, and Depthwise Cartilage Structure on Knee Joint Motion and Cartilage Response During Gait

In finite-element (FE) models of the knee joint, patella is often omitted. We investigated the importance of patella and quadriceps forces on the knee joint motion by creating an FE model of the subject's knee. In addition, depthwise strains and stresses in patellar cartilage with different tissue properties were determined. An FE model was created from subject's magnetic resonance images. Knee rotations, moments, and translational forces during gait were recorded in a motion laboratory and used as an input for the model. Three material models were implemented into the patellar cartilage: (1) homogeneous model, (2) inhomogeneous (arcadelike fibrils), and (3) random fibrils at the superficial zone, mimicking early stages of osteoarthritis (OA). Implementation of patella and quadriceps forces into the model substantially reduced the internal–external femoral rotations (versus without patella). The simulated rotations in the model with the patella matched the measured rotations at its best. In the inhomogeneous model, maximum principal stresses increased substantially in the middle zone of the cartilage. The early OA model showed increased compressive strains in the superficial and middle zones of the cartilage and decreased stresses and fibril strains especially in the middle zone. The results suggest that patella and quadriceps forces should be included in moment- and force-driven FE knee joint models. The results indicate that the middle zone has a major role in resisting shear forces in the patellar cartilage. Also, early degenerative changes in the collagen network substantially affect the cartilage depthwise response in the patella during walking.

In vivo deformation of thin cartilage layers: Feasibility and applicability of T2* mapping

The objectives of this study were as follows: (i) to assess segmentation consistency and scan precision of T2* mapping of human tibio-talar cartilage, and (ii) to monitor changes in T2* relaxation times of ankle cartilage immediately following a clinically relevant in vivo exercise and during recovery. Using multi-echo gradient recalled echo sequences, averaged T2* values were calculated for tibio-talar cartilage layers in 10 healthy volunteers. Segmentation consistency and scan precision were determined from two repeated segmentations and two repeated acquisitions with repositioning, respectively. Subsequently, acute in vivo cartilage loading responses were monitored by calculating averaged tibio-talar T2* values at rest, immediately after (i.e., deformation) and at 15 min (i.e., recovery) following a 30-repetition knee bending exercise. Precision errors attained 4-6% with excellent segmentation consistency point estimates (i.e., intra-rater ICC of 0.95) and acceptable limits of confidence. At deformation, T2* values were increased in both layers [+16.1 (10.7)%, p = 0.004 and +17.3 (15.3)%, p = 0.023, for the talus and tibia, respectively] whereas during recovery no significant changes could be established when comparing to baseline [talar cartilage: +5.2 (8.2)%, p = 0.26 and tibial cartilage: +6.6 (10.4)%, p = 0.23]. T2* mapping is a viable method to monitor deformational behavior in thin cartilage layers such as ankle cartilage. Longitudinal changes in T2* can be reliably appraised and require at least 4-6% differences to ascertain statistical significance. The ability to detect considerable change even after non-strenuous loading events, endorses T2* mapping as an innovative method to evaluate the effects of therapeutic exercise on thin cartilage layers. © 2015 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 34:771-778, 2016.

Probabilistic evaluation of the material properties of the in vivo subject-specific articular surface using a computational model

This article used probabilistic analysis to evaluate material properties of the in vivo subject-specific tibiofemoral (TF) joint model. Sensitivity analysis, based on a Monte Carlo (MC) method, was performed using a subject-specific finite element (FE) model generated from in vivo computed tomography (CT) and magnetic resonance imaging (MRI) data, subjected to two different loading conditions. Specifically, the effects of inherent uncertainty in ligament stiffness, horn attachment stiffness, and articular surface material properties were assessed using multifactorial global sensitivity analysis. The MRI images were taken before and after axial compression, and when the flexion condition had been maintained at up to 90 degree flexion in the subject-specific knee joint. The loading conditions of the probabilistic subject-specific FE model (axial compression and 90 degree flexion) were similar to the MRI acquisition setup. We were able to detect the influence of material parameters while maintaining the potential effect of parametric interactions. Throughout the in silico property optimization, a subject-specific FE model was used and less sensitive parameters were eliminated in the global sensitivity method. Soft tissue material properties were estimated using an optimization procedure that involved the minimization of the differences between the kinematics predicted by the subject-specific model and those obtained through in vivo subject-specific data. The results of this approach suggest that the articular surface mechanical properties could be found by using in vivo measurements, which clarifies the valuable tool for future subject-specific studies related to TF joint scaffolds, allografts and biologics. © 2016 Wiley Periodicals, Inc. J Biomed Mater Res Part B: Appl Biomater, 105B: 1390-1400, 2017.

Mechanical properties of normal and osteoarthritic human articular cartilage

Isotropic hyperelastic models have been used to determine the material properties of normal human cartilage, but there remains an incomplete understanding of how these properties may be altered by osteoarthritis. The aims of this study were to (1) measure the material constants of normal and osteoarthritic human knee cartilage using isotropic hyperelastic models; (2) determine whether the material constants correlate with histological measures of structure and/or cartilage tissue damage; and (3) quantify the abilities of two common isotropic hyperelastic material models, the neo-Hookean and Yeoh models, to describe articular cartilage contact force, area, and pressure. Small osteochondral specimens of normal and osteoarthritic condition were retrieved from human cadaveric knees and from the knees of patients undergoing total knee arthroplasty and tested in unconfined compression at loading rates and large strains representative of weight-bearing activity. Articular surface contact area and lateral deformation were measured concurrently and specimen-specific finite element models then were used to determine the hyperelastic material constants. Structural parameters were measured using histological techniques while the severity of cartilage damage was quantified using the OARSI grading scale. The hyperelastic material constants correlated significantly with OARSI grade, indicating that the mechanical properties of cartilage for large strains change with tissue damage. The measurements of contact area described anisotropy of the tissue constituting the superficial zone. The Yeoh model described contact force and pressure more accurately than the neo-Hookean model, whereas both models under-predicted contact area and poorly described the anisotropy of cartilage within the superficial zone. These results identify the limits by which isotropic hyperelastic material models may be used to describe cartilage contact variables. This study provides novel data for the mechanical properties of normal and osteoarthritic human articular cartilage and enhances our ability to model this tissue using simple isotropic hyperelastic materials.

In vivo articular cartilage deformation: noninvasive quantification of intratissue strain during joint contact in the human knee

The in vivo measurement of articular cartilage deformation is essential to understand how mechanical forces distribute throughout the healthy tissue and change over time in the pathologic joint. Displacements or strain may serve as a functional imaging biomarker for healthy, diseased, and repaired tissues, but unfortunately intratissue cartilage deformation in vivo is largely unknown. Here, we directly quantified for the first time deformation patterns through the thickness of tibiofemoral articular cartilage in healthy human volunteers. Magnetic resonance imaging acquisitions were synchronized with physiologically relevant compressive loading and used to visualize and measure regional displacement and strain of tibiofemoral articular cartilage in a sagittal plane. We found that compression (of 1/2 body weight) applied at the foot produced a sliding, rigid-body displacement at the tibiofemoral cartilage interface, that loading generated subject- and gender-specific and regionally complex patterns of intratissue strains, and that dominant cartilage strains (approaching 12%) were in shear. Maximum principle and shear strain measures in the tibia were correlated with body mass index. Our MRI-based approach may accelerate the development of regenerative therapies for diseased or damaged cartilage, which is currently limited by the lack of reliable in vivo methods for noninvasive assessment of functional changes following treatment.

Validation of a method for combining biplanar radiography and magnetic resonance imaging to estimate knee cartilage contact

Combining accurate bone kinematics data from biplane radiography with cartilage models from magnetic resonance imaging, it is possible to estimate tibiofemoral cartilage contact area and centroid location. Proper validation of such estimates, however, has not been performed under loading conditions approximating functional tasks, such as gait, squatting, and stair descent. The goal of this study was to perform an in vitro validation to resolve the accuracy of cartilage contact estimations in comparison to a laser scanning gold standard. Results demonstrated acceptable reliability and accuracy for both contact area and centroid location estimates. Root mean square errors in contact area averaged 8.4% and 4.4% of the medial and lateral compartmental areas, respectively. Modified Sorensen-Dice agreement scores of contact regions averaged 0.81 ± 0.07 for medial and 0.83 ± 0.07 for lateral compartments. These validated methods have applications for in vivo assessment of a variety of patient populations and physical activities, and may lead to greater understanding of the relationships between knee cartilage function, effects of joint injury and treatment, and the development of osteoarthritis.

Investigation of the mechanical behavior of kangaroo humeral head cartilage tissue by a porohyperelastic model based on the strain-rate-dependent permeability

Solid-interstitial fluid interaction, which depends on tissue permeability, is significant to the strain-rate-dependent mechanical behavior of humeral head (shoulder) cartilage. Due to anatomical and biomechanical similarities to that of the human shoulder, kangaroos present a suitable animal model. Therefore, indentation experiments were conducted on kangaroo shoulder cartilage tissues from low (10(-4)/s) to moderately high (10(-2)/s) strain-rates. A porohyperelastic model was developed based on the experimental characterization; and a permeability function that takes into account the effect of strain-rate on permeability (strain-rate-dependent permeability) was introduced into the model to investigate the effect of rate-dependent fluid flow on tissue response. The prediction of the model with the strain-rate-dependent permeability was compared with those of the models using constant permeability and strain-dependent permeability. Compared to the model with constant permeability, the models with strain-dependent and strain-rate-dependent permeability were able to better capture the experimental variation at all strain-rates (p < 0.05). Significant differences were not identified between models with strain-dependent and strain-rate-dependent permeability at strain-rate of 5 × 10(-3)/s (p = 0.179). However, at strain-rate of 10(-2)/s, the model with strain-rate-dependent permeability was significantly better at capturing the experimental results (p < 0.005). The findings thus revealed the significance of rate-dependent fluid flow on tissue behavior at large strain-rates, which provides insights into the mechanical deformation mechanisms of cartilage tissues.

An MRI-compatible loading device to assess knee joint cartilage deformation: Effect of preloading and inter-test repeatability

It has been suggested that the extent and location of cartilage deformation within a joint under compressive loading may be predictive of predisposition to further degeneration. To explore this relationship in detail requires the quantification of cartilage deformation under controlled loads on a per-patient basis in a longitudinal manner. Our objectives were (1) to design a device capable of applying controllable axial loads while ensuring repeatable within-patient tibiofemoral positioning during magnetic resonance imaging (MRI) scans and (2) to determine the duration for which load should be maintained prior to the image acquisition, for a reproducible measurement of cartilage deformation, within the restraints of a clinical setting. A displacement control loading device was manufactured from MRI-compatible materials and tested on four volunteers for the following five scans: an unloaded scan, two repeat immediate scans which were started immediately after the application of 50% body weight, and two repeat delayed scans started 12 min after load application. Outcome measures included within-patient changes in tibiofemoral position and cartilage deformation between repeat loaded scans. The differences in tibiofemoral position between repeat loaded scans were <1mm in translation and <2° in rotation. Cartilage deformations were more consistent in the delayed scans compared to the immediate scans. We conclude that our loading device can ensure repeatable tibiofemoral positioning to allow for longitudinal studies, and the delayed scan may enable us to obtain more reproducible measurements of cartilage deformation in a clinical setting.

In vivo cartilage strain increases following medial meniscal tear and correlates with synovial fluid matrix metalloproteinase activity

Meniscal tears are common injuries, and while partial meniscectomy is a frequent treatment option, general meniscus loss is a risk factor for the development of osteoarthritis. The goal of this study was to measure the in vivo tibiofemoral cartilage contact patterns in patients with meniscus tears in relation to biomarkers of cartilage catabolism in the synovial fluid of these joints. A combination of magnetic resonance imaging and biplanar fluoroscopy was used to determine the in vivo motion and cartilage contact mechanics of the knee. Subjects with isolated medial meniscus tears were analyzed while performing a quasi-static lunge, and the contralateral uninjured knee was used as a control. Synovial fluid was collected from the injured knee and matrix metalloproteinase (MMP) activity, sulfated glycosaminoglycan, cartilage oligomeric matrix protein, prostaglandin E2, and the collagen type II cleavage biomarker C2C were measured. Contact strain in the medial compartment increased significantly in the injured knees compared to contralateral control knees. In the lateral compartment, the contact strain in the injured knee was significantly increased only at the maximum flexion angle (105°). The average cartilage strain at maximum flexion positively correlated with total MMP activity in the synovial fluid. These findings show that meniscal injury leads to loss of normal joint function and increased strain of the articular cartilage, which correlated to elevated total MMP activity in the synovial fluid. The increased strain and total MMP activity may reflect, or potentially contribute to, the early development of osteoarthritis that is observed following meniscal injury.

In vivo measurement of localized tibiofemoral cartilage strains in response to dynamic activity

BACKGROUND

Altered local mechanical loading may disrupt normal cartilage homeostasis and play a role in the progression of osteoarthritis. Currently, there are limited data quantifying local cartilage strains in response to dynamic activity in normal or injured knees.

PURPOSE/HYPOTHESIS

To directly measure local tibiofemoral cartilage strains in response to a dynamic hopping activity in normal healthy knees. We hypothesized that local regions of cartilage will exhibit significant compressive strains in response to hopping, while overall compartmental averages may not.

STUDY DESIGN

Controlled laboratory study.

METHODS

Both knees of 8 healthy subjects underwent magnetic resonance imaging before and immediately after a dynamic hopping activity. Images were segmented and then used to create 3-dimensional surface models of bone and cartilage. These pre- and postactivity models were then registered using an iterative closest point technique to enable site-specific measurements of cartilage strain (defined as the normalized change in cartilage thickness before and after activity) on the femur and tibia.

RESULTS

Significant strains were observed in both the medial and lateral tibial cartilage, with each compartment averaging a decrease of 5%. However, these strains varied with location within each compartment, reaching a maximum compressive strain of 8% on the medial plateau and 7% on the lateral plateau. No significant averaged compartmental strains were observed in the medial or lateral femoral cartilage. However, local regions of the medial and lateral femoral cartilage experienced significant compressive strains, reaching maximums of 6% and 3%, respectively.

CONCLUSION

Local regions of both the femur and tibia experienced significant cartilage strains as a result of dynamic activity. An understanding of changes in cartilage strain distributions may help to elucidate the biomechanical factors contributing to cartilage degeneration after joint injury.

CLINICAL RELEVANCE

Site-specific measurements of in vivo cartilage strains are important because altered loading is believed to be a factor contributing to the development and progression of osteoarthritis. Specifically, this methodology and data could be used to evaluate the effects of soft tissue injuries (such as ligament or meniscus tears) on cartilage strains in response to dynamic activities of daily living.

Fiducial marker-based correction for involuntary motion in weight-bearing C-arm CT scanning of knees. II. Experiment

PURPOSE

A C-arm CT system has been shown to be capable of scanning a single cadaver leg under loaded conditions by virtue of its highly flexible acquisition trajectories. In Part I of this study, using the 4D XCAT-based numerical simulation, the authors predicted that the involuntary motion in the lower body of subjects in weight-bearing positions would seriously degrade image quality and the authors suggested three motion compensation methods by which the reconstructions could be corrected to provide diagnostic image quality. Here, the authors demonstrate that a flat-panel angiography system is appropriate for scanning both legs of subjects in vivo under weight-bearing conditions and further evaluate the three motion-correction algorithms using in vivo data.

METHODS

The geometry of a C-arm CT system for a horizontal scan trajectory was calibrated using the PDS-2 phantom. The authors acquired images of two healthy volunteers while lying supine on a table, standing, and squatting at several knee flexion angles. In order to identify the involuntary motion of the lower body, nine 1-mm-diameter tantalum fiducial markers were attached around the knee. The static mean marker position in 3D, a reference for motion compensation, was estimated by back-projecting detected markers in multiple projections using calibrated projection matrices and identifying the intersection points in 3D of the back-projected rays. Motion was corrected using three different methods (described in detail previously): (1) 2D projection shifting, (2) 2D deformable projection warping, and (3) 3D rigid body warping. For quantitative image quality analysis, SSIM indices for the three methods were compared using the supine data as a ground truth.

RESULTS

A 2D Euclidean distance-based metric of subjects' motion ranged from 0.85 mm (±0.49 mm) to 3.82 mm (±2.91 mm) (corresponding to 2.76 to 12.41 pixels) resulting in severe motion artifacts in 3D reconstructions. Shifting in 2D, 2D warping, and 3D warping improved the SSIM in the central slice by 20.22%, 16.83%, and 25.77% in the data with the largest motion among the five datasets (SCAN5); improvement in off-center slices was 18.94%, 29.14%, and 36.08%, respectively.

CONCLUSIONS

The authors showed that C-arm CT control can be implemented for nonstandard horizontal trajectories which enabled us to scan and successfully reconstruct both legs of volunteers in weight-bearing positions. As predicted using theoretical models, the proposed motion correction methods improved image quality by reducing motion artifacts in reconstructions; 3D warping performed better than the 2D methods, especially in off-center slices.

Subject-specific finite element modeling of the tibiofemoral joint based on CT, magnetic resonance imaging and dynamic stereo-radiography data in vivo

In this paper, we present a new methodology for subject-specific finite element modeling of the tibiofemoral joint based on in vivo computed tomography (CT), magnetic resonance imaging (MRI), and dynamic stereo-radiography (DSX) data. We implemented and compared two techniques to incorporate in vivo skeletal kinematics as boundary conditions: one used MRI-measured tibiofemoral kinematics in a nonweight-bearing supine position and allowed five degrees of freedom (excluding flexion-extension) at the joint in response to an axially applied force; the other used DSX-measured tibiofemoral kinematics in a weight-bearing standing position and permitted only axial translation in response to the same force. Verification and comparison of the model predictions employed data from a meniscus transplantation study subject with a meniscectomized and an intact knee. The model-predicted cartilage-cartilage contact areas were examined against "benchmarks" from a novel in situ contact area analysis (ISCAA) in which the intersection volume between nondeformed femoral and tibial cartilage was characterized to determine the contact. The results showed that the DSX-based model predicted contact areas in close alignment with the benchmarks, and outperformed the MRI-based model: the contact centroid predicted by the former was on average 85% closer to the benchmark location. The DSX-based FE model predictions also indicated that the (lateral) meniscectomy increased the contact area in the lateral compartment and increased the maximum contact pressure and maximum compressive stress in both compartments. We discuss the importance of accurate, task-specific skeletal kinematics in subject-specific FE modeling, along with the effects of simplifying assumptions and limitations.

Subject-specific patterns of femur-labrum contact are complex and vary in asymptomatic hips and hips with femoroacetabular impingement

BACKGROUND

Femoroacetabular impingement (FAI) may constrain hip articulation and cause chondrolabral damage, but to our knowledge, in vivo articulation and femur-labrum contact patterns have not been quantified.

PURPOSE

In this exploratory study, we describe the use of high-speed dual-fluoroscopy and model-based tracking to dynamically measure in vivo hip articulation and estimate the location of femur-labrum contact in six asymptomatic hips and three hips with FAI during the impingement examination. We asked: (1) Does femur-labrum contact occur at the terminal position of impingement? (2) Could range of motion (ROM) during the impingement examination appear decreased in hips with FAI? (3) Does the location of femur-labrum contact coincide with that of minimum bone-to-bone distance? (4) In the patients with FAI, does the location of femur-labrum contact qualitatively correspond to the location of damage observed intraoperatively?

METHODS

High-speed dual-fluoroscopy images were acquired continuously as the impingement examination was performed. CT arthrogram images of all subjects were segmented to generate three-dimensional (3-D) surfaces for the pelvis, femur, and labrum. Model-based tracking of the fluoroscopy images enabled dynamic kinematic observation of the 3-D surfaces. At the terminal position of the examination, the region of minimal bone-to-bone distance was compared with the estimated location of femur-labrum contact. Each patient with FAI underwent hip arthroscopy; the location of femur-labrum contact was compared qualitatively with damage found during surgery. As an exploratory study, statistics were not performed.

RESULTS

Femur-labrum contact was observed in both groups, but patterns of contact were subject-specific. At the terminal position of the impingement examination, internal rotation and adduction angles for each of the patients with FAI were less than the 95% confidence intervals (CIs) for the asymptomatic control subjects. The location of minimum bone-to-bone distance agreed with the region of femur-labrum contact in two of nine hips. The locations of chondrolabral damage identified during surgery qualitatively coincided with the region of femur-labrum contact.

CONCLUSIONS

Dual-fluoroscopy and model-based tracking provided the ability to assess hip kinematics in vivo during the entire impingement examination. The high variability in observed labrum-femur contact patterns at the terminal position of the examination provides evidence that subtle anatomic features could dictate underlying hip biomechanics. Although femur-labrum contact occurs in asymptomatic and symptomatic hips at the terminal position of the impingement examination, contact may occur at reduced adduction and internal rotation in patients with FAI. Use of minimum bone-to-bone distance may not appropriately identify the region of femur-labrum contact. Additional research, using a larger cohort and appropriate statistical tests, is required to confirm the findings of this exploratory study.