Biomechanics is an essential form of measurement in the understanding of the development and progression of osteoarthritis (OA). However, the number of participants in biomechanical studies are often small and there is limited ways to share or combine data from across institutions or studies. This is essential for applying modern machine learning methods, where large, complex datasets can be used to identify patterns in the data. Using these data-driven approaches, it could be possible to better predict the optimal interventions for patients at an early stage, potentially avoiding pain and inappropriate surgery or rehabilitation. In this project we developed a prototype database platform for combining and sharing biomechanics datasets. The database includes methods for importing and standardising data and associated variables, to create a seamless, searchable combined dataset of both healthy and knee OA biomechanics. Data was curated through calls to members of the OATech Network+ (Abstract
Objectives
Methods
Spinal disorders such as back pain incur a substantial societal and economic burden. Unfortunately, there is lack of understanding and treatment of these disorders are further impeded by the inability to assess spinal forces in vivo. The aim of this project is to address this challenge by developing and testing a novel image-driven approach that will assess the forces in an individual's spine in vivo by incorporating information acquired from multimodal imaging (magnetic resonance imaging (MRI) and biplane X-rays) in a subject-specific model. Magnetic resonance and biplane X-ray imaging are used to capture information about the anatomy, tissues, and motion of an individual's spine as they perform a range of everyday activities. This information is then utilised in a subject-specific computational model based on the finite element method to predict the forces in their spine. The project is also utilising novel machine learning algorithms and in vitro, six-axis mechanical testing on human, porcine and bovine samples to develop and test the modelling methods rigorously.Abstract
Objectives
Methods
Knee varus malalignment increases medial knee compartment loading and is associated with knee osteoarthritis (OA) progression and severity1. Altered biomechanical loading and dysregulation of joint tissue biology drive OA progression, but mechanistic links between these factors are lacking. Subchondral bone structural changes are biomechanically driven, involve bone resorption, immune cell influx, angiogenesis, and sensory nerve invasion, and contribute to joint destruction and pain2. We have investigated mechanisms underlying this involving RANKL and alkaline phosphatase (ALP), which reflect bone resorption and mineralisation respectively3 and the axonal guidance factor Sema3A. Sema3A is osteotropic, expressed by mechanically sensitive osteocytes, and an inhibitor of sensory nerve, blood vessel and immune cell invasion4. Sema3A is also differentially expressed in human OA bone5.HYPOTHESIS: Medial knee compartment overloading in varus knee malalignment patients causes dysregulation of bone derived Sema3A signalling directly linking joint biomechanics to pathology and pain. Synovial fluid obtained from 30 subjects with medial knee OA (KL grade II-IV) undergoing high tibial osteotomy surgery (HTO) was analysed by mesoscale discovery and ELISA analysis for inflammatory, neural and bone turnover markers. 11 of these patients had been previously analysed in a published patient-specific musculoskeletal model6 of gait estimating joint contact location, pressure, forces, and medial-lateral condyle load distribution in a published data set included in analyses. Data analysis was performed using Pearson's correlation matrices and principal component analyses. Principal Components (PCs) with eigenvalues greater than 1 were analysed.Abstract
OBJECTIVE
METHODS
Investigate Magnetic Resonance Imaging (MRI) as an alternative to Computerised Tomography (CT) when calculating kinematics using Biplane Video X-ray (BVX) by quantifying the accuracy of a combined MRI-BVX methodology by comparing with results from a gold-standard bead-based method. Written informed consent was given by one participant who had four tantalum beads implanted into their distal femur and proximal tibia from a previous study. Three-dimensional (3D) models of the femur and tibia were segmented (Simpleware Scan IP, Synopsis) from an MRI scan (Magnetom 3T Prisma, Siemens). Anatomical Coordinate Systems (ACS) were applied to the bone models using automated algorithms1. The beads were segmented from a previous CT and co-registered with the MRI bone models to calculate their positions. BVX (60 FPS, 1.25 ms pulse width) was recorded whilst the participant performed a lunge. The beads were tracked, and the ACS position of the femur and tibia were calculated at each frame (DSX Suite, C-Motion Inc.). The beads were digitally removed from the X-rays (MATLAB, MathWorks) allowing for blinded image-registration of the MRI models to the radiographs. The mean difference and standard deviation (STD) between bead-generated and image-registered bone poses were calculated for all degrees of freedom (DOF) for both bones. Using the principles defined by Grood and Suntay2, 6 DOF kinematics of the tibiofemoral joint were calculated (MATLAB, MathWorks). The mean difference and STD between these two sets of kinematics were calculated.Abstract
Objectives
Methods
Mechanical loading of joints with osteoarthritis (OA) results in pain-related functional impairment, altered joint mechanics and physiological nociceptor interactions leading to an experience of pain. However, the current tools to measure this are largely patient reported subjective impressions of a nociceptive impact. A direct measure of nociception may offer a more objective indicator. Specifically, movement-induced physiological responses to nociception may offer a useful way to monitor knee OA. In this study, we gathered preliminary data on healthy volunteers to analyse whether integrated biomechanical and physiological sensor datasets could display linked and quantifiable information to a nociceptive stimulus. Following ethical approval, 15 healthy volunteers completed 5 movement and stationary activities in 2 conditions; a control setting and then repeated with an applied quantified thermal pain stimulus to their right knee. An inertial measurement unit (IMU) and an electromyography (EMG) lower body marker set were tested and integrated with ground reaction force (GRF) data collection. Galvanic skin response electrodes for skin temperature and conductivity and photoplethysmography (PPG) sensors were manually timestamped to the integrated system. Pilot data showed EMG, GRF and IMU fluctuations within 0.5 seconds of each other in response to a thermal trigger. Preliminary analysis on the 15 participants tested has shown skin conductance, PPG, EMG, GRFs, joint angles and kinematics with varying increases and fluctuations during the thermal condition in comparison to the control condition. Preliminary results suggest physiological and biomechanical data outputs can be linked and identified in response to a defined nociceptive stimulus. Study data is currently founded on healthy volunteers as a proof-of-concept. Further exploratory statistical and sensor readout pattern analysis, alongside early and late-stage OA patient data collection, can provide the information for potential development of wearable nociceptive sensors to measure disease progression and treatment effectiveness.
Biplane video X-ray (BVX) – with models segmented from magnetic resonance imaging (MRI) – is used to directly track bones during dynamic activities. Investigating tibiofemoral kinematics helps to understand effects of disease, injury, and possible interventions. Develop a protocol and compare in-vivo kinematics during loaded dynamic activities using BVX and MRI. BVX (60 FPS) was captured whilst three healthy volunteers performed three repeats of lunge, stair ascent and gait. MRI scans were performed (Magnetom 3T Prisma, Siemens). 3D bone models of the tibia and femur were segmented (Simpleware Scan IP, Synopsis). Bone poses were obtained by manually matching bone models to X-rays (DSX Suite, C-Motion Inc.). Mean range of motion (ROM) of the contact points on the medial and lateral tibial plateau were calculated using custom MATLAB code (MathWorks). Results were filtered using an adaptive low pass Butterworth filter (Frequency range: 5-29Hz). Gait and Stair ascent activities from one participant's data showed increased ROM for medial-lateral (ML) translation in the medial compartment but decreased ROM in anterior-posterior (AP) translation when comparing against the same translations on the lateral compartment of the tibial plateau. Lunge activity showed increased ROM for both ML and AP translation in the medial compartment when compared with the lateral compartment. These results highlight the variability in condylar translations between different activities. Understanding healthy in-vivo kinematics across different activities allows the determination of suitable activities to best investigate the kinematic changes due to disease or injury and assess the efficacy of different interventions. Acknowledgements: This research was supported by the Engineering and Physical Sciences Research Council (EPSRC) doctoral training grant (EP/T517951/1).
Principal Component Analysis (PCA) is a useful method for analysing human motion data. The objective of this study was to use PCA to quantify the biggest variance in knee kinematics waveforms between a Non-Pathological (NP) group and individuals awaiting High Tibial Osteotomy (HTO) surgery. Thirty knees (29 participants) who were scheduled for HTO surgery were included in this study. Twenty-eight NP volunteers were recruited into the study. Human motion analysis was performed during level gait using a modified Cleveland marker set. Subjects walked at their self-selected speed for a minimum of 6 successful trials. Knee kinematics were calculated within Visual3D (C-Motion). The first three Principal Components (PCs) of each input variable were selected. Single-component reconstruction was performed alongside representative extremes of each PC to aid interpretation of the biomechanical feature reconstructed by each component.Abstract
Objectives
Methods
Current tools to measure pain are broadly subjective impressions of the impact of the nociceptive impulse felt by the patient. A direct measure of nociception may offer a more objective indicator. Specifically, movement-induced physiological responses to nociception may offer a useful way to monitor knee OA. In this proof-of-concept study, we evaluated whether integrated biomechanical and physiological sensor datasets could display linked and quantifiable information to a nociceptive stimulus. Following ethical approval, we applied a quantified thermal pain stimulus to a volunteer during stationary standing in a gait lab setting. An inertial measurement unit (IMU) and an electromyography (EMG) lower body marker set were tested and integrated with ground reaction force (GRF) data collection. Galvanic skin response electrodes and skin thermal sensors were manually timestamp linked to the integrated system.Abstract
Objectives
Method
Optical motion capture (OMC) is the current gold standard for motion analysis, however measuring patellofemoral kinematics is not possible using the technique. One approach to measuring in-vivo kinematics is to use biplane video X-ray (BVX) and 3D models generated from MRI to track the movement of the patellar. Understanding how the patellar is moving during different loaded dynamic activities can help with understanding the effects of different interventions when treating disease or injury. To develop a protocol and compare patellofemoral kinematics for different activities using biplane video X-ray (BVX) Two healthy volunteers performed level walk, lunge, and stair ascent activities while simultaneous capturing BVX and synchronised OMC. Participants undertook MR imaging (Magnetom 3T Prisma, Siemens) which was manually segmented into 3D bone models (Simpleware Scan IP, Synopsis). Bone position and orientation for the patellar and femur were calculated by manual matching of 3D Bone models to X-Rays (DSX Suite, C-Motion, Inc.). Patellofemoral kinematics were calculated using Visual 3D (C-Motion, Inc.).Abstract
Objective
Methods
Valgus high tibial osteotomy (HTO) represents an effective treatment for patients with medial compartment osteoarthritis (OA) in a varus knee. However, the mechanisms which cause this clinical improvement are unclear. Previous studies suggest a wider stance gait can reduce medial compartment loading via reduction in the external knee adduction moment (KAM); a measure implicated in progression of medial compartment OA. This study aimed to measure whether valgus HTO is associated with a postoperative increase in static stance width. 32 patients, recruited in the Biomechanics and Bioengineering Centre Versus Arthritis HTO study, underwent valgus (medial opening wedge) HTO. Weightbearing pre- and post- operative radiographs were taken showing both lower limbs. The horizontal distance, measured from a fixed point on the right talus to the corresponding point on the left, was divided by the talus width to give a standardised “stance width” for each radiograph. The difference between pre- and post- operative stance width was compared for each patient using a paired sample t-test.Abstract
OBJECTIVES
METHODS
Explore whether high tibial osteotomy (HTO) changes knee contact forces and to explore the relationship between the external knee adduction moment (EKAM) pre and 12 months post HTO. Three-dimensional gait analysis was performed on 17 patients pre and 12-months post HTO using a modified Cleveland marker-set. Tibiofemoral contact forces were calculated in SIMM. The scaled musculoskeletal model integrated an extended knee model allowing for 6 degrees of freedom in the tibiofemoral and patellofemoral joint. Joint angles were calculated using inverse kinematics then muscle and contact forces and secondary knee kinematics were estimated using the COMAC algorithm. Paired samples t-test were performed using SPSS version 25 (SPSS Inc., USA). Testing for normality was undertaken with Shapiro-Wilk. Pearson correlations established the relationships between EKAM1 to medial KCF1, and EKAM2 to medial KCF2, pre and post HTO.Abstract
Objective
Methods
Skeletal kinematics are traditionally measured by motion analysis methods such as optical motion capture (OMC). While easy to carry out and clinically relevant for certain applications, it is not suitable for analysing the ankle joint due to its anatomical complexity. A greater understanding of the function of healthy ankle joints could lead to an improvement in the success of ankle-replacement surgeries. Biplane video X-ray (BVX) is a technique that allows direct measurement of individual bones using highspeed, dynamic X-Rays. To develop a protocol to quantify in-vivo foot and ankle kinematics using a bespoke High-speed Dynamic Biplane X-ray system combined with OMC. Two healthy volunteers performed five level walks and step-down trials while simultaneous capturing BVX and synchronised OMC. participants undertook MR imaging (Magnetom 3T Prisma, Siemens) which was manually segmented into 3D bone models (Simpleware Scan IP, Synopsis). Bone position and orientation for the Talus, Tibia and Calcaneus were calculated by manual matching of 3D Bone models to X-Rays (DSX Suite, C-Motion, Inc.). OMC markers were tracked (QTM, Qualisys) and processed using Visual 3D (C-motion, Inc.).Abstract
Objective
Methods
One of the main surgical goals when performing a total knee replacement (TKR) is to ensure the implants are properly aligned and correctly sized; however, understanding the effect of alignment and rotation on the biomechanics of the knee during functional activities is limited. Cardiff University has unique access to a group of local patients who have relatively high frequency of poor alignment, and early failure. This provides a rare insight into how malalignment of TKR's can affect patients from a clinical and biomechanical point of view to determine how to best align a TKR. This study aims to explore relationship clinical surgical measurements of Implant alignment with in-vivo joint kinematics. 28 patient volunteers (with 32 Kinemax (Stryker) TKR's were recruited. Patients undertook single plane video fluoroscopy of the knee during a step-up and step-down task to determine TKR in-vivo kinematics and centre of rotation (COR). Joint Track image registration software (University of Florida, USA) was used to match CAD models of the implant to the x-ray images. Hip-Knee-Ankle (HKA) was measured using long-leg radiographs to determine frontal plane alignment. Posterior tibial slope angle was calculated using radiographs. An independent sample t-test was used to explore differences between neutral (HKA:-2° to 2°), varus (≥2°) and valgus alignment (≤-2°) groups. Other measures were explored across the whole cohort using Pearson's correlations (SPSS V23). There was found to be no statistical difference between groups or correlations for HKA. The exploratory analysis found that tibial slope correlated with Superior/Inferior translation ROM during step up (r=−0.601, p<0.001) and step down (r=−.512, p=0.03) the position of the COR heading towards the lateral (r=−.479, p=0.006) during step down. Initial results suggest no relationship between frontal plane alignment and in-vivo. Exploratory analyses have found other relationships that are worthy of further research and may be important in optimizing function.
Whilst home-based exercise rehabilitation plays a key role in determining patient outcomes following orthopaedic intervention (e.g. total knee replacement), it is very challenging for clinicians to objectively monitor patient progress, attribute functional improvement (or lack of) to adherence/non-adherence and ultimately prescribe personalised interventions. This research aimed to identify whether 4 knee rehabilitation exercises could be objectively distinguished from each other using lower body inertial measurement units (IMUs) and principle components analysis (PCA) in the hope to facilitate objective home monitoring of exercise rehabilitation. 5 healthy participants performed 4 repetitions of 4 exercises (knee flexion in sitting, knee extension, single leg step down and sit to stand) whilst wearing lower body IMU sensors (Xsens, Holland; sampling at 60 Hz). Anthropometric measurements and a static calibration were combined to create the biomechanical model, with 3D hip, knee and ankle angles computed using the Euler sequence ZXY. PCA was performed on time normalised (101 points) 3D joint angle data which reduced all joint angle waveforms into new uncorrelated PCs via an orthogonal transformation. Scatterplots of PC1 versus PC2 were used to visually inspect for clustering between the PC values for the 4 exercises. A one-way ANOVA was performed on the first 3 PC values for the 9 variables under analysis. Games-Howell post hoc tests identified variables that were significantly different between exercises. All exercises were clearly distinguishable using the PC scatterplot representing hip flexion-extension waveforms. ANOVA results revealed that PC1 for the knee flexion angle waveform was the only PC value statistically different across all exercises. Findings demonstrate clear potential to objectively distinguish between different knee rehabilitation exercises using IMU sensors and PCA. Flexion-extension angles at the hip and knee appear most suited for accurate separation, which will be further investigated on patient data and additional exercises.
Healthy cartilage is essential for optimal joint function. Although, articular cartilage defects are highly prevalent in the active population and might hamper joint function, the effect of articular cartilage defects on knee contact forces and pressures is not yet documented. Therefore, the present study compared knee contact forces and pressures between patients with a tibiofemoral cartilage defect and healthy controls. This might provide additional insights in movement adaptations and the role of altered loading in the progression from defect to OA. Experimental gait data was collected in 15 patients with isolated articular cartilage defects (8 medial-affected, 7 lateral-affected) and 19 healthy asymptomatic controls and was processed using a musculoskeletal model to calculate contact forces and pressures. Differences between medial-affected, lateral-affected and controls were evaluated using Kruskal-Wallis tests and individually compared using Mann-Whitney-U tests (alpha <0.05). The lateral-affected group walked significantly slower compared to the healthy controls. No adaptations in the movement pattern that resulted in decreased loading on the injured condyle were observed. Additionally, the location of loading was not significantly affected. The current results suggest that isolated cartilage defects do not induce changes in the knee joint loading pattern. Consequently, the involved condyle will be equally loaded, indicating that a similar amount of force should be distributed over the remaining cartilage surrounding the articular cartilage defect and may cause local degenerative changes in the cartilage. This in combination with inflammatory responses might play a key role in the progression from articular cartilage defect to a more severe OA phenotype.
The aim of this study was to perform a comprehensive evaluation of the changes in function from pre- to post-surgery in total and unilateral knee arthroplasty (UKA/TKA) patients. Twenty healthy (age 62.4 ±5.9, 11 male), 14 UKA (age 60.9 ±10.1, 8 male) and 17 TKA (age 67.2 ±8.1, 9 male) patients were studied. KA patients were assessed four weeks pre- and six months post-operation. Measures of perceived pain and function were collected using Oxford Knee Score (OKS) questionnaire. Tests of objective function included joint range of motion (RoM), ultrasound imaging, and 3-D motion analysis/inverse modelling from gait and sit-stand. An optimal set of variables was used to classify KA function using the Cardiff DST method. Pre-KA and healthy individuals were accurately classified (96%). Post-operation questionnaire measures of function improved for both UKA and TKA groups. However, observed measures of RoM, muscle atrophy and gait had only limited gains. This resulted in 57% of UKA and only 27% of TKA patients being classified as healthy post-operation. The results of this study show that 6 months post-surgery UKA patients had higher function than TKA. Using statistical approaches to combine functional assessments has provided an accurate platform to classify function and estimate changes from pre- to post-surgery. The clinical application of this tool requires further investigation and comparison to commonly used clinical techniques.
Assessing the efficacy of cervical orthoses in restricting spinal motion has historically proved challenging due to a poor understanding of spinal kinematics and the difficulty in accurately measuring spinal motion. This study is the first to use an 8 camera optoelectronic, passive marker, motion analysis system with a novel marker protocol to compare the effectiveness of the Aspen, Aspen Vista, Philadelphia, Miami-J and Miami-J Advanced collars. Restriction of cervical spine motion was assessed for physiological and functional range of motion (ROM). Nineteen healthy volunteers (12 female, 7 male) were fitted with collars by an approved physiotherapist. ProReflex (Qualisys, Sweden) infra-red cameras were used to track the movement of retro-reflective marker clusters attached to the head and trunk. 3-D kinematic data was collected from uncollared and collared subjects during forward flexion, extension, lateral bending and axial rotation for physiological ROM and during five activities of daily living (ADLs). ROM in the three clinical planes was analysed using the Qualisys Track Manager (Qualisys, Sweden) 6 Degree of Freedom calculation to determine head orientation relative to the trunk. For physiological ROM, the Aspen and Philadelphia were more effective at restricting flexion/extension than the Vista (p<0.001), Miami-J (p<0.001 and p<0.01) and Miami-J Advanced (p<0.01 and p<0.05). The Aspen was more effective at restricting rotation compared to the Vista (p<0.001) and Miami-J (p<0.05). The Vista was least effective at restricting lateral bending (p<0.001). Through functional ROM, the Vista was less effective than the Aspen (p<0.001) and other collars (p<0.01) at restricting flexion/extension. The Aspen and Miami-J Advanced were more effective at restricting rotation than the Vista (p<0.01 and p<0.05) and Miami-J (p<0.05). All the collars were comparable when restricting lateral bending. The Aspen is superior to, and the Aspen Vista inferior to, the other collars at restricting cervical spine motion through physiological ROM. Functional ROM observed during ADLs are less than those observed through physiological ROM. The Aspen Vista is inferior to the other collars at restricting motion through functional ROM. The Aspen collar again performs well, particularly at restricting rotation, but is otherwise comparable to the other collars at restricting motion through functional ranges.
High tibial Osteotomy (HTO) realigns the forces in the knee to slow the progression of osteoarthritis. This study relates the changes in knee joint biomechanics during level gait to glutamate signalling in the subchondral bone of patients pre and post HTO. Glutamate transmits mechanical signals in bone and activates glutamate receptors to influence inflammation, degeneration and nociception in arthritic joints. Thus glutamate signalling is a mechanism whereby mechanical load can directly modulate joint pathology and pain. 3D motion analysis was used to assess level gait prior to HTO (n=5) and postoperatively (n=2). A biomechanical model of each subject was created in Visual3D (C-motion. Inc) and used for biomechanical analysis. Gene expression was analysed by RT-PCR from bone cores from anterior and posterior drill holes, subdivided according to medial or lateral proximal tibia from HTO patients (n=5).BACKGROUND
METHODS
Motion analysis is routinely used in the clinical and research sectors to quantify joint biomechanics. It plays an important role in clinical assessments by aiding the physician to distinguish between primary movement abnormalities and any secondary compensatory mechanisms that may overshadow the cause of the problem. During a data collection session, a wealth of biomechanical data regarding joint and segment kinematics and kinetics are collected from patients performing daily activities. Objective classification can be used to automate a diagnosis from this data and has been used previously to analyse measurements of level gait [1]. It is of interest to assess the knee during stair-gait as this activity involves greater range of motion (ROM) of the lower limbs, larger forces and moments acting at the knee. The aim of the current study is to explore the use of an objective classifier [1] to characterise knee osteoarthritis (OA) and monitor functional recovery following a total knee replacement (TKR) using measurements from stair-gait.INTRODUCTION
AIM
Patient specific knee modelling has the potential to help understand the development of the mechanically induced degenerative disease, Osteoarthritis. A full joint contact model of the knee involves modelling the bones, ligaments, articular cartilage (AC) and meniscus, as well as, the kinematics and geometry of real joints. These finite element models will inevitably require great computational resource to run and it is desirable to find resource effective material model formulations which can accurately describe the mechanical behaviour of the soft tissues. Biphasic models (BIMs) have long been established as an effective formulation for modelling AC. However, the swelling behaviour caused by changes in the ionic phase is a major recovery mechanism and is neglected in the BIMs. It is therefore believed that BIMs alone are insufficient to fully describe the mechanical behaviour of AC. Instead, a thermal analogy method which is generically a BIM that includes the swelling behaviour has been thought to be suitable and has been validated against literature data using material parameters optimized to match the numerical and experimental results. To ensure the model is suitable for patient specific modelling where it will have the ability to reflect the individual AC material properties of the patients in the mechanical behaviour it predicts, two experiments have been planned and are currently being carried out using bovine AC. The first experiment is to investigate the diffusivity of the tissue in solutions of different molarity by measuring the change in tissue weight over time. Eleven explants are taken from the same bovine articular joint using a 6mm biopsy punch and are left in 10mM of PBS overnight to ensure ionic equilibrium has been reached before experiments are carried out. The explants are then placed in PBS solutions of molarities ranging from 0mM to 10mM and weighed at regular time intervals. In the final stage, the explants are then lyophilized and weighed for determining the volume of water in the tissues. Using Archimedes principle, the change in porosity of the tissue is found. A preliminary study has shown that explants submerged in a solution of 5mM has an approximately 4% change in weight after the first 24h and a further 1.73% change in the following 24h. Control specimens left in a solution of 10mM had a 0% change in weight. The second experiment is to carry out mechanical loading on the AC specimens while submerged in a solution of different ion concentrations. Experiments with various loading conditions are being investigated to explore their efficacy for validation. Preliminary compression tests have been carried out where steps of 1% strain was applied, giving a total of 10% strain. Between each step, strain was held constant until full relaxation has been achieved. The reaction force measured from the second experiment in conjunction with data collected from the first experiment will be compared to results predicted in the numerical model. This will allow the determination of whether thermal analogy is adequate or whether more complex triphasic models need to be considered. Furthermore, the development of these experimental methods will contribute to the validation of other AC material models in the future.
