Skip to main content

ORIGINAL RESEARCH article

Front. Public Health, 10 August 2022
Sec. Aging and Public Health
This article is part of the Research Topic Insights in Aging and Public Health: 2022 View all 21 articles

Identification of radiographic characteristics associated with pain in hallux valgus patients: A preliminary machine learning study

\nChenGuang Wang*&#x;ChenGuang Wang1*Chao Li&#x;Chao Li2Rui Zhang&#x;Rui Zhang1ZhiJun LiZhiJun Li1HuaFeng ZhangHuaFeng Zhang1Yuan ZhangYuan Zhang1Shen LiuShen Liu1XiaoYue ChiXiaoYue Chi1Rui Zhao
Rui Zhao1*
  • 1Department of Orthopedics Surgery, Tianjin Medical University General Hospital, Tianjin, China
  • 2Department of Orthopedics, Xiangyang Central Hospital, Affiliated Hospital of Hubei University of Arts and Science, Xiangyang, China

Objective: To investigate the association between the structural deformity and foot pain in hallux valgus (HV) patients using a multi-variate pattern analysis (MVPA) approach.

Methods: Plain radiographic metrics were calculated from 36 painful and 36 pain-free HV feet. In analysis 1, univariate analyses were performed to investigate the clinical and radiographic differences between painful and pain-free HV. In analysis 2, we investigated the pattern differences for radiographic metrics between these two groups using a MVPA approach utilizing a support vector machine. In analysis 3, sequential backward selection and exhaustive search were performed as a feature-selection procedure to identify an optimal feature subtype. In analysis 4, hierarchical clustering analysis was used to identify the optimal radiographic HV subtype associated with pain in HV.

Results: We found that: (1) relative to feet with pain-free HV, the painful ones exhibited a higher hallux valgus angle, i.e., the magnitude of distal metatarsal and phalangeal deviation; (2) painful HV could be accurately differentiated from pain-free HV via MVPA. Using sequential backward selection and exhaustive search, a 5-feature subset was identified with optimal performance for classifying HV as either painful or pain-free; and (3) by applying hierarchical clustering analysis, a radiographic subtype with an 80% pain incidence was identified.

Conclusion: The pain in HV is multifactorial and associated with a radiographic pattern measured by various angles on plain radiographs. The combination of hallux valgus angle, inter-phalangeal angle, distal metatarsal articular angle, metatarsal cuneiform angle and metatarsal protrusion distance showed the optimal classification performance between painful and pain-free HV.

Introduction

Hallux valgus (HV) is characterized as a combined deformity with mal-positioning of the first metatarsophalangeal joint caused by a lateral deviation of the great toe and a medial deviation of the first metatarsal bone. This can progress to significant metatarsophalangeal (MTP) joint subluxation, which is frequently encountered in foot and ankle clinic with a prevalence exceeding 20% in adults over 50 years of age (1). HV is often associated with the pain (2), impaired physical function (3) with consequent poorer general health (4) and diminished quality of life (5). Pain is the primary complaint related to HV (6), and typically the severity of the deformity directly correlates with the degree of pain. A prior publication noted a positive correlation between self-reported severity of HV (Short Form 36 or Foot Pain and Disability Index) and increasing levels of pain (5). However, there have also been conflicting reports of only a modest correlation between pain and structural HV deformity (i.e., HV angle measured from a radiographic image). Moreover, there have been contradictory reports concerning the link between foot pain and HV (5, 79). Discrepancies have been noted such as patients with a large HV angle yet without any pain, and others with only a mild HV angle that experience severe foot pain. The precise mechanism behind such phenomena has yet to be identified. In recent years, some authors have postulated alternative explanations that HV has affected the entire foot rather than just the MTP joint (1012). When pain could not be explained by the degree of the HV angle alone, which has been traditionally used to evaluate the deformity of MTP joint, other radiographic angles, such as the intermetatarsal angle, and the metatarsal cuneiform angle have been implicated as potentially associated with pain in HV patients.

