- 1Clinic for Small Animal Surgery, Department for Small Animals, Vetsuisse Faculty University of Zurich, Zürich, Switzerland
- 2Department of Veterinary Clinical Sciences, Colorado State University, Fort Collins, CO, United States
- 3Department of Biomedical Engineering, National Taiwan University, Taipei, Taiwan
Bone morphometry varies among dogs of different sizes and breeds. Studying these differences may help understand the predisposition of certain breeds for specific orthopedic pathologies. This study aimed to develop a statistical shape model (SSM) of the femur, patella, and tibia of dogs without any clinical orthopeadic abnormalities to analyze and compare morphological variations based on body weight and breed. A total of 97 CT scans were collected from different facilities and divided based on breed and body weight. The 3D models of the bones were obtained and aligned to a coordinate system. The SSM was created using principal component analysis (PCA) to analyze shape variations. The study found that the first few modes of variation accounted for a significant percentage of the total variation, with size/scale being the most prominent factor. The results provide valuable insights into normal anatomical variations and can be used for future research in understanding pathological bone morphologies and developing 3D imaging algorithms in veterinary medicine.
1 Introduction
Dogs have specific anatomical and physiological variations among breeds, body weight, and size (1, 2), that may predispose to pathologies. For example, brachycephalic dogs are subject to respiratory difficulties partially due to the shape of their skull (3). While some research has been done in small groups of dogs (4–7), no studies have provided a large-scale dataset to investigate bone anatomical variations among breeds or body weight. Since bone morphology has been shown to influence the development of certain diseases such as medial patellar luxation (MPL), large-scale methodologies to investigate bone geometry have a high potential for the prevention and treatment of these diseases (8).
Geometric morphology is the study of shape assessment using landmarks. It looks at the shape as a whole, rather than using a restricted set of anatomical parameters per se. There is limited information about the geometric morphology of the femur, tibia, and patella in dogs. Palierne et al. (9) measured different anatomical parameters for the normal femur. Savio et al. (10) and Soparat et al. (11) reported similar femoral varus angle (FVA) for large breed dogs (5.5°) and Pomeranian (5.85°). Similar studies have reported data for proximal tibia, but only for a single or two specific breeds (6, 12, 13). The relationship between geometric bone morphology and different orthopedic diseases is yet to be extensively investigated.
The canine hind limbs are subject to very common diseases such as hip dysplasia (HD), cranial cruciate ligament disease (CCLD), medial patella luxation (MPL), and osteochondrosis (OC) (14–16). Some studies have shown the correlation between the shape of the bones and the prevalence of these diseases by analyzing various anatomical bone measurements (17–20). The main regions susceptible to pathological variation are the extremities of the bones such as the femoral head and condyles, tibia plateau and the small bones such as the patella due to the previously mentioned disease (14–16). However, these anatomical variations can be very specific for certain breeds. The previously mentioned studies would suggest that most of the variations are in size, thickness and curvature of the bones but the specific localization and type of geometrical variation has not been investigated.
Statistical shape model (SSM) is a computational technique used in human medicine to study geometric bone morphology, and disease pathogenesis and perform medical imaging analysis, automatic and semi-automatic segmentation, and electrocardiogram simulation among other applications (21–31). Using principal component analysis (PCA), it can provide a list of the main geometrical properties (or variation parameters) of a set of 3D bone models, also known as modes, the main ones often being the size and thickness (22). PCA allows to analyze and visualize the complete 3D models by describing their shape with a mathematical algorithm, instead of quantifying each length, angle and size manually (10–13). This method provides a more global understanding of bone geometry and does not limit our understanding to local changes (21). SSM does not provide any information on the internal morphology (i.e., cortical thickness), the bone axis or joint angles which are parameters often used by clinicians. It only describes the shape of the bone with a mathematical equation. Based on our knowledge, SSM has never been used to research bone morphology in dogs.
This study aims to develop an SSM of the femur, patella, and tibia of dogs without any clinical orthopeadic abnormalities and to analyze and compare the morphological variations based on body weight and breeds. These data will provide a better understanding of the normal anatomical variations and can be used as a baseline for future studies, investigating morphologies of pathological bones, and developing and validating SSM-based 3D image reconstruction algorithms (32, 33).
2 Materials and methods
2.1 Data sample
A total of 97 computed tomography scans (CT) were collected retrospectively from three different facilities (University of Zürich, Switzerland, and University of Wisconsin, United States, University of Florida, United States). Inclusion criteria was the presence of at least the femur, patella, and proximal part of the tibia. Exclusion criteria was the presence of any pathologic abnormalities seen in the clinic or on the CT by experts, such as osteoarthritis, HD, MPL, osteochondrosis or previous surgery. This study included CT scans of both cadaveric limbs and live dogs.
2.2 3D model reconstruction
All hindlimbs were scanned with 0.8-mm slice thickness with 512×512 matrix, pixel size ranging between 0.1927 mm and 0.473 mm (University of Florida—Model Toshiba, University Zürich—Brilliance CT, Philips AG, Zurich, Switzerland, and University of Wisconin-Model GE Medical Systems) The images were manually segmented using an open-source software (3D Slicer). Initial segmentation was performed using a threshold value of 350, followed by performing a smoothing using a built-in function for closing (fill holes) with kernel size of 3.0 mm. Each bone (femur, patella, tibia, and fabella) was segmented separately to create an individual model.
2.3 Coordinate system application and meshing
The models were meshed using the engineering software Geomagic WRAP (Geomagic, Inc., Research Triangle Park, NC) to obtain 3D models of reduced size and uniformly sampled vertices (target edge length 2 mm). Complete bones and partial models (proximal distal part of the femur and proximal part of the tibia) are created for evaluations. All femur models were cut to isolate the region of interest with uniform length. Due to large variation of length and size, the length of the partial femur model was normalized to be two times the length of the condyles (partial distal model) and femoral head (partial proximal model). The length was measured to be the diameter of a sphere fitted in the condyles and the diameter of a sphere fitted in the femoral head. The partial proximal tibia model was cut at the most distal part of the tibial tuberosity (34–37).
The models were aligned to an anatomical coordinate system based on anatomic landmarks, similar to previous joint kinematics studies (38). The origin was centered on the point equally distant to the center of the two condyles. Where the center was determined by registering a sphere to the geometry of the femoral condyle. The Y-axis was set to be parallel to the anatomical axis of the femur in a proximal-distal direction, with a positive direction going proximal. The Z-axis was set to be in a latero-medial direction, with the positive direction going medial. The X-axis was then calculated to be perpendicular to the other two axes. For each group and set of bones (left and right femur, tibia and patella, and the different groups), reference/template models were created. The template models were created from the median bone model, remeshed to have uniform vertices of 2 mm (39). All the training models were registered to the template model using the nearest neighbor algorithm (MATLAB and Statistics Toolbox Release R2020a, The MathWorks, Inc., Natick, MA.) to unify the number of faces and vertices.
2.4 Statistical shape model
Once all the training models were registered and aligned to the template model, PCA was performed. This operation generated the results as the shape variation and Modes for individual bones. A Mode was defined as the modification of the mean model due to the variation of only one component of the PCA. This allowed individual analysis of each component. The Modes were numbered 1 to 5 based on the magnitude of variation that was detected in the bones. For analysis, the data for the femur were divided between body weight groups, breed, and bones (Table 1). The specific bone regions analyzed were the femoral condyle, femoral head, and the tibial plateau. A realignment of the femur model was performed with Geomagic WRAP to obtain more detailed results in some cases, such as for interpretation of the Mode 2. Observation of the models at different alignment was used to obtain more insight in localization of the deformity in the different Modes. They were aligned based on the shape of their extremities (for example: the second Mode of the complete femur). A total of 33 SSMs (mean, ±3 standard deviations) representing the first 5 Modes from the PCA were generated for visual evaluations (Figures 1–6) (25, 29).
Figure 1. First two modes for the PCA of the left femur from Group 1 [(A) Mode 1 viewed from cranial, (B) Mode 2 viewed from cranial]. They are for each mode 3 models to be seen: Mean model and ±3 SD.
Figure 2. Modes 2, 3, 4, and 5 for the PCA of the left femur [(A) mode 2 viewed from cranial after realignment, (B) mode 3 viewed from lateral, (C) mode 3 viewed from crainal, (D) mode 4 viewed form lateral, (E) mode 5 viewed from cranioproximal, (F) mode 5 viewed from lateral]. They are for each modes 3 Models to be seen: Mean model and ±3 SD (only the most external model is described, then comes the mean model and the opposite model).
Figure 3. First SSM mode for the left femur from Group 2 (viewed from cranial). They are 3 Models to be seen: Mean model and ±3 SD.
Figure 4. Modes 1 and 3 from the PCA of the left femur condyles [(A) mode 1 viewed from cranial, (B) mode 3 viewed from lateral] and modes 1 and 2 from the PCA of the left femoral head [(C) mode 1 viewed from cranial, (D) mode 2 viewed from cranial]. They are for each mode 3 Models to be seen: Mean model and ±3 SD.
Figure 5. Modes 1, 4, and 5 from the PCA of the left tibia [(A) mode 1 viewed from cranial, (B) mode 4 viewed from lateral, (C) mode 5 viewed from lateral]. Modes 1 and 4 from the PCA of the left tibia plateau [(D) mode 1 viewed from cranial, (E) mode 4 viewed from lateral]. They are for each mode 3 Models to be seen: Mean model and ±3 SD.
Figure 6. Modes 1, 2, and 3 from the PCA of the left patella [(A) mode 1 viewed from lateral, (B) mode 2 viewed from lateral, (C) mode 3 viewed from cranial]. They are for each mode 3 Models to be seen: Mean model and +/− 3 SD.
3 Results
From 97 samples, 61 models without any presence of clinical orthopeadic abnormalities of complete femur with head and condyles, 30 complete tibia, 72 tibia plateaus and 85 patella were collected after exclusion of pathological conditions (n = 11), from 41 different type of breeds (Tables 1, 2). Groups 1 through 6 are separated by weight class and breed for complete femur. The different weight classes were set as follows: dogs from 0 to 10 kg, dogs from 10 to 25 kg and dogs greater than 25 kg body weight. The breed groups consisted of 2 groups: one for the retrievers (Golden and Labrador) and one for the chondrodystrophic breeds (French bulldog, Maltese, Dachshund, Bichon Frise, Beagles). Groups 7 and 8 were separated to evaluate partial femurs (femoral condyles and femoral head). Due to limited sample collection for complete tibia, separate evaluation for breed and weight class were not performed. Group 9 and 10 were created for complete tibia and partial tibia (proximal end of the tibia) for analysis. The last SSM group is the patella. The dog population incorporated in this study has the media age of 10 years (range 0.5–17), the median weight of 27.4 Kg (range 3–65) with 30 different with type of breeds (Table 2).
Table 2. Every breed and number of dogs included in the study (some dogs only had one type of bones).
For the analysis of the complete bones, the first 5 Modes of variations accounted for around 99% of the total variation (Table 3). For the first group of the femur, the first Mode accounted for 93.5% of the total variation with only needing the first 2 Modes to add up to 99.1% of total variation (Table 3). The mean and the standard deviation models for Group 1 and Group 2 are shown in Figures 1–3. For all complete femur groups, size/scale seems to best represent the variation for Mode 1.
For the analysis of the head and condyles of the femur, the first 5 Modes of variations accounts for around 93.9%–98.0% of the total variation (Table 4). The tables also show the results for the Tibia, Tibia plateau and Patella. For all Groups the first 5 Modes describe more than 93% of the total variance and each impact of following Modes reduces significantly. The results from left and right side are similar. The most visually relevant Modes of the regions of femur, tibia, complete tibia, and patella are shown in Figures 4–6.
4 Discussion
This is the first study to utilize SSM to evaluate and compare the geometric morphology of the femur, tibia, and patella, from dogs without any clinical orthopeadic abnormalities of different sizes and breeds. Eleven groups of analysis were performed using SSM to explore the factors responsible for morphological variability. We found that regardless of size, breed, and bone type, the first Mode showed that the scale factor accounted for most variation (>90% of total variation). Also, it only takes 2–3 Modes to describe more than 99% of total shape variation. In clinical terms, this result suggests that the shape of the bones from a small dog are similar to bones from a larger dog. Another important finding for the femur and tibia was that the second most important source of variation was varus-valgus. Mode 3, 4 and 5 were more difficult to interpret, because the deformation of the shape may depend on multiple factors including shaft diameter and shape, torsion, and epiphysis shape. These results provide a new perspective on how to approach bone morphologic studies on a large scale.
4.1 Femur
The second Mode showed a strong varus/valgus deformation of the femur, confirming previous studies that showed strong (up to 4°) variation in normal dogs (9, 40, 41). This result was found in both chondrodystrophic and retriever groups suggesting that varus-valgus is an important anatomical parameter, irrespective of breed. To get further insight on the location of the deformation we realigned the models based on the two extremities and at the mid-diaphysis to observe the curvature of the bone. This study discovered that there is effectively some curvature of the femoral shaft, although it is less than expected. The condyles are also not superimposing, which could also be interpreted as a difference in femoral condyle size.
We found that Modes 3, 4 and 5 were influenced by several factors, which made them more difficult to interpret. It is also important to realize that collectively these 3 Modes accounted only for 0.62% of variation. The best interpretation of the third Mode is a combination of three different components: width of the condyles, femoral head size and craniocaudal thickness of the whole femur. Mode 4 also shows 2 components: procurvation and torsion of the femoral head. Mode 5 could also have thickness and malignment components.
The magnitude of the deformation dependent on Modes 3, 4 and 5 is small compared to the Mode 1 and 2. This is expected because the majority of the variations are from Mode 1. This suggests that when the size factor is eliminated, the normal femur is very similar across different breads. Even though for Groups 2 through Group 6, categorized by weight class and specific breeds, the variability between individual dogs is still dominated by the size/scale factor. This can be identified by the changes in the length of the femur. Group 4 shows the most similar bone length (see Table 1) resulting in less variability in size between Models from large dogs. This reduction of size variability is reflected with the lower percentage represented by Mode 1 compared to other groups. For Group 5, typical specificities of the chondrodystrophic breeds can be seen on the mean model, such as a smaller major trochanter, thicker condyles, or femoral head, and stronger procurvation, which can also be seen on the variant Modes.
To reduce the effect of the femoral length and to isolate the proximal and distal parts of the femur, all the femur models were cut. The only clear interpretation could be made from the first Mode of the femoral condyles SSM, which is the size factor. The rest of the Modes showed a mixture of different combined factors of condyle length, thickness, length of the femoral shaft above the patella grove, and malalignment. Previous studies investigated the femoral condyle size or other condyle variation as a morphological factor that may predispose to CCLD in dogs (42, 43) as well as humans (44, 45). Future studies could use SSM to compare these parameters in dogs with CCLD and normal dogs. For the femoral head SSM, size was the most dominant factor. Interestingly the second Mode represented the size of the femoral head independently of the major trochanter. Modes 3 to 5 represented some differences in alignment as well as a small difference in the anteversion angle of the femoral head. However, the relevance is to be questioned, as the differences are minor. It should be noted that no pathological deformation has been observed such as patellar groove depth differences typically seen in MPL disease, strong varus/valgus of the femoral shaft or deformation of the femoral head typically seen in hip dysplastic cases (18, 20). This indicates a good set of inclusion and exclusion criteria for our data sets.
4.2 Tibia
The first 2 Modes from the complete tibia SSM are similar to the femur SSM, representing the size and a combination of varus/valgus variability. While the third Mode showed a combination of different factors, the fourth Mode showed a difference in the craniocaudal thickness of the tibial shaft. The fifth Mode represented a difference in the thickness of the tibia plateau and distal extremity. Interestingly, no Modes show a difference in the position of the tibial tuberosity in relation to the tibial shaft or a torsion like the femur. The lower variability in morphology of the tibia relative to the femur is consistent with previous studies that show that most variation is situated in the femur (46). The Modes from the PCA focused on the tibial plateau follow a similar pattern, with Mode 1 being the size factor. The only other clear factor was observed in Mode 4, procurvation. Other Modes did not show any clear identification of variability.
4.3 Patella
The first Mode of the PCA of the patella shows size as the most important factor. Interestingly the second Mode shows an elongation or a more sharply cornered proximal extremity of the patella. Modes three to five most likely represent malalignment.
4.4 Interpretation
It should be noted that while some Modes show certain deformations, they can also be interpreted in different ways depending on which Modes are considered. For example, an elongation-type deformation (with a variation on the shaft of the bone only and none on the thickness) can be interpreted as either elongation or as the thickness of the bone, as an elongated bone would be a narrower bone after size adjustment (as the main focus is the proportional variance). Another type of biased interpretation is best shown with the varus/valgus example. If a shaft is seen to be curved while aligned on the condyle, the curvature is not necessarily based on the condyle. It could be situated anywhere on the shaft. Only further analysis with realignment could address this problem.
While the current study does not report measurement of any length or angles, some comparisons can be made with the literature. When reviewing the measurements from different studies on geometrical bone morphology, they show high variability. Additionally, when comparing breeds, the different values overlap between small and large breeds and do not show a clear difference. For example, when looking at the anteversion angle of the femur, some studies found for small breeds angles between 20° and 27°, while on large dogs the values were between 16° and 25° (7, 10, 47–49). While precise angles were not measured in this study, variations on procurvation and torsion which both influence those angles were observed. Stronger procurvation on smaller breeds was observed when comparing different breeds.
Since the number of studies on the tibial angles is limited, the results show more correlation. Smaller breeds have a greater tibia plateau angle (30° vs. 25° in large breeds), Z angle (69.2° to 58.8°), distal tibial axis/proximal tibial axis angle (10.8° to 4.5°), greater relative tibial tuberosity width (0.86° to 0.74°).
We could not see that in our study since we found an overall similarity between breeds. However, the overlapping from those studies is great when looking at the range of the values. Nevertheless, we could show a great variability which could explain the differences seen in those studies. Differences in varus/valgus, procurvation and thickness were observed from the tibia SSM, which all influence those angles (50, 51).
To the best of our knowledge, this study is the first to present SSM for dogs’ hind limbs. This method of measurement considers the whole bone shape rather than single measurements, which is a crucial step toward a better understanding of the anatomy and further research on diseased bones. The different results in measurements in the literature reinforce the need for SSM. One of the tools of SSM is classification, the possibility of classifying if a stifle from a specific case is normal or pathologic bone, without having to rely on individual clinicians’ measurements. With the development of SSMs from bones without any clinical orthopeadic abnormalities, the process can be extended to diseased dogs to observe qualitative and quantitative changes in the Modes for certain diseases such as MPL or HD. This could allow the description of the pathologic changes and the correction needed without the need for manual analysis. Another useful application is using this SSM model to help with the surgery on certain bones by facilitating the finding of landmarks which are needed to perform the operation.
4.5 Limitations
This study has several limitations. The most significant factor is the limited sample size. The complete set of femurs (n = 61) could be considered a large enough data set, but when divided into groups, the sample sizes for the individual groups are greatly reduced. PCA analysis is a powerful tool that can help represent variations of geometric morphology, but it does not identify a single factor for the significance of each Mode it generates. Therefore, each Mode must be visually analyzed and interpreted. While some of the Modes can be easily identified (Mode 1 for scale), other Modes can be more difficult to interpret due to the significantly smaller magnitude of shape variation presented. Other methods like independent component analysis can be used as an alternative, but they all have their strengths and weaknesses. Also, despite consistent methods of model preparation, errors from segmentation, pre-registration, and alignment cannot be eliminated. Lastly, due to the significant effect of the size variation, other clinically interesting variations of deformation could not be extensively identified. However, this could also be because the current training data sets were from dogs without any clinical orthopeadic abnormalities. And from the chondrodystrophic dog’s half of them were beagle which shows only a small amount of the typical differences seen in those breeds. Perhaps, when data from pathological cases are evaluated, these clinically interesting variations will be observed.
5 Conclusion
While SSM is well-known in human medicine, it has only been applied a few times in veterinary orthopedics (52, 53). The goal of this study was to take a further step in this research domain and to establish baseline statistical shape models and values for normal bones. Starting with the stifle (femur, tibia, and patella), SSM was used to evaluate and compare the geometric variation between breed and weight classes for dogs. The results from this study will be used as a reference for further research on investigating geometric morphology for clinically pathological bones and support the advancement of the usage of digital technologies in the veterinary field.
Data availability statement
The data analyzed in this study is subject to the following licenses/restrictions: all our CTs data are controlled by the University of Zürich. Requests to access these datasets should be directed to bpark@vetclinics.uzh.ch.
Ethics statement
Ethical approval was not required for the studies involving animals in accordance with the local legislation and institutional requirements because this was a retrospective study. Written informed consent was not obtained from the owners for the participation of their animals in this study because this was a retrospective study.
Author contributions
JH: Writing – original draft, Writing – review & editing, Conceptualization, Investigation, Methodology. AP: Conceptualization, Project administration, Supervision, Validation, Visualization, Writing – review & editing. JB: Writing – review & editing. T-WL: Writing – review & editing. SK: Conceptualization, Methodology, Supervision, Writing – review & editing. BP: Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Project administration, Software, Supervision, Validation, Visualization, Writing – review & editing.
Funding
The author(s) declare that no financial support was received for the research, authorship, and/or publication of this article.
Acknowledgments
I would like thank Federico Longo and Pavlos Natsios for their clinical insight and advice during the project.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1. Boyko, AR, Quignon, P, Li, L, Schoenebeck, JJ, Degenhardt, JD, Lohmueller, KE, et al. A simple genetic architecture underlies morphological variation in dogs. PLoS Biol. (2010) 8:e1000451. doi: 10.1371/journal.pbio.1000451
2. Oswald, H, Sharkey, M, Pade, D, and Martinez, MN. Canine gastrointestinal physiology: breeds variations that can influence drug absorption. Eur J Pharm Biopharm. (2015) 97:192–203. doi: 10.1016/j.ejpb.2015.09.009
3. Ekenstedt, KJ, Crosse, KR, and Risselada, M. Canine Brachycephaly: anatomy, pathology, genetics and welfare. J Comp Pathol. (2020) 176:109–15. doi: 10.1016/j.jcpa.2020.02.008
4. Matchwick, A, Bridges, JP, Mielke, B, Pead, MJ, Phillips, A, and Meeson, RL. Computed tomographic Mesurement of trochlear depth in three breeds of brachycephalic dog. Vet Comp Orthop Traumatol. (2021) 34:124–9. doi: 10.1055/s-0040-1719165
5. Al Aiyan, A, Richardson, K, Manchi, G, Plendl, J, and Brunnberg, L. Measurement of the femoral neck angle in medium and large breeds using computed tomography. Acta Vet Hung. (2019) 67:22–33. doi: 10.1556/004.2019.003
6. Dismukes, DI, Tomlinson, JL, Fox, DB, Cook, JL, and Song, KJ. Radiographic measurement of the proximal and distal mechanical joint angles in the canine tibia. Vet Surg. (2007) 36:699–704. doi: 10.1111/j.1532-950X.2007.00323.x
7. Tomlinson, J, Fox, D, Cook, JL, and Keller, GG. Measurement of femoral angles in four dog breed. Vet Surg. (2007) 36:593–8. doi: 10.1111/j.1532-950X.2007.00309.x
8. Vanden Berg-Foels, WS, Schwager, SJ, Todhunter, RJ, and Reeves, AP. Femoral head shape differences during development may identify hips at risk of degeneration. Ann Biomed Eng. (2011) 39:2955–63. doi: 10.1007/s10439-011-0393-3
9. Palierne, S, Asimus, E, Mathon, D, Meynaud-Collard, P, and Autefage, A. Geometric analysis of the proximal femur in diverse sample of dogs. Res Vet Sci. (2006) 80:243–52. doi: 10.1016/j.rvsc.2005.07.010
10. Savio, G, Baroni, T, Concheri, G, Baroni, E, Meneghello, R, and Longo, F. Computation of femoral canine morphometric parameters in three-dimensional geometrical models. Vet Surg. (2016) 45:987–95. doi: 10.1111/vsu.12550
11. Soparat, C, Wangdee, C, Chuthatep, S, and Kalpravidh, M. Radiographic measurments for femoral varus in Pomeranian dogs with and without medial patellar luxation. Vet Comp Orthop Traumatol. (2012) 25:197–201. doi: 10.3415/VCOT-11-04-0057
12. Comerford, EJ, Tarlton, JF, Avery, NC, Bailey, AJ, and Innes, JF. Distal femoral intercondylar notch dimensions and their relationship to composition and metabolism of the canine anterior cruciate ligament. Osteoarthr Cartil. (2006) 14:273–8. doi: 10.1016/j.joca.2005.09.001
13. Baroni, E, Matthias, RR, Marcellin-Little, DJ, Vezzoni, A, and Stebbins, ME. Comparison of radiographic assessments of the tibial plateau slope in dogs. Am J Vet Res. (2003) 64:586–9. doi: 10.2460/ajvr.2003.64.586
14. King, MD. Etipathogenesis of canine hip dysplasia, prevalence, and genetics. Vet Clin Small Anim. (2017) 47:753–67. doi: 10.1364/OE.520251
15. Harari, J. Osteochondrosis of the femur. Vet Clin North Am Small Anim Pract. (1998) 28:87–94. doi: 10.1016/S0195-5616(98)50006-X
16. Roush, JK. Canine Patellar Luxation. Vet Clin North Am Small Anim Pract. (1993) 23:855–68. doi: 10.1016/S0195-5616(93)50087-6
17. Butler, JR, and Gambino, J. Canine hip dysplasia diagnostic imaging. Vet Clin Small Anim. (2017) 47:777–93. doi: 10.1016/j.cvsm.2017.02.002
18. Lusetti, F, Bonardi, A, Eid, C, Brandstetter de Belesini, A, and Martini, FM. Pelvic limb alignment measured by computed tomography in purebred English bulldogs with medial patellar luxation. Vet Comp Orthop Traumatol. (2017) 30:200–8. doi: 10.3415/VCOT-16-07-0116
19. Reagan, JK. Canine hip dysplasia Creening within the United States. Vet Clin Small Anim. (2017) 47:795–805. doi: 10.1016/j.cvsm.2017.02.003
20. Newman, M, and Voss, K. Computeed tomographic evalutation of femoral and tibial conformation in English Staffordshire bull terriers with and without congenital medial patellar luxation. Vet Comp Orthop Traumatol. (2017) 30:191–9. doi: 10.3415/VCOT-16-12-0162
21. Ding, AS, Lu, A, Li, Z, Galaiya, D, Ishii, M, Siewerdsen, JH, et al. Statistical shape model of the temporal bone using segmentation propagation. Otol Neurotol. (2022) 43:679–87. doi: 10.1097/MAO.0000000000003554
22. Audenaert, EA, Pattyn, C, Steenackers, G, De Roeck, J, Vandermeulen, D, and Claes, P. Statistical shape modeling of skeletal anatomy for sex discrimination: their training size, sexual dimorphism, and asymmetry. Front Bioeng Biotechnol. (2019) 7:302. doi: 10.3389/fbioe.2019.00302
23. Heimann, T, and Meinzer, HP. Statistical shape models for 3D medical image segmentation: A review. Med Image Anal. (2009) 13:543–63. doi: 10.1016/j.media.2009.05.004
24. Bruse, JL, Zuluaga, MA, Khushnood, A, Mcleod, K, Ntsinjana, HN, Hsia, TY, et al. Detecting clinically meaningful shape clusters in medical image data: metrics analysis for hierarchical clustering applied to healthy and pathological aortic arches. IEEE Trans Biomed Eng. (2017) 64:2373–83. doi: 10.1109/TBME.2017.2655364
25. Cootes, TF, Hill, A, Taylor, CJ, and Haslam, J. The use of active shape models for locating structures in medical images. Image Vis Comput. (1994) 12:355–66. doi: 10.1016/0262-8856(94)90060-4
26. Lu, HY, Shih, KS, Lin, CC, Lu, TW, Li, SY, and Kuo, HW. Three-dimensional subject-specific knee shape reconstruction with asynchronous fluoroscopy images using statistical shape modeling. Front Bioeng Biotechnol. (2021) 9:736420. doi: 10.3389/fbioe.2021.736420
27. Lu, HY, Lin, CC, Shih, KS, Lu, TW, Kuo, MY, and Li, SY. Integration of statistical shape modeling and alternating interpolation-based model tracking technique for measuring knee kinematics in vivo using clinical interleaved bi-plane fluoroscopy. PeerJ. (2023) 11:e15371. doi: 10.7717/peerj.15371
28. Thodberg, HH, Kreiborg, S, Juul, A, and Pedersen, KD. The BoneXpert method for automated determination of skeletal maturity. IEEE Trans Med Imaging. (2009) 28:52–66. doi: 10.1109/TMI.2008.926067
29. Neumann, A, and Lorenz, C. Statistical shape model based segmentation of medical image. Comput Med Imaging Graph. (1998) 22:133–43. doi: 10.1016/S0895-6111(98)00015-9
30. Ambellan, F, Lamecker, H, von Tycowicz, C, and Zachow, S. Statistical shape models: understanding and mastering variation in anatomy. Adv Exp Med Biol. (2019) 1156:67–84. doi: 10.1007/978-3-030-19385-0_5
31. Nagel, C, Schuler, S, Dössel, O, and Loewe, A. A bi-atrial statistical shape model for large-scale in silico studies of human atria: Model development and application to ECG simulations. Med Image Anal. (2021) 74:102210. doi: 10.1016/j.media.2021.102210
32. Khanduja, V, Baelde, N, Dobbelaere, A, Van Houcke, J, Li, H, Pattyn, C, et al. Patient-specific assessment of dysmorphism of the femoral head–neck junction: a statistical shape model approach. Int J Med Robot. (2016) 12:765–72. doi: 10.1002/rcs.1726
33. Gerber, N, Reyes, M, Barazzetti, L, Kjer, HM, Vera, S, Stauber, M, et al. A multiscale imaging and modelling dataset of the human inner ear. Scientificdata. (2017) 4:132. doi: 10.1038/sdata.2017.132
34. Grood, ES, and Suntay, WJ. A joint coordinate system for the clinical description of three-dimensional motions: application to the knee. J Biomech Eng. (1983) 105:136–44.
35. Park, BH, Leffler, J, Franz, A, Dunbar, NJ, and Banks, SA. Kinematics of monoblock bicompartmental knee arthroplasty during weight-bearing activities. Knee Surg Sports Traumatol Arthrosc. (2015) 23:1756–62. doi: 10.1007/s00167-014-3427-1
36. 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 treated by tibial plateau leveling osteotomy. Am J Vet Res. (2014) 75:486–93. doi: 10.2460/ajvr.75.5.477
37. Kim, SE, Jones, SC, Lewis, DD, Banks, SA, Conrad, BP, Tremolada, G, et al. In-vivo three-dimensional knee kinematics during daily activities in dogs. J Orthop Res. (2015) 33:1603–10. doi: 10.1002/jor.22927
38. Tinga, S, Kim, SE, Banks, SA, Jones, SC, Park, BH, Pozzi, A, et al. Femorotibial kinematics in dogs with cranial cruciate ligament insufficiency: a three-dimensional in-vivo fluoroscopic analysis during walking. BMC Vet Res. (2018) 14:85. doi: 10.1186/s12917-018-1395-2
39. Smoger, LM, Fitzpatrick, CK, Clary, CW, Cyr, AJ, Maletsky, LP, Rullkoetter, PJ, et al. Statistical modeling to characterize relationships between knee anatomy and kinematics. J Orthop Res. (2015) 33:1620–30. doi: 10.1002/jor.22948
40. Swiderski, JK, Radecki, SV, Park, RD, and Palmer, RH. Comparison of radiographic and anatomic femoral varus angle measurements in normal dogs. Vet Surg. (2008) 37:43–8. doi: 10.1111/j.1532-950X.2007.00347.x
41. Brühschwein, A, Burg-Personnaz, J, Zöllner, M, Reese, S, and Meyer-Lindenberg, A. Comparison of CT-measured angles of pelvic limbs without patellar luxation of six canine breeds. Front Vet Sci. (2023) 10:1194167. doi: 10.3389/fvets.2023.1194167
42. Cunningham, DP, Mostafa, AA, Gordan-Evans, WJ, Boudrieau, RJ, and Griffon, DJ. Factors contributing to the variability of a predictive score for cranial cruciate ligament deficiency in Labrador retrievers. BMC Vet Res. (2017) 13:235. doi: 10.1186/s12917-017-1154-9
43. Ševčík, K, Hluchý, M, Ševčíková, M, Domaniža, M, and Ledecký, V. Inter- and intra-observer variations in radiographic evaluation of pelvic limbs in Yorkshire terriers with cranial cruciate ligament rupture and patellar luxation. Vet Sci. (2022) 9:179. doi: 10.3390/vetsci9040179
44. Pfeiffer, TR, Burnham, JM, Kanakamedala, AC, Hughes, JD, Zlotnicki, J, and Popchak, A. Distal femur morphology affects rotatory knee instability in patients with anterior cruciate ligament ruptures. Knee Surg Sports Traumatol Arthrosc. (2019) 27:1514. doi: 10.1007/s00167-018-5269-8
45. He, M, and Li, J. Increased lateral femoral condyle ratio measured by MRI is associated with higher risk of noncontact anterior cruciate ligament injury. BMC Musculoskelet Disord. (2022) 23:190. doi: 10.1186/s12891-022-05134-x
46. Yasukawa, S, Edamura, K, Tanegashima, K, Seki, M, Teshima, K, and Asano, K. Evaluation of bone deformities of the femur, tibia, and patella in toy poodles with medial patellar luxation using computed tomography. Vet Comp Orthop Traumatol. (2016) 29:29–38. doi: 10.3415/VCOT-15-05-0089
47. Olimpo, M, Piras, LA, and Peirone, B. Pelvic limb alignment in small breed dogs:a comparison between affected and free subjects from medial patellar luxation. Vet Ital. (2015) 52:45–50. doi: 10.12834/VetIt.71.206.3
48. Sasaki, A, Hidaka, Y, Mochizuki, M, and Honnami, M. Computed tomographic measurements of the sulcus angle of the femoral trochlea in small-breed dogs with and without medial patellar luxation. Vet Comp Orthop Traumatol. (2022) 35:314–20. doi: 10.1055/s-0042-1749151
49. Kara, ME, Sevil-Kilimci, F, Dilek, ÖG, and Onar, V. Proximal and distal alignment of normal canine femurs: A Morphometric analysis. Ann Anat. (2018) 217:125–8. doi: 10.1016/j.aanat.2018.02.006
50. Vedrine, B, Guillemot, A, Fontaine, D, Ragetly, GR, and Etchepareborde, S. Comparative anatomy of the proximal tibia in healthy Labrador retrievers and Yorkshire terriers. Vet Comp Orthop Traumatol. (2013) 26:266–70. doi: 10.3415/VCOT-12-02-0018
51. Aertsens, A, Rincon Alvarez, J, Poncet, CM, Beaufrère, H, and Ragetly, GR. Comparison of the tibia plateau angle between small and large dogs with cranial cruciate ligament disease. Vet Comp Orthop Traumatol. (2015) 28:385–90. doi: 10.3415/VCOT-14-12-0180
52. Van Houtte, J, Gao, X, Sijbers, J, and Zheng, G. 2D/3D registration with a statistical deformation model prior using deep learning. in: 2021 IEEE EMBS International Conference on Biomedical and Health Informatics (BHI). (2021) 1–4.
Keywords: statistical shape model, SSM, morphology, canine, femur, tibia
Citation: Huart J, Pozzi A, Bleedorn J, Lu T-W, Knell S and Park B (2024) Statistical shape modeling of the geometric morphology of the canine femur, tibia, and patella. Front. Vet. Sci. 11:1366827. doi: 10.3389/fvets.2024.1366827
Edited by:
Carmel T. Mooney, University College Dublin, IrelandReviewed by:
Adam Biedrzycki, University of Florida, United StatesMatthew Allen, University of Cambridge, United Kingdom
Copyright © 2024 Huart, Pozzi, Bleedorn, Lu, Knell and Park. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Brian Park, bpark@vetclinics.uzh.ch