Abstract
Adaptive variation in the craniofacial skeleton is a key component of resource specialization and habitat divergence in vertebrates, but the proximate genetic mechanisms that underlie complex patterns of craniofacial variation are largely unknown. Here we demonstrate that the Hedgehog (Hh) signaling pathway mediates widespread variation across a complex functional system that affects the kinematics of lower jaw depression--the opercular four-bar linkage apparatus--among Lake Malawi cichlids. By using a combined quantitative trait locus mapping and population genetics approach, we show that allelic variation in the Hh receptor, ptch1, affects the development of distinct bony elements in the head that represent two of three movable links in this functional system. The evolutionarily derived allele is found in species that feed from the water column, and is associated with shifts in anatomy that translate to a four-bar system capable of faster jaw rotation. Alternatively, the ancestral allele is found in species that feed on attached algae, and is associated with the development of a four-bar system that predicts slower jaw movement. Experimental manipulation of the Hh pathway during cichlid development recapitulates functionally salient natural variation in craniofacial geometry. In all, these results significantly extend our understanding of the mechanisms that fine-tune the craniofacial skeletal complex during adaptation to new foraging niches.