Abstract
Frogs have a highly conserved body plan, yet they employ a diverse array of locomotor modes, making them ideal organisms for investigating the relationships between morphology and locomotor function, in particular whether anatomical complexity is a prerequisite for functional complexity. We use diffusible iodine contrast-enhanced microCT (diceCT) imaging to digitally dissect the gross muscle anatomy of the pelvis and hindlimbs for 30 species of frogs representing five primary locomotor modes, including the first known detailed dissection for some of the world's smallest frogs, forming the largest digital comparative analysis of musculoskeletal structure in any vertebrate clade to date. By linking musculoskeletal dissections and phylogenetic comparative methods, we then quantify and compare relationships between anatomy and function across over 160 million years of anuran evolution. In summary, we have found that bone lengths and pelvic crest sizes are generally not reliable predictors of muscle sizes, which highlights important implications for future palaeontological studies. Our investigation also presents previously unreported differences in muscle anatomy between frogs specialising in different locomotor modes, including several of the smallest frog hindlimb muscles, which are extremely difficult to extract and measure using traditional approaches. Furthermore, we find evidence of many-to-one and one-to-many mapping of form to function across the phylogeny. Additionally, we perform the first quantitative analysis of how the degree of muscle separation can differ between frogs. We find evidence that phylogenetic history is the key contributing factor to muscle separation in the pelvis and thigh, while the separation of shank muscles is influenced more strongly by locomotor mode. Finally, our anatomical 3D reconstructions are published alongside this manuscript to contribute towards future research and serve as educational materials.