In recent decades, with the development of artificial intelligence and machine learning algorithms, multi-variate pattern analysis (MVPA) has been widely applied in the medical imaging field (1315). MVPA methods provide a rich characterization of imaging data, often in a data-driven manner. This technique has enabled researchers to identify the key, significantly important parameters of a complex disease in a multi-variate manner rather than exploring univariate differences between the disease and healthy controls (15). Prior studies exploring the correlates of foot pain in HV patients have relied only on univariate analyses to determine the association between pain and HV deformity (e.g., two sample T-test or univariate correlation analyses between radiographic findings and degree of pain). Although the results of univariate analysis were straight-forward and easy to interpret, only linear-relations between the radiographic findings and clinical assessments were considered. Pattern information consisting of different radiographic metrics to evaluate the deformity of the entire foot was largely overlooked. More importantly, HV affected the structure of the entire foot rather than just the MTP joint (1012). It has been shown that hallux valgus could affect the transverse arch structure and its force loading patterns (12). Furthermore, HV not only alters the geometry structure of forefoot but also causes biomechanics function changes of the foot (16, 17). A more comprehensive understanding of the deformity pattern associated with pain in HV would likely provide the basis for foot surgeons to develop novel surgical techniques to correct the HV deformity and resolve the pain.

Therefore, in this study, we first (1) conducted a MVPA to test whether the various radiographic angles could be used to classify painful HV patients from pain-free HV patients in a supervised manner; (2) performed a sequential backward selection and exhaustive search as a feature selection procedure to identify the optimal feature subset of radiographic features; and (3) carried out hierarchical clustering in an unsupervised manner to identify the radiographic HV subtype which was associated with pain in HV patients.

Materials and methods

Subjects

Written informed consent was obtained from each participant prior to each procedure. Ethical approval was granted by the Local Institutional Ethics Committee. The study comprised a total of 103 feet with HV, from 81 patients, recruited from January 2021 to January 2022 at both inpatient and outpatient departments of foot and ankle surgery. To control the degree of overall HV deformity, we used a HV angle above 20° as a cut-off value for patient recruitment rather than a traditional criterion for HV angle of >15° determined by standard radiographic images.

The exclusion criteria were as following: (1) history of foot or ankle surgery or trauma; (2) history of neurological diseases; (3) ankle joint or foot sprain within 3 months; (4) presence of concomitant hallux limitus (i.e., a minimum of 50 degrees of passive dorsiflexion at the first MTP joint was required) (2).

Clinical assessment

The demographic data of the participants including sex, age, education years, ethnicity was obtained by questionnaire. The total years of education was determined from the patient's self-report of completed years of education. The Visual Analog Scale (VAS) was implemented to assess the degree of foot pain in each HV foot. Feet with a VAS score of 0 were identified as pain-free. Feet with VAS scores >3 extending more than 3 months were classified as painful (18, 19). Painless feet were selected from the dataset and excluded those that meet the exclusion criteria. Therefore, 36 painless HV feet were obtained. To further performed machine learning analyses and to avoid class-imbalance problem, 36 painful HV feet from the dataset to match the painless HV feet with age, gender, education years and ethnicity. Finally, 72 HV feet from 63 patients, 36 painful and 36 pain-free, were included in our study.

Radiographic assessment

Bilateral weight bearing radiographs were obtained for all participants using a standardized procedure (tube to film distance 100 cm, angled 15° from vertical). All radiographic assessments were carried out in ImageJ toolbox. The assessment procedure of radiographic metrics we used were as follows. (1) Hallux Valgus angle (HVA): the hallux valgus angle is formed by the longitudinal axis of the first proximal phalanx and the longitudinal axis of the first metatarsus (20, 21). (2) Inter-Metatarsal Angle (IMA): the IMA, or metatarsus primus adducts angle, is the angle between the longitudinal axes of the first and second metatarsals (22). (3) Inter-Phalangeal Angle (IPA): the IPA is the angle between longitudinal axes of the proximal and distal phalanges (22). (4) Distal Metatarsal Articular Angle (DMAA): the DMAA is the angle between the perpendicular to the effective articular surface of the first metatarsal head and the longitudinal axis of the first metatarsal bone (22). (5) Distal Articular Set Angle (DASA): the DASA is the angle between the perpendicular to the effective articular surface of the 1st proximal phalanx and the longitudinal axis of the 1st proximal phalanx (22, 23). (6) Metatarsal Cuneiform Angle (MCA): the MCA is the measured angle between the longitudinal axis of the first metatarsal bone and the longitudinal axis of the medial cuneiform bone. (7) Metatarsal Adducent Angle (MAA): we used the method devised by Engel et al. (24) to measure the MAA. We measured the angle between the axis of the second metatarsal and the axis perpendicular to the transverse axis of the tarsus using the most lateral and most distal points of the joints of the cuboid, with the fifth metatarsal as a reference. (8) Metatarsal Protrusion Distance (MPD): located at the most distal extent of the second metatarsal bone, perpendicular to the longitudinal axis, this line reflects the distal-most protrusion of the second metatarsal. A line parallel to this is now constructed coursing through the distal-most point of the first metatarsal. The difference between these 2 parallel lines is the MPD (25). (9) Sesamoid Position (SP): The position of the tibial sesamoid was graded I, II, III, IV, V, VI, or VII according to its position relative to the functional longitudinal axis of the first metatarsal (26).

