Morphologic analysis of the subtalar joint using statistical shape modeling

Influence of Talar and Calcaneal Morphology on Subtalar Kinematics During Walking

BACKGROUND

Cadaver biomechanical testing suggests that the morphology of articulating bones contributes to the stability of the joints and determines their kinematics; however, there are no studies examining the correlation between bone morphology and kinematics of the subtalar joint. The purpose of this study was to investigate the influence of talar and calcaneal morphology on subtalar kinematics during walking in healthy individuals.

METHODS

Forty ankles (20 healthy subjects, 10 women/10 men) were included. Participants walked at a self-selected pace while synchronized biplane radiographs of the hindfoot were acquired at 100 images per second during stance. Motion of the talus and calcaneus was tracked using a validated volumetric model-based tracking process, and subtalar kinematics were calculated. Talar and calcaneal morphology were evaluated using statistical shape modeling. Pearson correlation coefficients were used to assess the relationship between subtalar kinematics and the morphology features of the talus and calcaneus.

RESULTS

This study found that a shallower posterior facet of the talus was correlated with the subtalar joint being in more dorsiflexion, more inversion, and more internal rotation, and higher curvature in the posterior facet was correlated with more inversion and eversion range of motion during stance. In the calcaneus, a gentler slope of the middle facet was correlated with greater subtalar inversion.

CONCLUSION

The morphology of the posterior facet of the talus was found to a primary factor driving multiplanar subtalar joint kinematics during the stance phase of gait.

CLINICAL RELEVANCE

This new knowledge relating form and function in the hindfoot may assist in identifying individuals susceptible to subtalar instability and in improving implant design to achieve desired kinematics after surgery.

Human in vivo midtarsal and subtalar joint kinematics during walking, running and hopping

The interaction among joints of the midtarsal complex and subtalar joint is important for locomotor function; however, its complexity poses substantial challenges in quantifying the joints' motions. We determine the mobility of these joints across locomotion tasks and investigate the influence of individual talus morphology on their motion. Using highly accurate biplanar videoradiography, three-dimensional bone kinematics were captured during walking, running and hopping. We calculated the axis of rotation of the midtarsal complex and subtalar joint for the landing and push-off phases. A comparison was made between these rotation axes and the morphological subtalar axis. Measurement included total rotation about and the orientation of the rotation axes in the direction of the subtalar joint and its deviation via spatial angles for both phases. The rotation axes of all three bones relative to the talus closely align with the morphological subtalar axis. This suggests that the midtarsal and subtalar joints' motions might be described by one commonly oriented axis. Despite having such an axis, the location of the axes and ranges of motion differed among the bones. Our results provide a novel perspective of healthy foot function across different sagittal plane-dominant locomotion tasks underscoring the importance of quantifying midtarsal complex and subtalar motion while accounting for an individual's talus morphology.

A Correspondence-Based Network Approach for Groupwise Analysis of Patient-Specific Spatiotemporal Data

Methods for statistically analyzing patient-specific data that vary both spatially and over time are currently either limited to summary statistics or require elaborate surface registration. We propose a new method, called correspondence-based network analysis, which leverages particle-based shape modeling to establish correspondence across a population and preserve patient-specific measurements and predictions through statistical analysis. Herein, we evaluated this method using three published datasets of the hip describing cortical bone thickness of the proximal femur, cartilage contact stress, and dynamic joint space between control and patient cohorts to evaluate activity- and group-based differences, as applicable, using traditional statistical parametric mapping (SPM) and our proposed spatially considerate correspondence-based network analysis approach. The network approach was insensitive to correspondence density, while the traditional application of SPM showed decreasing area of the region of significance with increasing correspondence density. In comparison to SPM, the network approach identified broader and more connected regions of significance for all three datasets. The correspondence-based network analysis approach identified differences between groups and activities without loss of subject and spatial specificity which could improve clinical interpretation of results.

Regional plantar forces and surface geometry variations of a chronic ankle instability population described by statistical shape modelling

BACKGROUND

Understanding detailed foot morphology as well as regional plantar forces could provide insight into foot function and provide recommendation for footwear design for chronic ankle instability (CAI) people.

RESEARCH QUESTION

