Heteroblasty refers to the changes in leaf shape and size (allometry) along stems. Although evolutionary changes involving heteroblasty might contribute to leaf diversity, little is known of the extent to which heteroblasty differs between species or how it might relate to other aspects of allometry or other developmental transitions. Here, we develop a computational model that can quantify differences in leaf allometry between Antirrhinum (snapdragon) species, including variation in heteroblasty. It allows the underlying genes to be mapped in inter-species hybrids, and their effects to be studied in similar genetic backgrounds. Heteroblasty correlates with overall variation in leaf allometry, so species with smaller, rounder leaves produce their largest leaves earlier in development. This involves genes that affect both characters together and is exaggerated by additional genes with multiplicative effects on leaf size. A further heteroblasty gene also alters leaf spacing, but none affect other developmental transitions, including flowering. We suggest that differences in heteroblasty have co-evolved with overall leaf shape and size in Antirrhinum because these characters are constrained by common underlying genes. By contrast, heteroblasty is not correlated with other developmental transitions, with the exception of internode length, suggesting independent genetic control and evolution.