Analysis 1: Univariate analyses for exploring radiographic differences between painful and pain-free HV feet

Two sample t-tests were performed for continuous variables and chi-square test was used for categorical variables to compare the baseline demographic data between two groups. Subsequently, Shapro-Wilk tests were performed to check the normal distribution of the radiographic measurements and QQ plots were also generated for illustrating normality of our data. All statistical analyses and figures were generated using GraphPad prism 9. Two sample t-tests were performed for continuous variables and chi-square test was used for categorical variables to compare the radiographic metrics between groups. Moreover, parametric test was performed for between group comparison for radiographic metrics which passed the Shapiro-Wilk test, otherwise nonparametric test was performed (i.e., Mann-Whitney test). The significant level was set to P < 0.05. Moreover, to elucidate the associations among clinical measures. Pearson correlation analyses were performed to demonstrate the association among radiographic metrics.

Analysis 2: Multi-variate pattern analyses (MVPA) using a support vector machine (SVM) for classifying painful HV feet from pain-free ones

The SVM was use to distinguish painful HV feet from pain-free ones using radiographic features according to the homemade MATLAB script based on LibSVM's (http://www.csie.ntu.edu.tw/~cjlin/libsvm/) implementation of linear SVM using the default parameters: We used linear kernel, the penalty coefficient (c) was set to 1, other parameters were not adjustable therefore we made no further adjustment. To overcome the loss of generalization due to the relative sample size, the Leave-one-out-cross-validation (LOOCV) technique was used in this study. One of the available data-points in LOOCV was held-out, and the model was trained using the remainder of the data in the dataset and tested against the held-out data. This procedure was repeated until all data-points had been held-out once as the testing data. Because nearly the entire dataset was used for training, and the trained model was close to the real one, we expected the LOOCV error bias to be small. The corresponding P-value for the classification for painful HV vs. pain-free HV was calculated from the null distribution obtained from 10,000 random permutation tests by randomly shuffling the labels of samples in the dataset. The p-values were calculated as a proportion of the number of permutations generated that were greater than or equal to actual classification accuracy, and the total number of permutations. If none of 10,000 permutations reached the actual accuracy, the P-value was P < 0.0001.

Analysis 3: Sequential backward selection and exhaustive search for identifying an optimal feature subset via MVPA

The sequential backward selection procedure was initiated from the full set of features; at each iteration we compared the performances of different models built by sequentially removing each of the features from the current set of candidate features. We then excluded the feature whose subtraction resulted in the most increase in classification accuracy from the next iteration the feature. This procedure was repeated until the classification accuracy showed no increase by excluding features.

However, using sequential backward selection could end up with a global optimum rather than a local optimal feature subset. The exhaustive search was also performed to enumerate all possible models for classifying painful HV vs. pain-free HV. The P-value for the classification was also generated using permutation test as described in analysis 2.

Analysis 4: Hierarchical clustering analysis for identifying radiographic HV subtype associated with pain in HV