This study presented 3-dimensional statistical shape models of feet from three different populations including CAI, copers and healthy individuals, with regional plantar forces also acquired.

METHODS

Sixty-six males (22 participants per group) were included in this study to capture 3-dimensional foot shapes under a standing condition and regional plantar forces during a cutting maneuver. Principal component analysis was performed to generate a mean foot shape of each group as well as modes of variations. A generalized procrustes analysis was used to achieve rapid registration of mean shapes. Besides, regional plantar forces and contact duration among these three populations were compared.

RESULTS

For 3-dimensional foot shapes, although no significant differences of the average distance between each mode and mean shape were found among three populations, there were subtle variations in mean shapes. The CAI population presented a more bulging of the lateral malleolus; copers were characterized by the flexion of the lesser toes, a more bulging of the medial foot in the sagittal plane; and healthy individuals showed a greater heel width and a more bulging of the heel in the sagittal plane. In terms of plantar forces, healthy individuals had significantly greater summated plantar forces and greater plantar forces in the lateral heel area during the early contact phase compared to copers and CAI participants.

SIGNIFICANCE

Overall, this study suggested that repetitive ankle sprains may lead to the bulging of the lateral malleolus. Further, CAI and copers seem to stabilize the ankle joint by medially shifting the center of pressure compared to healthy individuals under the static and less challenging dynamic conditions.

Unique shape variations of hind and midfoot bones in flatfoot subjects-A statistical shape modeling approach

Flatfoot deformity is a prevalent hind- and midfoot disorder. Given its complexity, single-plane radiological measurements omit case-specific joint interaction and bone shape variations. Three-dimensional medical imaging assessment using statistical shape models provides a complete approach in characterizing bone shape variations unique to flatfoot condition. This study used statistical shape models to define specific bone shape variations of the subtalar, talonavicular, and calcaneocuboid joints that characterize flatfoot deformity, that differentiate them from healthy controls. Bones of the aforementioned joints were segmented from computed tomography scans of 40 feet. The three-dimensional hindfoot alignment angle categorized the population into 18 flatfoot subjects (≥7° valgus) and 22 controls. Statistical shape models for each joint were defined using the entire study cohort. For each joint, an average weighted shape parameter was calculated for each mode of variation, and then compared between flatfoot and controls. Significance was set at p < 0.05, with values between 0.05 ≤ p < 0.1 considered trending towards significance. The flatfoot population showed a more adducted talar head, inferiorly inclined talar neck, and posteriorly orientated medial subtalar articulation compare to controls, coupled with more navicular eversion, shallower navicular cup, and more prominent navicular tuberosity. The calcaneocuboid joint presented trends of a more adducted calcaneus, more abducted cuboid, narrower calcaneal roof, and less prominent cuboid beak compared to controls. Statistical shape model analysis identified unique shape variations which may enhance understanding and computer-aided models of the intricacies of flatfoot, leading to better diagnosis and, ultimately, surgical treatment.

Biomechanical Insights Afforded by Shape Modeling in the Foot and Ankle

Advancements in volumetric imaging makes it possible to generate high-resolution three-dimensional reconstructions of bones in throughout the foot and ankle. The use of weightbearing computed tomography allows for the analysis of joint relationships in a consistent natural position that can be used for statistical shape modeling. Using statistical shape modeling, a population-based statistical model is created that can be used to compare mean bone shape morphology and identify anatomical modes of variation. A review is presented to highlight the current work using statistical shape modeling in the foot and ankle with a future view of the impact on clinical care.

Statistical multi-level shape models for scalable modeling of multi-organ anatomies

Statistical shape modeling is an indispensable tool in the quantitative analysis of anatomies. Particle-based shape modeling (PSM) is a state-of-the-art approach that enables the learning of population-level shape representation from medical imaging data (e.g., CT, MRI) and the associated 3D models of anatomy generated from them. PSM optimizes the placement of a dense set of landmarks (i.e., correspondence points) on a given shape cohort. PSM supports multi-organ modeling as a particular case of the conventional single-organ framework via a global statistical model, where multi-structure anatomy is considered as a single structure. However, global multi-organ models are not scalable for many organs, induce anatomical inconsistencies, and result in entangled shape statistics where modes of shape variation reflect both within- and between-organ variations. Hence, there is a need for an efficient modeling approach that can capture the inter-organ relations (i.e., pose variations) of the complex anatomy while simultaneously optimizing the morphological changes of each organ and capturing the population-level statistics. This paper leverages the PSM approach and proposes a new approach for correspondence-point optimization of multiple organs that overcomes these limitations. The central idea of multilevel component analysis, is that the shape statistics consists of two mutually orthogonal subspaces: the within-organ subspace and the between-organ subspace. We formulate the correspondence optimization objective using this generative model. We evaluate the proposed method using synthetic shape data and clinical data for articulated joint structures of the spine, foot and ankle, and hip joint.

