Tibiofemoral Joint Forces during the Stance Phase of Gait after ACL Reconstruction ()
1. Introduction
The main goals of anterior cruciate ligament reconstruction (ACLR) are to restore knee function and prevent development of osteoarthritis (OA). However, the incidence of early-onset OA remains higher in patients following ACLR [1]. It has been reported that abnormal motion and/or moments at the knee could lead to the development of early OA following ACLR [2,3]. To our knowledge, the tibiofemoral joint (TFJ) force during gait has not been studied as a possible contributory factor for OA development after ACLR.
Previous studies have reported TFJ loads for normal healthy individuals during gait [4,5]. Gardinier et al. [6] reported altered TFJ loading during gait in subjects who had sustained a complete ACL rupture within the past seven months. They found that patients walked with decreased force on their ACL deficient knee as compared to their contralateral knee. These results are contrary to a recent report that females who have undergone ACLR had increased tibiofemoral compressive forces during a single-leg drop-land task [7]. To our knowledge, there are no reports describing TFJ loading during gait in subjects following ACLR.
We use principal component analysis (PCA) to explore differences between ACLR subjects’ and control subjects’ TFJ loads during the stance phase of gait. PCA is a dimension reduction method that transforms a set of observations of chosen variables that may be correlated with one another into a set of uncorrelated variables, called principal components [8,9]. PCA has previously been applied to gait waveform data to characterize differences among patients with a specific pathology and normal control subjects [10-13]. PCA allows consideration of the entire waveforms in the analysis as opposed to extracting arbitrary parameters.
The purpose of this study was to compare the TFJ forces and muscle forces of ACL reconstructed knees to those of BMI-matched control subjects during the stance phase of gait. We used the AnyBody Modeling System (Version 5.2, AnyBody Technology A/S, Aalborg, Denmark) to scale a musculoskeletal model to each individual and estimate tibiofemoral and muscle forces based on each individual’s gait data. We hypothesized that the use of PCA would allow us to characterize alterations in three-dimensional TFJ loads and muscle forces after ACL reconstruction as compared to a healthy control population.
2. Methods
Sixteen subjects, eight (three male, five female) ACLR subjects (mean BMI 25.5, SD 4.1 kg/m2) and eight BMImatched (five male, three female) control subjects with no history of lower extremity trauma were recruited (mean BMI 24.3, SD 4.3 kg/m2) following institutional review board approval and informed consent. The ACLR subjects had sustained a unilateral ACL tear and seven were surgically reconstructed with patellar tendon grafts and one with a hamstring tendon graft more than seven months prior to testing (average 93 months). All subjects had a normal contralateral knee.
Prior to testing, retro-reflective markers were placed over bony landmarks including the sacrum, C7, and bilateral acromion, ASIS, PSIS, medial and lateral epicondyles, and medial and lateral malleoli, calcaneous, head of the fifth metatarsal, and dorsum of the foot. Arrays of markers were attached to the thighs and shanks using elastic wrap. A nine camera video-based opto-electronic system (Qualisys AB, Gothenburg, Sweden) was used for 3D motion capture as subjects walked barefoot at a selfselected speed on a 10 meter walkway instrumented with three force plates (AMTI, Watertown, MA, USA). Data were recorded for multiple trials to ensure that each foot made at least three clean footfalls on a force plate. Electromyogram (EMG) electrodes (Trigno, Delsys Inc., Boston, MA, USA) were placed bilaterally over the muscle belly of the rectus femoris, vastus medialis, biceps femoris, semitendinosus, and medial gastrocnemius. EMG data were collected at 2000 Hz, rectified, low-pass filtered, and normalized by the maximum value during the stance phase.
Motion capture and force plate data were imported into the AnyBody Modeling System to estimate the TFJ forces for each subject [14] using the Twente Lower Extremity Model [15]. The model consists of 12 body segments: HAT (head, arms, and trunk), pelvis, and right and left femur, patella, tibia, talus, and foot. The model contains 11 joints: L5S1 and left and right hip, knee, patella/femur, talocrural and subtalar. The L5S1 and hip joints were modeled as a ball-and-socket, and the knee, talocrural, and subtalar joints were defined as a hinge. The patella could rotate with respect to the femur, but the orientation and position of the patella was dependent upon the knee flexion angle. The orientation and position of the pelvis with respect to the 3D global reference frame along with the joint rotations resulted in a model with 21 degrees of freedom. Each leg contained 56 muscles whose mechanical effect was modeled by 159 simple muscle slips, each consisting of a contractile element [15].
The model was scaled in order to match each subject’s anthropometry using a static, standing reference trial. The model was morphed using radial basis functions to match the assumed bony landmarks based on the marker positions [16]. An anthropometric data set [15] was used to model mass, inertia points, and muscle sites/geometry for all segments. The muscle attachment sites and geometries were scaled using a linear geometry scaling law. The muscle strength was scaled according to a length, mass, fat scaling which takes body mass index into account [17].
Inverse dynamics was performed and muscle forces were distributed by using a cubic polynomial optimization scheme that minimizes the sum of the cubes of muscle activations (force/maximum force) at each time step. The compression-distraction (CD), antero-posterior (AP), and medio-lateral (ML) TFJ force (in the tibial reference frame) during stance was averaged for each subject for three trials of gait for the ACLR knee. The same procedure was performed for the corresponding knees (right/ left) of BMI-matched control subjects. All gait waveforms were resampled to 101 values corresponding to 100% of the stance phase of gait (approximately 0% - 60% of the gait cycle). For each subject, force data were averaged from three trials for analysis.
We applied PCA to the averaged gait waveforms of control data to develop principal component (PC) models [10] for each measure, CD force, AP force, ML force, and gastrocnemius, vastii, rectus femoris, and hamstring muscle forces. Each PC model can be considered a projection of the data from the p-dimensional space defined by the original variables (101) to a k-dimensional hyperplane defined by the principal components, where k < p [9,10]. The principal components are optimal in that they explain the maximal amount of variance in the original variables.
The principal component models were then used to reconstruct the original waveform of each control subject’s knee, and the residual was calculated as the difference between the original and reconstructed waveforms. The sum of squares of the residuals, Q, was then calculated as a measure of the perpendicular distance of each knee from the hyper-plane defined by the PC model [10]. Similarly, the Mahalanobis distance, T2, is a measure of the distance of each knee from the center of the hyperplane. The Mahalanobis distance is a weighted sum of squares of the PC scores [10]. Upper limits for each of these measures, T2 and Q, were derived from the normal subjects’ data and used as a reference for comparing the ACLR data. Lower limits were not derived since both measures are squared quantities. Finally, the PC models developed from the control data were applied to the corresponding ACLR subject data, and values of T2 and Q were calculated for each ACLR knee for each gait measure. These values were compared to the control limits to determine deviation from normal.
3. Results
PC models were developed for the CD, AP, ML, gastrocnemius, vastii, rectus femoris, and hamstring forces (Table 1). The number of PCs used in each model (k) was chosen through an 85% trace criterion [9]. Of the eight ACLR knees analyzed, five showed deviation in at least one of the TFJ forces (Table 2). In four of these five subjects, the abnormal TFJ force was accompanied by an alteration in one of the major muscles crossing the knee joint (hamstrings, quadriceps, and/or gastrocnemius).
The average tibiofemoral compressive force for both the ACLR and control groups was biphasic in shape, with the largest peak occurring during terminal stance and a lesser peak occurring at the end of the loading response (Figure 1(a)). Four ACLR subjects were identified as falling outside of the normal range for the Q value. The corresponding waveforms for these subjects show deviations from normal in magnitude (both increased and decreased), shape, and phase shift (Figure 1(b)).
The ML TFJ force was also biphasic in shape with peaks near the end of loading response and during terminal stance (Figure 2(a)). Again, four of the ACLR subjects had Q values which fell outside of the range defined by control subjects. As seen in Figure 2(b), the ACLR subjects deviated from normal in a variety of ways with no noticeable trend.
The mean AP TFJ force for both groups is shown in Figure 3. None of the ACL subjects had T2 or Q values above the normal range of control subjects.
Of the eight ACLR knees, four displayed an abnormal TFJ compressive force during stance. In three of these four ACLR knees that displayed abnormal compressive forces, one of the major muscles/muscle groups crossing the knee also deviated from the control group. Subject ACL4 displayed a TFJ compression force with a reduced magnitude from approximately 10% - 70% of the stance phase with a slightly delayed second peak (Figure 4(e)). When analyzing the corresponding muscle forces, the force in the vastii muscle group had a decreased magni-
Table 1. Principal component models.
(a)(b)
Figure 1. (a) Mean TFJ compression force for the control (solid) and ACLR (dashed) groups; (b) Four ACLR subjects with abnormal compression force (dashed) and mean for control group (solid).
(a)(b)
Figure 2. (a) Mean TFJ mediolateral force for the control (solid) and ACLR (dashed) groups; (b) Four ACLR subjects with abnormal mediolateral force (dashed) and mean for control group (solid).
tude from approximately 10% - 65% of the stance phase (Figure 4(b)) which was likely the cause of the reduced magnitude in the compression force. Also, the gastrocnemius force was delayed in time and slightly decreased in magnitude compared to controls (Figure 4(c)). The delay in gastrocnemius activation explains the delayed
Figure 3. Mean TFJ antero-posterior force for the control (solid) and ACLR (dashed) groups.
second peak in the compression force.
Subject ACL7 displayed an abnormal TFJ compressive force which had increased first and second peaks during stance (Figure 5(e)). The rectus femoris muscle force was also found to lie outside the bounds of the normal control subjects and appears increased in magnitude from approximately 10% - 85% of stance (Figure 5(d)).
Subject ACL8 displayed a gastrocnemius force increased in magnitude at approximately 70% stance (Figure 6(c)) which is likely the cause for the increased TFJ compressive force which also occurs at 70% stance (Figure 6(e)). This subject’s hamstring force was also considered abnormal when compared to control subjects (Figure 6(a)).
Finally, subject ACL9 had a compressive tibiofemoral force which deviated from that of the control group (Figure 7(e)), but none of this subject’s individual muscle
forces could explain the deviation in the compressive force (Figures 6(a)-(d)).
For the subjects who displayed abnormal TFJ forces, EMG data is plotted along with the muscle force to assess if the muscles in the model are active in the time that would be expected as determined from EMG. In cases where a muscle group force is shown, a representative muscle from that group was chosen from the available EMG data. EMG from the semitendinosus is plotted along with the hamstring muscle force, and the vastusmedialis EMG is plotted for the vastii muscle group. For the gastrocnemius force, EMG data was collected from the medial gastrocnemius. The EMG signal was not plotted if the signal to noise ratio was not at least four. Control subject EMG data for the same muscles is plotted below each graph to show the timing of normal muscle activity during the stance phase [18].
For the four ACL subjects shown, the predicted hamstrings muscle force activation coincides well with the timing of the semitendinosus EMG. The predicted vastii muscle force timing coincides well with EMG during the first half of stance phase, but some discrepancies were noted during the latter portion of stance in subjects ACL4 (Figure 4(b)) and ACL8 (Figure 6(b)). Subject ACL7 displayed rectus femoris EMG activity during the early part of stance, but the predicted muscle activity occurred predominantly later in stance (Figure 5(d)). Subject ACL8 displayed rectus femoris EMG activity both during early stance (0% - 35%) and later in stance (70% - 90%), but the predicted muscle force only occurred during the latter portion of the stance phase (Figure 6(d)). Subject ACL4 showed a delay in both the predicted muscle activity and EMG of the gastrocnemius (Figure 4(c)). In the other three subjects, the timing of the predicted muscle force coincided with that of the control subjects but the gastrocnemius EMG was delayed slightly behind the predicted muscle force activity.
4. Discussion
We hypothesized that alterations in three-dimensional knee forces during the stance phase of gait could be detected as outliers following ACL reconstruction using PCA. The results show that abnormalities in gait waveforms can be detected in subjects who have undergone ACL reconstruction. Of the eight ACLR subjects analyzed, five were found to have an abnormality in either or both the CD and ML TFJ force. Four of these five subjects also had predicted muscle forces which deviated from those of the control population. The joint loading abnormalities appeared to be subject-specific, as there was variability across subject outcomes. This subjectspecific outcome is similar to what we found in our earlier work studying the knee joint moments during gait [13].
To our knowledge, this is the first report of knee joint forces and muscle forces during gait in subjects after ACL reconstruction. Tsai et al. [7] reported the peak TFJ compressive force in females who had undergone ACLR during a single-leg drop-land task. They found that the ACLR subjects had an increased peak compressive force when compared to a control group using an MRI-based EMG-driven model. We believe our results may differ due to the nature of the two different tasks, drop-land versus gait, and the modeling approach. Their model was EMG-driven, whereas ours was an inverse dynamics model. Also, their model only contained 10 muscles (hamstrings, quadriceps, and gastrocnemius). They did not incorporate the other muscles crossing the knee joint (popliteus, plantaris, gracilis, and sartorius) or other muscles which have been shown to affect tibiofemoral joint loading (i.e. soleus and gluteusmedius) [5]. Also, they do not show individual subject outcomes, so it is unclear how their subjects varied from normal on an individual basis.
It has been reported that acutely injured patients with ACL deficiency (ACLD) walk with decreased force on their injured knee compared to their uninjured knee [6]. This report was based on an EMG-driven model, and each subjects’ ACLD knee was compared to the contralateral knee. Our study differs in that we compared the ACLR knees to knees of a healthy control population. Also, the subjects in the prior study [6] had all been injured within the past seven months before testing, whereas our subjects were an average of 93 months post-surgery. Whereas ACLD subjects in the prior study [6] walked with decreased force on the ACLD knee, we found that outcomes for our ACLR subjects were subject-specific in whether their TFJ forces were increased, decreased, or in the range of healthy control subjects. However, both studies report similar TFJ contact force patterns and magnitudes. The timing of muscle activation in our ACLR and control subjects were also similar to the ACL deficient subjects, although magnitudes could not be compared since they did not report the body mass for the representative set of muscle forces shown.
There are several limitations that need to be acknowledged. Firstly, the knee was modeled as a hinge joint which did not require muscles to balance the adductionabduction and internal-external rotation torques in the model. Secondly, the influence of ligaments, cartilage, fluid, and other soft tissues was not included in this model. However, the passive force contributions by the ligaments, etc. have been shown to be small [19] and would have likely had little influence on our results. Also, the muscle physiological cross-sectional areas and moment arms were scaled from cadaveric data based on each subject’s mass and height since these values are difficult to measure in vivo. Predicted muscle forces have been shown to be sensitive to these parameters [20]. We also chose to distribute the muscle forces based on minimizing the sum of cubed muscle activations. Predictions of muscle force have also been shown to be sensitive to the optimality principle being applied [21]. The percentage of total muscle fibers activated (muscle activation) is related to muscle fatigue, and experimental evidence has shown this relationship to be approximately cubic [22]. Thus, while imperfect, we believe minimizing the sum of cubed muscle activation is a reasonable choice for the optimization criterion.
Finally, the sample sizes for both groups were relatively small. The results presented in this work were based on the assumption that the control dataset is truly representative of the normal population. Due to the small sample size we used to develop normal models and determine tolerance limits for comparison, appropriate caution must be exercised in interpreting the results.
5. Conclusion
Principal component analysis allowed us to identify TFJ loading waveforms and predicted muscle force waveforms which deviated from normal as determined using a control population. We found that of the eight ACLR knees, five knees displayed abnormal TFJ loading during the stance phase of gait. In four of these cases, the abnormal joint loading could be attributed to abnormal muscle forces of the major muscles crossing the knee joint. We believe that each subject has a unique response to their injury, reconstructive surgery, and rehabilitation. Further work is required to understand if this is due to surgical (e.g. graft orientation, placement, and tension) or patient variables (e.g. neuromuscular function, proprioception, and muscle strength) or both.
6. Acknowledgements
This work was funded in part by a grant from the Fed-Ex Institute of Technology to the University of Memphis.
NOTES