- Research article
- Open Access
Normal patellofemoral kinematic patterns during daily activities in dogs
BMC Veterinary Researchvolume 12, Article number: 262 (2016)
Patellar abnormalities are a common cause of pain and lameness in dogs; however, in vivo the relative motion between the femur and patella in dogs is not well described. The objective of this study was to define normal in vivo sagittal plane patellofemoral kinematics in three axes of motion using non-invasive methods. We hypothesized patellofemoral alignment in the sagittal plane would tightly correlate with the femorotibial flexion angle. Six healthy dogs without orthopedic disease underwent computed tomography (CT) of their hind limbs to create 3-D models of the patella and femur. Normal stifle joint motion was captured via flat-panel imaging while each dog performed a series of routine activities, including sitting, walking, and trotting. The 3-D models of the patella and femur were digitally superimposed over the radiographic images with shape-matching software and the precise movement of the patella relative to the femur was calculated.
As the femorotibial joint flexed, the patellofemoral joint also flexed and the patella moved caudally and distally within the femoral trochlea during each activity. Patellar flexion and distal translation during walk and sit were linearly coupled with the femorotibial flexion angle. Offset was evident while trotting, where patella poses were significantly different between early and late swing phase (p ≤ 0.003). Patellar flexion ranged from 51 to 6° while trotting. The largest flexion angle (92°) occurred during sit. The patella traversed the entire proximodistal length of the femoral trochlea during these daily activities.
Using single-plane flat-panel imaging, we demonstrated normal in vivo patellofemoral kinematics is tightly coupled with femorotibial kinematics; however, trot kinematic patterns did not follow the path defined by walking and stand-to-sit motions. Our normal data can be used in future studies to help define patellofemoral joint kinematics in dogs with stifle abnormalities.
The stifle is a complex synovial joint consisting of the femorotibial and patellofemoral joints. Stifle motion occurs in three planes resulting directly from the intricate anatomical relationship between the distal femur, proximal tibia and fibula, patella, pelvic limb musculature, and the stifle joint capsule and its associated ligaments . Patellofemoral abnormalities are common in dogs . Patellar luxation is a highly prevalent disorder affecting the stifles of dogs and alters stifle mechanics . Cranial cruciate ligament insufficiency is also common in dogs and has been shown to disrupt normal patellofemoral joint motion [3, 4].
Abnormal motion or kinematics of joints can lead to cartilage degradation, inflammation, pain, lameness, and progressive osteoarthritis . The traditional method utilized for analyzing canine stifle kinematics involves attaching reflective markers to the overlying skin of tissue landmarks with movement captured by digital cameras [5–8]. These studies, while non-invasive, do not provide precise data regarding movement of the underlying bones; instead these studies only yield general information on the angles and velocities of femorotibial flexion-extension movement. Imprecise placement of the reflective markers and skin motion can introduce variability in femorotibial kinematic data during gait analysis . Other methods for evaluating stifle kinematics include radiography and goniometry; however, these methods are not able to track the position of the patella during motion in dogs in vivo [9, 10].
Patellar kinematics is well described in humans and is typically referenced in relation to the femoral trochlea; of which tilt, rotation, flexion, and shift are of particular importance [11, 12]. Despite the high prevalence of patellofemoral problems in dogs, a similar description of normal patellar motion in this species has not been well described. Research evaluating the patellar kinematics following cranial cruciate ligament transection using a 2-D digital technique in a cadaveric model found patellar flexion angle was altered, suggesting abnormal patellofemoral biomechanics may play a role in the development of patellofemoral osteoarthritis in dogs with cranial cruciate insufficiency . Cadaveric studies, however, have several limitations and cannot fully simulate complex in vivo stifle biomechanics.
Joint kinematics can be quantified in a precise manner in-vivo by using fluoroscopic or flat-panel methods, where the region of interest is imaged while the subject performs tasks such as walking. In-vivo femorotibial kinematics in normal and cranial-cruciate ligament deficient stifles of dogs have been characterized using continuous horizontal-beam imaging [13–15]. To the investigators’ knowledge, in vivo dynamic patellofemoral kinematics has not been studied in dogs. The objective of this study was to define the normal in vivo sagittal plane patellofemoral kinematics in dogs during daily activities using horizontal-beam flat-panel imaging. We hypothesized patellofemoral alignment in the sagittal plane (patellar flexion-extension angle, craniocaudal translation, and proximodistal translation) would tightly correlate with the femorotibial flexion angle. We also hypothesized patellar movement would be restricted to the length of the femoral trochlea.
Six healthy Labrador Retrievers (four males, two females) were studied. All procedures were approved by the University of Florida Institution’s Animal Care and Use Committee. Dogs underwent complete physical and orthopedic examinations prior to data collection. The mean age was 4 years (range 1–7 years) and mean weight was 28 kg (range 26–32 kg). These dogs were confirmed via computed tomographic (CT) analysis to be free of pelvic limb orthopedic abnormalities.
Computed tomography analysesFootnote 1 were conducted to obtain stationary images of both hind limbs from the coxofemoral to the tarsocrural joint. Transverse images were obtained with a slice thickness of 0.5 mm.
Horizontal-beam lateral projection flat-panel imagesFootnote 2 of stifles were acquired while each dog walked at a velocity of 1.1 m/s (2.5 mph) and trotted at a velocity of 2 m/s (4.5 mph) on a treadmill for three separate trials with three to five strides per trial, and while the dog underwent a stand-to-sit motion for two separate trials. Images were acquired using a pulse width of 1 ms at 30 frames per second and an image area of 400 × 300 mm, producing a pixel size of 0.39 mm × 0.39 mm and image resolution of 1,024 × 1,024 pixels. Radiographic configurations supplied a 72 kV beam with a 50 mA beam current.
Three-dimensional model creation and coordinate assignation
Three-dimensional models of the patella and femur were constructed from the CT analyses of the subjects using the Digital Imaging and Communication in Medicine (DICOM) images and an open source 3-D-segmentation program.Footnote 3 The bone models were imported into reverse-engineering softwareFootnote 4 and coordinate systems were assigned based on local anatomical landmarks of the patella and femur [13, 16]. Femoral coordinates were applied such that the mediolateral axis (x-axis) passed through the center of the lateral and medial femoral condyles with the femoral origin located at the mid-point between the condyles (Fig. 1). The proximodistal axis (y-axis) passed proximally along the femoral shaft, perpendicular to the mediolateral axis in the plane common to the center of both femoral condyles and the femoral head. Patellar coordinates were applied such that the mediolateral axis (x-axis) passed through the most medial and lateral points on the circumference of the bone with the patellar origin defined as the mid-point of the axis. The proximodistal axis (y-axis) was defined as a line that passed through the most proximal and distal aspects of the patella, perpendicular to the mediolateral axis. The craniocaudal axes (z-axis) for the femur and patella were created from the cross product of the mediolateral and proximodistal axes, thus creating a Cartesian coordinate system.
Three-dimensional to Two-dimensional shape-matching
The 3-D bone models of the patella and femur and the single-plane flat-panel images were imported into open-source shape-matching softwareFootnote 5 (JointTrack). The 3-D models were digitally superimposed over the lateral flat-panel images with JointTrack and the model’s projected silhouettes were manipulated via translation and rotation such that the anatomic contours of the models precisely overlapped the corresponding contours of the 2-D flat-panel images (Fig. 2). The patella was positioned centrally within the trochlea groove, such that the center of the patella’s articulating surface remained as congruent as possible with the trochlear groove in the axial plane. This process was performed by one individual (EM) and was repeated for each frame of the various activities and cycles performed. The orientation of the shape-matched bone models was used to calculate the relative alignment between the patella and femur for patellar flexion-extension angle (degrees), proximodistal translation (mm), and mediolateral translation (mm) using the custom-written MatlabFootnote 6 program. The custom Matlab program transformed the data into clinically relevant patellofemoral poses in six degrees of freedom . The program was utilized such that the degrees of freedom were reported as patella relative to femur. Kinematic poses were normalized and interpolated to allow for comparative analysis of each cycle shape-matched. The kinematic properties were described in relation to femorotibial flexion angle, which have been reported for this exact data set in a separate study .
A statistical packageFootnote 7 was used for all analyses. Differences between single points during the gait cycle at equivalent flexion angles (offset) were determined by a paired t-test. A Pearson’s correlation coefficient was performed on the interpolated patellofemoral and femorotibial data to determine the extent of a linear relationship. The data was considered as strongly correlated if the r > 0.8 and weakly correlated in r < 0.3. For all statistical analyses performed, p < 0.05 was considered statistically significant.
As the femorotibial joint flexed, the patella also flexed (Fig. 3) and translated caudally (Fig. 4) and distally (Fig. 5) within the femoral trochlea for each activity. Patellar flexion, caudal translation, and distal translation during sit were linearly correlated with femorotibial flexion angle (Figs. 3, 4 and 5). All three measured degrees of freedom were correlated during walk and sit, however no correlation was found for these parameters during trot (Table 1). Offset, defined as significant differences in patellofemoral alignment at identical femorotibial flexion angles, was evident during trot, where the patella poses differed in early and late swing phase in each of the three measured degrees of freedom (Fig. 6; flexion-extension, p = 0.002; craniocaudal, p = 0.003; proximodistal, p = 0.002). Patellar flexion ranged from a mean of 51 to 6° during trot, and the largest patellar flexion angle of 92° was evident during sit (Table 1). Patellar flexion during the swing phase of walk ranged from 10° to a maximum of 33°, while a range of 7 to 13° was observed during stance (Table 2). During trot, patellar flexion angles ranged from 13 to 51° during swing and from 7 to 17° during stance (Table 2). The patella also traversed the length of the femoral trochlea during these daily activities, with nearly the entire patella positioned distal to the trochlear groove in deep flexion during sit (Fig. 7). The patella was confined within the proximal aspect of the trochlear groove during extension of all activities (Fig. 7).
Using single-plane flat-panel imaging, we demonstrated non-invasively that normal in vivo patellofemoral kinematics are tightly coupled with femorotibial kinematics. When the daily activities were assessed collectively, we observed patellar flexion angle increased and the patella translated distally as the femorotibial joint flexed. The highest demand activity, trotting, produced kinematic patterns that did not follow the same pathways seen during walking and stand-to-sit motions. Although a portion of the patella was positioned distal to the trochlear groove during deep flexion while sitting, the base of the patella always remained positioned within the femoral trochlea during the activities analyzed in these dogs.
Our results are similar to the findings from a canine cadaveric investigation, in which passive patellofemoral joint motion was described over a smaller range of femorotibial motion . Both studies demonstrated that all three kinematic parameters were linearly related to the femorotibial flexion angle. Our results, however, differ slightly with respect to the magnitude of observed change; for example, changes in femorotibial flexion angle from 90 to 150° induced proximodistal patella translation of approximately 13 mm in the cadaver study versus approximately 18 mm in our in vivo study . The change in patellar flexion of approximately 25° was reported in the cadaver study, whereas we found a change of 44° over the equivalent femorotibial range of motion . There are several explanations for the discrepancies between the studies. Most obviously, our investigation was an in vivo dynamic analysis, accounting for all the complex forces acting on the patella in vivo. Anatomic differences between breeds may have been a factor, as our study used Labrador Retrievers while the cadaver study used mix-breed dogs. Equally importantly, variations in coordinate assignation and reference points have been shown to dramatically affect patellar tracking patterns in humans, even within the same individual . To the author’s knowledge, this is only the second report characterizing patellar kinematics in dogs. Standardized coordinate systems and reference points for future studies would allow more meaningful comparison between investigations.
The most significant finding of the study was that the relationship between patellar poses and femorotibial flexion angle varied according to the phase of the gait cycle during trotting; the patella was positioned more proximal, more cranial, and more flexed in early swing phase when compared to late swing phase at the identical femorotibial flexion angle (Fig. 6). This offset likely contributed to the lack of a statistical linear correlation with the femorotibial flexion angle. The cause of this offset upon entering and exiting swing during trotting remains to be clarified. Potential causes include the varying magnitude at which the pelvic limb musculature, particularly the quadriceps muscle group, are contracting and acting on the patella, as well as the secondary motions of the tibia such as internal-external rotation and craniocaudal tibial translation. While an increased force of quadriceps contraction did not have significant effect on patellar movement in a human cadaveric study, an in vivo MRI study demonstrated the resting patellar position in the trochlea groove could be altered by muscular contraction during either open- or closed-chain exercises [18, 19]. Surface electromyographic studies in dogs have demonstrated widely varying patterns of vastus lateralis contraction according to differing activities, which could alter patellofemoral poses during different phases of the gait cycle as observed in our study . Electromyographic studies performed concurrently with radiographic kinematic analyses might improve our understanding of how muscular contraction alters patellar kinematics in dogs.
The single-plane flat-panel imaging and shape-matching methodology utilized in this study was previously validated for femorotibial kinematics in dogs [15, 21]. The precision of this methodology for this joint was determined to be within 1.28 mm and 1.58° with an intraobserver variability of less than 0.52 mm and 0.91° for translation and rotations, respectively . These values reported for the femorotibial joint are unlikely to be directly applicable for the patellofemoral joint due to differences in bone geometry. Studies performed in humans using similar methodology to evaluate the patellofemoral joint kinematics reported accuracy within 0.6° for in-plane rotations and 1.5 mm for in-plane translations . Analysis in this study was confined to sagittal plane translations and rotation due to the uniplanar nature of the technique and the patellofemoral joint anatomy. Although studies have not been performed to specifically validate the accuracy of our methodology for determining patellofemoral joint kinematics in dogs, a pilot series assessing repeatability of 3-D to 2-D image registration for this joint appeared to consistent with the results observed in other studies [23, 24].
The 3-D patella model was approximated to remain central within the trochlear groove such that the center of the patella’s articulating surface remained as congruent as possible with the trochlear groove in the axial plane. We utilized the 3-D geometry of the trochlear groove to define medial-lateral translation as well as patellar tilt. The kinematic parameters in the sagittal plane may have been different if the assumptions were not made. Studies in humans have shown the axial plane topography of the trochlear groove could be used to predict axial plane patellar kinematics, supporting our methodology . However, coronal plane rotation, conventionally known as patellar rotation, was not readily predictable by trochlear anatomy in humans . Bi-plane radiographic analysis would be required to gain a more thorough intricate understanding of canine patellofemoral joint motion in vivo.
Our data should be useful as a baseline of normal patellofemoral kinematics in dogs, against which comparisons can be made in future studies. The main goal of this study was to characterize normal patellofemoral motion in order to define the change in kinematics caused by cranial cruciate ligament rupture in a future study by our group. Because the most profound femorotibial kinematic abnormalities with cranial cruciate ligament rupture occur in the sagittal plane, it is logical to expect that the patellofemoral kinematics we reported, including patellar flexion, proximodistal translation, and craniocaudal translation, may be disrupted by the condition. Indeed, cadaveric studies have found the cranial cruciate ligament rupture can alter patella alignment and patellofemoral contact mechanics . Patellofemoral mechanics are also of particular interest with the surgical treatment of cranial cruciate ligament rupture as stifle extensor mechanism abnormalities frequently occur following various procedures used to address the disease . Our results are likely less applicable to patellar luxation, where the major abnormalities in motion are in the coronal plane.
The study used CT based models and single-plane flat-panel imaging to quantify patellofemoral motion during commonly performed daily activity in dogs. There was significant linear correlation between femorotibial flexion angle and each of the three degrees of freedom in the sagittal plane during sit and walk. The patella utilizes the full length of the trochlear groove during daily activity. Trotting induced offset, where patellar poses differed according to the phase of the gait cycle, despite identical femorotibial flexion angles. By defining normal in vivo kinematic patterns non-invasively, new knowledge of the normal canine stifle has been generated which may lead to the improvement and development of new surgical correction methods for stifle abnormalities.
Toshiba Aquilon 8, Toshiba American Medical Systems Inc, Tustin, CA.
Toshiba Infinix-I flat panel C-arm fluoroscope, Toshiba American Medical Systems Inc, Tustin, CA.
ITK-SNAP, Paul A. Yushkevich et al., http://www.itksnap.org.
Geomagic Studio, Geomagic Inc., Research Triangle Park, NC.
JointTrack, University of Florida, FL, http://sourceforge.net/projects/jointtrack/.
Matlab, MathWorks, Natick, MA.
g SigmaPlot 12, Systat Software, Inc. San Jose, CA.
Digital imaging and communication in medicine
Carpenter DH, Cooper RC. Mini review of canine stifle joint anatomy. Anat Histol Embryol. 2000;29:321–9.
Johnson JA, Austin C, Breur GJ. Incidence of canine appendicular musculoskeletal disorders in 16 veterinary teaching hospitals from 980 through 1989. Vet Comp Orthopaed. 1994;7:56–69.
Pozzi A, Dunbar NJ, Kim SE. Effect of tibial plateau leveling osteotomy on patellofemoral alignment: a study using canine cadavers. Vet J. 2013;198:98–102.
Guerrero TG, Pozzi A, Dunbar N, Kipfer N, Haessig M, Horodyski MB, et al. Effect of tibial tuberosity advancement on the contact mechanics and the alignment of the patellofemoral and femorotibial joints. Vet Surg. 2011;40:839–48.
Durant AM, Millis DL, Headrick JF. Kinematics of stair ascent in healthy dogs. Vet Comp Orthopaed. 2011;2:99–105.
Holler PJ, Brazda V, Dal-Bianco B, Lewy E, Mueller MC, Peham C, et al. Kinematic motion analysis of the joints of the forelimbs and hind limbs of dogs during walking exercise regimens. Am J Vet Res. 2010;71:734–40.
Agostinho FS, Rahal SC, Miqueleto NSML, Verdugo MR, Inamassu LR, El-Warrak AO. Kinematic analysis of Labrador retrievers and rottweilers trotting on a treadmill. Vet Comp Orthopaed. 2011;3:185–91.
Torres BT, Whitlock D, Reynolds LR, Fu YC, Navik JA, Speas AL, et al. The effect of marker location variability on noninvasive canine stifle kinematics. Vet Surg. 2011;40:715–9.
Thomas TM, Marcellin-Little DJ, Roe SC, Lascelles BDX, Brosey BP. Comparison of measurements obtained by use of an electrogoniometer and a universal plastic goniometer for the assessment of joint motion in dogs. Am J Vet Res. 2006;67:1974–9.
Torres BT, Punke JP, Fu YC, Navik JA, Speas AL, Sornborger A, et al. Comparison of canine stifle kinematic data collected with three different targeting models. Vet Surg. 2010;39:504–12.
Katchburian MV, Bull AMJ, Shih YF, Heatley FW, Amis AA. Measurement of patellar tracking: assessment and analysis of the literature. Clin Orthop Relat R. 2003;412:241–59.
Bull AMJ, Katchburian MV, Shih YF, Amis AA. Standardisation of the description of patellofemoral motion and comparison between different techniques. Knee Surg Sport Tr A. 2002;10:184–93.
Tashman S, Anderst W. In vivo measurement of dynamic joint motion using high speed biplane radiography and CT: application to canine ACL deficiency. J Biomech Eng. 2003;125:238–45.
Rey J, Fischer MS, Böttcher P. Sagittal joint instability in the cranial cruciate ligament insufficient canine stifle. Caudal slippage of the femur and not cranial tibial subluxation. Tierarztl Prax Ausg K Kleintiere Heimtiere. 2014;42:151–6.
Kim SE, Jones SC, Lewis DD, Banks SA, Conrad BP, Tremolada G, et al. In-vivo three-dimensional stifle kinematics during daily activities in dogs. J Orthop Res. 2015;33:1603–10.
Goh JCH, Lee PYC, Bose K. A cadaver study of the function of the oblique part of vastus medialis. J Bone Joint Surg. 1995;77:225–31.
Grood ES, Suntay WJ. A joint coordinate system for clinical description of three-dimenstional motions. Application to the knee. J Biomech Eng. 1983;105:136–44.
Nagamine R, Otani T, White SE, McCarthy DS, Whiteside LA. Patellar tracking measurement in the normal knee. J Orthop Res. 1995;13:115–22.
Felicio LR, Saad MC, Liporaci RF, Baffa Ado P, dos Santos AC, Bevilaqua-Grossi D. Evaluating patellar kinematics through magnetic resonance imaging during open- and closed-kinetic-chain exercises. J Sport Rehabil. 2010;19:1–11.
Breitfuss K, Franz M, Peham C, Bockstahler B. Surface electromyography of the vastus lateralis, biceps femoris, and gluteus medius muscle in sound dogs during walking and specific physiotherapeutic exercises. Vet Surg. 2014. doi:10.1111/j.1532-950X.2014.12302.x.
Jones SC, Kim SE, Banks SA, Conrad BP, Abbasi AZ, Tremolada G, et al. Accuracy of noninvasive, single-plane fluoroscopic analysis for measurement of three-dimensional femorotibial joint poses in dogs. Am J Vet Res. 2014;75:477–85.
Akbarshahi M, Schache AG, Fernandez JW, Baker R, Banks S, Pandy MG. Non-invasive assessment of soft-tissue artifact and its effect on knee joint kinematics during functional activity. J Biomech. 2010;43:1292–301.
Fregly BJ, Rahman HA, Banks SA. Theoretical accuracy of model-based shape matching for measuring natural knee kinematics with single-plane fluoroscopy. J Biomech Eng. 2005;127:692–9.
Mahfouz MR, Hoff WA, Komistek RD, Dennis DA. A robust method for registration of three-dimensional knee implant models to two-dimensional fluoroscopy images. IEEE Trans Med Imaging. 2003;22:1561–74.
Varadarajan KM, Freiberg AA, Gill TJ, Rubash HE, Li G. Relationship between three-dimensional geometry of the trochlear groove and in vivo patellar tracking during weight-bearing knee flexion. J Biomech Eng. 2010. doi:10.1115/1.4001360.
Mattern KL, Berry CR, Peck JN, De Han JJ. Radiographic and ultrasonographic evaluation of the patellar ligament following tibial plateau leveling osteotomy. Vet Radiol Ultrasoun. 2006;2:185–91.
The authors thank Mr Tim Vojt for assisting with producing the graph figures.
This study was funded by grants for the Merial Veterinary Student Scholars Program and Departmental Funds at the University of Florida College of Veterinary Medicine.
Availability of data and materials
The dataset supporting the conclusions of this article is available in Open Science Framework repository, https://osf.io/2j4bz/
EJM shape-matched the 3-D bone models of the femur and patella onto the flat-panel images, analyzed the data, and drafted the manuscript. SEK participated in the design of the study, assisted in data analysis, and helped draft the manuscript. JC created the 3-D bone models and trained EJM in shape-matching. SCJ provided the shape-matched data for femorotibial kinematics. SAB, DDL and AP participated in the design of the study, as well as drafting of the manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
All procedures were approved by the University of Florida Institution’s Animal Care and Use Committee. Privately owned dogs were recruited to this study. Informed consent was obtained from owners of the animals in this study.