Comprehensive personalized ankle joint shape analysis of children with cerebral palsy from pediatric MRI

Cerebral palsy, a common physical disability in childhood, often causes abnormal patterns of movement and posture. To better understand the pathology and improve rehabilitation of patients, a comprehensive bone shape analysis approach is proposed in this article. First, a group analysis is performed on a clinical MRI dataset using two state-of-the-art shape analysis methods: ShapeWorks and a voxel-based method relying on Advanced Normalization Tools (ANTs) registration. Second, an analysis of three bones of the ankle is done to provide a complete view of the ankle joint. Third, a bone shape analysis is carried out at subject level to highlight variability patterns for personnalized understanding of deformities.

Coverage maps demonstrate 3D Chopart joint subluxation in weightbearing CT of progressive collapsing foot deformity

A key element of the peritalar subluxation (PTS) seen in progressive collapsing foot deformity (PCFD) occurs through the transverse tarsal joint complex. However, the normal and pathological relations of these joints are not well understood. The objective of this study to compare Chopart articular coverages between PCFD patients and controls using weight-bearing computed tomography (WBCT). In this retrospective case control study, 20 patients with PCFD and 20 matched controls were evaluated. Distance and coverage mapping techniques were used to evaluate the talonavicular and calcaneocuboid interfaces. Principal axes were used to divide the talar head into 6 regions (medial/central/lateral and plantar/dorsal) and the calcaneocuboid interface into 4 regions. Repeated selections were performed to evaluate reliability of joint interface identification. Surface selections had high reliability with an ICC > 0.99. Talar head coverage decreases in plantarmedial and dorsalmedial (- 79%, p = 0.003 and - 77%, p = 0.00004) regions were seen with corresponding increases in plantarlateral and dorsolateral regions (30%, p = 0.0003 and 21%, p = 0.002) in PCFD. Calcaneocuboid coverage decreased in plantar and medial regions (- 12%, p = 0.006 and - 9%, p = 0.037) and increased in the lateral region (13%, p = 0.002). Significant subluxation occurs across the medial regions of the talar head and the plantar medial regions of the calcaneocuboid joint. Coverage and distance mapping provide a baseline for understanding Chopart joint changes in PCFD under full weightbearing conditions.

Summary Report of the Arthritis Foundation and the American Orthopaedic Foot & Ankle Society's Symposium on Targets for Osteoarthritis Research: Part 2: Treatment Options

This second of a 2-part series of articles recounts the key points presented in a collaborative symposium sponsored jointly by the Arthritis Foundation and the American Orthopaedic Foot & Ankle Society with the intent to survey current treatment options for osteoarthritis (OA) of the foot and ankle. A meeting was held virtually on December 10, 2021. A group of experts were invited to present brief synopses of the current state of knowledge and research in this area. Topics were chosen by meeting organizers, who then identified and invited the expert speakers. Part 2 overviews the current treatment options, including orthotics, non-joint destructive procedures, as well as arthroscopies and arthroplasties in ankles and feet. Opportunities for future research are also discussed, such as developments in surgical options for ankle and the first metatarsophalangeal joint. The OA scientific community, including funding agencies, academia, industry, and regulatory agencies, must recognize the importance to patients of addressing the foot and ankle with improved basic, translational, and clinical research.

LEVEL OF EVIDENCE

Level V, review article/expert opinion.

Interaction of loading and ligament injuries in subtalar joint instability quantified by 3D weightbearing computed tomography

