Geographical patterns of morphological variation and environmental correlates in contact zones: a multi-scale approach using two Mediterranean vipers (Serpentes)

Adaptation to similar selective pressures can explain morphological convergence between closely related species in contact zones. Geostatistics and Geographical Information System were used to identify multi-scale patterns of morphological variability and test the hypothesis of morphological convergence due to local environmental pressures in the contact areas between Vipera aspis and Vipera latastei in the Iberian Peninsula. Nine morphological traits from 630 and 362 vipers at regional and local scale, respectively, were interpolated by Kriging to generate surfaces of morphological variation. Kriging is a geostatistical algorithm that allows investigating the spatial structure of data with statistical models. At both scales, a convergent north–south pattern in morphological variability was observed and the contact areas were identified as integration zones where intermediate vipers are found. Significant correlations were found between surfaces of univariate and multivariate traits, with precipitation and temperature seasonality. Thus, several morphological traits were apparently under local environmental selection. Nevertheless, the influence of biotic pressures and gene flow on morphological convergence of vipers in contact zones deserves further study.