RESUMEN
Pygopodids are elongate, functionally limbless geckos found throughout Australia. The clade presents low taxonomic diversity (â¼45 spp.), but a variety of cranial morphologies, habitat use, and locomotor abilities that vary between and within genera. In order to assess potential relationships between cranial morphology and ecology, computed tomography scans of 29 species were used for 3D geometric morphometric analysis. A combination of 24 static landmarks and 20 sliding semi-landmarks were subjected to Generalized Procrustes Alignment. Disparity in cranial shape was visualized through Principal Component Analysis, and a multivariate analysis of variance (MANOVA) was used to test for an association between shape, habitat, and diet. A subset of 27 species with well-resolved phylogenetic relationships was used to generate a phylomorphospace and conduct phylogeny-corrected MANOVA. Similar analyses were done solely on Aprasia taxa to explore species-level variation. Most of the variation across pygopodids was described by principal component (PC) 1(54%: cranial roof width, parabasisphenoid, and occipital length), PC2 (12%: snout elongation and braincase width), and PC3 (6%: elongation and shape of the palate and rostrum). Without phylogenetic correction, both habitat and diet were significant influencers of variation in cranial morphology. However, in the phylogeny-corrected MANOVA, habitat remained weakly significant, but not diet, which can be explained by generic-level differences in ecology rather than among species. Our results demonstrate that at higher levels, phylogeny has a strong effect on morphology, but that influence may be due to small sample size when comparing genera. However, because some closely related taxa occupy distant regions of morphospace, diverging diets, and use of fossorial habitats may contribute to variation seen in these geckos.