Hierarchical clustering of the 72 HV samples was performed based on homemade MATLAB script using Euclidean distance as distance measure. All other parameters were set to default values. The dendrogram along with silhouette index, Calinski-Harabasz index and Davies-Bouldin index were conducted to determine the optimal number of clusters. Subsequently, the mean of radiographic metrics of each subtype were illustrated in a radar map along with the X-ray of the most representative sample of each subtype (e.g., with the closest distance to the cluster-average). Last, the prevalence of foot pain in each subtype was compared using chi-square test between subtypes.

Results

Demographic data

The demographic data and the clinical assessments of all participants are summarized in Table 1. There were no significant inter-group differences with regards to age, sex, or education years (P > 0.05).

TABLE 1
www.frontiersin.org

Table 1. Demographic data of painful and pain-free hallux valgus (HV).

Analysis 1: Univariate analyses for exploring radiographic differences between painful and pain-free HV feet

Relative to the pain-free group, the painful group exhibited a significantly higher HVA (P = 0.001, Mann-Whitney U = 364), DMAA (P = 0.001, Mann-Whitney U = 369.5) and a higher level of sesamoid subluxation (P = 0.04, Mann-Whitney U = 477). No significant differences between groups were observed in the following metrics: IMA, IPA, DASA, MCA, MAA, MPD (Figure 1). Significant positive correlations were observed between the following metrics: HVA and DMAA (R = 0.84, P < 0.001), HVA and IMA (R = 0.53, P < 0.001), IMA and DMAA (R = 0.46, P < 0.005), HVA and SP (R = 0.53, P < 0.001), IMA and SP (R = 0.55, P < 0.001), and DMAA and SP (R = 0.61, P < 0.001). Significant negative correlation was observed between IMA and MCA (R = −0.32, P < 0.05) (Figure 2), however, this negative correlation did not survive the Bonferroni correction.

FIGURE 1
www.frontiersin.org

Figure 1. Univariate differences for radiographic metrics between painful and pain-free HV. HVp, Hallux Valgus with Pain; HVnp, Hallux Valgus without Pain; HVA, Hallux Valgus Angle; IMA, Inter-Metatarsal Angle; IPA, Inter-Phalangeal Angle; DMAA, Distal Metatarsal Articular Angle; DASA, Distal Articular Set Angle; MCA, Metatarsal Cuneiform Angle; MAA, Metatarsal Adducent Angle; MPD, Metatarsal Protrusion Distance; SP, Sesamoid Position. **P < 0.005; *P < 0.05. ns means non-significant.

FIGURE 2
www.frontiersin.org

Figure 2. Associations among radiographic metrics in HV. The correlation coefficients are shown in the center of each box. HVA, Hallux Valgus Angle; IMA, Inter-Metatarsal Angle; IPA, Inter-Phalangeal Angle; DMAA, Distal Metatarsal Articular Angle; DASA, Distal Articular Set Angle; MCA, Metatarsal Cuneiform Angle; MAA, Metatarsal Adducent Angle; MPD, Metatarsal Protrusion Distance; SP, Sesamoid Position.

Analysis 2: Multi-variate pattern analyses (MVPA) using a support vector machine (SVM) for classifying painful HV feet from pain-free ones

The MVPA was performed to determine whether the pattern of radiographic features could be used to distinguish painful HV feet from pain-free ones. The painful group could be accurately distinguished from the pain-free group using the entire feature set (i.e., all 9 radiographic features) with a 76.4% classification accuracy. The corresponding P-value was 0.0054 (Figure 3). This result indicated that a pattern of radiographic features of the foot might accurately predict pain.

FIGURE 3
www.frontiersin.org

Figure 3. The null distributions for Painful HV vs. Pain-free HV. The left panel illustrates the null distribution of 10,000 permutations for Painful HV vs. Pain-free HV using the entire feature set. The right panel illustrates the null distribution of 10,000 permutations for Painful HV vs. Pain-free HV using the optimal feature set including hallux valgus angle, inter-phalangeal angle, distal metatarsal articular angle, metatarsal adducent angle, metatarsal protrusion distance.

Analysis 3: Sequential backward selection and exhaustive search for identifying an optimal feature subset via MVPA

