Morphospace analysis leads to an evo-devo model of digit patterning

Biological forms occupy a constrained portion of theoretical morphospaces. Developmental models accounting for empirical morphospaces are necessary to achieve a better understanding of this phenomenon. We analyzed the phalangeal formulas (PFs) in lizards and relatives' hands by comparing them with a set of simulated PFs that compose a theoretical morphospace. We detected that: (1) the empirical morphospace is severely limited in size, (2) the PFs comply with two properties of phalangeal count per digit, namely the ordering rule (DI ≤ DII ≤ DIII ≤ DIV ≥ DV), and the contiguity relationship (neighbor digits differ on average in one phalanx), (3) the totality of the PFs can be categorized into four categories of hands aligned along a feasibility gradient. We also reconstructed the evolution of PFs and found a stepwise trajectory from the plesiomorphic PF towards reduced conditions. Finally, we propose a developmental model as the generative mechanism behind the PFs. It is consistent with the bulk of evidence managed and involves an ordered digit primordia initialization timed with periodic signals of joint formation coming from digit tips. Our approach is also useful to address the study of other meristic sequences in nature such as dental, floral, and branchial formulas.