Despite decades of research since its first description, subtalar joint instability remains a diagnostic enigma within the concept of hindfoot instability. This could be attributed to current imaging techniques, which are impeded by two-dimensional measurements. Therefore, we used weightbearing computed tomography imaging to quantify three-dimensional displacement associated with subtalar joint instability. Three-dimensional models were generated in seven paired cadaver specimens to compute talocalcaneal displacement after different patterns of axial load (85 kg) combined with torque in internal and external rotation (10 Nm). Sequential imaging was repeated in the subtalar joint containing intact ligaments to determine reference displacement. Afterward, the interosseus talocalcaneal ligament (ITCL) or calcaneofibular ligament (CFL) was sectioned, then the ITCL with CFL and after the ITCL, CFL with the deltoid ligament (DL). The highest translation could be detected in the dorsal direction and the highest rotation occurred in the internal direction when external torque was applied to the foot without load. These displacements differed significantly from the condition containing intact ligaments, with a mean difference of 1.6 mm (95% CI, 1.3 to 1.9) for dorsal translation and a mean of 12.4° (95% CI, 10.1 to 14.8) for internal rotation. Clinical relevance: Our study provides a novel and noninvasive analysis to quantify subtalar joint instability based on three-dimensional WBCT imaging. This approach overcomes former studies using trans-osseous fixation to determine three-dimensional subtalar joint displacement and implements an imaging device and software modalities that are readily available. Based on our findings, we recommend applying torque in external rotation to the foot to optimize the detection of subtalar joint instability.

Assessment of Progressive Collapsing Foot Deformity Using Semiautomated 3D Measurements Derived From Weightbearing CT Scans

BACKGROUND

In progressive collapsing foot deformity (PCFD), hind- and midfoot deformities can be hard to characterize based on weightbearing plain radiography. Semiautomated 3-dimensional (3D) measurements derived from weightbearing computed tomography (WBCT) scans may provide a more accurate deformity assessment. In the present study, automated 3D measurements based on WBCT were used to compare hindfoot alignment of healthy individuals to patients with PCFD.

METHODS

The WBCT scans of 20 patients treated at our institution with either a flexible (N = 10) or rigid (N = 10) PCFD were compared with the WBCT scans of a control group of 30 healthy individuals. Using semiautomated image analysis software, from each set of 3D voxel images, we measured the talar tilt (TT), hindfoot moment arm (HMA), talocalcaneal angle (TCA; axial/lateral), talonavicular coverage (TNC), and talocalcaneal overlap (TCO). The presence of medial facet subluxation as well as sinus tarsi/subfibular impingement was additionally assessed.

RESULTS

With the exception of the TCA (axial/lateral), the analyzed measurements differed between healthy individuals and patients with PCFD. The TCA axial correlated with the TNC in patients with PCFD. An increased TCO combined with sinus tarsi impingement raised the probability of predicting a deformity as rigid.

CONCLUSION

Using 3D measurements, in this relatively small cohort of patients, we identified relevant variables associated with a clinical presentation of flexible or rigid PCFD. An increased TCO combined with sinus tarsi impingement raised the probability of predicting a deformity as rigid. Such WBCT-based markers possibly can help the surgeon in decision-making regarding the appropriate surgical strategy (eg, osteotomies vs realignment arthrodesis). However, prospective studies are necessary to confirm the utility of the proposed parameters in the treatment of PCFD.

LEVEL OF EVIDENCE

Level III, case-control study.

The Assessment of Ankle Osteoarthritis with Weight-Bearing Computed Tomography

The standard for diagnostic radiographic imaging in foot and ankle surgery was until 2012 radiographs with full weight-bearing without any useful alternative. Weight-bearing cone-beam computed tomography (WBCT) was introduced 2012 for foot and ankle use as a new technology that allows 3D imaging with full weight-bearing which should be not influenced by projection and/or foot orientation. The assessment of ankle osteoarthritis with WBCT including the description of healthy status, effect of alignment and7or (in)stability is extensively illustrated in this review article.

Multi-level multi-domain statistical shape model of the subtalar, talonavicular, and calcaneocuboid joints