Because Analysis 2 indicated that the painful HV could be accurately distinguished from pain-free HV using the entire feature set of radiographic measures, the next step was to define an optimal subset of radiographic features that would correlate with HV-related pain. Sequential backward selection procedure was performed, and we found a 5-feature subset with an optimum performance of 80.1% classification accuracy (P = 0.0014; Figure 3). Moreover, the exhaustive search also identified the same subset. The 5-feature subset consisted of the following 5 metrics: HVA, IPA, DMAA, MAA and MPD.

Analysis 4: Hierarchical clustering analysis for identifying radiographic HV subtype associated with pain in HV

The dendrogram showed that the optimal number of clusters was 3 (Figure 4A). The results of Silhouette index, Calinski-Harabasz index and Davies-Bouldin index analyses also showed the optimal number of clusters was 3 (Figure 4B). After hierarchical clustering, 3 subtypes were identified: subtype 1 had the largest sample size (e.g., 38 HV), subtype 3 had the smallest sample size (e.g., 9 HV), and the sample size of subtype 2 was 25. The characteristics of each subtype were illustrated in the radar map and the most representative images of each subtype were also illustrated (e.g., with the closest distance to the cluster-average; Figure 5). Furthermore, the incidence of foot pain was highest in subtype 2 (e.g., 80%; Figure 6).

FIGURE 4
www.frontiersin.org

Figure 4. In (A), the dendrogram of the hierarchical clustering is illustrated. The y-axis shows the study-specific distances between the clusters as height h. When h = 0, each volume forms its own cluster, and h = 2 corresponds to splitting the volumes into 13 clusters. The x-axis shows the number for each subject. The results of Silhouette index, Calinski-Harabasz index and Davies-Bouldin index analyses are shown in (B).

FIGURE 5
www.frontiersin.org

Figure 5. The characteristics of each subtype are illustrated in the radar map and the most representative images of each subtype are illustrated (e.g., with the closest distance to the cluster-average). HVA, Hallux Valgus Angle; IMA, Inter-Metatarsal Angle; IPA, Inter-Phalangeal Angle; DMAA, Distal Metatarsal Articular Angle; DASA, Distal Articular Set Angle; MCA, Metatarsal Cuneiform Angle; MAA, Metatarsal Adducent Angle; MPD, Metatarsal Protrusion Distance; SP, Sesamoid Position.

FIGURE 6
www.frontiersin.org

Figure 6. The incidence of foot pain in each radiographic subtype. ***P < 0.001.

Discussion

The three main findings of the present study were: (1) relative to the pain-free group, the painful group exhibited higher HVA, DMAA, and these two metrics were positively correlated; (2) the painful group could be accurately distinguished from the pain-free group using a MVPA. Using sequential backward selection and exhaustive search, an optimal subset was developed that contained 5 radiographic metrics: HVA, IPA, DMAA, MAA and MPD that had optimal performance for classifying painful and pain-free HV; and (3) by applying hierarchical clustering analysis, a radiographic subtype associated with foot pain was also identified.

Radiographic correlates of pain in HV patients illustrated by univariate analysis

Previous studies had correlated higher HVA with foot pain in HV patients, to which the present study provided further evidence that increased degrees of HV deformity were associated with worse pain in HV. This finding corroborated prior research that had shown pain intensity to be directly correlated with higher degrees of HV deformity. However, whether HV deformity is the major cause of foot pain remains controversial (2, 6, 7, 9, 10). Foot pain in patients with HV has also been associated with patient characteristics such as general health status and occupational physical activity levels rather than severity of hallux deviation (2). Our results partly support with the previous reports that pain intensity is correlated with HV deformity in HV patients. The reasons for contradiction reports were two folds: (1) the inclusion criteria vary from different studies. In our current study, we included HV patients with moderate to severe deformity (i.e., hallux valgus angle above 20°). Several past studies have a relatively wide inclusion criteria that using a HVA above 15° as an inclusion criterion (2, 6, 7). The biomechanical compensatory changes of the foot could partly neutralize the pain and might affect the association between pain intensity and degree of deformity; (2) past studies only explore the univariate association between HV deformity and pain intensity, ignored subtypes which determined by the combination of various angles measured from the HV foot. In our current study, we identified a radiographic painful subtype and two painless subtypes. Interestingly, one of the painless subtypes has a HV angle close to painful subtype, indicating that pain in HV patients is multifactorial.

