Sexual dimorphism and sex-specific features of pediatric hand bone shapes via statistical shape modeling
Authors: Edwin Tay, Amir A. Zadpoor, Nazli Tümer
Organizations: Department of Biomechanical Engineering, Faculty of Mechanical Engineering, Delft University of Technology (TU Delft), Mekelweg 2, 2628CD, Delft, The Netherlands
Understanding the envelope of shape variation encompassing both sexual dimorphism and growth-related changes in living tissue is a non-trivial task. In this work, we use statistical shape models (SSMs) and analyze demographic-stratified shape parameters to quantify sexual dimorphism and sex-specific shape changes. To this end, we developed a two-stage detection-segmentation framework that accurately segmented n=27664 unique bone masks from a radiographic dataset of n=1456 subjects with diverse ages and sexes (HD95 = 1.512±0.005, ASSD 0.449±001, Dice = 0.987±0.00004, IoU = 0.975±0.0001). Following SSM fitting, our follow-up analyses of stratified shape parameters highlighted the morphology of sexual dimorphisms for each bone and demonstrated sex-specific differences in shape development. Our results were consistent with existing literature and provided an unambiguous, data-driven methodology for characterizing tissue shape and its population-scale variation. This approach provides a potential methodological foundation for developing population-based models of anatomical variation and could be extended to investigate tissue growth and shape changes over time. Future work should evaluate the applicability of this methodology on alternative anatomical structures and imaging modalities, with the goal of validating these growth trajectories on longitudinal data.
Table A1: Age term elected per shape mode for bones of the third ray. Each cell names the lowest-BIC candidate for that bone and mode, and its BIC margin over the runner-up. cub denotes cubic, quad denotes quadratic, and lin denotes linear. ∼ marks a margin below 2, which is no meaningful preference.
Figure 1: Comparison of the different fits of varying forms for the age terms. Each subfigure title denotes the difference in the Bayesian information criterion (BIC); a lower BIC indicates a better fit.
Bone
modes
T2
F
df
D
p
Distal phalanx 1
PC1–8
61.712
7.677
8, 1447
0.4118
p≤0.001
Distal phalanx 2
PC1–7
65.67
9.343
7, 1448
0.4248
p≤0.001
Distal phalanx 4
PC1–7
51.661
7.35
7, 1448
0.3767
p≤0.001
Distal phalanx 5
PC1–7
69.628
9.906
7, 1448
0.4374
p≤0.001
Medial phalanx 1
PC1–7
108.65
15.46
7, 1448
0.5464
p≤0.001
Medial phalanx 3
PC1–7
97.191
13.83
7, 1448
0.5167
p≤0.001
Appendix
Table B1: Hotelling’s T2 over sex-stratified shape parameters encompassing c shape modes retained for analysis. F and p -values denote statistical significance whilst T2 and D quantifies the distance between males and females. Results indicate that all bones are sexually dimorphic.
Bone
Pillai
s
Pillai /s
F
df
p ( F )
p (Freedman–Lane)
Distal phalange 1
0.07335
3
0.02445
4.521
24, 4329
p≤0.001
p≤0.001
Distal phalange 2
0.07228
2
0.03614
7.74
14, 2890
p≤0.001
p≤0.001
Distal phalange 4
0.1063
3
0.03542
7.576
21, 4332
p≤0.001
p≤0.001
Distal phalange 5
0.1161
3
0.03869
8.302
21, 4332
p≤0.001
p≤0.001
Medial phalange 1
0.1503
3
0.05008
10.88
21, 4332
p≤0.001
p≤0.001
Medial phalange 3
0.2033
3
0.06777
15
21, 4332
p≤0.001
p≤0.001
Appendix
Table B2: Results from Pillai’s trace on sex × age interaction. Raw Pillai’s trace magnitudes are normalized with respect to s degrees of freedom, giving Pillai /s . We report raw p -values based on analytic F , but because residuals are heteroskedastic, we used permutation tests via Freedman-Lane sampling with Rademacher sign flipping to obtain a more robust indicator of significance. Results indicate that all bones have sex-specific ontogeny.
Dept. of Computer Science and Engineering Özyeğgin University İstanbul, Türkiye · Dept. of Artificial Intelligence and Data Engineering Özyeğgin University İstanbul, Türkiye · Osteoid Inc. İstanbul, Türkiye