Our current research aims to develop technologies to predict spinal loads in vivo using a combination of imaging and modelling methods. To ensure the project's success and inform future applications of the technology, we sought to understand the opinions and perspectives of patients and the public. A 90-minute public and patient involvement event was developed in collaboration with Exeter Science Centre and held on World Spine Day 2023. The event involved a brief introduction to the project goals followed by an interactive questionnaire to gauge the participants’ background knowledge and interest. The participants then discussed five topics: communication, future directions of the research, concerns about the research protocol, concerns about data, and interest in the project team and research process. A final questionnaire was used to determine their thoughts about the event.Background
Methods
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).
This study aimed to examine the effect of high tibial osteotomy (HTO) on the ankle and subtalar joints via analysis of static radiographic alignment. We hypothesised that surgical alteration of the alignment of the proximal tibia would result in compensatory distal changes. 35 patients recruited as part of the wider Biomechanics and Bioengineering Centre Versus Arthritis HTO study between 2011 and 2018 had pre- and postoperative full-length weightbearing radiographs taken of their lower limbs. In addition to standard alignment measures of the limb and knee (mechanical tibiofemoral angle, Mikulicz point, medial proximal tibial angle), additional measures were taken of the ankle/subtalar joints (lateral distal tibial angle, ground-talus angle, joint line convergence angle of the ankle) as well as a novel measure of stance width. Results were compared using a paired T-test and Pearson's correlation coefficient. Following HTO, there was a significant (5.4°) change in subtalar alignment. Ground-talus angle appeared related both to the level of malalignment preoperatively and the magnitude of the alignment change caused by the HTO surgery; suggesting subtalar positioning as a key adaptive mechanism. In addition to compensatory changes within the subtalar joints, the patients on average had a 31% wider stance following HTO. These two mechanisms do not appear to be correlated but the morphology of the tibial plafond may influence which compensatory mechanisms are employed by different subgroups of HTO patients. These findings are of vital importance in clinical practice both to anticipate potential changes to the ankle and subtalar joints following HTO but it could also open up wider indications for HTO in the treatment of ankle malalignment and osteoarthritis.
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
Medial knee OA effects approximately 4.1 million people in England. Non-surgical strategies to lower knee joint loading is commonly researched in the knee OA literature as a method to alleviate pain and discomfort. Medial knee OA is much more prevalent than lateral knee OA due to the weight bearing line passing medial to the knee causing an external knee adduction moment (KAM). Numerous potential gait retraining strategies have been proposed to reduce either the first and/or the second peak KAM, including: toe-in gait, toe-out gait, lateral trunk lean and medial thrust gait. Gait retraining has been researched with little regard to the biomechanical consequences at the hip and ankle joints. This systematic review aimed to establish whether gait retraining can reduce medial knee loading as assessed by first and second peak KAMs, establish what are the biomechanical effects a reduced KAM has on other lower limb joint biomechanics and outline patient/participant reported outcomes on how easy the gait retraining style was to implement. The protocol for this systematic review was registered with PROSPERO on the 23rd January 2018 (registration ID: CRD42018085738). 13 databases were searched by one author (J.B.B). Additionally, PROSPERO was searched for ongoing or recently completed systematic reviews. Risk of bias was assessed using the Downs and Black quality index. Search: Group one consisted of keywords “walk” OR “gait”. Keywords “knee” OR “adduction moment” built up the second group. Group three consisted “osteoarthriti” OR “arthriti” OR “osteo arthriti”, OR “OA”. Group four included “hip” OR “ankle”. the searched results of each group were combined with conjunction “AND” in all fields. Out of the eight different gait retraining strategies identified, trunk lean reduced first peak KAM the most, which was evaluated in 3 studies, reducing first peak KAM by 20%-65%. There was a lack of collective pelvic, hip and/or ankle joint biomechanical variables reported across all 11 studies. Of eight gait retraining styles identified, the strategy that reduced first peak KAM the most was an increased lateral trunk lean, which was evaluated in 3 different studies. This is the first systematic review that has highlighted that there is limited evidence of the biomechanical consequences of a reduced knee joint load has on the pelvic, hip and/or ankle joints when undertaking gait retraining protocols. Future studies assessing gait retraining strategies should provide biomechanical outputs for other lower limb joints other than the knee joint, as well as providing participant perceptions on the level of difficulty the gait style is to perform.
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.
Risk factors for poor outcomes after total knee replacement (TKR) have been identified, but the underlying causes are not fully understood. The aim of this research was to establish the relationship between measurable gait parameters and patients' subjective function, pre and post total knee replacement. 25 subjects underwent gait analysis, before and one year following total knee replacement. Patient reported function was investigated using the Activities of Daily Living Scale of the Knee Outcome Survey (KOS). Gait analysis was performed using infrared cameras and reflective marker clusters. Correlation between motion analysis data and patient reported function was investigate. Whilst multiple gait parameters correlated with KOS score preoperatively, there was no correlation after TKR. Three preoperative measurements correlated with the improvement in score a subject achieved following surgery: These were preoperative rate of extension in swing, total range of flexion from heel strike and time point of maximum stance extension. Our results suggest that whilst preoperatively there is a close relationship between knee biomechanical function and patient reported function, after TKR factors other than biomechanical function determine patient outcomes.
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.