The pain of HV is multifactorial, determined by multiple radiographic measurements

In past studies, a large proportion of HV patients with variability of pain intensity remained unexplained. These controversial results were mainly based on univariate analyses, with no attention to associations of other radiographic variables. In the current study, the results of clustering analyses showed that subtypes 1 and 2 exhibited no differences in HVA yet there was significant difference in pain incidence. This finding further supported the idea that pattern information of radiographic variables was crucial in correlation with HV-related pain. The direct comparison of HVA between pain and non-pain participants and liner regression are straight-forward procedures that reveal the association between pain and HVA. However, it only characterizes the quantitative liner association between HVA and pain intensity, which neglects the relationships of HVA and other radiographic measures. Using the MVPA approach and pattern information, the structural morphology of the entire foot was characterized. Recently, studies have shown that HV affected not only the first metatarsophalangeal joint but also the structure of the entire foot including the transverse arch structure (12) and sesamoid platform that led to the biomechanical changes in the foot (11). These findings indicated that the hallux valgus deformity was a multiplanar deformity that affected a wide range of foot structures and could not be adequately evaluated by HVA alone. Our current conclusion that the pattern of radiographic measures could be used to classify painful and pain-free HV also supported this concept. By applying feature selection procedures such as sequential backward selection and exhaustive search, an optimal feature subset which included HVA, IPA, DMAA, MAA and MPD was identified with the best classification performance. Therefore, the pattern of these metrics appeared to have strong correlation to foot pain in patients with HV.

Clinical significance

In our current study, we found that pain of HV was multifactorial and determined by multiple radiographic measurements. By applying sequential backward selection and exhaustive search, 5 measurements were identified as the optimal feature subset which showed optimal performance classifying painful and pain-free HV feet. It is worth mentioning that the importance of these metrics could not be ranked by our current procedure, and ranking the importance of these metrics would be partial considering our current data failed to include other pain-related factors such as soft tissue of the forefoot, which is also crucial to the pathology of HV. In clinical practice, surgical approaches were designed to correct HVA, DMAA and IMA. By applying distal soft tissue (e.g., McBride procedure) surgery and osteotomies (e.g., Chevron osteotomy, Scarf osteotomy, etc.), these deformities could be significantly corrected (i.e., HVA, IMA) and were effective to alleviate symptoms in HV patients (21). Therefore, from clinical aspect, these metrics would be more of importance than other measurements (3). The deformities measured by IPA, MAA and MPD were corrected only in patients with specific symptoms (e.g., kiss corn, metatarsus adductus and metastatic metatarsalgia) or with residual interphalangeus following metatarsal osteotomy (23). In general, our current results were consistent with clinical experience that the corrective surgery for HV should be holistic rather than only correct the HVA, DMAA and IMA. Our current results provided theoretic basis for this clinical consensus. More importantly, our results suggested that there is a pain-free HV subtype which exhibited similar HVA, DMAA and IMA to painful HV subtype. The differences between these two subtypes were IPA, MCA, indicating that when metatarsal osteotomy was unsatisfactory, other procedures for correcting IPA and MCA [e.g., Akin procedure (23) and Lapidus procedure] would be of potential efficacy for reducing clinical symptoms in HV patients. To our knowledge, the present study is the first conducted using MVPA to identify radiographic correlates of pain in HV patients. Future studies are still needed to validate our results.

Limitations

Regarding limitations in our current study: (1) we had no external validation set to verify our current results, so multi-center studies with larger sample sizes are needed; (2) we only included the axial images of plane radiographs, future studies are needed to encompass sagittal or even 3-dimensional imaging; (3) our current study is retrospective, lacking assessments of psychological factors which are also considered influential in HV patients with foot pain; (4) our current sample size is relatively small according to the Vapnik-Chervonenkis Dimension of the linear classifier. SVM is well-suited for modeling small samples with powerful predictability, and we also conducted leave-one-out-cross-validation (LOOCV) to include as much data as possible in the training set for model training to minimize the possible effect of small sample size. Still, future study is needed to enlarge the sample size.

Conclusion

The pain in HV is multifactorial and associated with a radiographic pattern measured by the various angles on plain radiographs. The combination of HVA, IPA, DMAA, MAA and MPD showed the optimal classification performance between painful and pain-free HV.

