Species delimitation is often based on a single or very few genetic or phenetic traits, something which leads to misinterpretations and often does not provide information about evolutionary processes. Here, we investigated the diversity pattern of multiple phenetic traits of the two extant species of Salamandrina, a genus split only after molecular traits had been studied but the two species of which are phenetically very similar. The phenetic traits we studied are size, external body shape and head colour pattern, in a model comparison framework using non-linear mixed models and unsupervised and supervised clustering. Overall, we found high levels of intra-specific variability for body size and shape, depending on population belonging and habitat, while differences between species were generally lower. The habitat the salamanders dwell in also seems important for colour pattern. Basing on our findings, from the methodological point of view, we suggest (i) to take into account the variability at population level when testing for higher level variability, and (ii) a semi-supervised learning approach to high dimensional data. We also showed that different phenotypic traits of the same organism could result from different evolutionary routes. Local adaptation is likely responsible for body size and shape variability, with selective pressures more similar across species than within them. Head colour pattern also depends on habitat, differently from ventral colour pattern (not studied in this paper) which likely evolved under genetic drift.
Keywords: evolutionary processes; image analysis; local adaptation; multivariate clustering; phenetic characters.