Evolution of a complex phenotype with biphasic ontogeny: contribution of development versus function and climatic variation to skull modularity in toads

The theory of morphological integration and modularity predicts that if functional correlations among traits are relevant to mean population fitness, the genetic basis of development will be molded by stabilizing selection to match functional patterns. Yet, how much functional interactions actually shape the fitness landscape is still an open question. We used the anuran skull as a model of a complex phenotype for which we can separate developmental and functional modularity. We hypothesized that functional modularity associated to functional demands of the adult skull would overcome developmental modularity associated to bone origin at the larval phase because metamorphosis would erase the developmental signal. We tested this hypothesis in toad species of the Rhinella granulosa complex using species phenotypic correlation pattern (P-matrices). Given that the toad species are distributed in very distinct habitats and the skull has important functions related to climatic conditions, we a...