Data availability statement

The original contributions presented in the study are included in the article/Supplementary materials, further inquiries can be directed to the corresponding author/s.

Ethics statement

The studies involving human participants were reviewed and approved by Tianjin Medical Univerisity General Hospital. The patients/participants provided their written informed consent to participate in this study.

Author contributions

CW designed the study. RZhan, CL, and RZhao analyzed the data and wrote the manuscript. ZL, YZ, SL, and HZ collected the data. XC help analyzed the data. RZhao and CW revised the manuscript. All authors contributed to the article and approved the submitted version.

Acknowledgments

The authors would like to express their gratitude to EditSprings (https://www.editsprings.cn) for the expert linguistic services provided.

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.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpubh.2022.943026/full#supplementary-material

References

1. Menz HB, Marshall M, Thomas MJ, et al. Incidence and Progression of Hallux Valgus: a Prospective Cohort Study. Arthritis Care Res (Hoboken). (2021). doi: 10.1002/acr.24754. [Epub ahead of print].

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Hurn SE, Vicenzino BT, Smith MD. Correlates of foot pain severity in adults with hallux valgus: a cross-sectional study. J Foot Ankle Res. (2014) 7:32. doi: 10.1186/1757-1146-7-32

PubMed Abstract | CrossRef Full Text | Google Scholar

3. D'Arcangelo PR, Landorf KB, Munteanu SE, et al. Radiographic correlates of hallux valgus severity in older people. J Foot Ankle Res. (2010) 3:20. doi: 10.1186/1757-1146-3-20

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Mickle KJ, Munro BJ, Lord SR, et al. Cross-sectional analysis of foot function, functional ability, and health-related quality of life in older people with disabling foot pain. Arthritis Care Res (Hoboken). (2011) 63:1592–8. doi: 10.1002/acr.20578

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Menz HB, Roddy E, Thomas E, Croft PR. Impact of hallux valgus severity on general and foot-specific health-related quality of life. Arthritis Care Res (Hoboken). (2011) 63:396–404. doi: 10.1002/acr.20396

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Coşkun G, Talu B, Bek N, Bayramlar KY. Effects of hallux valgus deformity on rear foot position, pain, function, and quality of life of women. J Phys Ther Sci. (2016) 28:781–7. doi: 10.1589/jpts.28.781

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Roddy E, Zhang W, Doherty M. Prevalence and associations of hallux valgus in a primary care population. Arthritis Rheum. (2008) 59:857–62. doi: 10.1002/art.23709

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Cho NH, Kim S, Kwon DJ, Kim HA. The prevalence of hallux valgus and its association with foot pain and function in a rural Korean community. J Bone Joint Surg Br. (2009) 91:494–8. doi: 10.1302/0301-620X.91B4.21925

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Badlissi F, Dunn JE, Link CL, et al. Foot musculoskeletal disorders, pain, and foot-related functional limitation in older persons. J Am Geriatr Soc. (2005) 53:1029–33. doi: 10.1111/j.1532-5415.2005.53315.x

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Wagner E, Wagner P. Metatarsal pronation in hallux valgus deformity: a review. J Am Acad Orthop Surg Glob Res Rev. (2020) 4:e20.00091. doi: 10.5435/JAAOSGlobal-D-20-00091

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Wong DW-C, Wang Y, Chen TL-W, et al. Finite element analysis of generalized ligament laxity on the deterioration of hallux valgus deformity (Bunion). Front Bioeng Biotechnol. (2020) 8:571192–571192. doi: 10.3389/fbioe.2020.571192

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Zeidan H, Ryo E, Suzuki Y, et al. Detailed analysis of the transverse arch of hallux valgus feet with and without pain using weightbearing ultrasound imaging and precise force sensors. PLoS ONE. (2020) 15:e0226914–e0226914. doi: 10.1371/journal.pone.0226914

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Su Q, Zhao R, Wang S, et al. Identification and therapeutic outcome prediction of cervical spondylotic myelopathy based on the functional connectivity from resting-state functional MRI data: a preliminary machine learning study. Front Neurol. (2021) 12:711880. doi: 10.3389/fneur.2021.711880

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Zhao R, Guo X, Wang Y, et al. Functional MRI evidence for primary motor cortex plasticity contributes to the disease's severity and prognosis of cervical spondylotic myelopathy patients. Eur Radiol. (2022). doi: 10.1007/s00330-021-08488-3

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Zhao R, Song Y, Guo X, et al. Enhanced information flow from cerebellum to secondary visual cortices leads to better surgery outcome in degenerative cervical myelopathy patients: a stochastic dynamic causal modeling study with functional magnetic resonance imaging. Front Hum Neurosci. (2021) 15:632829. doi: 10.3389/fnhum.2021.632829

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Xiang L, Mei Q, Wang A, et al. Evaluating function in the hallux valgus foot following a 12-week minimalist footwear intervention: a pilot computational analysis. J Biomech. (2022) 132:110941. doi: 10.1016/j.jbiomech.2022.110941

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Schneider W. Distal soft tissue procedure in hallux valgus surgery: biomechanical background and technique. Int Orthop. (2013) 37:1669–75. doi: 10.1007/s00264-013-1959-5

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Fillingim RB, Loeser JD, Baron R, Edwards RR. Assessment of chronic pain: domains, methods, and mechanisms. J Pain. (2016) 17:T10–20. doi: 10.1016/j.jpain.2015.08.010

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Treede RD, Rief W, Barke A, et al. Chronic pain as a symptom or a disease: the IASP classification of chronic pain for the International Classification of Diseases (ICD-11). Pain. (2019) 160:19–27. doi: 10.1097/j.pain.0000000000001384

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Shima H, Okuda R, Yasuda T, et al. Radiographic measurements in patients with hallux valgus before and after proximal crescentic osteotomy. J Bone Joint Surg Am. (2009) 91:1369–76. doi: 10.2106/JBJS.H.00483

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Welck MJ, Al-Khudairi N. Imaging of hallux valgus: how to approach the deformity. Foot Ankle Clin. (2018) 23:183–92. doi: 10.1016/j.fcl.2018.01.002

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Gentili A, Masih S, Yao L, Seeger LL. Pictorial review: foot axes and angles. Br J Radiol. (1996) 69:968–74. doi: 10.1259/0007-1285-69-826-968

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Dellonder-Frigolé J, Matas-Pareja M, Calle-García JA, et al. Revisiting the original Akin technique: a comparative analysis of radiographic outcome. Foot Ankle Surg. (2021). doi: 10.1016/j.fas.2021.05.004

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Engel E, Erlick N, Krems I. A simplified metatarsus adductus angle. J Am Podiatry Assoc. (1983) 73:620–8. doi: 10.7547/87507315-73-12-620

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Osher L, Blazer MM, Buck S, Biernacki T. Accurate determination of relative metatarsal protrusion with a small intermetatarsal angle: a novel simplified method. J Foot Ankle Surg. (2014) 53:548–56. doi: 10.1053/j.jfas.2014.04.021

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Shi GG, Henning P, Marks RM. Correlation of postoperative position of the sesamoids after chevron osteotomy with outcome. Foot Ankle Int. (2016) 37:274–80. doi: 10.1177/1071100715624147

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: hallux valgus, multi-variate pattern analysis, pain, support vector machine, hierarchical clustering

Citation: Wang C, Li C, Zhang R, Li Z, Zhang H, Zhang Y, Liu S, Chi X and Zhao R (2022) Identification of radiographic characteristics associated with pain in hallux valgus patients: A preliminary machine learning study. Front. Public Health 10:943026. doi: 10.3389/fpubh.2022.943026

Received: 13 May 2022; Accepted: 22 July 2022;
Published: 10 August 2022.

Edited by:

Yaodong Gu, Ningbo University, China

Reviewed by:

Xuanzhen Cen, Óbuda University, Hungary
Liangliang Xiang, University of Auckland, New Zealand

Copyright © 2022 Wang, Li, Zhang, Li, Zhang, Zhang, Liu, Chi and Zhao. 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: Rui Zhao, zhaorui@tmu.edu.cn; ChenGuang Wang, WangCG_TJGH@163.com

These authors have contributed equally to this work

Disclaimer: 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.