Traditionally, two-dimensional conventional radiographs have been the primary tool to measure the complex morphology of the foot and ankle. However, the subtalar, talonavicular, and calcaneocuboid joints are challenging to assess due to their bone morphology and locations within the ankle. Weightbearing computed tomography is a novel high-resolution volumetric imaging mechanism that allows detailed generation of 3D bone reconstructions. This study aimed to develop a multi-domain statistical shape model to assess morphologic and alignment variation of the subtalar, talonavicular, and calcaneocuboid joints across an asymptomatic population and calculate 3D joint measurements in a consistent weightbearing position. Specific joint measurements included joint space distance, congruence, and coverage. Noteworthy anatomical variation predominantly included the talus and calcaneus, specifically an inverse relationship regarding talar dome heightening and calcaneal shortening. While there was minimal navicular and cuboid shape variation, there were alignment variations within these joints; the most notable is the rotational aspect about the anterior-posterior axis. This study also found that multi-domain modeling may be able to predict joint space distance measurements within a population. Additionally, variation across a population of these four bones may be driven far more by morphology than by alignment variation based on all three joint measurements. These data are beneficial in furthering our understanding of joint-level morphology and alignment variants to guide advancements in ankle joint pathological care and operative treatments.

Assessment of Hindfoot Alignment Comparing Weightbearing Radiography to Weightbearing Computed Tomography

BACKGROUND

Hindfoot alignment view (HAV) radiographs are widely utilized for 2-dimensional (2D) radiographic assessment of hindfoot alignment; however, the development of weightbearing computed tomography (WBCT) may provide more accurate methods of quantifying 3-dimensional (3D) hindfoot alignment. The aim of this study was to compare the 2D calcaneal moment arm measurements on HAV radiographs with WBCT.

METHODS

This retrospective cohort study included 375 consecutive patients with both HAV radiographs and WBCT imaging. Measurement of the 2D hindfoot alignment moment arm was compared between both imaging modalities. The potential confounding influence of valgus/varus/neutral alignment, presence of hardware, and motion artifact were further analyzed.

RESULTS

The intraclass correlation coefficients (ICCs) of interobserver and intraobserver reliability for measurements with both imaging modalities were excellent. Both modalities were highly correlated (Spearman coefficient, 0.930; P < .001). HAV radiographs exhibited a mean calcaneal moment arm difference of 3.9 mm in the varus direction compared with WBCT (95% CI, -4.9 to 12.8). The difference of hindfoot alignment between both modalities was comparable in subgroups with neutral/valgus/varus alignment, presence of hardware, and motion artifact.

CONCLUSION

Both HAV radiographs and WBCT are highly reliable and highly correlated imaging methods for assessing hindfoot alignment. Measurements were not influenced by severe malalignment, the presence of hardware, or motion artifact on WBCT. On average, HAV radiographs overestimated 3.9 mm of varus alignment as compared with WBCT.

LEVEL OF EVIDENCE

Level III, retrospective comparative study.

Assignment of local coordinate systems and methods to calculate tibiotalar and subtalar kinematics: A systematic review

The introduction of biplane fluoroscopy has created the ability to evaluate in vivo motion, enabling six degree-of-freedom measurement of the tibiotalar and subtalar joints. Although the International Society of Biomechanics defines a standard method of assigning local coordinate systems for the ankle joint complex, standards for the tibiotalar and subtalar joints are lacking. The objective of this systematic review was to summarize and appraise the existing literature that (1) defined coordinate systems for the tibia, talus, and/or calcaneus or (2) assigned kinematic definitions for the tibiotalar and/or subtalar joints. A systematic literature search was developed with search results limited to English Language from 2006 through 2020. Articles were screened by two independent reviewers based on title and abstract. Methodological quality was evaluated using a modified assessment tool. Following screening, 52 articles were identified as having met inclusion criteria. Methodological assessment of these articles varied in quality from 61 to 97. Included articles adopted primary methods for defining coordinate systems that included: (1) anatomical coordinate system (ACS) based on individual bone landmarks and/or geometric shapes, (2) orthogonal principal axes, and (3) interactive closest point (ICP) registration. Common methods for calculating kinematics included: (1) joint coordinate system (JCS) to calculate rotation and translation, (2) Cardan/Euler sequences, and (3) inclination and deviation angles for helical angles. The methods each have strengths and weaknesses. This summarized knowledge should provide the basis for the foot and ankle biomechanics community to create an accepted standard for calculating and reporting tibiotalar and subtalar kinematics.